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

Phase Change Heat Conduction

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 13, 2011, 10:47
Default Phase Change Heat Conduction
  #1
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Hi all

I am a new user of OPEN FOAM and I want to model a "Phase Change Heat Conduction Problem" with it.

Would you please let me know if there is a tutorial for phase change heat conduction simulation with OPEN FOAM? There are some tutorials for OPENFOAM on the web but most of them are related to flow problems and not the problem that I want to model!

Regards!
hawkeye321 is offline   Reply With Quote

Old   April 14, 2011, 09:07
Default Phase Change in OF
  #2
Senior Member
 
fabian_roesler's Avatar
 
Fabian Roesler
Join Date: Mar 2009
Location: Bad Friedrichshall, Germany
Posts: 154
Rep Power: 8
fabian_roesler is on a distinguished road
Hi mahdiiowa

there is no solver in OF like you describe it. You will have to code it by yourself. The solver which just handles scalar transport is the scalarTransportFoam http://openfoamwiki.net/index.php/ScalarTransportFoam. There you have to delete the convective term and add a source term for the enthalpy jump during the phase change. There is much literature around describing the enthalpy method.

Regards

Fabian
fabian_roesler is offline   Reply With Quote

Old   April 16, 2011, 14:29
Default
  #3
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Dear Fabian
Thanks for you reply!
hawkeye321 is offline   Reply With Quote

Old   April 16, 2011, 15:14
Default
  #4
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Hi
I am now trying to modify scalarTransportFOAM solver! But I cannot where the source file for this solver is!
My directory

Home/OpenFOAM/magma-1.7.1/Run/tutorials/basic/scalarTransportFOAM/pitzDaily

Just contains files

0 constant system
!
I would be so thankful if someone lets me know how I can find the main solver for scalarTransportFOAM!
hawkeye321 is offline   Reply With Quote

Old   April 16, 2011, 15:48
Default
  #5
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Hi
I am now trying to modify scalarTransportFOAM solver! But I cannot where the source file for this solver is!
My directory

Home/OpenFOAM/magma-1.7.1/Run/tutorials/basic/scalarTransportFOAM/pitzDaily

Just contains files

0 constant system
!
I would be so thankful if someone lets me know how I can find the main solver for scalarTransportFOAM!
hawkeye321 is offline   Reply With Quote

Old   April 17, 2011, 04:12
Default
  #6
Senior Member
 
Fabian Braennstroem
Join Date: Mar 2009
Posts: 407
Rep Power: 10
braennstroem is on a distinguished road
Hello,

