Skip to content

Navigation Menu

Sign in
Appearance settings

Search code, repositories, users, issues, pull requests...

Provide feedback

We read every piece of feedback, and take your input very seriously.

Saved searches

Use saved searches to filter your results more quickly

Appearance settings

Update initialization for turbulent mixing layer and add validation results #879

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

Open
wants to merge 26 commits into
base: master
Choose a base branch
Loading
from

Conversation

hyeoksu-lee
Copy link
Contributor

@hyeoksu-lee hyeoksu-lee commented Jun 11, 2025

Description

This PR includes an update on initialization method for turbulent mixing layer. Also, Matlab scripts for computing turbulence statistics are added to examples/3D_turb_mixing/turbulence_stat as well as validation results.

Type of change

Please delete options that are not relevant.

  • Bug fix (non-breaking change which fixes an issue)
  • New feature (non-breaking change which adds functionality)
  • Something else

Scope

  • This PR comprises a set of related changes with a common goal

If you cannot check the above box, please split your PR into multiple PRs that each have a common goal.

How Has This Been Tested?

Ran examples/3D_turb_mixing on Carpenter

Streamwise velocity field at final time
vel1

Vorticity magnitude at final time
omega_mag

Evolution of Reynolds stress
tstep_0

Evolution of TKE budget
tstep_0

TKE budget averaged over self-similar period
avg_self_similar

Momentum thickness growth rate (see run_turbulence.m for references in legend)
growth_rate

Checklist

  • I have added comments for the new code
  • I added Doxygen docstrings to the new code
  • I have made corresponding changes to the documentation (docs/)
  • I have added regression tests to the test suite so that people can verify in the future that the feature is behaving as expected
  • I have added example cases in examples/ that demonstrate my new feature performing as expected.
    They run to completion and demonstrate "interesting physics"
  • I ran ./mfc.sh format before committing my code
  • New and existing tests pass locally with my changes, including with GPU capability enabled (both NVIDIA hardware with NVHPC compilers and AMD hardware with CRAY compilers) and disabled
  • This PR does not introduce any repeated code (it follows the DRY principle)
  • I cannot think of a way to condense this code and reduce any introduced additional line count

@hyeoksu-lee hyeoksu-lee requested review from a team and sbryngelson as code owners June 11, 2025 19:32
Copy link

codecov bot commented Jun 11, 2025

Codecov Report

Attention: Patch coverage is 76.05634% with 17 lines in your changes missing coverage. Please review.

Project coverage is 43.71%. Comparing base (8e94d0a) to head (a29daed).
Report is 1 commits behind head on master.

Files with missing lines Patch % Lines
src/pre_process/m_perturbation.fpp 75.00% 14 Missing and 3 partials ⚠️
Additional details and impacted files
@@            Coverage Diff             @@
##           master     #879      +/-   ##
==========================================
- Coverage   45.64%   43.71%   -1.93%     
==========================================
  Files          68       68              
  Lines       18646    18450     -196     
  Branches     2249     2238      -11     
==========================================
- Hits         8511     8066     -445     
- Misses       8775     9071     +296     
+ Partials     1360     1313      -47     

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

🚀 New features to boost your workflow:
  • ❄️ Test Analytics: Detect flaky tests, report on failures, and find test suite problems.

@@ -38,7 +38,7 @@ module m_start_up

use m_helper

#ifdef MFC_MPI
#ifdef MFC_MP
Copy link
Member

Choose a reason for hiding this comment

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

bug

Copy link
Contributor Author

Choose a reason for hiding this comment

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

fixed

@@ -10,7 +10,7 @@ module m_start_up

use m_derived_types !< Definitions of the derived types

use m_global_parameters !< Global parameters for the code
use m_global_parameters !< Global parameters for the codeI
Copy link
Member

Choose a reason for hiding this comment

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

?

Copy link
Contributor Author

Choose a reason for hiding this comment

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

fixed

set(gca,'TickLabelInterpreter','latex');

% Pantano & Sarkar (2002)
load reference_data/Pantano_Sarkar_2002/production.dat;
Copy link
Member

Choose a reason for hiding this comment

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

we can't use this script because these data are unavailable. can you put the relevant data into one combined .dat (or .mat) file in this directory in version control? i assume the file will be small

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Sorry I thought I added them. I've added them now.

Copy link
Contributor Author

Choose a reason for hiding this comment

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

I can put all of them in a single file. I will update this soon

Copy link
Contributor Author

Choose a reason for hiding this comment

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

Done

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Development

Successfully merging this pull request may close these issues.

2 participants
Morty Proxy This is a proxified and sanitized view of the page, visit original site.