# temperature equation for one phase in multiphase solver

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

 April 13, 2010, 05:06 temperature equation for one phase in multiphase solver #1 Senior Member     Roman Thiele Join Date: Aug 2009 Location: London, UK Posts: 368 Rep Power: 13 Sponsored Links Hej, I have been working with the interPhaseChangeFoam solver for a while now, and have started writing a phase change model based on temperature. Right now the temperature is still given via the dictionary input. I wanted to add the energy equation for the fluid phase to this solver. The governing equation is which I tried implementing like this Code: ```{ // const volScalarField rhoAlpha = cp * rho1 * alpha1; // for time derivative volScalarField rhoAlpha ( "rhoAlpha", cp * rho1 * alpha1 ); surfaceScalarField phiAlpha ( "phiAlpha", phi * rho1 * cp * fvc::interpolate(alpha1) //check with rhoPhi from the alphaEqn.H file // rhoPhi * cp //* fvc::interpolate(alpha1) ); volScalarField lambdaFAlpha ( "lambdaFAlpha", lambdaF * alpha1 ); Pair > mDotAlphal = twoPhaseProperties->mDotAlphal(); const volScalarField mDotAlphac = mDotAlphal[0](); fvScalarMatrix TEqn ( fvm::ddt(rhoAlpha,T) + fvm::div(phiAlpha, T) - fvm::laplacian(lambdaFAlpha, T) == mDotAlphac * ifg ); TEqn.relax(); TEqn.solve(); Info << "Max(T): " << max(T).value() << " K " << "Min(T): " << min(T).value() << " K" << endl; }``` Unfortunately, I have had no success in actually using this within the solver, the results are completely bogus. One thing that is obvious is that there might be a problem with the alpha1 attachment to the equation. My first thought was to maybe be able to extract the mesh that contains water alpha1>0, and impose boundary conditions onto the interface cells. Is there a way to do it, or does somebody else see what could be wrong with the implementation of this equation? __________________ ~roman

 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 ameya OpenFOAM Running, Solving & CFD 0 December 5, 2008 02:16 james FLUENT 5 July 18, 2006 02:01 james FLUENT 0 June 19, 2006 03:54 Kerem FLUENT 4 March 27, 2006 08:20 Nicolas Main CFD Forum 0 December 13, 2003 15:19