Solver blow up due to negative psi values on several cells -- CFD Online Discussion Forums

February 26, 2020, 11:28

Solver blow up due to negative psi values on several cells

Soll

New Member

Zeng

Join Date: Jul 2017

Location: London, UK

Posts: 16

Rep Power: 9

I'm currently using rhoCentralFoam to simulate a nozzle flow with body-force.
The solver blows up after several iterations. I did some research and found the problem occurs when the solver excutes

Code:

volScalarField c("c", sqrt(thermo.Cp()/thermo.Cv()*rPsi));

I write out the related field. It turns out that the rPsi field (1/Psi) has negative values (see attached fig), but the other part of the field is rather continuous.
Anyone has idea on this problem and how to debug? Any suggestions are appreciated.

Attached Images

rPsi_field.png (23.7 KB, 10 views)

« Previous Thread | Next Thread »

Posting Rules
You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is Off Trackbacks are Off Pingbacks are On Refbacks are On Forum Rules

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Other] mesh airfoil NACA0012	anand_30	OpenFOAM Meshing & Mesh Conversion	13	March 7, 2022 18:22
"unknown psi type" error in reactingFoam	Speno93	OpenFOAM Running, Solving & CFD	2	May 25, 2021 10:15
New solver errors, wmake	Sedullo	OpenFOAM Programming & Development	4	January 24, 2020 07:43
[snappyHexMesh] sHM layer process keeps getting killed	MBttR	OpenFOAM Meshing & Mesh Conversion	4	August 15, 2016 04:21
[blockMesh] Errors during blockMesh meshing	Madeleine P. Vincent	OpenFOAM Meshing & Mesh Conversion	51	May 30, 2016 11:51

All times are GMT -4. The time now is 03:19.