|February 5, 2014, 13:05||
Source Term due to evaporation in energy transport equation
Join Date: May 2012
Posts: 55Rep Power: 7
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:
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:
fvScalarMatrix TEqn ( fvm::ddt(rhoCp,T) + fvm::div(rhoCpPhi,T) - fvm::laplacian(k12,T) == S //sharp explicit source term (only in interface cells) );
|March 24, 2014, 04:16||
Join Date: Jul 2012
Posts: 34Rep Power: 6
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||
Join Date: Aug 2014
Location: Orlando, Fl
Posts: 38Rep Power: 4
Would you please post your code for your pEqn and alphaEqn? I am working on this problem and am having similar difficulties.
|Thread||Thread Starter||Forum||Replies||Last Post|
|add source term to energy equation by udf||tomas||Fluent UDF and Scheme Programming||7||September 10, 2013 02:08|
|build problem swak4Foam OF 2.2.0||mcathela||OpenFOAM Installation||14||April 23, 2013 13:59|
|funkySetFields compilation error||tayo||OpenFOAM||39||December 3, 2012 06:18|
|Is it possible to use divergence in the source term of an scalar transport equation?||jannnesss||CFX||0||January 8, 2010 20:53|
|UDFs for Scalar Eqn - Fluid/Solid HT||Greg Perkins||FLUENT||0||October 11, 2000 03:43|