|
[Sponsors] |
temperature equation for one phase in multiphase solver |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 21 ![]() |
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<tmp<volScalarField> > 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; } 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 | Thread Starter | Forum | Replies | Last Post |
Transport Equation for temperature variance | ameya | OpenFOAM Running, Solving & CFD | 0 | December 5, 2008 02:16 |
how to active Temperature equation in FLUENT | james | FLUENT | 5 | July 18, 2006 03:01 |
how to active temperature equation | james | FLUENT | 0 | June 19, 2006 04:54 |
multiphase UDF, variables for phase velocities | Kerem | FLUENT | 4 | March 27, 2006 09:20 |
2 phase flow equation in 1D | Nicolas | Main CFD Forum | 0 | December 13, 2003 15:19 |