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

Problem with validation of regionCouple solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 22, 2014, 19:09
Default Problem with validation of regionCouple solver
  #1
Member
 
Tony Ladd
Join Date: Aug 2013
Posts: 32
Rep Power: 3
tladd is on a distinguished road
I am working on a solver for Stefan problems - i.e. a sharp interface (typically moving) separating two different regions. I started with conjugateHeatFoam from foam-extend-3.0. I solve a Laplace equation for the pressure in two regions with different permeability and then solve for a concentration field (convection-diffusion) in the inner region. Solver is attached (stefanFoam.tar.gz).

In the test problem the interface between the regions has a cosine profile x=1+A*cos(2*pi*y) with -0.5<y<0.5, with A constant. For clarity I choose both perms to be unity so that on the boundary between the domains the pressure should be simply -x if the pressure at the inlet is zero and the pressure gradient at the outlet is -1. The plots show the pressure from OpenFOAM (blue) vs the analytic solution (red). If the amplitude A is 0.02 (or larger) the solution is correct but for A=0.01 the variation in pressure is much too small (by more than an order of magnitude). In fact there is a more or less discontinuous change in behavior; A=0.014 is fine and A=0.013 is not. I have played with grid resolution, domain size, fvSchemes and fvSolution but to no avail.

The illustrated problems can be run from the attached tarballs. The python functions in scripts/tools/py can be used to make new blockMeshDict files if needed.

Does anyone have any suggestions?

Regards

Tony
Attached Images
File Type: png p_01.png (28.5 KB, 12 views)
File Type: png p_02.png (38.7 KB, 10 views)
Attached Files
File Type: gz stefanFoam.tar.gz (2.8 KB, 4 views)
File Type: gz stefan2D_01.tar.gz (31.9 KB, 1 views)
File Type: gz stefan2D_02.tar.gz (41.3 KB, 2 views)

Last edited by tladd; May 23, 2014 at 09:12. Reason: Correct x to -x in paragraph 2 (after simply)
tladd 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
Is there a problem in the Euler solver? Combas SU2 4 March 28, 2014 17:48
[ANSYS Meshing] Problem with axisymmetric solver Rasel ANSYS Meshing & Geometry 0 June 20, 2012 05:35
CFX-Pre problem, pls help!!! cth_yao CFX 0 February 17, 2012 01:52
mesh.update problem in a new FSI solver ICL OpenFOAM 0 October 8, 2011 15:16
patching problem unsteady solver yellow-stuff Main CFD Forum 0 September 25, 2009 02:26


All times are GMT -4. The time now is 07:19.