CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

rhoSimpleFoam

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

Like Tree2Likes
  • 1 Post By FloK
  • 1 Post By eugene

Reply
 
LinkBack Thread Tools Display Modes
Old   March 31, 2010, 04:03
Default rhoSimpleFoam
  #1
Member
 
Claudio
Join Date: Mar 2010
Posts: 57
Rep Power: 7
claco is on a distinguished road
Dear Sirs,
my name is Claudio Comis, and I am a PhD candidate in Energetic at the Department of Mechanical Engineering of University di Padova. I deal with aerodynamic shape optimization, and recently Iíve heard about OPENFOAM as CFD solver.
Bearing in mind that I am not "well experienced", anyhow, I'm starting using Openfoam (1.6 release) for my calculations.
I chose rhoSimpleFoam as solver for my calculations (external, compressible aerodynamics, i.e. fuselage aerodynamics).
I used to use Fluent, where I imposed pressureinlet and pressureoutlet as b.c.s.: what are exactly the equivalent conditions in Openfoam (that is, I have to impose total pressure and static pressure at inlet, while static pressure at outlet of my domain)?
Moreover, in Fluent I used to set an operating pressure (equal to 0): where can I set this value in Openfoam?
Finally, when I launch the solver, the following message appears, but I can't understand where is the problem:

Build : 1.6-f802ff2d6c5a
Exec : rhoSimpleFoam
Date : Mar 30 2010
Time : 14:40:33
Host : claudio-laptop
PID : 14087
Case : /home/claudio/OpenFOAM/claudio-1.6/run/tutorials/caso_mio
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kOmegaSST
#0 Foam::error:rintStack(Foam::Ostream&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#6 Foam::compressible::RASModels::kOmegaSST::F2() const in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#7 Foam::compressible::RASModels::kOmegaSST::kOmegaSS T(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#8 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<Foam::compressible::RASModels::kOmegaS ST>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#9 Foam::compressible::RASModel::New(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#10 main in "/home/claudio/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam"
#11 __libc_start_main in "/lib/libc.so.6"
#12 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception
************************************************** **
Can you please help me?
I want to make it clear first that non-zero values were assigned to both k and omega.
Moreover, I took fvSchemes and fvSolution from others tutorials, but I'm not sure they are correct for my case.
I think it is not a mesh problem, since I have already performed a calculation with potentialFoam solver.
I hope I didn't bore you.

Yours sincerely,

Claudio Comis
--
Ing. Claudio Comis Da Ronco
University of Padova
Department of Mechanical Engineering
Via Venezia, 1 - 35131 Padova, Italy
Phone 3494552408
claudio.comis@unipd.it
claudio.comis@alice.it
claco is offline   Reply With Quote

Old   March 31, 2010, 10:37
Default
  #2
Member
 
MSR CHANDRA MURTHY
Join Date: Mar 2009
Posts: 32
Rep Power: 8
chandramurthy is on a distinguished road
what is order of Mach number in the free stream region? if it is supersonic, dont use rhosimplefoam, better to use sonicFoam. fix all outlet variables with zerogradient bc (extrapolate from upstream). if u can explain your problem setup with all case files, i can suggest better.
chandramurthy is offline   Reply With Quote

Old   March 31, 2010, 11:07
Default
  #3
Member
 
Claudio
Join Date: Mar 2010
Posts: 57
Rep Power: 7
claco is on a distinguished road
Thank you,
Mach number is below 1, since this is a subsonic case.
I would like to simulate a flux of 50 m/s impinging my fuselage
I post the U,p,T, k,and omega files.
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 101325;

boundaryField
{
symmetry
{
type symmetryPlane;
}

coda
{
type zeroGradient;
}
cabina
{
type zeroGradient;
}

fusoliera
{
type zeroGradient;
}

muso
{
type zeroGradient;
}

inlet
{
type totalPressure;//zeroGradient;
p0 uniform 1531.25;
gamma 1.4;
U U;
phi phi;
psi none;//si;
rho none;//rho;
//value uniform 1531.25;

}

outlet
{
type fixedValue;
value uniform 101325;
}
}


FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (50 0 0);

boundaryField
muso
{
type fixedValue;
value uniform (0 0 0);
}
inlet
{
type outletInlet;
// phi phi;
//rho rho;
outletValue uniform (0 0 0);//(50 0 0);
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
//value uniform (0 0 0);
}
}{
symmetry
{
type symmetryPlane;
}
coda
{
type fixedValue;
value uniform (0 0 0);
}
fusoliera
{
type fixedValue;
value uniform (0 0 0);
}
cabina
{
type fixedValue;
value uniform (0 0 0);
}
muso
{
type fixedValue;
value uniform (0 0 0);
}
inlet
{
type outletInlet;
// phi phi;
//rho rho;
outletValue uniform (0 0 0);
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
//value uniform (0 0 0);
}
}

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 1 0 0 0];

internalField uniform 289;

boundaryField
{
symmetry
{
type symmetryPlane;
}

inlet
{
//type fixedValue;
type totalTemperature;
T0 uniform 289;
U U;
phi phi;
psi psi;
gamma 1.4;
value uniform 289;
}

outlet
{
type inletOutletTotalTemperature;
T0 uniform 289;
U U;
phi phi;
psi psi;
gamma 1.4;
//type fixedValue;
//value uniform 289;
}
muso
{
type zeroGradient;//inletOutlet;
//inletValue uniform 289;
//value uniform 289;
}

cabina
{
type zeroGradient;//inletOutlet;
//inletValue uniform 289;
// value uniform 289;

}
fusoliera
{
type zeroGradient;//inletOutlet;
// inletValue uniform 289;
//value uniform 289;
}
coda
{
type zeroGradient;//inletOutlet;
// inletValue uniform 289;
// value uniform 289;
}
}


FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.37;

boundaryField
{
symmetry
{
type symmetryPlane;
}
inlet
{
type turbulentIntensityKineticEnergyInlet;
intensity 0.01;
value uniform 0.37;
// type fixedValue;
// value uniform 0.37;
}
outlet
{
type inletOutlet;
inletValue uniform 0.37;
value uniform 0.37;
}
coda
{
type compressible::kqRWallFunction;
value uniform 1e-20;
}
fusoliera
{
type compressible::kqRWallFunction;
value uniform 1e-20;
}
cabina
{
type compressible::kqRWallFunction;
value uniform 1e-20;
}
muso
{
type compressible::kqRWallFunction;
value uniform 1e-20;
}
}


FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 -1 0 0 0 0];

internalField uniform 1.11;

boundaryField
{
symmetry
{
type symmetryPlane;
}
inlet
{
type fixedValue;
value uniform 1.11;
}
outlet
{
// type zeroGradient;
type inletOutlet;
inletValue uniform 1.11;
value uniform 1.11;
}
coda
{
type compressible:megaWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 1.11;
}
fusoliera
{
type compressible:megaWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 1.11;
}
cabina
{
type compressible:megaWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 1.11;
}
muso
{
type compressible:megaWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 1.11;
}
}


I'm not sure that the 0 files are correct. I have taken the values of p, U, k, T and the omega from the Fluent case, just to be sure. (the operating pressure in Fluent is 1012325, but I don't know where to set this value in OF)
Finally, the question remains:
I chose rhoSimpleFoam as solver for my calculations (external, compressible aerodynamics, i.e. fuselage aerodynamics).
I used to use Fluent, where I imposed pressureinlet (po=1531 [Pa] and an initial gauge pressure of 0 [Pa]) and pressureoutlet (p=0 [Pa]
) as b.c.s.: what are exactly the equivalent conditions in Openfoam (that is, I have to impose total pressure and static pressure at inlet, while static pressure at outlet of my domain)?
Thank You in Advance.
Claudio
claco is offline   Reply With Quote

Old   April 6, 2010, 12:49
Default
  #4
