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

Always crash when solve a C-D nozzle flow field using rhoCentralFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 9, 2011, 07:13
Question Always crash when solve a C-D nozzle flow field using rhoCentralFoam
  #1
New Member
 
Hawk
Join Date: Jan 2011
Posts: 5
Rep Power: 6
hawklion is on a distinguished road
dear all:
I have used the rhocentralFoam in OpenFoam to solve an CD nozzle with subsonic inlet and supersonic outlet.
But there is always an error after serval iteration steps even if I have tried varied boundary conditions for the case (you can find them in the notes of the P、U、T files in the 0 folder).
Fluent mesh file and source files can be found in the attachments.
Can anyone give me some advice?
thanks!

source code"0"
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 1013250;

boundaryField
{
symmetry.5
{
type symmetryPlane;

}
wall.2
{
type zeroGradient;
}
/* pressure_inlet.3
{
type totalPressure;
p0 uniform 1.01325e+06;//_in
U U;
phi phi;
rho none;
psi none;
gamma 1.4;
value uniform 1.01325e+06;
}
*/
/*
pressure_outlet.4
{
type waveTransmissive;
value uniform 100;
field p;
gamma 1.4;
phi phi;
rho rho;
psi psi;
lInf 1;
fieldInf 100;
}*/
pressure_inlet.3
{
type fixedValue;
value uniform 1.01325e+06;
}
pressure_outlet.4
{
type zeroGradient;
}
frontAndBackPlanes
{
type empty;

}
}
// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 2000;
boundaryField
{
symmetry.5
{
type symmetryPlane;

}
wall.2
{
type zeroGradient;
}
pressure_outlet.4
{
type zeroGradient;
}
pressure_inlet.3
{
type fixedValue;
value uniform 2000;
}

frontAndBackPlanes
{
type empty;

}
}
// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (30 0 0);
boundaryField
{

symmetry.5
{
type symmetryPlane;

}
wall.2
{
type fixedValue;
value uniform (0 0 0);
}
pressure_outlet.4
{
type zeroGradient;
}
pressure_inlet.3
{
type pressureInletVelocity;
// type pressureDirectedInletOutletVelocity;
value uniform (30 0 0);
}
frontAndBackPlanes
{
type empty;

}
}
// ************************************************** *********************** //
source code "constant"

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
thermoType ePsiThermo<pureMixture<constTransport<specieThermo <eConstThermo<perfectGas>>>>>;
mixture air 1 28.9 717.5 0 1.82e-05 0.73;
//thermoType hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>;
//mixture air 1 28.9 1007 0 1.84e-05 0.7;
// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
simulationType laminar;

// ************************************************** *********************** //
source code"system"
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application rhoCentralFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 100;
deltaT 1e-4;
writeControl adjustableRunTime;
writeInterval 1e-1;
cycleWrite 0;
writeFormat ascii;
writePrecision 15;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
adjustTimeStep yes;
maxCo 0.2;
maxDeltaT 1;

// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
fluxScheme Kurganov;
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(tauMC) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
reconstruct(rho) MUSCL;
reconstruct(U) MUSCLV;
reconstruct(T) MUSCL;
}
snGradSchemes
{
default corrected;
}

// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
rho
{
solver diagonal;
}
rhoU
{
solver diagonal;
}
rhoE
{
solver diagonal;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-09;
relTol 0.01;
}
e
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-10 relTol 0;
}

}

// ************************************************** *********************** //
Attached Files
File Type: zip CDnozzle_OF.zip (72.2 KB, 25 views)
File Type: zip mesh.zip (30.5 KB, 46 views)
hawklion 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
Can 'shock waves' occur in viscous fluid flows? diaw Main CFD Forum 105 November 20, 2009 05:19
Stability Problem with sonicFoam for Nozzle Flow Julian K. OpenFOAM 2 June 29, 2009 19:46
Problem with rhoSimpleFoam matteo_gautero OpenFOAM Running, Solving & CFD 0 February 28, 2008 07:51
Nonphysical flow field while using coodles solver ankgupta8um OpenFOAM Running, Solving & CFD 5 January 26, 2008 17:54
Inviscid Drag at subsonic, subcritical Mach # Axel Rohde Main CFD Forum 1 November 19, 2001 13:19


All times are GMT -4. The time now is 08:57.