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

Add a Physics for THM compressible Euler flow model #28528

Draft
wants to merge 95 commits into
base: next
Choose a base branch
from

Conversation

GiudGiud
Copy link
Contributor

@GiudGiud GiudGiud commented Sep 3, 2024

I hope to figure out at least with VACE:

  • channels
  • inlets/outlets
  • simple junctions

with FV

  • channels
  • inlets/outlets
  • simple junctions
  • orienting pipes
  • heat
  • scalars

I ll remove the WCNSFV part from the PR and put it into another PR. I m just trying things here to make sure the design will suit our needs

@GiudGiud GiudGiud self-assigned this Sep 3, 2024
Copy link
Contributor

@joshuahansel joshuahansel left a comment

Choose a reason for hiding this comment

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

Just some initial thoughts to get started.

[]

[pipe]
type = PhysicsFlowChannel
Copy link
Contributor

Choose a reason for hiding this comment

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

Is it possible to have physics = myflowphysics here?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

yes, I started differently but I ll correct course

[inlet]
type = PhysicsInletMassFlowRateTemperature
input = 'pipe:in'
physics = 'all'
Copy link
Contributor

Choose a reason for hiding this comment

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

Can we get rid of this parameter? The BC component should be able to query the connected component.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

good point

/**
* Sets up the single-phase flow equations using Euler's equations and a RDG discretization
*/
class THMWCNSFVFlowPhysics : public virtual ThermalHydraulicsFlowPhysics, public virtual WCNSFVFlowPhysics
Copy link
Contributor

Choose a reason for hiding this comment

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

This will need a test of course

Copy link
Contributor Author

Choose a reason for hiding this comment

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

oh this is like super not ready. will do

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 realized that later in the review 😄


// boundary flux user object
// we add them in addBCs for convenience
if (boundary_type == InletTypeEnum::MdotTemperature)
Copy link
Contributor

Choose a reason for hiding this comment

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

Would someone be able to add a new BC formulation in their application without modifying THM?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

yes, by using a XYZFlowChannel instead of the Physics one
or by creating a new physics

}

void
THMVACESinglePhaseFlowPhysics::addInletBoundaries()
Copy link
Contributor

Choose a reason for hiding this comment

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

Do we need to distinguish inlet/outlet, or can we just have a generic addFlowBCs() or whatever?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

modeling after WCNSFV here. in wcnsfv we do need to distinguish

Copy link
Contributor

@joshuahansel joshuahansel Sep 3, 2024

Choose a reason for hiding this comment

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

Sometimes the distinction can be blurry, like with reversible, and then there are walls.

@GiudGiud
Copy link
Contributor Author

GiudGiud commented Sep 3, 2024

thanks for the early feedback!

Suppress boundary parameters
This is better, if someone adds addUO in the TH physics base class, it wont compile until they clear the non-trivial override
Rename channel test inputs for convenience

We could make this more general, as in ICs for any variable
…table

Maybe it s time we make BlockRestrictableDynamic
… a2) type terms

as volumetric sources in 1D flow channels
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants