# problem in solving such an equation

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 May 26, 2005, 11:44 problem in solving such an equation #1 R Guest   Posts: n/a Sponsored Links 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
 Sponsored Links

 May 27, 2005, 00:48 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

 May 27, 2005, 03:25 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.

 May 27, 2005, 05:01 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).

 May 27, 2005, 07:04 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

 May 27, 2005, 07:49 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post cheng1988sjtu OpenFOAM Bugs 15 May 1, 2016 16:12 lordvon OpenFOAM Running, Solving & CFD 15 October 19, 2015 13:52 nileshjrane OpenFOAM Running, Solving & CFD 4 February 25, 2013 05:13 nileshjrane OpenFOAM Running, Solving & CFD 8 August 26, 2010 12:50 nedved OpenFOAM Running, Solving & CFD 1 November 25, 2008 21:21

 Sponsored Links

All times are GMT -4. The time now is 09:41.

 Contact Us - CFD Online - Top