/ ( )(H17/9 )
ADVENTURECluster http://www.alde.co.jp
2 OpenFOAM for beginner
LES ( H17/3)3D (in-house)
CFDOpenFOAM/FrontFlow/UPACS/SPHYSICS/FreeCFD
CFD(intel Core i7 920 , 12GB, geforce GTS250)(intel Core i7 920 , 12GB, geforce GTS250)
(H22/2)( )
OpenFOAM( )
2 OpenFOAM for beginner
( )1.1 reactingFoam1.2 dieselFoam1.3 XiFoam
example
OpenFOAM extension or development versionOpenFOAM extension or development version
OpenFOAM
※ OpenFOAM tutorial
2 OpenFOAM for beginner
/①XiFoam(H.Weller b-Ξ )②engineFoam( XiFoam)②engineFoam( XiFoam)
③reactingFoam(Golovitchev@chalmers PaSR )④rhoReacingFoam(reactingFoam density based solver version)
dieselFoam( PaSR )dieselEngineFoam( )fireFoam( )
coalChemistryFoam( )coalChemistryFoam( )porousExplicitSourceReactingParcelFoam( )
EBU/EDM SimpleFoam ?Ref : AALBORG C. Andersen , N. E. L. Nielsen
Numerical investigation of a BFR using OpenFOAM
2 OpenFOAM for beginner
1.1 reactingFoamchalmers Golovitchev PaSRchalmers Golovitchev PaSR
→eddy disspation / eddy break up
0
U
constant system
CH4( )
U
T
p
Ydefault
k( )
polyMesh
thermophysicalProperties
turbulentProperties(RAS/LES)
chemistryproperties
reactions
thermo.compressibleGas
controlDictpolyMeshpolyMesh
fvSolution
fvSchemes
2 OpenFOAM for beginner
1.1 reactingFoam/ constant/ constant
①thermophysicalPropertiesFoamFile{}// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //thermoType hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;
inertSpecie N2;inertSpecie N2;
chemistryReader foamChemistryReader;
foamChemistryFile "$FOAM_CASE/constant/reactions";
foamChemistryThermoFile "$FOAM_CASE/constant/thermo.compressibleGas";
2 OpenFOAM for beginner
1.1 reactingFoam/ constant
★reactingFoam thermoType
thermoType hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;or
thermoType hsPsiMixtureThermo<multiComponentMixture<gasThermoPhysics>>;
★ → N2InertSpecie N2 ;
★ / → Of or chemkinchemistryReader foamChemistryReader ; // OF
orchemistryReader chemkinReader ; // chemkin
2 OpenFOAM for beginner
1.1 reactingFoam/ constant
★ /foamChemistryFile "$FOAM_CASE/constant/reactions";
orCHEMKINFile "$FOAM_CASE/constant/chem.inp";
ELEMENTSH O C NENDSPECIE C7H16 O2 N2 CO2 H2O
species(
O2H2OCH4CO2N2 C7H16 O2 N2 CO2 H2O
ENDREACTIONSC7H16 + 11O2 => 7CO2 + 8H2O 5.00E+8 0.0 15780.0! 1
FORD / C7H16 0.25 /FORD / O2 1.5 /
END
N2);reactions(
irreversibleArrheniusReactionCH4 + 2O2 = CO2 + 2H2O(5.2e16 0 14906)
);
2 OpenFOAM for beginner
1.1 reactingFoam/ constant
★ foamChemistryThermo "thermo.compressibleGas";or
CHEMKINThermoFile "therm.dat";(O2 O2 1 31.9988
200 5000 10003.69758 0.00061352 -1.25884e-07 1.77528e-11 -1.13644e-15 -1233.93 3.189173.21294 0.00112749 -5.75615e-07 1.31388e-09 -8.76855e-13 -1005.25 6.034741.67212e-06 170.672
H2O H2O 1 18.0153cp = R((((a4T + a3)T + a2)T + a1)T + a0R
≒1000Kjanaf a0 a4 a5 a6janaf a0 a4 a5 a6
) H = RT(((((a4/5T + a3/4)T + a2/3)T + a1/2)T+a0)+a5R
S = R((((a4/4T + a3/3)T + a2/2)T + a1)T + a0RlnT+a6R
2 OpenFOAM for beginner
1.1 reactingFoam/ constant/ constant
②chemistryPropertiespsiChemistryModel ODEChemistryModel<gasThermoPhysics>;
chemistry on;
chemistrySolver ode;
initialChemicalTimeStep 1e-07;initialChemicalTimeStep 1e-07;
turbulentReaction on;
Cmix Cmix [ 0 0 0 0 0 0 0 ] 0.1;
2 OpenFOAM for beginner
1.1 reactingFoam/ constant/ constant
★ reactingFoam 2psiChemistryModel ODEChemistryModel<gasThermoPhysics>;
orpsiChemistryModel ODEChemistryModel<icoPoly8ThermoPhysics>;
★chemistry on;chemistry on;
turbulentReaction on;
Cmix Cmix [ 0 0 0 0 0 0 0 ] 0.1;
2 OpenFOAM for beginner
1.1 reactingFoam/ constant
★initialChemicalTimeStep 1e-07;
★Cmix Cmix [ 0 0 0 0 0 0 0 ] 0.1;
※PaSRtk = Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon());tk = Cmix*sqrt(turbulence->muEff()/rho/turbulence->epsilon());
Turbulentreaction onkappa = (runTime.deltaT() + tc)/(runTime.deltaT() + tc + tk);Turbulentreaction offkappa = 1.0;
chemistrySh = kappa*chemistry.Sh()(); →
2 OpenFOAM for beginner
1.1 reactingFoam/ constant
★Sequential or EulerImplicit or ode
chemistrySolver ode;
sequentialCoeffs{
cTauChem 0.001;
EulerImplicitCoeffs{
cTauChem 0.05;equilibriumRateLimiter off;
}
tc
cTauChem 0.001;}
odeCoeffs{
ODESolver SIBS;eps 0.05;scale 1;
}
}
-SIBS : Semi-Implicit Bulirsh-Stoer
KRR4 : Kaps-Rentrop
RK : Runge-Kutta
2 OpenFOAM for beginner
1.2 dieselFoamreactingFoam ( )reactingFoam ( )
(ref M06005 )
0
U
T
constantsystem
controlDictthermophysicalProperties
CH4( )
T
p
Ydefault
k( )
polyMesh
combustionProperties
turbulentProperties(RAS/LES)
chemistryproperties
sprayProperties
injectorProperties
polyMesh fvSolution
fvSchemesft / fu
chemkin
chem.inp
therm.dat
chemkin
2 OpenFOAM for beginner
1.2 dieselFoam/ constant
①thermophysicalProperties
thermoType hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>;
CHEMKINFile "$FOAM_CASE/chemkin/chem.inp";
CHEMKINThermoFile "~OpenFOAM/thermoData/therm.dat";
reactingFoam
inertSpecie N2;
liquidComponents ( C7H16 ); ←liquidProperties{
C7H16 C7H16 defaultCoeffs; ← ( )}
2 OpenFOAM for beginner
1.2 dieselFoam/ constant
①thermophysicalProperties
★ ( CHEMKIN File )
LiquidComponents(C7H16C6H14C6H14);
liquidProperties{
C7H16 C7H16 defaultCoeffs;C6H14 C6H14 defaultCoeffs;
}
2 OpenFOAM for beginner
1.2 dieselFoam/ constant/ constant
②combustionProperties
Cmix Cmix [ 0 0 0 0 0 0 0 ] 1; ← PaSR ( )
ignitionProperties1{
ignite off; ← on/off
ignitionPoint ignitionPoint [ 0 1 0 0 0 0 0 ] ( 0.2 0 0.02 ); ←
timing timing [ 0 0 1 0 0 0 0 ] 0; ←
duration duration [ 0 0 1 0 0 0 0 ] 1; ←}
2 OpenFOAM for beginner
1.2 dieselFoam/ constant/ constant
③injectorProperties
({
injectorType unitInjector; ←
unitInjectorProps{
commonRailInjectormultiHoleInjectorswirlInjector
{Position (0 0.0995 0); ← injector ( 0.5mm)Direction (0 -1 0); ← (-y )Diameter 0.00019; ← nozzleCd 0.9; ← mass 6e-06; ← NParcels 5000; ←
2 OpenFOAM for beginner
1.2 dieselFoam/ constant
③injectorProperties 2X(
1.0 ← ( 1.0));massFlowRateProfile ← mass flow(
(0 0.1272)
liquidProperties
(0 0.1272)(4.16667e-05 6.1634)
temperatureProfile ←(
(0.0 320.0)(0.00125 320.0)
);
2 OpenFOAM for beginner
1.2 dieselFoam/ constant
④sprayProperties 1
InterpolationSchemes {
U cellPointFace;rho cell;p cell;T cell;
}}
subCycles 2; ← minimum number of Lagrangian sub cycles
atomizationModel off; ← How atomization is treated( LISA/blobsSheetAtomization)
includeOscillation yes; ← droplet deformation; will effect droplet drag coef
2 OpenFOAM for beginner
1.2 dieselFoam/ constant
④sprayProperties 2
breakupModel ReitzKHRT; ← If secondary break up is used( reitzDiwakar/ETAB/TAB/SHF/noBreakup)
injectorModel hollowConeInjector; ← Which injector model to use( blobsSwirl/Chomiak/pressureSwirl/constant)
collisionModel off; ← particle – particle interactioncollisionModel off; ← particle – particle interaction( ORourke/trajectoryModel)
evaporationModel standardEvaporationModel; ← which evaporation model to use ( RutlandFlashBoil/saturateEvaporationModel)
2 OpenFOAM for beginner
1.2 dieselFoam/ constant
④sprayProperties 3
heatTransferModel RanzMarshall; ← Particle heat transfer model(or no heat transfer model)
dispersionModel off; ← If turbulent dispersion is used or not
dragModel standardDragModel; ← Particle drag modeldragModel standardDragModel; ← Particle drag model(or no drag model)
wallModel reflect; ← What happens to particles hitting the walls(or remove )
2 OpenFOAM for beginner
1.3 XiFoamImperialCollege Gosman H.Weller( OpenCFD)
Ref 24th Symposium (international) on Combution 1998/pp899-907Ref 24th Symposium (international) on Combution 1998/pp899-907
flame surface density model
→
0
U
constant system
muSgs / k ( )
U
T / Tu
p
b / Xi
Su
polyMesh
thermophysicalProperties
turbulentProperties(RAS/LES)
combustionProperties
controlDictpolyMeshpolyMesh
fvSolution
fvSchemes
2 OpenFOAM for beginner
1.3 XiFoam/ constant
①thermophysicalProperties ThermoType (1)< (2) < (3) < (4) < (5) < (6) > > > > >
hhuMixtureThermo<homogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>;
( )hhuMixtureThermo<inhomogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>;anafThermo<perfectGas>>>>>;
hhuMixtureThermo<veryInhomogeneousMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>;
hhuMixtureThermo<egrMixture<sutherlandTransport<specieThermo<janafThermo<perfectGas>>>>>;
janaf +
2 OpenFOAM for beginner
1.3 XiFoam/ constant/ constant
①thermophysicalProperties 2ThermoType (1)< (2) < (3) < (4) < (5) < (6) > > > > >
( )hhuMixtureThermo<inhomogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>;
hhuMixtureThermo<veryInhomogeneousMixture<constTransport<specieThermo<hhhuMixtureThermo<veryInhomogeneousMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>;
hhuMixtureThermo<egrMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>;
+
2 OpenFOAM for beginner
1.3 XiFoam/ constant/ constant
②combustionProperties 1
laminarFlameSpeedCorrelation constant; ← Gulders
fuel Propane; ←
Su Su [0 1 -1 0 0 0 0] 0.135; ←Su Su [0 1 -1 0 0 0 0] 0.135; ←
SuModel transport; ← (algebraic or fixed )
equivalenceRatio equivalenceRatio [0 0 0 0 0 0 0] 0.6;
sigmaExt sigmaExt [0 0 -1 0 0 0 0] 338;
2 OpenFOAM for beginner
1.3 XiFoam/ constant/ constant
②combustionProperties 2
XiModel transport; ← Xi algebraic or fixed
XiCoef XiCoef [0 0 0 0 0 0 0] 0.62;
XiShapeCoef XiShapeCoef [0 0 0 0 0 0 0] 1;XiModel
uPrimeCoef uPrimeCoef [0 0 0 0 0 0 0] 1;
2 OpenFOAM for beginner
1.3 XiFoam/ constant
②combustionProperties 3
GuldersCoeffs ← Guilders{
3
Methane
{W 0.422;
Propane
{W 0.446;
IsoOctane
{W 0.4658;W 0.422;
eta 0.15;xi 5.18;alpha 2;beta -0.5;f 2.3;}
W 0.446;eta 0.12;xi 4.95;alpha 1.77;beta -0.2;f 2.3;}
W 0.4658;eta -0.326;xi 4.48;alpha 1.56;beta -0.22;f 2.3;}
Gulder.O.L 23th Symp(int) on Combust ,pp743-
2 OpenFOAM for beginner
1.3 XiFoam/ constant/ constant
②combustionProperties 4
ignite yes; ← on/off
IgnitionSites ←(
{location (0.005 -0.02 0); ←location (0.005 -0.02 0); ←Diameter 0.003; ←start 0; ←duration 0.05; ←strength 20; ← ( )
});
2 OpenFOAM for beginner
1.3 XiFoam/ constant/ constant
②combustionProperties 5
ignitionSphereFraction 0; ← ( ) 3D
ignitionThickness ignitionThickness [0 1 0 0 0 0 0] 0.001; ←
ignitionCircleFraction 1; ← ( ) 2D
ignitionKernelArea ignitionKernelArea [0 2 0 0 0 0 0] 0; ←
2 OpenFOAM for beginner
example
/ XiFoamLES(Weller pitzDaily 2D RANS or 3D LES )RAS(moriyoshiHomogeneous→RAS(moriyoshiHomogeneous→
reactingFoam2D RANS(k-e model) 1
dieselFoam3D aachenBomb RANS
chemkin 35 115 3954 290
engineFoam(LANL Kiva ?)fireFoam(poolFire Fire Dynamics Simulator ?)
2 OpenFOAM for beginner
OpenFOAM extension or development version
rostock LTTExtensionsubversion 1.6.x
http://janus.fms.uni-rostock.de/viewvc/repository/OpenFOAM/branches/
CERFACS (colin,Poinsot et al)thickenedFlame model XiFoam
PDF (ILDM, S.Pope)
freiberg OF+Cantera(alternate(steady)ReactingFoam)OpenFoamWiki 1.5.xOpenFoamWiki 1.5.x
http://openfoamwiki.net/index.php/Contrib_alternateReactingFoam
Cantera interface
reactingFoam (ex GRI-3.0)
2 OpenFOAM for beginner
OpenFOAM extension or development version
Sweden chalmersOpenFOAM ph.D
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2009/
(Håkan Nilsson )
reactingFoam dieselFoamph.D
5th openFoam Workshop
IHIShape optimization of valve channel with incompressible to
compressible simulation for pressure Loss calculation
2 OpenFOAM for beginner
OpenFOAM
★ 7/30 CAE
OpenFOAMOpenFOAM
OpenFOAM
(LES Flamelet DT )
. B 75(759), 2299-2307, 2009-11-25. B 75(759), 2299-2307, 2009-11-25
→OpenFOAM( )
※ 1 FrontFlow/red
FFR ( )
OpenFoam
2 OpenFOAM for beginner
OpenFOAM
ex) reactingFoam +
5th OpenFOAM workshopdev extension
2 OpenFOAM for beginner