CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

rhoPimpleFoam fails on an Euler vortex

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 29, 2019, 10:08
Post rhoPimpleFoam fails on an Euler vortex
  #1
New Member
 
Thomas
Join Date: Jul 2018
Posts: 10
Rep Power: 4
Roxxor91 is on a distinguished road
Hi all,
I would like to reproduce the Euler vortex in that paper with the same mesh motion and eventually incorporate my own mesh motion.

My problem has two steps:
  1. Simulate the Euler vortex
  2. Incorporate a mesh motion

I followed the description of the vortex problem (Initial values for pressure, velocity and density, "periodic" boundary conditions, same geometry) but it keeps failing to advect the vortex in a correct manner. Either it crashes because of a temperature that goes below 0 or the simulation goes on with crazy values and has obviously no physical sense.

I attached the case where simulation go on but fields are non sense (fvOption does not seem to have an effect though). To see how I run it, check the script. I gave the 0 folder because you would have to compile and run makeInitialConditions binary otherwise.

My opinion is that the problem comes from BCs. I set them as CYCLIC because it is what looked the most liek periodic (prescribed in the paper). I run OF6 on Linux Mint 18.2 x64.

How can I make rhoPimpleFoam run my problem ?



EDIT: This is my checkMesh

Quote:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

--> FOAM Warning :
From function static Foam::instantList Foam::timeSelector::select0(Foam::Time&, const Foam::argList&)
in file db/Time/timeSelector.C at line 274
No time specified or available, selecting 'constant'
Create polyMesh for time = constant

Time = constant

Mesh stats
points: 2676
internal points: 0
faces: 17886
internal faces: 12538
cells: 7606
faces per cell: 4
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 1

Overall number of cells of each type:
hexahedra: 0
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 7606
polyhedra: 0

Checking topology...
Boundary definition OK.
***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
emptyPatches_empt 4800 2542 ok (non-closed singly connected)
top_cyc 158 121 ok (non-closed singly connected)
bottom_cyc 158 121 ok (non-closed singly connected)
inlet_cyc 116 90 ok (non-closed singly connected)
outlet_cyc 116 90 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-10 -7.5 0) (10 7.5 0.875)
Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
Mesh has 2 solution (non-empty) directions (1 1 0)
***Number of edges not aligned with or perpendicular to non-empty directions: 4346
<<Writing 2492 points on non-aligned edges to set nonAlignedEdges
Boundary openness (-4.02954e-18 1.07454e-17 -1.25923e-18) OK.
Max cell openness = 7.40149e-17 OK.
Max aspect ratio = 3.08972 OK.
Minimum face area = 0.0947931. Maximum face area = 0.333659. Face area magnitudes OK.
Min volume = 0.0130641. Max volume = 0.0537569. Total volume = 262.5. Cell volumes OK.
Mesh non-orthogonality Max: 52.7007 average: 33.1246
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.591481 OK.
**Error in coupled point location: 202 faces have their 0th or consecutive vertex not opposite their coupled equivalent. Average mismatch 0.163664.
<<Writing 202 faces with incorrectly matched 0th (or consecutive) vertex to set coupledFaces

Failed 2 mesh checks.

End


EDIT2: I ran rhoCentralFoam on the same test case, just add to specify div(tauMC) "Gauss linear;" in fvSchemes, and got the famous "printStack" error:

Quote:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.00457721, Final residual = 2.35364e-07, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00456667, Final residual = 2.31317e-07, No Iterations 1
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 0.00407362, Final residual = 3.6188e-07, No Iterations 1
ExecutionTime = 7.49 s ClockTime = 7 s

Mean and max Courant Numbers = 0.215573 0.414158
deltaT = 0.0001
Time = 0.0106

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.0045694, Final residual = 1.93503e-07, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00453489, Final residual = 1.9462e-07, No Iterations 1
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 0.00402989, Final residual = 3.32675e-07, No Iterations 1
ExecutionTime = 7.55 s ClockTime = 7 s

