# Source Term due to evaporation in energy transport equation

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

 February 5, 2014, 13:05 Source Term due to evaporation in energy transport equation #1 Member   Join Date: May 2012 Posts: 55 Rep Power: 7 Hello foamers, I'm struggling with the implementation of an evaporation model in interFoam. So far I was able to implement a interface reconstruction algorithm, which selects interface cells and their neighboring cells. With this algorithm I'm able to calculate the local evaporating mass flux: The interface reconstruction provides the geometric data needed to calculate the temperature gradient of the liquid and vapor side. The calculation of the temperature gradient is not performed directly in the cells that contain a part of the interface but in the neighboring cells that contain pure liquid or gas. where d_int is the distance from the interface to the centre of the neighboring cell. The mass flux (mDot) through the liquid-vapor interface can be transferred to a energy source term which is included in the energy transport equation: The source term for the continuity equation accounts for the divergence of the velocity field in cells where mass is removed or added: Due to the fact that the velocity field is not free of divergence in cells which contain source terms, an additional source term is added to the transport equation: The source terms in the continuity and transport equation are smeared over a band of cells to avoid numerical errors as Hardt&Wondra mentioned in their paper: http://www.sciencedirect.com/science...21999108001228 The implementation of the source terms in the transport and continuity equation seems to be qualitatively correct. I've tested the code with the stefan problem test case. The liquid is evaporating, if I turn off the source term S in the energy equation. But if source term S is turned on, the temperature of the interface is oscillating, which causes an oscillating evaporation flux. Depending on the time step, the fluid is evaporating or condensating. The corresponding code for solving the energy equation is: Code:  fvScalarMatrix TEqn ( fvm::ddt(rhoCp,T) + fvm::div(rhoCpPhi,T) - fvm::laplacian(k12,T) == S //sharp explicit source term (only in interface cells) ); I assume the source of the oscillations is the explicit source term. How can I solve the energy equation so that the interface temperature is at saturation temperature?

 March 24, 2014, 04:16 #2 Member   Mahdi S Join Date: Jul 2012 Posts: 34 Rep Power: 6 Hi Would you tell me how you implement the transport and continuity source terms? I think you added the source term into pEqu and then use the implicit MULE algorithm for VOF transport equation. If so, I think it is important to check the order solver reads. I mean the heat transfer source term is in Tequ while it might be called before pEqu.

 September 3, 2014, 21:24 #3 Member   james wilson Join Date: Aug 2014 Location: Orlando, Fl Posts: 38 Rep Power: 4 Would you please post your code for your pEqn and alphaEqn? I am working on this problem and am having similar difficulties. Thanks, James

 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 tomas Fluent UDF and Scheme Programming 7 September 10, 2013 02:08 mcathela OpenFOAM Installation 14 April 23, 2013 13:59 tayo OpenFOAM 39 December 3, 2012 06:18 jannnesss CFX 0 January 8, 2010 20:53 Greg Perkins FLUENT 0 October 11, 2000 03:43

All times are GMT -4. The time now is 23:57.