Running a high Re low mach rhoPimpleFoam

CaptCrash · November 15, 2022, 13:56

Hello everyone!
I am a long time reader, first time poster. I am trying to run a rhoPimpleFoam for a steady state high reynolds number low mach airfoil but running into a supersonic flow solution. Here are the pertinent simulation parameters
Chord (m): 10.85838
Re: 95,624,744
Angle of Attack (deg): 5

Uxyz (m/s): [128.1220, 11.20922, 0]
Uinf (m/s): 128.611
Minf (Mach): 0.378
Tinf (K): 288.15
Pinf (Pa): 101325
rho (kg/m^3): 1.225
Pr: 0.7323
Hf (j/kg): 2.8833272E5

Y+ target: 30
Turbulence Properties: RAS, kOmegaSST, I = 0.01%

k: 2.48E-4
nut: 1.2E-2
omega: 2.07E-2

I've attached as many of the files as possible to this string to help. I used gmsh for the mesh generation and those files are included as well to generate a full mesh which, after generating and saving the mesh in gmsh, can be converted using gmshToFoam <file_name.msh>. I considered running a normalized chord (1 meter) in water at 107.404m/s. The issue being this will be part of a series of test to generate CL/CD/CM plots up to the trans-sonic regime so I would like to capture those effects. Thank you very much in advance for your assistance.

CaptCrash · November 16, 2022, 18:07

The following information is what I calculated for the inflation layers. The grid is a C-H grid. I've included an image of the C-H grid, the computations for the inflation layers, and the checkMesh output. The C-H grid has has the following properties.

Screenshot from 2022-11-16 14-56-18.jpg
BAC310 AIR.jpg
checkmeshreport.txt

CaptCrash · November 19, 2022, 13:45

I think one issue stems from the high aspect ratio mesh cells. I am using GMSH but I can not figure out how to generate an inflation layer and maintain a structured CH mesh. Every topic or example I've seen has been an inflation layer connected to an unstructured mesh. I used this video (slowed down to 25% playback):https://www.youtube.com/watch?v=bekRbU7rtZE as the basis for generating the mesh, except for the wake spacing. I tried using snappyHexMesh earlier in the process but the blender plugin doesnt seem to work on linux for some reason.

Thanks in advance for any help.

November 16, 2022, 18:07	Follow up information	#2
CaptCrash New Member Join Date: Nov 2022 Posts: 7 Rep Power: 3	The following information is what I calculated for the inflation layers. The grid is a C-H grid. I've included an image of the C-H grid, the computations for the inflation layers, and the checkMesh output. The C-H grid has has the following properties. Screenshot from 2022-11-16 14-56-18.jpg BAC310 AIR.jpg checkmeshreport.txt

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
high Temperature shoot up and high mach flow.	skdubey	ANSYS	1	February 19, 2022 06:35
demarcation between high Re number and Low Re number	pdewcfd	Main CFD Forum	0	March 31, 2021 03:24
High Mach number flow	Taendyr	Main CFD Forum	8	November 29, 2020 07:15
Multiphase, high Mach simulations	alexvaleije	FLUENT	0	June 6, 2019 06:53
replacing of shock tube high pressure part with a boundary condition for low pressure	immortality	Main CFD Forum	0	May 2, 2013 13:30

November 19, 2022, 13:45		#3
CaptCrash New Member Join Date: Nov 2022 Posts: 7 Rep Power: 3	I think one issue stems from the high aspect ratio mesh cells. I am using GMSH but I can not figure out how to generate an inflation layer and maintain a structured CH mesh. Every topic or example I've seen has been an inflation layer connected to an unstructured mesh. I used this video (slowed down to 25% playback):https://www.youtube.com/watch?v=bekRbU7rtZE as the basis for generating the mesh, except for the wake spacing. I tried using snappyHexMesh earlier in the process but the blender plugin doesnt seem to work on linux for some reason. Thanks in advance for any help.