Turbulent forcing#535
Turbulent forcing#535baperry2 merged 32 commits intoAMReX-Combustion:developmentfrom efhunt:TurbulentForcing
Conversation
…art of a simulation. (Essential for Maintained HIT workflow)
|
Thanks for the PR and for including documentation and test cases to demonstrate the new capability! It may be a bit before we can review this thoroughly, but there are a few things you can start on now to get things moving along:
|
…me doing this so may not work)
…t case. Not added FlameSheet turbulent forcing as I do not think its testing anything additional than HIT.
|
Thank you for your response and suggestions. The amendments have been made and I have integrated the test cases into HITDecay (now renamed HIT) and FlameSheet. |
baperry2
left a comment
baperry2
left a comment
There was a problem hiding this comment.
Thanks again for this contribution. Looks pretty good, mostly there are just some extraneous files that can be deleted.
…uggested by dmontgomeryNREL)
…equested by baperry2)
|
I am not quite sure why there is a branch conflict and how to resolve it? Any advice or help would be appreciated. |
|
Git merges can be a bit weird with the submodules. I resolved on the command line by doing a merge with development than updating and committing the PelePhysics submodule to the latest version from that PR. I'll merge both PRs once the tests all pass. |
|
Fantastic, thanks! |
3 participants
This PR adds the required hooks to use the turbulent forcing functions in PelePhysics (requires PR#588). The turbulent forcing is used to generate and maintain homogeneous isotropic turbulence. The method used is the same as that used in PeleLM and IAMR (Andrew Aspden's forcing, see Aspden et al (CAMCoS 2008)).
Additionally, it is an essential component of the workflow to start a flame sheet with already developed HIT. Therefore functionality to read in the velocities from an existing plotfile and add them to the state data is provided.
Example configurations have been provided in Exec/RegTests/TurbForce_Incompressible for generating the inital turbulent velocity plotfile and TurbForce_FlameSheet for a turbulent flame sheet.
The code has been tested using PeleLMeX on CPUs (by myself gcc (Debian 12.2.0-14+deb12u1) 12.2.0)) and Max Schneider on GPUs (ml Stages/2024 NVHPC/25.1-CUDA-12 OpenMPI/4.1.5 CMake/3.26.3)