the attached solver does something similar, I hope it helps... it is based on Zaki Saldi's work. There is a validation case of Solomon's master thesis (I think... it's a long time ago).
Best Regards!
Fabian
Attached Files
File Type: gz meltFoam.tar.gz (5.8 KB, 146 views)
braennstroem is offline   Reply With Quote

Old   April 18, 2011, 14:18
Default
  #7
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Dear Braennstroem
Thank you very much for the file and your kind considerations!
hawkeye321 is offline   Reply With Quote

Old   April 18, 2011, 15:16
Default
  #8
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Dear Braennstroem
Thank you again for the file! I tried to run it in a similar way to running icoFoam!
I set the current directory to
magma@r-soll34:~/OpenFOAM/magma-1.7.1/run/TestMahdi/meltFoam/testCase>
I run the blockMesh!
But when I want to run the meltFoam by typing meltFoam I get the following error

If 'meltFoam' is not a typo you can use command-not-found to lookup the package that contains it, like this:
cnf meltFoam
!
Would you please help me how to resolve this error!
Regards!
hawkeye321 is offline   Reply With Quote

Old   April 19, 2011, 09:36
Default
  #9
New Member
 
hossein
Join Date: Dec 2010
Posts: 14
Rep Power: 6
hossein.mashhadi is on a distinguished road
Dear Mahdi

i think you have not compiled the solver MeltFoam so first do it and after that you can create a case like IcoFoam tutorial and run MeltFoam.

hope it help you.
hossein.mashhadi is offline   Reply With Quote

Old   April 26, 2011, 17:31
Default
  #10
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Dear Hossein
I want to compile pisoFoam.C, butI get the following error

%%%%%
mkdir: cannot create directory `linux64GccDPOpt': Permission denied
/bin/sh: linux64GccDPOpt/options: No such file or directory
make: *** [linux64GccDPOpt/options] Error 1
/opt/openfoam171/wmake/MakefileFiles:39: linux64GccDPOpt/options: No such file or directory
make: *** No rule to make target `linux64GccDPOpt/options'. Stop.
wmake error: file 'Make/linux64GccDPOpt/objectFiles' could not be created
%%%%%%
Would you please help me hove to fix it?

Regards
hawkeye321 is offline   Reply With Quote

Old   April 27, 2011, 00:25
Default
  #11
New Member
 
hossein
Join Date: Dec 2010
Posts: 14
Rep Power: 6
hossein.mashhadi is on a distinguished road
Quote:
Originally Posted by mahdiiowa View Post
Dear Hossein
I want to compile pisoFoam.C, butI get the following error

%%%%%
mkdir: cannot create directory `linux64GccDPOpt': Permission denied
/bin/sh: linux64GccDPOpt/options: No such file or directory
make: *** [linux64GccDPOpt/options] Error 1
/opt/openfoam171/wmake/MakefileFiles:39: linux64GccDPOpt/options: No such file or directory
make: *** No rule to make target `linux64GccDPOpt/options'. Stop.
wmake error: file 'Make/linux64GccDPOpt/objectFiles' could not be created
%%%%%%
Would you please help me hove to fix it?

Regards
Hi Mahdi

i think you should log in as root user (by sudo bash command you can do it) or in file "file" in Make folder change the directory to "FOAM_USER_APPBIN".

i hope it works.
hossein.mashhadi is offline   Reply With Quote

Old   April 28, 2011, 14:35
Default
  #12
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Dear Hossein
Thank you very much for your help!
I logged in as a root user but now I have another problem even in compiling pisoFoam,(standard solver in Users Guide, Page 70)! the error is this
%%%%%%%%%%%%%%%%%%%
g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-40 -I/opt/openfoam171/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam171/src/transportModels -I/opt/openfoam171/src/transportModels/incompressible/singlePhaseTransportModel -I/opt/openfoam171/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam171/src/OpenFOAM/lnInclude -I/opt/openfoam171/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed Make/linux64GccDPOpt/pisoFoam.o -L/opt/openfoam171/lib/linux64GccDPOpt \
-lincompressibleTurbulenceModel -lincompressibleRASModels -lincompressibleLESModels -lincompressibleTransportModels -lfiniteVolume -lmeshTools -lOpenFOAM -liberty -ldl -lm -o /opt/openfoam171/applications/bin/linux64GccDPOpt/pisoFoam
/usr/lib64/gcc/x86_64-suse-linux/4.5/../../../../x86_64-suse-linux/bin/ld: cannot find -liberty
collect2: ld returned 1 exit status
make: *** [/opt/openfoam171/applications/bin/linux64GccDPOpt/pisoFoam] Error 1
magma@r-soll34:/opt/openfoam171/applications/solvers/incompressible/pisoFoam>
%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%
I do not know why it says "cannot find liberty"!!!
hawkeye321 is offline   Reply With Quote

Old   April 29, 2011, 03:53
Default
  #13
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 919
Rep Power: 17
akidess will become famous soon enough
Installation of OpenFOAM 1.7.1 to OprnSUSE 11.2
akidess is offline   Reply With Quote

Old   April 30, 2011, 14:20
Default
  #14
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
Thank you very much akidess!!!
The problem has been fixed and now I can compile successfully!!!
hawkeye321 is offline   Reply With Quote

Old   April 30, 2011, 14:24
Default
  #15
Member
 
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 5
hawkeye321 is an unknown quantity at this point
I want to edit scalarTransportFoam.C by adding another if loop to it!!! in such a way that for T > 5 ; another term be added to the equation!
I did it this way
%%%%%
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
{
fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi, T)
- fvm::laplacian(DT, T)
);
if (TEqn> 5)
{
TEqn += fvm::laplacian(DT, T);
};

TEqn.solve();

}
but when I compile it I get the following error
%%%%%%
ConductionMT.C: In function ‘int main(int, char**)’:
ConductionMT.C:66:23: error: no match for ‘operator>’ in ‘TEqn > 5’
/opt/openfoam171/src/finiteVolume/lnInclude/readSIMPLEControls.H:6:10: warning: unused variable ‘momentumPredictor’
/opt/openfoam171/src/finiteVolume/lnInclude/readSIMPLEControls.H:9:10: warning: unused variable ‘transonic’
make: *** [Make/linux64GccDPOpt/ConductionMT.o] Error 1
%%%%%%

Would you please help me fix the problem??
hawkeye321 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
Two phase flow with phase change Ahmad Al-Zoubi CFX 1 November 26, 2008 04:59
phase change process shen qingqing FLUENT 0 September 7, 2008 09:29
simulation of phase change material melting rosol FLUENT 0 March 25, 2008 01:54
Conduction heat not balanced? Nepal CFX 7 May 31, 2005 11:51
Phase change process Sang-jin Lee CD-adapco 5 December 5, 2001 06:17


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