Re: DASKR: Initial Conditions for orifice flow equation
- From: Fred <fred.krogh@xxxxxxxxx>
- Date: Sun, 6 Apr 2008 15:18:10 -0700 (PDT)
On Apr 5, 2:22 am, gurry <frederick.the.f...@xxxxxxxxx> wrote:
Hi,There is enough information to know, but perhaps you could get around
I am trying to solve the equation for flow through an orifice using
DASKR (the newest DASSL code).
My equations are:
equation[0] = y[0] - sqrt(y[1])
equation[1] = yPrime[1] - K;
y[1] here is the pressure drop across the orifice and y[0] is the flow
through it. and K is the rate of change of pressure drop (and it is a
positive number).
My initial flow and pressure drop at time = 0 are:
y[0] = 0
y[1] = 0
I have a problem, however, with the initial value of derivatives.
While the value of derivative of y[1] is alwayws going to be K, the
derivative of y[0] at time = 0 must be inifinity. On paper this is
okay because yPrime[0] is going to decay as y[1] increases linearly.
However, in the numerical solver I can't in anyway supply and initial
value of infinity to yPrime[0].
I would like to know if anyone else has faced similar problems with
this kinda equation. If so, how did you handle it?
\
the problem by introducing arc length as the independent variable.
You then have one more differential equation to solve for the
independent variable.
Fred
.
- Follow-Ups:
- References:
- Prev by Date: multidimensional scaling
- Next by Date: Elements of Propulsion Gas Turbines and Rockets - Mattingly
- Previous by thread: DASKR: Initial Conditions for orifice flow equation
- Next by thread: Re: DASKR: Initial Conditions for orifice flow equation
- Index(es):
Relevant Pages
|
