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

Probable bug with laplacian steady state

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 23, 2009, 09:53
Default Probable bug with laplacian steady state
  #1
New Member
 
Valerio Novaresio
Join Date: Mar 2009
Location: Polonghera, Cuneo, Italy
Posts: 27
Rep Power: 8
novyno is on a distinguished road
Hi foamers...

I noticed this stange behaviour with OpenFOAM.

I've built a solver that only solve the laplacian (steady state version of laplacianFoam). I've tested this solver with the cavity geometry (icoFoam folder) where I set the boundaries value of my field in this way:

FixedWalls: zeroGradient
movingWall: fixedValue uniform 1
InternalField: fixedValue uniform 1

The analitic solution of this problem is that values of the field converged to the value 1, so in this case the initial solution correspond to the exact solution. The expected behaviour is that at first "time" step the solver of the linear system only does one (or zero) step and the convergence is achived. But, even if the correct solution is achived, the solvers needs many iterations and the initialResidual value doesn't decrese.

I've tried many experiment with different initial values and my oppinion is that ther is a converging problem to the uniform solution.

Someone can say me why?
__________________
...The best way to acquire new knowledge is to share it...
novyno is offline   Reply With Quote

Old   November 23, 2009, 20:31
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Ciao Valerio,

did you try to run the case with the original laplacianFoam? It should give you a set of lines with initial and final residuals equal to zero.

The solution to that problem is always identical to the initial condition since all the fluxes are zero (on fixedWalls by default, and due to the equal value of T on movingWall and inside the system), as a consequence the initial residual should be zero also in the unsteady case.

How did you code the steady solver?

P.S. I'm Italian too
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
mass flow in is not equal to mass flow out saii CFX 2 September 18, 2009 08:07
Mass Diffusion: Transient and Steady State BC rval CFX 3 November 19, 2008 01:52
plz help,urgent, vof model steady state Garima Chaudhary FLUENT 2 May 30, 2007 04:38
steady state, laminar vof_model Garima Chaudhary FLUENT 0 May 24, 2007 03:11


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