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

Multiphase hybrid #915

Merged
merged 49 commits into from
Oct 6, 2023
Merged

Multiphase hybrid #915

merged 49 commits into from
Oct 6, 2023

Conversation

mbkuhn
Copy link
Contributor

@mbkuhn mbkuhn commented Oct 5, 2023

lots of optional features, some will be removed in the future. Best working version is default.

-- commenting or removing extra stuff
* also, found 1.5 * i_th works for pseudo-CFL
- corresponding version of levelset2vof
-huge bug fix for this method
@mbkuhn mbkuhn requested a review from psakievich October 5, 2023 19:37
@mbkuhn mbkuhn merged commit b0b00f3 into Exawind:multiphase_dev Oct 6, 2023
mbkuhn added a commit to mbkuhn/amr-wind that referenced this pull request Oct 11, 2023
* option to sharpen nalu vof field or not

* putting previous commit into practice

* important options for sharpening and masking

* able to turn off mac projection

* getting toward a working sharpening
-- sloppy commit

* provides stability for very small epsilon (interface thickness)

* removing unnecessary things, clean up

* more unnecessary

* default to a smaller interface thickness

* masking fluxes is sufficient, no need to mask whole term

* density fluxes and update (to make sure it works)

* best upwinding approach so far

* begin to look at the rest of the iteration

* ability to turn off strong overset coupling to nodal projection

* initialize pressure in sloshing tank case

* modify pressure by hydrostatic source term
-- commenting or removing extra stuff

* momentum sharpening as included in paper
* also, found 1.5 * i_th works for pseudo-CFL

* to go with last change

* ability to specify tolerance for convergence, not just steps

* to go with previous commit

* keep vof bounded

* initialize pressure field for other cases (like waves)

* specify relative length scale by input argument

* move copy to old, misc

* Neumann version of youngs normal
- corresponding version of levelset2vof

* discrete sharpening, much improved from other one

* likely unused stuff, add for history

* more reasonable asdf function, better for small nalu domains

* more flexible flux computation and commented debugging stuff

* populate ghost cells for iblank fields
-huge bug fix for this method

* making vof margin for flux upwinding/downwinding an input argument

* changing defaults of sharpening tool

* commented print statements for debugging

* include boundedness of vof and recalc of density

* new FuzzyInterface case

* calculate face velocity in the upwinding/downwinding form like the vof fluxes

* changing the sharpen margin default
* essential to better results with flat interface

* do density adjustment as intended

* options to ignore nalu pressure and to turn off pressure source term in sharpening

* cleaning up the code and putting in more options

* replacing the pressure gradient with hydrostatic
-- should make this optional

* reapply hydrostatic guess of pressure gradient

* turning latest commits into input options

* switching to full upwinding for velocity sharpening

* ramping the pseudo-time seems to work better
- less noisy velocity field

* ability to sharpen pressure gradient

* 3 big changes
- calculate pseudo dt based on pseudo CFL
- make routine suited to multiple levels
- remove unnecessary parts of  disable_pressure_from_nalu

* set current working version of options as default
mbkuhn added a commit to mbkuhn/amr-wind that referenced this pull request Oct 11, 2023
* option to sharpen nalu vof field or not

* putting previous commit into practice

* important options for sharpening and masking

* able to turn off mac projection

* getting toward a working sharpening
-- sloppy commit

* provides stability for very small epsilon (interface thickness)

* removing unnecessary things, clean up

* more unnecessary

* default to a smaller interface thickness

* masking fluxes is sufficient, no need to mask whole term

* density fluxes and update (to make sure it works)

* best upwinding approach so far

* begin to look at the rest of the iteration

* ability to turn off strong overset coupling to nodal projection

* initialize pressure in sloshing tank case

* modify pressure by hydrostatic source term
-- commenting or removing extra stuff

* momentum sharpening as included in paper
* also, found 1.5 * i_th works for pseudo-CFL

* to go with last change

* ability to specify tolerance for convergence, not just steps

* to go with previous commit

* keep vof bounded

* initialize pressure field for other cases (like waves)

* specify relative length scale by input argument

* move copy to old, misc

* Neumann version of youngs normal
- corresponding version of levelset2vof

* discrete sharpening, much improved from other one

* likely unused stuff, add for history

* more reasonable asdf function, better for small nalu domains

* more flexible flux computation and commented debugging stuff

* populate ghost cells for iblank fields
-huge bug fix for this method

* making vof margin for flux upwinding/downwinding an input argument

* changing defaults of sharpening tool

* commented print statements for debugging

* include boundedness of vof and recalc of density

* new FuzzyInterface case

* calculate face velocity in the upwinding/downwinding form like the vof fluxes

* changing the sharpen margin default
* essential to better results with flat interface

* do density adjustment as intended

* options to ignore nalu pressure and to turn off pressure source term in sharpening

* cleaning up the code and putting in more options

* replacing the pressure gradient with hydrostatic
-- should make this optional

* reapply hydrostatic guess of pressure gradient

* turning latest commits into input options

* switching to full upwinding for velocity sharpening

* ramping the pseudo-time seems to work better
- less noisy velocity field

* ability to sharpen pressure gradient

* 3 big changes
- calculate pseudo dt based on pseudo CFL
- make routine suited to multiple levels
- remove unnecessary parts of  disable_pressure_from_nalu

* set current working version of options as default
mbkuhn added a commit to mbkuhn/amr-wind that referenced this pull request Oct 11, 2023
* option to sharpen nalu vof field or not

* putting previous commit into practice

* important options for sharpening and masking

* able to turn off mac projection

* getting toward a working sharpening
-- sloppy commit

* provides stability for very small epsilon (interface thickness)

* removing unnecessary things, clean up

* more unnecessary

* default to a smaller interface thickness

* masking fluxes is sufficient, no need to mask whole term

* density fluxes and update (to make sure it works)

* best upwinding approach so far

* begin to look at the rest of the iteration

* ability to turn off strong overset coupling to nodal projection

* initialize pressure in sloshing tank case

* modify pressure by hydrostatic source term
-- commenting or removing extra stuff

* momentum sharpening as included in paper
* also, found 1.5 * i_th works for pseudo-CFL

* to go with last change

* ability to specify tolerance for convergence, not just steps

* to go with previous commit

* keep vof bounded

* initialize pressure field for other cases (like waves)

* specify relative length scale by input argument

* move copy to old, misc

* Neumann version of youngs normal
- corresponding version of levelset2vof

* discrete sharpening, much improved from other one

* likely unused stuff, add for history

* more reasonable asdf function, better for small nalu domains

* more flexible flux computation and commented debugging stuff

* populate ghost cells for iblank fields
-huge bug fix for this method

* making vof margin for flux upwinding/downwinding an input argument

* changing defaults of sharpening tool

* commented print statements for debugging

* include boundedness of vof and recalc of density

* new FuzzyInterface case

* calculate face velocity in the upwinding/downwinding form like the vof fluxes

* changing the sharpen margin default
* essential to better results with flat interface

* do density adjustment as intended

* options to ignore nalu pressure and to turn off pressure source term in sharpening

* cleaning up the code and putting in more options

* replacing the pressure gradient with hydrostatic
-- should make this optional

* reapply hydrostatic guess of pressure gradient

* turning latest commits into input options

* switching to full upwinding for velocity sharpening

* ramping the pseudo-time seems to work better
- less noisy velocity field

* ability to sharpen pressure gradient

* 3 big changes
- calculate pseudo dt based on pseudo CFL
- make routine suited to multiple levels
- remove unnecessary parts of  disable_pressure_from_nalu

* set current working version of options as default
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.

1 participant