Commit 935e5a8f authored by Henry Weller's avatar Henry Weller
Browse files

EulerDdtScheme: evaluate dimensioned internal field expressions rather than...

EulerDdtScheme: evaluate dimensioned internal field expressions rather than primitive field expressions
Ensures dimension consistency and simplifies the expressions
parent bc2fd6c3
......@@ -122,12 +122,10 @@ EulerDdtScheme<Type>::fvcDdt
new GeometricField<Type, fvPatchField, volMesh>
(
ddtIOobject,
mesh(),
rDeltaT.dimensions()*vf.dimensions(),
rDeltaT.value()*
rDeltaT*
(
vf.primitiveField()
- vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
vf()
- vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
),
rDeltaT.value()*
(
......@@ -174,12 +172,10 @@ EulerDdtScheme<Type>::fvcDdt
new GeometricField<Type, fvPatchField, volMesh>
(
ddtIOobject,
mesh(),
rDeltaT.dimensions()*rho.dimensions()*vf.dimensions(),
rDeltaT.value()*rho.value()*
rDeltaT*rho*
(
vf.primitiveField()
- vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
vf()
- vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
),
rDeltaT.value()*rho.value()*
(
......@@ -226,13 +222,11 @@ EulerDdtScheme<Type>::fvcDdt
new GeometricField<Type, fvPatchField, volMesh>
(
ddtIOobject,
mesh(),
rDeltaT.dimensions()*rho.dimensions()*vf.dimensions(),
rDeltaT.value()*
rDeltaT*
(
rho.primitiveField()*vf.primitiveField()
- rho.oldTime().primitiveField()
*vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
rho()*vf()
- rho.oldTime()()
*vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
),
rDeltaT.value()*
(
......@@ -282,18 +276,15 @@ EulerDdtScheme<Type>::fvcDdt
new GeometricField<Type, fvPatchField, volMesh>
(
ddtIOobject,
mesh(),
rDeltaT.dimensions()
*alpha.dimensions()*rho.dimensions()*vf.dimensions(),
rDeltaT.value()*
(
alpha.primitiveField()
*rho.primitiveField()
*vf.primitiveField()
alpha()
*rho()
*vf()
- alpha.oldTime().primitiveField()
*rho.oldTime().primitiveField()
*vf.oldTime().primitiveField()*mesh().Vsc0()/mesh().Vsc()
- alpha.oldTime()()
*rho.oldTime()()
*vf.oldTime()()*mesh().Vsc0()/mesh().Vsc()
),
rDeltaT.value()*
(
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment