CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Solving non-conservative equations

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

LinkBack Thread Tools Display Modes
Old   July 14, 2010, 06:15
Default Solving non-conservative equations
Senior Member
Join Date: Apr 2009
Posts: 118
Rep Power: 10
lost.identity is on a distinguished road
I'm trying to analytically solve a simple diffusion problem (written in non-conservative form) with that of numerical simulation that essentially solves the equations in conservative form.

The transport equation which I'm solving numerically is
\frac{\partial\rho{c}}{\partial{t}} + \nabla\cdot\rho\boldsymbol{u}c = \nabla\cdot(D\nabla{c}) ------------------------(1)

where c is the scalar variable I'm interested in looking at, D is the diffusivity (a constant), \rho is the density and \boldsymbol{u} is the velocity vector. Since it's a diffusion problem I set the velocities to zero, i.e. the second (convection) term on the LHS of the above equation becomes zero. So essentially I'm solving the equation

\frac{\partial\rho{c}}{\partial{t}}  = \nabla\cdot(D\nabla{c})

When I solve the above equation analytically assume that density is constant (which is also a constant in my numerical calculations), and then I obtain the following equation

\frac{\partial{c}}{\partial{t}}  = \frac{D}{\rho}\nabla^2{c}--------------(2)

which is just the diffusion equation. The problem is my solutions are a bit different (i.e. numerical and analytic solutions. This might well be a problem in my discretisation etc. But first of all I just wanted to see if what I'm doing make sense, i.e comparing the numerical results from equation (1) with the analytical result of equation (2). Note that in solving equation (1) in my code even though I set velocities to zero, I still solve the momentum equations.
lost.identity is offline   Reply With Quote


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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 15:56.