tutorials: Fixes to TDAC cases
Properties have been removed that are set in the standard TDAC ".cfg" file, and ".orig" files have been used to better ensure that cleanCase restores the original state. Sandia has also had it's TDAC parameters slightly tweaked for stability.
This commit is contained in:
@ -36,8 +36,6 @@ odeCoeffs
|
||||
|
||||
reduction
|
||||
{
|
||||
tolerance 1e-4;
|
||||
|
||||
// Search initiating set (SIS) of species, needed for most methods
|
||||
initialSet
|
||||
(
|
||||
|
||||
@ -1,11 +0,0 @@
|
||||
#!/bin/sh
|
||||
cd ${0%/*} || exit 1 # Run from this directory
|
||||
|
||||
# Source tutorial clean functions
|
||||
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
|
||||
|
||||
cleanCase
|
||||
|
||||
rm -f 0/T
|
||||
|
||||
#------------------------------------------------------------------------------
|
||||
@ -15,121 +15,32 @@ FoamFile
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#includeEtc "caseDicts/solvers/chemistry/TDAC/chemistryPropertiesFlame.cfg"
|
||||
|
||||
chemistryType
|
||||
{
|
||||
solver ode;
|
||||
method TDAC;
|
||||
}
|
||||
|
||||
chemistry on;
|
||||
chemistry on;
|
||||
|
||||
initialChemicalTimeStep 1e-07;
|
||||
initialChemicalTimeStep 1e-7;
|
||||
|
||||
odeCoeffs
|
||||
{
|
||||
solver seulex;
|
||||
absTol 1e-08;
|
||||
relTol 0.1;
|
||||
}
|
||||
|
||||
reduction
|
||||
{
|
||||
// Activate reduction
|
||||
active on;
|
||||
|
||||
// Switch logging of the reduction statistics and performance
|
||||
log on;
|
||||
|
||||
// Tolerance depends on the reduction method, see details for each method
|
||||
tolerance 1e-4;
|
||||
|
||||
// Available methods: DRG, DAC, DRGEP, PFA, EFA
|
||||
method DAC;
|
||||
|
||||
// Search initiating set (SIS) of species, needed for most methods
|
||||
initialSet
|
||||
(
|
||||
CO
|
||||
CH4
|
||||
HO2
|
||||
);
|
||||
|
||||
// For DAC, option to automatically change the SIS switch from HO2 to H2O
|
||||
// and CO to CO2, + disable fuel
|
||||
automaticSIS off;
|
||||
|
||||
// When automaticSIS, the method needs to know the fuel
|
||||
fuelSpecies
|
||||
(
|
||||
(CH4 1)
|
||||
);
|
||||
absTol 1e-8;
|
||||
relTol 1e-1;
|
||||
}
|
||||
|
||||
tabulation
|
||||
{
|
||||
// Activate tabulation
|
||||
active on;
|
||||
|
||||
// Switch logging of the tabulation statistics and performance
|
||||
log on;
|
||||
|
||||
printProportion off;
|
||||
|
||||
printNumRetrieve off;
|
||||
|
||||
// Tolerance used for retrieve and grow
|
||||
tolerance 0.003;
|
||||
|
||||
// ISAT is the only method currently available
|
||||
method ISAT;
|
||||
|
||||
// Scale factors used in the definition of the ellipsoid of accuracy
|
||||
scaleFactor
|
||||
{
|
||||
otherSpecies 1;
|
||||
Temperature 1000;
|
||||
Pressure 1e15;
|
||||
deltaT 1;
|
||||
}
|
||||
|
||||
// Maximum number of leafs stored in the binary tree
|
||||
maxNLeafs 5000;
|
||||
|
||||
// Maximum life time of the leafs (in time steps) used in unsteady
|
||||
// simulations to force renewal of the stored chemPoints and keep the tree
|
||||
// small
|
||||
chPMaxLifeTime 1000;
|
||||
|
||||
// Maximum number of growth allowed on a chemPoint to avoid distorted
|
||||
// chemPoints
|
||||
maxGrowth 100;
|
||||
|
||||
// Number of time steps between analysis of the tree to remove old
|
||||
// chemPoints or try to balance it
|
||||
checkEntireTreeInterval 500;
|
||||
|
||||
// Parameters used to decide whether to balance or not if the tree's depth
|
||||
// is larger than maxDepthFactor*log2(nLeafs) then balance the tree
|
||||
maxDepthFactor 2;
|
||||
|
||||
// Try to balance the tree only if the size of the tree is greater
|
||||
minBalanceThreshold 30;
|
||||
|
||||
// Activate the use of a MRU (most recently used) list
|
||||
MRURetrieve false;
|
||||
|
||||
// Maximum size of the MRU list
|
||||
maxMRUSize 0;
|
||||
|
||||
// Allow to grow points
|
||||
growPoints true;
|
||||
|
||||
// When mechanism reduction is used, new dimensions might be added
|
||||
// maxNumNewDim set the maximum number of new dimensions added during a
|
||||
// growth
|
||||
maxNumNewDim 10;
|
||||
|
||||
variableTimeStep true;
|
||||
}
|
||||
|
||||
#include "reactionsGRI"
|
||||
|
||||
@ -20,7 +20,6 @@ then
|
||||
runApplication -a foamDictionary -entry "startTime" -set "0" system/controlDict
|
||||
runApplication -a foamDictionary -entry "writeInterval" -set "1500" system/controlDict
|
||||
runApplication -a foamDictionary -entry "endTime" -set "1500" system/controlDict
|
||||
rm -rf constant/chemistryProperties > /dev/null 2>&1
|
||||
runApplication -a foamDictionary -entry "chemistry" -set "off" constant/chemistryProperties
|
||||
|
||||
runApplication $application
|
||||
@ -29,12 +28,9 @@ then
|
||||
runApplication -a foamDictionary -entry "startTime" -set "1500" system/controlDict
|
||||
runApplication -a foamDictionary -entry "writeInterval" -set "100" system/controlDict
|
||||
runApplication -a foamDictionary -entry "endTime" -set "5000" system/controlDict
|
||||
rm -rf constant/chemistryProperties > /dev/null 2>&1
|
||||
runApplication -a foamDictionary -entry "chemistry" -set "on" constant/chemistryProperties
|
||||
fi
|
||||
|
||||
runApplication -o $application
|
||||
|
||||
rm -rf constant/chemistryProperties > /dev/null 2>&1
|
||||
|
||||
#------------------------------------------------------------------------------
|
||||
|
||||
@ -26,7 +26,6 @@ chemistryType
|
||||
chemistry on;
|
||||
|
||||
initialChemicalTimeStep 1e-7;
|
||||
//maxChemicalTimeStep 1e-3;
|
||||
|
||||
odeCoeffs
|
||||
{
|
||||
@ -35,14 +34,13 @@ odeCoeffs
|
||||
relTol 1e-1;
|
||||
}
|
||||
|
||||
reduction
|
||||
{
|
||||
tolerance 1e-4;
|
||||
}
|
||||
|
||||
tabulation
|
||||
{
|
||||
tolerance 3e-3;
|
||||
// Scale factors used in the definition of the ellipsoid of accuracy
|
||||
scaleFactor
|
||||
{
|
||||
Temperature 1000;
|
||||
}
|
||||
}
|
||||
|
||||
#include "reactionsGRI"
|
||||
|
||||
Reference in New Issue
Block a user