# Problem with validation of regionCouple solver

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

May 22, 2014, 19:09
Problem with validation of regionCouple solver
#1
Member

Join Date: Aug 2013
Posts: 47
Rep Power: 11
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
 p_01.png (28.5 KB, 22 views) p_02.png (38.7 KB, 19 views)
Attached Files
 stefanFoam.tar.gz (2.8 KB, 8 views) stefan2D_01.tar.gz (31.9 KB, 3 views) stefan2D_02.tar.gz (41.3 KB, 3 views)

Last edited by tladd; May 23, 2014 at 09:12. Reason: Correct x to -x in paragraph 2 (after simply)