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

Help for the initial field 4 lid-driven flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 21, 2008, 02:42
Default Help for the initial field 4 lid-driven flow
  #1
John
Guest
 
Posts: n/a
Dear All,

I want to set the initial condition for a lid-driven square cavity flow, which has a u=1 top BC and static BC at other boundaries. I use a Chebychev spectral collocation solver.

Suppose initially the fluid in the cavity is at rest, then there will be serious gradient of u with respect y , which cause the following computation diverge.

Did anyone do similar work? Or could you give any hint?

many thanks in advance

John
  Reply With Quote

Old   October 21, 2008, 08:13
Default Re: Help for the initial field 4 lid-driven flow
  #2
mettler
Guest
 
Posts: n/a
are you using the stream-function/vorticity equations to solve this? That is the best way to go about it,in my opinion.
  Reply With Quote

Old   October 21, 2008, 10:26
Default Re: Help for the initial field 4 lid-driven flow
  #3
otd
Guest
 
Posts: n/a
I solved this with a transient technique (algorithm based on the sola half of sola-vof).

At least for low Reynolds numbers (laminar flows), it converges quickly starting with u = v = 0 as an initial condition.
  Reply With Quote

Old   October 21, 2008, 22:24
Default Re: Help for the initial field 4 lid-driven flow
  #4
John
Guest
 
Posts: n/a
I use the primitive variables for a SPECTRAL (not FVM, FEM ,etc !) collocation solver, say, velocity and pressure. If at t=0, I give a static flowfield as the initial solution except the top boundary (u=1), then the discontinuity in velocity near that boundary will cause an extremely large gradient of du/dy and then ruin the following calculation.

However, this seems not a problem for FVM or FEM.

So I just want ask this question to anyone who has the experience of spectral solver ?

Thanks any way
  Reply With Quote

Old   October 22, 2008, 05:18
Default Re: Help for the initial field 4 lid-driven flow
  #5
Jed
Guest
 
Posts: n/a
I've used spectral collocation for this sort of problem. What sort of linear and nonlinear solvers are you using that don't converge? In what way do they not converge? One method is to use continuation in Re (start by solving the linear problem with Re=0) to solve steady N-S (this will fail as you approach transition to turbulence). If you are doing unsteady flow, use the steady solution or the linear Stokes solution as initial condition (you only need a smooth initial condition which is compatible with the boundary conditions).
  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
alphaEqn.H in twoPhaseEulerFoam cheng1988sjtu OpenFOAM Bugs 14 September 18, 2011 13:46
Forces in OF15 richard OpenFOAM Running, Solving & CFD 166 September 17, 2011 03:24
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM 6 April 12, 2011 11:24
Modeling in micron scale using icoFoam m9819348 OpenFOAM Running, Solving & CFD 7 October 27, 2007 00:36


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