New Reply

reduced navier stokes in mathematica urgent help please

 
Share Thread Thread Tools
Jan7-13, 10:57 AM   #1
 

reduced navier stokes in mathematica urgent help please


ok I am modelling airflow in the upper airway for application i obstructive sleep apnoea, but I have hit a brick wall with mathematica. I have a system of 3 differential equations with boundary conditions, and I need to solve to find 3 unknown functions numerically so that they may be plotted in various graphs.

The equations are as follows:

D[a[x]*u[x], x] == 0,
u[x] u'[x] == -p'[x],
p[x] - 1 == 2 (1 - ((a[x])^(-3/2))) - 50 (a''[x]).

with boundary conditions:

u[0] == 0.1, a[0] == 1, a[10] == 1, p[10] == 1.

so initially I tried to use NDSolve like so..
Code:
NDSolve[{D[a[x]*u[x], x] == 0, u[x] u'[x] == -p'[x], 
  p[x] - 1 == 2 (1 - ((a[x])^(-3/2))) - 50 (a''[x]), u[0] == 0.1, 
  a[0] == 1, a[10] == 1, p[10] == 1}, {a}, {x, 0, 10}]
but mathematica does this:
Code:
Power::infy: "Infinite expression 1/0.^(3/2) encountered. "
Infinity::indet: Indeterminate expression 0. ComplexInfinity encountered. >>
Infinity::indet: Indeterminate expression 0. ComplexInfinity encountered. >>
General::stop: Further output of Infinity::indet will be suppressed during this calculation. >>
NDSolve::ndnum: Encountered non-numerical value for a derivative at x == 0.`. >>
which is super annoying, any pointers as to where I'm going wrong would be great. I'm not even sure if I should be using NDSolve so let me know what you think.
thanks in advance
a.
PhysOrg.com
PhysOrg
science news on PhysOrg.com

>> Ants and carnivorous plants conspire for mutualistic feeding
>> Forecast for Titan: Wild weather could be ahead
>> Researchers stitch defects into the world's thinnest semiconductor
New Reply

Tags
mathematica, navier-stokes, urgent
Thread Tools


Similar Threads for: reduced navier stokes in mathematica urgent help please
Thread Forum Replies
Navier Stokes Eq Differential Equations 1
Navier stokes Differential Equations 6
Navier Stokes? Calculus & Beyond Homework 0
Why is there a prize of $1 million for showing that there are . . . ? Differential Equations 4
What are Navier-stokes equations and why are they difficult to solve? General Math 5