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

Convection-Diffution with implicit source

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 16, 2010, 05:40
Default Convection-Diffusion with implicit source
  #1
Senior Member
 
Ivan Flaminio Cozza
Join Date: Mar 2009
Location: Torino, Piemonte, Italia
Posts: 207
Rep Power: 9
ivan_cozza is on a distinguished road
Send a message via MSN to ivan_cozza
Hi Foamers,
I need to solve this convection-diffusion problem:

d(ut_i)/dt + d(U_j ut_i) / dx_j - nu d2(ut_i) / dx_j2 = - ut_j d(U_i) / dx_j

where ut is my variable, U is a velocity field known.

I implemented it like this:

fvVectorMatrix utEqn
(
fvm::ddt(ut)
+ fvm::div(phi, ut)
- fvm::laplacian(nu, ut)
==
- fvc::div(phiUt, U)
);

utEqn.relax();
utEqn.solve();

but, I prefer to introduce - fvc::div(phiUtf, U) in an implicit way.
When I try to do - fvm::div(phiUtf, U), it gives me error.
What can be a good way to do it?
Thank you,
Ivan

Last edited by ivan_cozza; February 16, 2010 at 09:17. Reason: Error in typing
ivan_cozza 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
Compiling gmshFoam with OpenFOAM-1.5 BlGene Open Source Meshers: Gmsh, Netgen, CGNS, ... 10 August 6, 2009 04:26
OpenFOAM on MinGW crosscompiler hosted on Linux allenzhao OpenFOAM Installation 127 January 30, 2009 20:08
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 11, 2000 03:43


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