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

Turbulent inflow generation logic patch #1695

Merged
merged 57 commits into from
Jul 22, 2024

Conversation

dustinma324
Copy link
Contributor

@dustinma324 dustinma324 commented Jul 17, 2024

This patch modifies the logic in two functions, calc_tpi_update() and apply_tpi(). In this version, a MultiFab is created within DataStructs/TurbPertStruct.H to store the random perturbation amplitude per cell within the perturbation box. The randomization step is now pulled out of apply_tpi(), which was incorrectly being called at every iteration, and placed within calc_tpi_update().

With this change, procedures such as calculating mean box velocity, perturbation amplitude, and assigning random perturbations are only triggered when the update interval is surpassed by the local elapsed time per box. This ensures that the randomization of perturbation assignment within the perturbation box occurs only once per update interval. The Array4 that stores the random perturbation is then used in apply_tpi() per iteration as a direct addition to the source term of the desired field.

  • This patch does not alter the behavior of acoustic waves when applying perturbations to the source term of the rhoTheta field.
  • Incompressible box perturbation now working with movie provided
movie_incompressible_BPM.mp4

Ting-Hsuan (Dustin) Ma added 30 commits June 25, 2024 08:07
…ariables. Renamed BPM and CPM tags to source and direct to keep neutrality. Need to add Ec, tp, and cell_data equations for direction add.
…ithin TimeIntegration/ERF_Advance.cpp for CPM cases. Need to implement options to swap back and forth on use side between BPM and CPM. Switch gear and start runtime optimization.
…ded outputs for PB_amplitude/PB_updateTime in seperate file when erf.v=1. Need to find a good way to generalize amplitude. Currently breaks code when amplitude is too large. CONTINUE with CPM debug, fast_rhs_fun debug with PB, and documentation with .rst file.
…n debug with PB, and documentation with .rst file
… boxes. Added in perturbation inflow .rst documents. CONTINUE with CPM debug, fast_rhs_fun debug with PB.
… boxes. Added in perturbation inflow .rst documents. CONTINUE with CPM debug, fast_rhs_fun debug with PB.
… format into file + verbose tags. Fixed parallel issue of boxes. Added in perturbation inflow .rst documents. CONTINUE with CPM debug, fast_rhs_fun debug with PB.
…ose tags. Fixed parallel issue of boxes. Added in perturbation inflow .rst documents. CONTINUE with CPM debug, fast_rhs_fun debug with PB.
…Currently blows up after 1st update interval.
… (2015). DOI links return Erorror 403, actual link returns Error 404.
@dustinma324 dustinma324 changed the title Turbulence inflow generation logic patch Turbulent inflow generation logic patch Jul 17, 2024
@dustinma324 dustinma324 force-pushed the perturbation_box_method branch from 2acb16e to d9d7c86 Compare July 22, 2024 22:16
Ting-Hsuan (Dustin) Ma added 2 commits July 22, 2024 15:19
@AMLattanzi AMLattanzi merged commit c011b2e into erf-model:development Jul 22, 2024
17 of 18 checks passed
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