CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

SIMPLE for viscous NavierStokes

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   July 11, 2005, 12:44
Default I have modified the UEqn.H and
  #1
New Member
 
Juergen Almanstoetter
Join Date: Mar 2009
Posts: 10
Rep Power: 8
almanstoetter is on a distinguished road
I have modified the UEqn.H and hEqn.h in buoyantSimpleFoam to represent steady laminar viscous Navier-Stokes:

volTensorField gradU = fvc::grad(U);
tmp<fvvectormatrix> UEqn
(
fvm::div(phi, U)
- fvm::Sp(fvc::div(phi),U)
- fvm::laplacian(2*mu+lambda, U)
);

UEqn().relax();

solve
(
UEqn() ==
fvc::div
(
mu*gradU.T()+lambda*I*tr(gradU)-(mu+lambda)*gradU
)
- fvc::grad(pd) + fvc::grad(rho)*gh
);

fvScalarMatrix hEqn
(
fvm::div(phi,h)
- fvm::Sp(fvc::div(phi), h)
- fvm::laplacian(k/Cp, h)
);

hEqn.relax();

solve
(
hEqn ==
fvc::div(phi/fvc::interpolate(rho)*fvc::interpolate(p))
- p*fvc::div(phi/fvc::interpolate(rho))
+ ((lambda*I*tr(gradU)+mu*(gradU+gradU.T())) && gradU)


Do UEqn.H() and UEqn.A() in pEqn.H still give the required diagonal and off-diagonal entries ? Has the field phi to be extended to use the SIMPLE algorithm in pEqn.H ?


Juergen
almanstoetter is offline   Reply With Quote

Old   July 11, 2005, 13:15
Default UEqn.H() will not be correct a
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
UEqn.H() will not be correct as you have currently formulated UEqn. You need to move the fvc::div term from the UEqn solve statement to the UEqn construction:

tmp<fvvectormatrix> UEqn
(
fvm::div(phi, U)
- fvm::Sp(fvc::div(phi),U)
- fvm::laplacian(2*mu+lambda, U)
- fvc::div
(
mu*gradU.T()+lambda*I*tr(gradU)-(mu+lambda)*gradU
)
);

UEqn().relax();

solve
(
UEqn() == - fvc::grad(pd) + fvc::grad(rho)*gh
);
henry is offline   Reply With Quote

Old   March 10, 2006, 00:23
Default Hi friends, I also want to
  #3
Senior Member
 
Guoxiang
Join Date: Mar 2009
Posts: 109
Rep Power: 8
liugx212 is on a distinguished road
Hi friends,

I also want to use the buoyantFoam case to simulate the sequestration. It means that I need add the Darcy term and source term in solver. But I am not sure what do I do. Could you please give some advice?

Thanks advance.
Guoxiang
liugx212 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Viscous heating Pipin FLUENT 3 July 14, 2008 12:05
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm Yan Kai Main CFD Forum 0 April 18, 2007 03:48
About Phase Coupled SIMPLE (PC-SIMPLE) algorithm Yan Kai FLUENT 0 April 13, 2007 23:17
navierstokes equation raj Main CFD Forum 2 August 23, 2006 06:26
Navierstokes solver vivek Main CFD Forum 2 July 27, 2005 06:00


All times are GMT -4. The time now is 13:44.