[gaza]

Hi Foamers
I downloaded the solver for evaporation from this post (at the end is a link to git source with hints how to compile it provided by Wyldckat):

http://www.cfd-online.com/Forums/ope...come-join.html

I wanted to validate this solver on simplest case, namely 1D Stefan Problem. However
my results are in strong contradiction with analytical solution. For example: for thermophysical properties which I used, the interface position after 100 s should be about 7.62 mm (thermophysical properties are taken from Bitan Shu conference article). But numerical results are about 0.5 mm. Maybe I set up my test case in a wrong way?
Can anybody check it? Below is my test case. You can run it with Allrun script (after solver compilation according to Wyldckat hints).

[gaza]

Is this solver works well? Does anybody know the answer?

[nesca87]

Hello pblasiak,

I know it's an old thread, but did you happen to solve your problem?

Regards,

Helvines

[gaza]

Hi,
I did not solve the problem. I think the solver is good but only for evaporation and I cannot set it up correctly.

[vidyadhar]

Originally Posted by gaza

Hi,
I did not solve the problem. I think the solver is good but only for evaporation and I cannot set it up correctly.

Hello pblasiak,
I am also trying to validate the stefan problem as you described: To track the displacement of evaporating liquid.
I am trying by using Fluent. However I am facing difficulties. Can this be done using OpenFoam.
Any advise is appreciated.


Thanks in advance!
vidyadhar

[gaza]

Originally Posted by vidyadhar

Hello pblasiak,
I am also trying to validate the stefan problem as you described: To track the displacement of evaporating liquid.
I am trying by using Fluent. However I am facing difficulties. Can this be done using OpenFoam.
Any advise is appreciated.


Thanks in advance!
vidyadhar

Hello vidyadhar,
Please try Nima Sam solver:
http://nimasamkhaniani.ir/openfoam/phasechangeheatfoam/

[vidyadhar]

Originally Posted by gaza

Hello vidyadhar,
Please try Nima Sam solver:
http://nimasamkhaniani.ir/openfoam/phasechangeheatfoam/

Thank you pblasiak for suggesting the Nima Sam solver. I think the code works only on OpenFoam 2.2 version
However, there were errors when I tried to compile with openFoam 4 installed on my PC. The errors are as below:


vid:~/phaseChangeHeatFoam/Application$ ./Allwmake
+ wmake libso smoothInterfaceProperties
wmakeLnInclude: linking include files to ./lnInclude
Making dependency list for source file smoothInterfaceProperties.C
could not open file alphaContactAngleFvPatchScalarField.H for source file smoothInterfaceProperties.C due to No such file or directory
could not open file cyclicAMILduInterface.H for source file smoothInterfaceProperties.C due to No such file or directory
could not open file cyclicAMIPolyPatch.H for source file smoothInterfaceProperties.C due to No such file or directory
g++ -std=c++0x -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3 -DNoRepository -ftemplate-depth-100 -I/opt/openfoam4/src/transportModels/twoPhaseInterfaceProperties/alphaContactAngle/alphaContactAngle -I/opt/openfoam4/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam4/src/OpenFOAM/lnInclude -I/opt/openfoam4/src/OSspecific/POSIX/lnInclude -fPIC -c smoothInterfaceProperties.C -o Make/linux64GccDPInt32Opt/smoothInterfaceProperties.o
smoothInterfaceProperties.C:27:49: fatal error: alphaContactAngleFvPatchScalarField.H: No such file or directory
compilation terminated.
/opt/openfoam4/wmake/rules/General/transform:8: recipe for target 'Make/linux64GccDPInt32Opt/smoothInterfaceProperties.o' failed
make: *** [Make/linux64GccDPInt32Opt/smoothInterfaceProperties.o] Error 1
+ wmake phaseChangeHeatFoam
Making dependency list for source file phaseChangeTwoPhaseMixtures/Hardt/Hardt.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/Lee/Lee.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/newPhaseChangeTwoPhaseMixture.C
Making dependency list for source file phaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture/phaseChangeTwoPhaseMixture.C
Making dependency list for source file twoPhaseMixtureI/twoPhaseMixtureI.C
Making dependency list for source file phaseChangeHeatFoam.C
could not open file turbulenceModel.H for source file phaseChangeHeatFoam.C due to No such file or directory
could not open file fvIOoptionList.H for source file phaseChangeHeatFoam.C due to No such file or directory
g++ -std=c++0x -m64 -Dlinux64 -DWM_ARCH_OPTION=64 -DWM_DP -DWM_LABEL_SIZE=32 -Wall -Wextra -Wold-style-cast -Wnon-virtual-dtor -Wno-unused-parameter -Wno-invalid-offsetof -O3 -DNoRepository -ftemplate-depth-100 -I/opt/openfoam4/src/transportModels -I/opt/openfoam4/src/transportModels/incompressible/lnInclude -I/opt/openfoam4/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam4/src/turbulenceModels/incompressible/turbulenceModel -IphaseChangeTwoPhaseMixtures/phaseChangeTwoPhaseMixture -I/opt/openfoam4/src/finiteVolume/lnInclude -I/opt/openfoam4/src/meshTools/lnInclude -I/opt/openfoam4/src/fvOptions/lnInclude -I/opt/openfoam4/src/sampling/lnInclude -IlnInclude -I. -I/opt/openfoam4/src/OpenFOAM/lnInclude -I/opt/openfoam4/src/OSspecific/POSIX/lnInclude -fPIC -c phaseChangeHeatFoam.C -o Make/linux64GccDPInt32Opt/phaseChangeHeatFoam.o
phaseChangeHeatFoam.C:46:29: fatal error: turbulenceModel.H: No such file or directory
compilation terminated.
/opt/openfoam4/wmake/rules/General/transform:8: recipe for target 'Make/linux64GccDPInt32Opt/phaseChangeHeatFoam.o' failed
make: *** [Make/linux64GccDPInt32Opt/phaseChangeHeatFoam.o] Error 1



Can you let me know how to resolve these issues.
Thanks & Regards,
vidyadhar

[nimasam]

it works only at OpenFOAM-2.2, for other versions of OpenFOAM minor changes in source code is needed.


Best regards

[vidyadhar]

Originally Posted by nimasam

it works only at OpenFOAM-2.2, for other versions of OpenFOAM minor changes in source code is needed.


Best regards

Hello Nima Samkhaniani,


I have installed phaseChangeHeatFoam and have followed the example case given for validating with Stefan problem.


Query: Can I use this solver to predict heat and mass transfer rates due to evaporation from a static liquid-vapor interface? The interface position is maintained at the same location by balancing the mass lost by evaporation with mass supply from the liquid side.

I want to determine mass and heat transfer from the static meniscus in the figure attached herewith.


Thanks & Regards,
Vidyadhar