Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Reduced navier stokes in mathematica help please

  1. Jan 7, 2013 #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 (Text):

    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 (Text):

    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.
     
  2. jcsd
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

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



Similar Discussions: Reduced navier stokes in mathematica help please
Loading...