mirror of
https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git
synced 2025-12-08 06:37:44 +00:00
Update rCFD fluidized bed case to improve stability.
This commit is contained in:
Thomas Lichtenegger
@ -52,7 +52,7 @@ surfaceScalarFields
|
||||
|
||||
//verbose true;
|
||||
|
||||
couplingSubStep 3;
|
||||
couplingSubStep 7;
|
||||
|
||||
initialRecSteps 150;
|
||||
|
||||
|
||||
@ -1,49 +0,0 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
object thermophysicalProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
thermoType
|
||||
{
|
||||
type hePsiThermo;
|
||||
mixture pureMixture;
|
||||
transport const;
|
||||
thermo eConst;
|
||||
equationOfState perfectGas;
|
||||
specie specie;
|
||||
energy sensibleInternalEnergy;
|
||||
}
|
||||
|
||||
mixture
|
||||
{
|
||||
specie
|
||||
{
|
||||
nMoles 1;
|
||||
molWeight 42.1;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Cv 1310;
|
||||
Hf 0;
|
||||
}
|
||||
transport
|
||||
{
|
||||
mu 1e-05;
|
||||
Pr 0.72;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -27,4 +27,8 @@ Cv Cv [ 0 2 -2 -1 0 0 0 ] 1310;
|
||||
|
||||
molM molM [1 0 0 0 -1 0 0 ] 0.0421;
|
||||
|
||||
TMax TMax [0 0 0 1 0 0 0] 400;
|
||||
|
||||
TMin TMin [0 0 0 1 0 0 0] 250;
|
||||
|
||||
// ************************************************************************* //
|
||||
|
||||
20488
tutorials/rcfdemSolverCoupledHeatTransfer/FluidizedBed_HeatSources/rCFD/CFD/orig.0/T
Normal file → Executable file
20488
tutorials/rcfdemSolverCoupledHeatTransfer/FluidizedBed_HeatSources/rCFD/CFD/orig.0/T
Normal file → Executable file
File diff suppressed because it is too large
Load Diff
@ -51,6 +51,63 @@ libs (
|
||||
|
||||
functions
|
||||
{
|
||||
Tave
|
||||
{
|
||||
type volFieldValue;
|
||||
libs ("libfieldFunctionObjects.so");
|
||||
|
||||
log true;
|
||||
writeControl timeStep;
|
||||
writeInterval 10;
|
||||
writeFields false;
|
||||
|
||||
regionType all;
|
||||
operation volAverage;
|
||||
|
||||
weightField voidfractionRec;
|
||||
|
||||
fields
|
||||
(
|
||||
T
|
||||
);
|
||||
}
|
||||
|
||||
fluxprobe0
|
||||
{
|
||||
type surfaceFieldValue;
|
||||
libs ("libfieldFunctionObjects.so");
|
||||
writeControl timeStep;
|
||||
writeInterval 100;
|
||||
log true;
|
||||
writeFields false;
|
||||
regionType faceZone;
|
||||
name fzcenter;
|
||||
operation sum;
|
||||
|
||||
fields
|
||||
(
|
||||
phi
|
||||
);
|
||||
}
|
||||
|
||||
fluxprobe1
|
||||
{
|
||||
type surfaceFieldValue;
|
||||
libs ("libfieldFunctionObjects.so");
|
||||
writeControl timeStep;
|
||||
writeInterval 100;
|
||||
log true;
|
||||
writeFields false;
|
||||
regionType faceZone;
|
||||
name fz0;
|
||||
operation sum;
|
||||
|
||||
fields
|
||||
(
|
||||
phi
|
||||
);
|
||||
}
|
||||
|
||||
fieldAverage1
|
||||
{
|
||||
type fieldAverage;
|
||||
@ -125,13 +182,12 @@ functions
|
||||
);
|
||||
}
|
||||
|
||||
|
||||
fieldOutput
|
||||
{
|
||||
type writeObjects;
|
||||
functionObjectLibs ( "libutilityFunctionObjects.so" );
|
||||
exclusiveWriting true;
|
||||
objects ("rhoRecMean" "voidfraction" "voidfractionRec" "voidfractionMean" "voidfractionRecMean" "addSourceMean" "phiSMean" "TMean" "T" "particleTemp" "particleTempMean");
|
||||
objects ("rhoRecMean" "voidfraction" "voidfractionRec" "voidfractionMean" "voidfractionRecMean" "addSourceMean" "TMean" "T" "partTemp" "partTempMean");
|
||||
writeControl timeStep;
|
||||
writeInterval 40000;
|
||||
}
|
||||
@ -177,10 +233,10 @@ functions
|
||||
scalarList newWeights(2);
|
||||
newWeights[0] = w0;
|
||||
newWeights[1] = 1-w0;
|
||||
/*
|
||||
scalarList newWeights(1);
|
||||
newWeights[0] = 1.0;
|
||||
*/
|
||||
|
||||
// scalarList newWeights(1);
|
||||
// newWeights[0] = 1.0;
|
||||
|
||||
weightDict.set("weights",newWeights);
|
||||
#};
|
||||
}
|
||||
|
||||
@ -28,13 +28,6 @@ gradSchemes
|
||||
divSchemes
|
||||
{
|
||||
default Gauss linear;
|
||||
div(phi,U) Gauss limitedLinearV 1;
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
div(phi,epsilon) Gauss limitedLinear 1;
|
||||
div(phi,R) Gauss limitedLinear 1;
|
||||
div(R) Gauss linear;
|
||||
div(phi,nuTilda) Gauss limitedLinear 1;
|
||||
div((nuEff*dev(T(grad(U))))) Gauss linear;
|
||||
div(phi,T) Gauss limitedLinear 1;
|
||||
// div(phi,T) Gauss upwind;
|
||||
}
|
||||
@ -57,7 +50,6 @@ snGradSchemes
|
||||
fluxRequired
|
||||
{
|
||||
default no;
|
||||
p ;
|
||||
T ;
|
||||
}
|
||||
|
||||
|
||||
@ -17,14 +17,6 @@ FoamFile
|
||||
|
||||
solvers
|
||||
{
|
||||
"(p|rho)"
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
tolerance 1e-6;
|
||||
relTol 0.01;
|
||||
}
|
||||
|
||||
"(correctedField)"
|
||||
{
|
||||
solver PCG;
|
||||
@ -33,26 +25,6 @@ solvers
|
||||
relTol 0.05;
|
||||
}
|
||||
|
||||
"(p|rho)Final"
|
||||
{
|
||||
$p;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
"(U|k|e|epsilon|R|nuTilda|c)"
|
||||
{
|
||||
solver smoothSolver;
|
||||
smoother GaussSeidel;
|
||||
tolerance 1e-05;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
"(U|e|k|nuTilda)Final"
|
||||
{
|
||||
$U;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
T
|
||||
{
|
||||
solver PBiCG;
|
||||
|
||||
@ -39,7 +39,6 @@ variable rfPW2 equal 0.1
|
||||
variable Tpart equal 330
|
||||
|
||||
variable dt equal 0.0025
|
||||
#0.001
|
||||
variable skin equal 0.0005
|
||||
###############################################
|
||||
|
||||
@ -112,13 +111,9 @@ variable Tave equal v_sumT1/v_np
|
||||
|
||||
fix printheat all print 10 "${time} ${Tave}" file ../DEM/temp_ave.txt title "#time T_ave"
|
||||
|
||||
#fix tdist1 all ave/histo 1000 1 1000 324 380 56 f_Temp[0] mode vector file ../DEM/temp_histo.txt title1 "Particle temperatures [K]"
|
||||
#fix tdist2 all ave/histo 1000 1 1000 324 380 448 f_Temp[0] mode vector file ../DEM/temp_histo_fine.txt title1 "Particle temperatures [K]"
|
||||
|
||||
fix tdist1 all ave/histo 400 1 400 324 380 56 f_Temp[0] mode vector file ../DEM/temp_histo.txt title1 "Particle temperatures [K]"
|
||||
fix tdist2 all ave/histo 400 1 400 324 380 448 f_Temp[0] mode vector file ../DEM/temp_histo_fine.txt title1 "Particle temperatures [K]"
|
||||
|
||||
|
||||
fix integr all nve/sphere
|
||||
|
||||
#screen output
|
||||
@ -131,3 +126,6 @@ dump dmp all custom/vtk 10000 ../DEM/post/dump*.liggghts_coupled.vtk id type x
|
||||
|
||||
run 1
|
||||
set region total property/atom Temp ${Tpart}
|
||||
|
||||
# heat sources are already set in the restart file; if they would not have been - or if restarting in a different way -, set them now
|
||||
#set region insReg property/atom heatSource 3.3379e-04 # = V_p x \dot{q}
|
||||
|
||||
Reference in New Issue
Block a user