Member
 
Claudio
Join Date: Mar 2010
Posts: 57
Rep Power: 7
claco is on a distinguished road
Dear Sirs,

now I have slightly changed my boundary conditions, but I always get this message:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.6-f802ff2d6c5a
Exec : rhoSimpleFoam
Date : Apr 06 2010
Time : 18:34:31
Host : claudio-laptop
PID : 10703
Case : /home/claudio/OpenFOAM/claudio-1.6/run/tutorials/caso_mio
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
Prt 1;
}


Starting time loop

Time = 1



incompatible dimensions for operation
[U[0 1 -2 0 0 0 0] ] + [U[1 -2 -2 0 0 0 0] ]#0 Foam::error:rintStack(Foam::Ostream&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/claudio/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 void Foam::checkMethod<Foam::Vector<double> >(Foam::fvMatrix<Foam::Vector<double> > const&, Foam::fvMatrix<Foam::Vector<double> > const&, char const*) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam"
#3 Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam:perator+<Foam::Vector<double> >(Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > const&, Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > const&) in "/home/claudio/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam"
#4 main in "/home/claudio/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhoSimpleFoam"
#5 __libc_start_main in "/lib/libc.so.6"
#6 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116


From function checkMethod(const fvMatrix<Type>&, const fvMatrix<Type>&)
in file /home/dm2/henry/OpenFOAM/OpenFOAM-1.6/src/finiteVolume/lnInclude/fvMatrix.C at line 1195.

FOAM aborting

Aborted



I cannot find what is wrong with my settings.
Moreover I would like to know what is the meaning of rho,phi, U and psi in the totalPressure boundary condition type.
Thank You in advance.
Claudio
claco is offline   Reply With Quote

Old   April 8, 2010, 08:30
Default
  #5
New Member
 
Join Date: Dec 2009
Posts: 3
Rep Power: 7
FloK is on a distinguished road
Hi,
the "incompatible dimensions" error occurs at my simulations when i copy the p file of an incompressible case to a compressible one.
In incompressible cases, the dimensions of pressure are m^2/s^2, its "density normalized pressure".
In compressible cases the dimensions are kg/(m*s^2) = Pa.
Check your dimensions, if there is anything right.

Regards,
FloK
mm.abdollahzadeh likes this.
FloK is offline   Reply With Quote

Old   April 15, 2010, 10:39
Default
  #6
Member
 
Claudio
Join Date: Mar 2010
Posts: 57
Rep Power: 7
claco is on a distinguished road
Hi all,

I checked the dimensions of p, and they are correct ([1 -1 -2 0 0 0 0]).
However, an incompatible dimensions error occurs.
I check all the other dimensions, and they are correct.
claco is offline   Reply With Quote

Old   April 16, 2010, 02:44
Default
  #7
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Try to have your case running using the rhoPimpleFoam tutorial as reference (you'll probably benefit of pimpleFoam capability of performing iterations in each time step anyway).

Best,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   April 20, 2010, 04:32
Default
  #8
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 12
eugene is on a distinguished road
Hi,

There are two fields for which the dimensions change between compressible and incompressible: p and phi. It looks to me as if you are using the phi field from the incompressible run. Simply delete phi from your starting directory (it will be recreated automatically and correctly) and try again.
namsivag likes this.
eugene is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pressure instability with rhoSimpleFoam daniel_mills OpenFOAM Running, Solving & CFD 44 February 17, 2011 18:08
Problem with rhoSimpleFoam mecbe2002 OpenFOAM 3 April 11, 2010 00:54
TRANSONIC FLOW in RHOSIMPLEFOAM dinonettis OpenFOAM Running, Solving & CFD 9 October 21, 2009 12:02
RhoSimpleFoam FoamX spv24 OpenFOAM Running, Solving & CFD 1 July 21, 2008 10:29
Stability startup problems with rhoSimpleFoam olesen OpenFOAM Running, Solving & CFD 1 July 18, 2006 08:09


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