search

Exawind / nalu-wind

Solver for wind farm simulations targeting exascale computational platforms
https://nalu-wind.readthedocs.io
Other
122 stars 83 forks source link

Default value of SDRWallFactor: 1 -> 10 #1254

Closed BumseokLee closed 1 month ago

BumseokLee commented 5 months ago

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 commented 5 months ago

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 commented 5 months ago

What value are the NASA codes using?

BumseokLee commented 5 months ago

@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 .