|
[Sponsors] |
July 14, 2010, 06:15 |
Solving non-conservative equations
|
#1 |
Senior Member
Join Date: Apr 2009
Posts: 118
Rep Power: 17 |
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 ------------------------(1) where c is the scalar variable I'm interested in looking at, D is the diffusivity (a constant), is the density and 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 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 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. |
|
|
|
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 18:07 |