 |
 |
 |
|
|
|
[Sponsors] | ||||
February 17, 2017, 10:01
|
rhoSimpleFoam error
|
#1 |
|
Member
Jan Surwiło
Join Date: Feb 2017
Posts: 31
Rep Power: 9  |
Hello,
I'm new in openFoam. I did the icoFoam and simpleFoam simulation. But now I want to try rhoSimpleFoam with k-e turbulent model and i face problem. Mu case is simple 2,5D cylinder in channel. Velocity inlet v=2m/s on left, pressure outlet p=0Pa on right. Up and down as wall and front/back as empty. I copy the tutorial files from /tutorials/compressible/rhoSimpleFoam/squareBend/. I changed the /0 to my boundary condition and try to run simulate. Byt at the beginning of the calculation I obtain error. I pass my BC below as well as error. boundary Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
6
(
cylinder
{
type wall;
inGroups 1(wall);
nFaces 188;
startFace 33632;
}
up_wall
{
type wall;
inGroups 1(wall);
nFaces 64;
startFace 33820;
}
down_wall
{
type wall;
inGroups 1(wall);
nFaces 64;
startFace 33884;
}
inlet
{
type patch;
nFaces 31;
startFace 33948;
}
outlet
{
type patch;
nFaces 77;
startFace 33979;
}
front_back_empty
{
type empty;
// inGroups 1(wall);
nFaces 33844;
startFace 34056;
}
)
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
cylinder
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
up_wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
down_wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front_back_empty
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.0006;
boundaryField
{
cylinder
{
type kqRWallFunction;
value uniform 0;
}
up_wall
{
type kqRWallFunction;
value uniform 0;
}
down_wall
{
type kqRWallFunction;
value uniform 0;
}
inlet
{
type fixedValue;
value uniform 0.0006;
}
outlet
{
type fixedValue;
value uniform 0.0006;
}
front_back_empty
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 0.0001207;
boundaryField
{
cylinder
{
type epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0001207;
}
up_wall
{
type epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0001207;
}
down_wall
{
type epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0001207;
}
inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.02;
value uniform 0.0001207;
}
outlet
{
type inletOutlet;
inletValue uniform 0.0001207;
value uniform 0.0001207;
}
front_back_empty
{
type empty
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
cylinder
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0001207;
}
up_wall
{
type nutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0001207;
}
down_wall
{
type nutWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.0001207;
}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
front_back_empty
{
type empty
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
"(cylinder|up_wall|down_wall)"
{
type zeroGradient;
}
inlet
{
type zeroGradient;
//type mixed;
refValue uniform ;
refGradient uniform 0;
valueFraction uniform 0;
}
outlet
{
type fixedValue;
value uniform 110000;
//type mixed;
//refValue uniform 110000;
//refGradient uniform 0;
//valueFraction uniform 1;
//type transonicOutletPressure;
//U U;
//phi phi;
//gamma 1.4;
//psi psi;
//pInf uniform 110000;
}
front_back_empty
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 300;
boundaryField
{
"(cylinder|up_wall|down_wall)"
{
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform 300;
}
outlet
{
type inletOutlet;
//type zeroGradient;
value uniform 300;
inletValue uniform 300;
}
front_back_empty
{
type empty;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
"(cylinder|up_wall|down_wall)"
{
type noSlip;
}
inlet
{
type fixedValue;
value uniform (2 0 0);
}
outlet
{
type zeroGradient;
}
front_back_empty
{
type empty;
}
}
// ************************************************************************* //
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.1
Exec : rhoSimpleFoam
Date : Feb 17 2017
Time : 14:35:44
Host : "jan-VirtualBox"
PID : 11498
Case : /home/jan/open_test/cylinder_rho_simple
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0
SIMPLE: convergence criteria
field p tolerance 0.001
field U tolerance 0.0001
field e tolerance 0.001
field "(k|epsilon|omega)" tolerance 0.001
Reading thermophysical properties
Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport sutherland;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::heThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) at ??:?
#4 Foam::psiThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#5 Foam::autoPtr<Foam::psiThermo> Foam::basicThermo::New<Foam::psiThermo>(Foam::fvMesh const&, Foam::word const&) at ??:?
#6 Foam::psiThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7 ? at ??:?
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? at ??:?
Floating point exception (core dumped)
|
|
|
|
|
|
February 19, 2017, 18:31
|
|
#2 |
|
Senior Member
Agustín Villa
Join Date: Apr 2013
Location: Alcorcón
Posts: 314
Rep Power: 15 |
Hi,
I am not an expert of compressible simulations, but are you sure that your pressure field is ok? In compressible flow, I think having a 0 Pa pressure is not good... Can you show the thermophysical properties fiel |
|
|
|
|
|
|
February 20, 2017, 04:42
|
|
#3 |
|
Member
Jan Surwiło
Join Date: Feb 2017
Posts: 31
Rep Power: 9 |
Thank you for your interest in my topic. At the beginning I have to say that I'm not very good in theory of fluid flow. I used to work with Fluent and CFX and never wandering how it works deep inside. I just trust results it gives. When I start working with OpenFOAM i realized that it is necessary to understand the basics physics to calculate e.g. k and epsilon. So I have much to learn.
In Fluent and CFX it is possible to set pressure boundary condition with 0 Pa and velocity inlet for compressible (ideal gas) model. So I guess my settings are wrong. Here is the thermophysical properties: Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType
{
type hePsiThermo;
mixture pureMixture;
transport sutherland;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}
mixture
{
specie
{
nMoles 1;
molWeight 28.9;
}
thermodynamics
{
Cp 1007;
Hf 0;
}
transport
{
As 1.4792e-06;
Ts 116;
}
}
// ************************************************************************* //
I tried to find any explanation what As Ts and Hf is but I didn't find any even in FOAM guide. Is any instruction I can find all abbreviation? Thank you in advance |
|
|
|
|
|
|
October 15, 2019, 10:53
|
As, Ts Sutherland Model
|
#4 | |
|
New Member
Wan Li
Join Date: Jun 2019
Location: Edinburgh
Posts: 3
Rep Power: 7 |
Quote:
As and Ts relate to Sutherland model for viscosity calculation. Please check the following link : https://www.cfd-online.com/Wiki/Sutherland%27s_law Hf = heat of fusion : Check following link on thermodynamic models : https://www.openfoam.com/documentati...mophysical.php |
||
|
|
|
||
|
June 1, 2021, 13:23
|
|
#5 | |
|
Member
Join Date: Mar 2019
Posts: 81
Rep Power: 7 |
Quote:
I'm facing exactly the same issue with rhoSimpleFoam. Just wondering if you managed to solve the issue at the end. Please advise  Cheers, MJ |
||
|
|
|
||
 |
| Tags |
| rhosimplefoam error |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| [blockMesh] blockMesh with double grading. | spwater | OpenFOAM Meshing & Mesh Conversion | 92 | January 12, 2019 10:00 |
| DPM udf error | haghshenasfard | FLUENT | 0 | April 13, 2016 07:35 |
| [swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh | gschaider | OpenFOAM Community Contributions | 300 | October 29, 2014 19:00 |
| OpenFOAM without MPI | kokizzu | OpenFOAM Installation | 4 | May 26, 2014 10:17 |
| DecomposePar links against liblamso0 with OpenMPI | jens_klostermann | OpenFOAM Bugs | 11 | June 28, 2007 18:51 |
Linear Mode
