exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bug fix: Maxwell transport. M. Nobrega

Browse source
This commit is contained in:
Hrvoje Jasak 2015-04-29 13:59:22 +01:00 committed by Dominik Christ
parent c080289664
commit a6d1c4021f
1 changed files with 6 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
src/transportModels/viscoelastic/viscoelasticLaws/Maxwell/Maxwell.C
View file

@ -77,7 +77,6 @@ Foam::tmp<Foam::fvVectorMatrix> Foam::Maxwell::divTau(volVectorField& U) const
- fvc::laplacian(etaPEff/rho_, U, "laplacian(etaPEff,U)") - fvc::laplacian(etaPEff/rho_, U, "laplacian(etaPEff,U)")
+ fvm::laplacian( (etaPEff + etaS_)/rho_, U, "laplacian(etaPEff+etaS,U)") + fvm::laplacian( (etaPEff + etaS_)/rho_, U, "laplacian(etaPEff+etaS,U)")
); );
} }
@ -86,6 +85,9 @@ void Foam::Maxwell::correct()
// Velocity gradient tensor // Velocity gradient tensor
volTensorField L = fvc::grad(U()); volTensorField L = fvc::grad(U());
// Convected derivate term
volTensorField C = tau_ & L;
// Twice the rate of deformation tensor // Twice the rate of deformation tensor
volSymmTensorField twoD = twoSymm(L); volSymmTensorField twoD = twoSymm(L);
@ -93,9 +95,11 @@ void Foam::Maxwell::correct()
fvSymmTensorMatrix tauEqn fvSymmTensorMatrix tauEqn
( (
fvm::ddt(tau_) fvm::ddt(tau_)
+ fvm::div(phi(), tau_)
== ==
etaP_/lambda_*twoD etaP_/lambda_*twoD
- fvm::Sp( 1/lambda_, tau_ ) + twoSymm(C)
- fvm::Sp( 1/lambda_, tau_)
); );
tauEqn.relax(); tauEqn.relax();