1. Limited time only! Sign up for a free 30min personal tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Solving non-conservative equations

  1. Jul 14, 2010 #1
    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

    [tex]\frac{\partial\rho{c}}{\partial{t}} + \nabla\cdot\rho\boldsymbol{u}c = \nabla\cdot(D\nabla{c})[/tex] ---(1)

    where c is the scalar variable I'm interested in looking at, D is the diffusivity (a constant), [tex]\rho[/tex] is the density and [tex]\boldsymbol{u}[/tex] 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

    [tex]\frac{\partial\rho{c}}{\partial{t}} = \nabla\cdot(D\nabla{c})[/tex]

    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

    [tex]\frac{\partial{c}}{\partial{t}} = \frac{D}{\rho}\nabla^2{c}[/tex] --------------(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.
  2. jcsd
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook

Can you offer guidance or do you also need help?
Draft saved Draft deleted