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

Default value of SDRWallFactor: 1 -> 10 #1254

Merged
merged 1 commit into from
Jul 30, 2024
Merged

Conversation

BumseokLee
Copy link
Contributor

The SDRWallFactor is a constant for the wall boundary condition of the omega in the SST turbulence model.

1

However, there are some uncertainties about this in Nalu-Wind. The original SST model by Menter uses 10 for the SDRWallFactor. Currently, in Nalu-Wind, users can specify the SDRWallFactor in the input yaml file, and the default value is 1 if not explicitly specified. Many users run RANS simulations with the default value. Some reg test cases in the repository used 0.625. Marc orignally performed the validation of the NASA’s 2D Zero Pressure Gradient (ZPG) Flat Plate Verification Case with SDRWallFactor of 10.

The wall boundary condition plays a crucial role in transition simulations because the transition model triggers laminar-turbulent transition based on flow quantities close to the middle of the boundary layer, which is very near the wall.

For this reason, I tested different values of SDRWallFactor for the NASA TMR’s 2D ZPG flat plate and NACA 0012 airfoil cases and adjusted the default value to 10 based on the test results, which I believe is the correct value.

I've attached some results for the ZPG flat pate cases below

2 3 4

Please refer to the attached file for more detailed analysis. Any discussion about the is welcomed!

SDR_Wall_factor.pdf

@marchdf
Copy link
Contributor

marchdf commented Apr 5, 2024

Fine by me. I think the only reason we didn't do this way back in the day is that it would change everyone's cases/runs and we seemed averse to doing that. Seems like going to 10 helps drags but hurts lift (based on comparing to NASA codes).

@psakievich
Copy link
Contributor

What value are the NASA codes using?

@BumseokLee
Copy link
Contributor Author

@marchdf Did you mean the lift coefficient of NACA0012 case that I included in the file? The lift coefficients from 0.625 and 10 are almost identical. There are only minor differences in the drag coefficient at AoA=15deg, although the trends of the omega residual are significantly different.

@psakievich NASA's codes use 10. I've attached the SST model page from NASA Langley's Turbulence Modeling Resource Website here: https://turbmodels.larc.nasa.gov/sst.html .

@marchdf marchdf enabled auto-merge (rebase) July 30, 2024 21:50
@marchdf marchdf merged commit abe69db into Exawind:master Jul 30, 2024
3 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.

4 participants