Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

NEMO: Native Air-7 gas model implementation for inviscid flows #1555

Merged
merged 36 commits into from
Mar 10, 2022

Conversation

WallyMaier
Copy link
Contributor

@WallyMaier WallyMaier commented Mar 4, 2022

Proposed Changes

@jtneedels and I have worked on implementing the Air-7 gas model to the NEMO native fluid model. This allows for simulations of inviscid weakly ionized flows with the addition of NO+ and e-.

I am open to adding a test case for this, similar to the thermal bath case, since they run extremely fast, though I am open to idea.

Work is ongoing to extend this viscous simulations.
Updates from #1533 are added here as well.

PR Checklist

  • I am submitting my contribution to the develop branch.
  • My contribution generates no new compiler warnings (try with the '-Wall -Wextra -Wno-unused-parameter -Wno-empty-body' compiler flags, or simply --warnlevel=2 when using meson).
  • My contribution is commented and consistent with SU2 style.
  • I have added a test case that demonstrates my contribution, if necessary.
  • I have updated appropriate documentation (Tutorials, Docs Page, config_template.cpp) , if necessary.

@pcarruscag
Copy link
Member

Beyoncé's rule of programming, if you like it put a test on it.

@@ -137,9 +136,6 @@ CSU2TCLib::CSU2TCLib(const CConfig* config, unsigned short val_nDim, bool viscou
Tcb_b.resize(nReactions,0.0);

/*--- Assign gas properties ---*/
// Wall mass fractions for catalytic boundaries
Wall_Catalycity[0] = 0.999;
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

So did we make a decision on where to keep this data? I think in my very stale for the cat wall branch I had it as a config option, which I think is a good choice, but would it be better to store it here in fluid model? I think I remember getting push back for having it as a config option.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hmmm, I do recall the issue. I'm not sure what's the best way forward here. Either way, hardcoding values doesn't seem like the right move. (If I understand catalyticity correct....)

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Chatted about this, I think we agreed to use the free-stream values as the default, rather than having hard-coded values in fluid model, so I agree with removal here.


//TODO: Implement Collision Integral Data for AIR-7 (JN)
/*--- Collision integral data ---*/
//TODO: ELECTRONS HAvE NO OMEGAS....NEED TO ALTER LOOPS? or add dummy vector.....
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think the general reasoning is electrons are too small to contribute meaningfully to transport properties, so their collision integrals are effectively zero. A dummy vector makes sense I think.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

👍 Ill pull more from Scalabrin and park during the viscous implementation.

@@ -1627,7 +2051,7 @@ vector<su2double>& CSU2TCLib::ComputeTemperatures(vector<su2double>& val_rhos, s
void CSU2TCLib::GetChemistryEquilConstants(unsigned short iReaction){

if (gas_model == "O2"){

// THESE ARE UNUSED. SHOULD WE KEEP???? Good for future?
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yeah I'd keep.

@WallyMaier WallyMaier merged commit 4d2854a into develop Mar 10, 2022
@WallyMaier WallyMaier deleted the feature_NEMO_chem branch March 10, 2022 16:46
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants