noSlip BC for rhoCentralFoam

tum-student · July 7, 2022, 18:40

Hello,
I am simulating a diamond airfoil at supersonic flow at Mach 3. I am using rhoCentralFoam as a solver and k-w SST turbulence model to solve this case. When I use velocity boundary condition as "slip" for the airfoil boundary, my run is successful without any problems, but when I change the BC to "noSlip". I am getting a "Negative Initial Temperature T0" error and the run crashes.
I have also tried to limit the temperature using fvOptions, but it's not working. I have also checked by refining the mesh, but still, the problem exists.

Please help me.

Other boundary conditions for airfoil are as follows:
p - type zeroGradient
T - type zeroGradient
alphat - type compressible::alphatWallFunction, value uniform 8.88e-9
epsilon - type epsilonWallFunction, value uniform 365627.58
k - type kqRWallFunction, value uniform 901.31
nut - type nutkWallFunction, value uniform 0.199
omega - type omegaWallFunction, value uniform 4507.4

chc · January 23, 2024, 20:50

Hello,
Have you had any luck in solving this issue? I keep getting those negative Temperature errors, however, after a few tenths of a second in the simulation. Thank you.

July 7, 2022, 18:40	noSlip BC for rhoCentralFoam	#1
tum-student New Member Gurjot Ahuja Join Date: Jun 2022 Posts: 5 Rep Power: 4	Hello, I am simulating a diamond airfoil at supersonic flow at Mach 3. I am using rhoCentralFoam as a solver and k-w SST turbulence model to solve this case. When I use velocity boundary condition as "slip" for the airfoil boundary, my run is successful without any problems, but when I change the BC to "noSlip". I am getting a "Negative Initial Temperature T0" error and the run crashes. I have also tried to limit the temperature using fvOptions, but it's not working. I have also checked by refining the mesh, but still, the problem exists. Please help me. Other boundary conditions for airfoil are as follows: p - type zeroGradient T - type zeroGradient alphat - type compressible::alphatWallFunction, value uniform 8.88e-9 epsilon - type epsilonWallFunction, value uniform 365627.58 k - type kqRWallFunction, value uniform 901.31 nut - type nutkWallFunction, value uniform 0.199 omega - type omegaWallFunction, value uniform 4507.4 Last edited by tum-student; July 10, 2022 at 10:53.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Openfoam Turbulent Car Simulation	meganrw	OpenFOAM Running, Solving & CFD	0	July 5, 2021 16:50
Wall Conditions : fixedValue or noSlip	Magistrane	OpenFOAM Running, Solving & CFD	4	September 28, 2020 07:55
ACMI rotating zone: which frame is used for fixedValue, noSlip, movingWallVelocity?	zommuter	OpenFOAM Running, Solving & CFD	0	September 2, 2020 04:13
Problem with boundary layer - noSlip condition not being respected - OpenFOAM	heitorvitorc	OpenFOAM Post-Processing	1	August 14, 2019 06:14
multiphaseEulerFoam/tank discharge	ves	OpenFOAM Running, Solving & CFD	8	December 31, 2017 15:59

January 23, 2024, 20:50		#2
chc New Member Join Date: Nov 2023 Posts: 18 Rep Power: 3	Hello, Have you had any luck in solving this issue? I keep getting those negative Temperature errors, however, after a few tenths of a second in the simulation. Thank you.