CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General 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 Search this Thread Display Modes
Old   October 21, 2008, 03: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, 09: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, 11: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, 23: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, 06: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 Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Forces in OF15 richard OpenFOAM Running, Solving & CFD 180 July 9, 2018 11:54
alphaEqn.H in twoPhaseEulerFoam cheng1988sjtu OpenFOAM Bugs 15 May 1, 2016 17:12
Extrusion with OpenFoam problem No. Iterations 0 Lord Kelvin OpenFOAM Running, Solving & CFD 8 March 28, 2016 12:08
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 05:03
Modeling in micron scale using icoFoam m9819348 OpenFOAM Running, Solving & CFD 7 October 27, 2007 01:36


All times are GMT -4. The time now is 15:34.