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

How to solve 1st step of semi implicit mtd

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 10, 2006, 04:41
Default 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
  Reply With Quote

Old   April 10, 2006, 05:47
Default 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
  Reply With Quote

Old   April 15, 2006, 04:44
Default 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.
  Reply With Quote

Old   April 15, 2006, 06:51
Default 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

  Reply With Quote

Old   May 3, 2006, 08:46
Default 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.
  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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Modeling in micron scale using icoFoam m9819348 OpenFOAM Running, Solving & CFD 7 October 27, 2007 00:36
Doubt on Implicit Methods analyse In India Main CFD Forum 10 March 9, 2007 04:01
Calculating the 1st step of a semi-implicit scheme zonexo Main CFD Forum 0 September 29, 2005 10:05
Physical Reason for stability of Implicit Schemes? radhakrishnan Main CFD Forum 25 December 7, 2000 16:25


All times are GMT -4. The time now is 01:20.