# code in f90

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

 April 6, 2002, 11:01 code in f90 #1 dav Guest   Posts: n/a Dear All I need to code the following equation U_{t} + F(U)_{x} = S(U) i've read lots of books and ref. etc. but I couldn't code the methods in f90. i was able only to code the case when S(U) = 0 using upwind methods etc. so can somebody help me by giveing me a 90 code to solve any linear or nonlinear case. many thanks

 April 6, 2002, 22:28 Re: code in f90 #2 versi Guest   Posts: n/a Use Strang Splitting when S(U) really is a physcial source term. In each time step you first solve U_t + F_x=0(PDE), then solve U_t=S (ODE). Check some text books on CFD.

 April 8, 2002, 14:24 Re: code in f90 #3 Johnix Guest   Posts: n/a I'm not sure whether I have understood your question. But I have an idea... Maybe u can try it. since F=F(U)_{x}, you can also regard S(U) as S(U,x) i.e. S(U)_{x} then you can get: U_{t} + F(U)_{x} - S(U)_{x} = 0 => U_{t} + G(U)_{x} = 0 where G = F - S

 April 10, 2002, 07:34 Re: code in f90 #4 dav Guest   Posts: n/a S(U) is physical sourc term. so say we solve U_t + F_x=0(PDE) using upwind methods, what shall i use to solve U_t=S (ODE)!!!

 April 10, 2002, 08:12 Re: code in f90 #5 Upwind Guest   Posts: n/a If you are using an upwind scheme for the flux the it is necessary to incoporate source terms in the upwind formulation otherwise the results could be inaccurate depending upon the strength of the source terms. You may want to consult the following paper for the upwind procedure to incorporate source terms: Construction of Second-Order TVD Schemes for Non- homogeneous Hyperbolic Conservation Laws, L. Gascon and J.M. Corberan, J. of Computational Physics, Vol. 172, pp. 261-297(2001).

 April 10, 2002, 21:20 Re: code in f90 #6 versi Guest   Posts: n/a Generally, U_t=S(U) is very stiff, so you must use a stiff ODE solver.

 April 11, 2002, 09:48 Re: code in f90 #7 dav Guest   Posts: n/a what do you mean by: U_t=S(U) is very stiff thanks

 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 RajeshAero Main CFD Forum 5 November 10, 2011 06:48 John C. Chien Main CFD Forum 19 May 17, 2001 15:56 John C. Chien Main CFD Forum 54 April 23, 2001 08:10 Bernhard Mueck Main CFD Forum 10 February 16, 2000 11:07 Heinz Wilkening Main CFD Forum 38 March 5, 1999 12:44

All times are GMT -4. The time now is 02:45.