Mean and max Courant Numbers = 0.215788 0.426271
deltaT = 0.0001
Time = 0.0107

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.00457203, Final residual = 1.74982e-07, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00452333, Final residual = 1.78578e-07, No Iterations 1
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 0.0040209, Final residual = 3.12121e-07, No Iterations 1
ExecutionTime = 7.61 s ClockTime = 8 s

Mean and max Courant Numbers = 0.216055 0.437246
deltaT = 0.0001
Time = 0.0108

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.00457214, Final residual = 1.85512e-07, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00453382, Final residual = 1.8473e-07, No Iterations 1
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 0.004061, Final residual = 3.0108e-07, No Iterations 1
ExecutionTime = 7.7 s ClockTime = 8 s

Mean and max Courant Numbers = 0.216374 0.447017
deltaT = 0.0001
Time = 0.0109

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 0.00459904, Final residual = 2.19477e-07, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00456113, Final residual = 2.02625e-07, No Iterations 1
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 ? in "/home/thomas/OpenFOAM/OpenFOAM-6/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? in "/home/thomas/OpenFOAM/OpenFOAM-6/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam"
Floating point exception (core dumped)
Attached Files
File Type: zip eulerVortex.zip (143.4 KB, 1 views)
Roxxor91 is offline   Reply With Quote

Old   September 29, 2019, 11:27
Default
  #2
New Member
 
Thomas
Join Date: Jul 2018
Posts: 10
Rep Power: 4
Roxxor91 is on a distinguished road
EDIT3: I have been able to reproduce a negative initial temperature error: if I switch from sutherland to hConst in thermodynamicProperties, I an error quickly which I think might also be related to the coefficients ( I chose mu=1e-3 and Pr = 0.7):
Quote:


--> FOAM FATAL ERROR:
Negative initial temperature T0: -39.0873

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam:erfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
in file /home/thomas/OpenFOAM/OpenFOAM-6/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 54.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#4 ? in "/home/thomas/OpenFOAM/OpenFOAM-6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? in "/home/thomas/OpenFOAM/OpenFOAM-6/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
Aborted (core dumped)
Roxxor91 is offline   Reply With Quote

Old   September 29, 2019, 14:18
Default Changing to rhoCentralFoam
  #3
New Member
 
Thomas
Join Date: Jul 2018
Posts: 10
Rep Power: 4
Roxxor91 is on a distinguished road
I figured out rhoPimpleFoam was originally suited for transcient simulation and I decided to switch to rhoCentralFoam.
I have slightly better results (still not correct). I have attached the new case and 3 snapshots of the pressue fields that dissipates very quickly. I guess it is due to the boundary conditions.

What I care about in the simulation is the advection of the initial vortex, so if there is a way to do it without periodic BCs, I'll take it even if it contradicts the paper..

Anyway I am still keen to receive suggestions and remarks
Attached Images
File Type: png 0.png (15.8 KB, 8 views)
File Type: png 3.png (29.3 KB, 6 views)
File Type: png 9.png (26.0 KB, 7 views)
Attached Files
File Type: zip eulerVortex.zip (163.7 KB, 1 views)
Roxxor91 is offline   Reply With Quote

Reply

Tags
ale, dynamicmeshdict, openfoam 1806, rhocentralfoam, rhopimplefoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
2D simulation of vortex shedder with rhoPimpleFoam linkamp OpenFOAM Running, Solving & CFD 0 October 31, 2016 11:00
simulation of irrotational vortex fails matthias OpenFOAM 3 October 21, 2013 05:49
Presenation of vortex strength and velocity concept? fruitkiwi Main CFD Forum 0 September 26, 2012 23:08
vortex cause pressure gradient or pressure gradient induce vortex? fruitkiwi Main CFD Forum 4 June 12, 2012 02:12
rhoPimpleFoam and sort of vortex tube pcaste12 OpenFOAM 0 July 15, 2011 10:10


All times are GMT -4. The time now is 11:39.