CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

problem in solving such an equation

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   May 26, 2005, 11:44
Default problem in solving such an equation
  #1
R
Guest
 
Posts: n/a
Hi, I have problem in solving such an equation

dy/dt=u(t)-S(r)*[1+(dy/dr)^2] where u(t)=U[1+Vsin(w*t)]

We are solving for y(r) and U is an constant =37, V >1. S(r)=S*(r-a)^(1/3) and S is also constant =3.33.

Follwing time depedent boundary conditions is used.

dy/dt=0 if u(t)>=S(a) dy/dr=0 if u(t)<=S(a)

I do not have problem when I solve this equations when V<1. Also I do not have problem when S(r)=S even if V>1.

Please help me thanks

  Reply With Quote

Old   May 27, 2005, 00:48
Default Re: problem in solving such an equation
  #2
harish
Guest
 
Posts: n/a
What kind of numerical scheme are you using to solve the problem ?

-Harish
  Reply With Quote

Old   May 27, 2005, 03:25
Default Re: problem in solving such an equation
  #3
R
Guest
 
Posts: n/a
Actually I have used CN, Fully implicit and expilcit with back differencing and central differencing schemes.
  Reply With Quote

Old   May 27, 2005, 05:01
Default Re: problem in solving such an equation
  #4
Tom
Guest
 
Posts: n/a
Have you tried to write down the analytic solution - your pde is only first order so you can write down a general expression for the characteristics - from which you should see the source of your problem!

Also isn't S(a)=0 your boundary condition so that dy/dt = 0 on r=a (I assame r=a is the boundary from your notation).
  Reply With Quote

Old   May 27, 2005, 07:04
Default Re: problem in solving such an equation
  #5
R
Guest
 
Posts: n/a
I have tried the analytic solutions of such equation. The analytic solution is only possible for the case of V<<1. But the case of V>1 analytical solution is wrong.

Can you give me examples how can analytical solutions helps in such case.

Yes you are write about the boundary condition. R

  Reply With Quote

Old   May 27, 2005, 07:49
Default Re: problem in solving such an equation
  #6
Tom
Guest
 
Posts: n/a
I'm not sure what you mean by the analytic solution being wrong - unless you mean you tried to obtain a series solution in V<<1? Otherwise the (correct) analytic solution cannot be wrong since it is the solution to the equation!

What I meant by anaytic solution was that you can reduce your pde, using the method of characteristics, to a set of odes for

dt/ds, dr/ds, dy/ds dp/ds and dq/ds where p=dy/dt, q=dy/dr and s is the variable along the characteristic.

The analytic solution helps for two reasons:

(1) It gives something to test your numerics against and

(2) If there is a problem with the existence of the solution (i.e. finite time blowup) then you will be able to see why and in what parameter regions it happens.
  Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SLTS+rhoPisoFoam: what is rDeltaT??? nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13
alphaEqn.H in twoPhaseEulerFoam cheng1988sjtu OpenFOAM Bugs 14 September 18, 2011 13:46
Velocity blows up suddenly after 30,000+ iterations lordvon OpenFOAM 14 October 27, 2010 11:13
Error while running rhoPisoFoam.. nileshjrane OpenFOAM Running, Solving & CFD 8 August 26, 2010 12:50
SimpleFoam k and epsilon bounded nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21


All times are GMT -4. The time now is 11:43.