# help with Laplace equation

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

 February 28, 2013, 14:05 help with Laplace equation #1 New Member   Join Date: Feb 2013 Posts: 3 Rep Power: 6 Hi, 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 ( ) = 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 George PS: General comment about OpenFoam I feel I have wasted a lot of time by not find it earlier

 March 2, 2013, 15:32 What am I doing wrong?? (comsol - openfoam comparison) #2 New Member   Join Date: Feb 2013 Posts: 3 Rep Power: 6 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 = 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

 November 20, 2013, 10:22 #3 New Member   prithviyesudas Join Date: Jan 2013 Posts: 6 Rep Power: 6 i am trying the same problem with same boundary condition. but getting this error while running blockMesh. FOAM FATAL ERROR: Cannot open mesh description file "/home/guest/OpenFOAM/guest-2.1.1/run/tutorials/cfdproject/constant/polyMesh/blockMeshDict" 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.

 Tags diffusion coefficient, laplaciafoam

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Mihail CFX 7 September 7, 2014 06:27 fisch OpenFOAM Running, Solving & CFD 10 July 16, 2013 22:13 cuteapathy CFX 14 March 20, 2012 07:45 Sas CFX 15 July 13, 2010 08:56 saii CFX 2 September 18, 2009 08:07

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