Register to reply

Reduced navier stokes in mathematica help please

Share this thread:
01af
#1
Jan7-13, 10:57 AM
P: 1
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..
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:
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.
Phys.Org News Partner Science news on Phys.org
FIXD tells car drivers via smartphone what is wrong
Team pioneers strategy for creating new materials
Team defines new biodiversity metric

Register to reply

Related Discussions
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