CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

help with Laplace equation

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 28, 2013, 13:05
Default help with Laplace equation
New Member
Join Date: Feb 2013
Posts: 3
Rep Power: 13
giorgk is on a distinguished road
I started using openFoam few days ago and after I went through alot of documentation and tutorials I decided to try my own very simple example to solve the steady state diffusion equation \nabla (\Gamma \nabla \phi) = 0. According to my OpenFoam understanding so far this is equivalent to solving :
fvm::ddt(T) - fvm::laplacian(DT, T) using steadyState as ddtSchemes.

The domain I uses is also a simple one. A box with dirichlet Bcs left and right, fix gradient on the top and zero gradient everywhere else.
The first run looked pretty good as the solution field varied linearly between the dirichlet conditions.
Then as I increased the fixed Gradient on the top the solution field started forming a Hump, which is reasonable again.

Then while keeping the same fixed gradient value that gave me the Hump in my last run I started increasing the diffuse term DT. I would then expect the Hump to start reducing and actually after some large DT value I would expected to get a solution that varied linearly. However the term DT seems to have no impact whatsoever.
Can anyone tell me what I am doing wrong?

Thank you


PS: General comment about OpenFoam
I feel I have wasted a lot of time by not find it earlier
giorgk is offline   Reply With Quote

Old   March 2, 2013, 14:32
Default What am I doing wrong?? (comsol - openfoam comparison)
New Member
Join Date: Feb 2013
Posts: 3
Rep Power: 13
giorgk is on a distinguished road
I've tried the simplest thing ever and I cannot reproduce in OpenFoam the same solution as in Comsol. (And the one in Comsol is the correct one).
I'm solving laplacian steady state \nabla(K\nabla h) = 0.
My domain is a box with 5km x 5 km and 300 m height.
I assign dirichlet BCs on the left and right 20 m and 40 m respectively and neumann on the top equal to 4.6296e-06 m/s (or 0.4 m/day). The K value is 0.0011574 m/s (100 m/day), and solve it as steady state.
I'm attaching below the results from Comsol and OpenFoam and they are VERY different
This is Comsol solution

and this is OpenFoam solution

So am I using the wrong equation in OpenFoam?

Please someone help

Thank you
giorgk is offline   Reply With Quote

Old   November 20, 2013, 09:22
New Member
Join Date: Jan 2013
Posts: 6
Rep Power: 13
prithvi yesudas is on a distinguished road
i am trying the same problem with same boundary condition. but getting this error while running blockMesh.
Cannot open mesh description file

From function blockMesh
in file blockMeshApp.C at line 148.
Can u pls upload me your case files since i require to run this problem urgently.
prithvi yesudas is offline   Reply With Quote


diffusion coefficient, laplaciafoam

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
mass flow in is not equal to mass flow out saii CFX 12 March 19, 2018 05:21
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
pisoFoam laplace equation fisch OpenFOAM Running, Solving & CFD 10 July 16, 2013 22:13
error message cuteapathy CFX 14 March 20, 2012 06:45
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56

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