# How to solve 1st step of semi implicit mtd

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

 April 10, 2006, 04:41 How to solve 1st step of semi implicit mtd #1 newbie Guest   Posts: n/a I'm trying to write a ns code using fractional step. I've read that most people use semi-implicit method ie CN2+AB2. AB2 requires info at 2 time steps. But in the beginning, how is the answer to the 1st step obtained? Thank you

 April 10, 2006, 05:47 Re: How to solve 1st step of semi implicit mtd #2 Darren Guest   Posts: n/a In my very humble opinion... You can use a 1st order method (forward Euler) in step 1 and then switch over to CN/AB2 in the following steps. But, maybe a better way is to just use an RK method for the whole simulation, (like TVD RK2) which doesn't rely on old time step info, is 2nd order accurate and has better stability properties. Hope this helps

 April 15, 2006, 04:44 Re: How to solve 1st step of semi implicit mtd #3 versi Guest   Posts: n/a Crank-Nicolson (CN) is for diffusion terms and better numerical statbility. TVD RK2 has worse stability, I do not recommend it for diffusion problems.

 April 15, 2006, 06:51 Re: How to solve 1st step of semi implicit mtd #4 Darren Guest   Posts: n/a Sorry, I didn't make it very clear in my first post. You can do an operator split and treat the non-linear terms with TVD RK2 and the viscous terms with CN, because as Versi said you need to treat diffusion implicitly to get accetable stability. There are loads of ways to do the time stepping though, so I'm not trying to say this is the "best" or anything although I would be interested to hear what other people do. Darren

 May 3, 2006, 08:46 Re: How to solve 1st step of semi implicit mtd #5 Fuka Guest   Posts: n/a Forward Euler is sufficient for starting a second order multistep method, because local error is second order.

 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 xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33 m9819348 OpenFOAM Running, Solving & CFD 7 October 27, 2007 00:36 analyse In India Main CFD Forum 10 March 9, 2007 04:01 zonexo Main CFD Forum 0 September 29, 2005 10:05 radhakrishnan Main CFD Forum 25 December 7, 2000 16:25

All times are GMT -4. The time now is 00:56.