Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Menu
Open sidebar
Development
openfoam
Commits
94f952ab
Commit
94f952ab
authored
Oct 05, 2015
by
Henry Weller
Browse files
reactingMultiphaseEulerFoam: Completed LTS support in multuphaseSystem
parent
9577ad2d
Changes
1
Hide whitespace changes
Inline
Side-by-side
applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/multiphaseSystem/multiphaseSystem.C
View file @
94f952ab
...
...
@@ -155,19 +155,43 @@ void Foam::multiphaseSystem::solveAlphas()
}
}
MULES
::
limit
(
1
.
0
/
mesh_
.
time
().
deltaT
().
value
(),
// ***HGW add support for LTS
geometricOneField
(),
phase
,
phi_
,
alphaPhiCorr
,
zeroField
(),
zeroField
(),
phase
.
alphaMax
(),
0
,
true
);
if
(
fv
::
localEulerDdt
::
enabled
(
mesh_
))
{
const
volScalarField
&
rDeltaT
=
fv
::
localEulerDdt
::
localRDeltaT
(
mesh_
);
MULES
::
limit
(
rDeltaT
,
geometricOneField
(),
phase
,
phi_
,
alphaPhiCorr
,
zeroField
(),
zeroField
(),
phase
.
alphaMax
(),
0
,
true
);
}
else
{
const
scalar
rDeltaT
=
1
.
0
/
mesh_
.
time
().
deltaTValue
();
MULES
::
limit
(
rDeltaT
,
geometricOneField
(),
phase
,
phi_
,
alphaPhiCorr
,
zeroField
(),
zeroField
(),
phase
.
alphaMax
(),
0
,
true
);
}
}
MULES
::
limitSum
(
alphaPhiCorrs
);
...
...
@@ -481,7 +505,7 @@ Foam::multiphaseSystem::multiphaseSystem
IOobject
(
"alphas"
,
mesh
.
time
().
timeName
(),
mesh
_
.
time
().
timeName
(),
mesh
,
IOobject
::
NO_READ
,
IOobject
::
AUTO_WRITE
...
...
@@ -496,13 +520,13 @@ Foam::multiphaseSystem::multiphaseSystem
deltaN_
(
"deltaN"
,
1e-8
/
pow
(
average
(
mesh
.
V
()),
1
.
0
/
3
.
0
)
1e-8
/
pow
(
average
(
mesh
_
.
V
()),
1
.
0
/
3
.
0
)
)
{
forAll
(
phases
(),
phasei
)
{
volScalarField
&
alphai
=
phases
()[
phasei
];
mesh
.
setFluxRequired
(
alphai
.
name
());
mesh
_
.
setFluxRequired
(
alphai
.
name
());
}
}
...
...
@@ -599,13 +623,12 @@ Foam::multiphaseSystem::nearInterface() const
void
Foam
::
multiphaseSystem
::
solve
()
{
const
fvMesh
&
mesh
=
this
->
mesh
();
const
Time
&
runTime
=
mesh
.
time
();
const
Time
&
runTime
=
mesh_
.
time
();
const
dictionary
&
alphaControls
=
mesh_
.
solverDict
(
"alpha"
);
label
nAlphaSubCycles
(
readLabel
(
alphaControls
.
lookup
(
"nAlphaSubCycles"
)));
bool
LTS
=
fv
::
localEulerDdt
::
enabled
(
mesh
);
bool
LTS
=
fv
::
localEulerDdt
::
enabled
(
mesh
_
);
if
(
nAlphaSubCycles
>
1
)
{
...
...
@@ -614,7 +637,7 @@ void Foam::multiphaseSystem::solve()
if
(
LTS
)
{
trSubDeltaT
=
fv
::
localEulerDdt
::
localRSubDeltaT
(
mesh
,
nAlphaSubCycles
);
fv
::
localEulerDdt
::
localRSubDeltaT
(
mesh
_
,
nAlphaSubCycles
);
}
dimensionedScalar
totalDeltaT
=
runTime
.
deltaT
();
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment