SIMPLE for viscous NavierStokes

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 July 11, 2005, 12:44 I have modified the UEqn.H and #1 New Member   Juergen Almanstoetter Join Date: Mar 2009 Posts: 10 Rep Power: 8 I have modified the UEqn.H and hEqn.h in buoyantSimpleFoam to represent steady laminar viscous Navier-Stokes: volTensorField gradU = fvc::grad(U); tmp 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

 July 11, 2005, 13:15 UEqn.H() will not be correct a #2 Senior Member   Join Date: Mar 2009 Posts: 854 Rep Power: 13 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 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 );

 March 10, 2006, 00:23 Hi friends, I also want to #3 Senior Member   Guoxiang Join Date: Mar 2009 Posts: 109 Rep Power: 8 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Pipin FLUENT 3 July 14, 2008 12:05 Yan Kai Main CFD Forum 0 April 18, 2007 03:48 Yan Kai FLUENT 0 April 13, 2007 23:17 raj Main CFD Forum 2 August 23, 2006 06:26 vivek Main CFD Forum 2 July 27, 2005 06:00

All times are GMT -4. The time now is 19:24.

 Contact Us - CFD Online - Top