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

A source term in a vector equation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 15, 2010, 05:35
Default A source term in a vector equation
  #1
New Member
 
Join Date: Sep 2010
Posts: 3
Rep Power: 6
CpCaniggia is on a distinguished road
Hello everyone,
I am solving a passive vector equation, along with the Navier-Stokes equations. There is a source term related to the main variable (AU & fvc::grad(U)). How can I treat it in an implicit way? The code is as follows,

fvVectorMatrix AUEqn
(
fvm::ddt(AU)
+ fvm::div(phi, AU)
- fvm::laplacian(PrRa_12, AU)
==
- AU & fvc::grad(U)
);

AUEqn.relax();

if (momentumPredictor)
{
solve
(
AUEqn == Arhok*g - fvc::grad(Ap)
);
}

In the present code, "AU & fvc::grad(U)" is treated explicitly. However, the solution is unstable. So I wonder if there is any method to treat it implicitly.

From OpenFoam manual I know that the source term can be solved implicitly by fvm::SuSp(rho , AU). However, rho here must be volScalarField while fvc::grad(U)) in my equation is a volTensorField.

Thank you!
CpCaniggia 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
momentum source term zwdi FLUENT 13 December 5, 2013 18:58
DxFoam reader update hjasak OpenFOAM Post-Processing 69 April 24, 2008 01:24
DecomposePar links against liblamso0 with OpenMPI jens_klostermann OpenFOAM Bugs 11 June 28, 2007 17:51
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 13, 2000 23:03
UDFs for Scalar Eqn - Fluid/Solid HT Greg Perkins FLUENT 0 October 11, 2000 03:43


All times are GMT -4. The time now is 04:21.