research/laminardiffusionflame
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

added job.mpi sbatch script and system/ and constant/ update

Browse source
This commit is contained in:
ignis 2018-11-11 05:32:07 +09:00
parent a0759f2acc
commit 94a6f182b6
6 changed files with 26 additions and 56 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
constant/thermophysicalProperties
View file

@ -30,8 +30,8 @@ inertSpecie N2;
chemistryReader foamChemistryReader; chemistryReader foamChemistryReader;
foamChemistryFile "$FOAM_CASE/constant/prager/grimech30.foam"; foamChemistryFile "$FOAM_CASE/constant/mech/methane-air-2step-BFER/reactions";
foamChemistryThermoFile "$FOAM_CASE/constant/prager/thermo30.foam"; foamChemistryThermoFile "$FOAM_CASE/constant/mech/methane-air-2step-BFER/thermos";
// ************************************************************************* // // ************************************************************************* //

12
job.mpi Normal file
View file

@ -0,0 +1,12 @@
#!/bin/bash
#SBATCH -J testLaminar # Job name
#SBATCH -o job.%j.out # Name of stdout output file (%j expands to jobId)
#SBATCH -n 100 # Total number of mpi tasks requested
#SBATCH --ntasks-per-node=20
#SBATCH --mail-user=ignis@postech.ac.kr
#SBATCH --mail-type=ALL
# Launch MPI-based executable
prun reactingFoam -parallel

6
system/controlDict
View file

@ -29,7 +29,7 @@ deltaT 1.0e-6;
writeControl adjustableRunTime; // clockTime; // timeStep; // writeControl adjustableRunTime; // clockTime; // timeStep; //
writeInterval 1000; // 10000; // 0.05; // 0.000078; //3600; // writeInterval 0.05; // 0.000078; //3600; // 1000; // 10000; //
purgeWrite 0; purgeWrite 0;
@ -41,11 +41,13 @@ secondaryWriteInterval 501; // 3600; //
secondaryPurgeWrite 1; secondaryPurgeWrite 1;
*/ */
/*
secondaryWriteControl clockTime; // timeStep; // secondaryWriteControl clockTime; // timeStep; //
secondaryWriteInterval 3000; // 1; secondaryWriteInterval 3000; // 1;
secondaryPurgeWrite 2; secondaryPurgeWrite 2;
*/
writeFormat ascii; writeFormat ascii;
@ -59,7 +61,7 @@ timePrecision 8;
runTimeModifiable true; runTimeModifiable true;
adjustTimeStep no; adjustTimeStep yes;
maxCo 0.4; maxCo 0.4;

11
system/fvOptions
View file

@ -18,20 +18,21 @@ FoamFile
source1 source1
{ {
type fixedTemperatureConstraint; type fixedTemperatureConstraint;
active false; active true;
timeStart 0.0; timeStart 0.0;
duration 0.1; duration 0.1;
selectionMode cellSet;
cellSet ignitionCells;
fixedTemperatureConstraintCoeffs fixedTemperatureConstraintCoeffs
{ {
selectionMode cellSet;
cellSet ignitionCells;
mode uniform; mode uniform;
temperature table temperature table
2 2
( (
(0 300) (0 300)
(0.01 2000) (0.01 2250)
) )
; ;
} }

21
system/fvSchemes
View file

@ -29,26 +29,13 @@ divSchemes
{ {
default none; default none;
/*
div(phi,U) Gauss linearUpwindV cellLimited 1;
div(phi,Yi_h) Gauss linearUpwind cellLimited 1;
div(phi,K) Gauss linearUpwind cellLimited 1;
div(phid,p) Gauss linearUpwind cellLimited 1;
div(phi,epsilon) Gauss linearUpwind cellLimited 1;
div(phi,k) Gauss linearUpwind cellLimited 1;
div((muEff*dev2(T(grad(U))))) Gauss linear;
div((interpolate((((De|ng)|Te)*grad(Te)))&S),ne) Gauss linearUpwind cellLimited 1;
div(ve,ne) Gauss linearUpwind cellLimited 1;
*/
div(phi,U) Gauss limitedLinearV 1; div(phi,U) Gauss limitedLinearV 1;
div(phi,Yi_h) Gauss limitedLinear 1; div(phi,Yi_h) Gauss limitedLinear 1;
div(phi,K) Gauss limitedLinear 1; div(phi,K) Gauss limitedLinear 1;
div(phid,p) Gauss limitedLinear 1; div(phid,p) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1; div(phi,epsilon) Gauss limitedLinear 1;
div(phi,k) Gauss limitedLinear 1; div(phi,k) Gauss limitedLinear 1;
div((muEff*dev2(T(grad(U))))) Gauss linear; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div((interpolate((((De|ng)|Te)*grad(Te)))&S),ne) Gauss limitedLinear 1;
div(ve,ne) Gauss limitedLinear 1;
} }
laplacianSchemes laplacianSchemes
@ -66,11 +53,5 @@ snGradSchemes
default orthogonal; default orthogonal;
} }
fluxRequired
{
default no;
p;
}
// ************************************************************************* // // ************************************************************************* //

28
system/fvSolution
View file

@ -53,33 +53,7 @@ solvers
Yi Yi
{ {
solver smoothSolver; $hFinal;
smoother DILU;
tolerance 1e-6;
relTol 0;
// $hFinal;
}
ne
{
$Yi;
tolerance 1e-9;
nNeSubCycles 1000;
}
Phi
{
solver PCG;
preconditioner DIC;
tolerance 1e-08;
relTol 0.2;
}
PhiFinal
{
$Phi;
tolerance 1e-08;
relTol 0.0;
} }
} }