-
Notifications
You must be signed in to change notification settings - Fork 57
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 Surface Tension #387
Add Surface Tension #387
Conversation
…o Surface-Tension
This reverts commit 91acd10.
@wilfonba check s_rectangle subroutine in m_patches.fpp, it only defines smooth_patch_id when ib is not true (lines 893-894), and you added a further check at line 928 where that variable is no longer defined, which is causing a bug in the 3D IBM case |
I forgot I'd added that stuff. Thanks for taking a look @anandrdbz! |
@wilfonba sync w/ master |
…o Surface-Tension
…o Surface-Tension
should be able to debug these |
c29ed36
to
0092aeb
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
lgtm on my side.
- I fixed
./mfc.sh count_diff
and overall improved how it's written. - I reverted the changes you made in the
3D_weak_scaling
case where you hardcoded the time-stepping parameters.
Note: |
Let me fix the macOS no debug runners. Not sure why those would have failed, but I haven't ran the test suite with no debug recently |
Related to new issue #443 |
Can merge this PR without that runner working |
Description
This PR adds the surface tension model of Schmidmayer et. al 2017 to MFC.
Fixes #276
Type of change
Scope
How Has This Been Tested?
x[y]_beg
and ghost cell extrapolation forx[y]_end
. The case file is available inexamples/2D_equilibration/
. This simulation was run on M2 Apple Silicone.test.mp4
x[y]_beg
and ghost cell extrapolation forx[y]_end
. The case file is available inexamples/2D_recoverCircle
. This simulation was run on 2 V100 GPUs.test.mp4
x[y,z]_beg
and ghost cell extrapolation atx[yz]_end
. This simulation was ran using 4 V100 GPUs on Pace Phoenix.test.mp4
x[y,z]_beg
and ghost cell extrapolation atx[yz]_end
. This simulation was ran using 8 MI250 GPUs on OLCF Frontier. While equilibration is not reached, the relaxation shown is the asymmetric and the expected behavior. One octant of the cube is simulated and reflections used for visualization. This case is given as an example inexamples/3D_surface_tension
.test.mp4
Checklist
docs/
)examples/
that demonstrate my new feature performing as expected./mfc.sh format
before committing my codeIf your code changes any code source files (anything in
src/
)To make sure the code is performing as expected on GPU devices, I have:
nvtx
ranges so that they can be identified in profiles./mfc.sh run XXXX --gpu -t simulation --nsys
, and have attached the output file (.nsys-rep
) and plain text results to this PR nsysRaw.txtexamples/3D_weak_scaling
with 32GB per GPU with 8 MI250 GPUs on both the master branch and this fork with the following time/time-step:Master: 0.555s
This Fork: 0.539s
I'm not sure why this fork is faster, but these results were consistent across three runs.