I want to find proper solver. (multiphase & heat exchange)
 

Hello! I'm a newb of OF.
i want to find a proper solver that can handle multiphase and heat transfer.
What i want to do is shoot hot air into a pipe half full of water. I expect hot air make water to steam.
Is it possible?? And is there any solver to handle it?
I will appreciate it in advance.

In the ESI version of OpenFOAM there is the solver chtMultiRegionTwoPhaseEulerFoam. But I am not sure, if you can use different gases (air/steam) with it.

See: https://www.openfoam.com/releases/op...nd-physics.php

Thank you so much jherb.
i will try to do this solver.
if i will success, i will compare the result.
Thank you again!

Hi,

Here you have the list and description of OpenFOAM solvers : https://openfoam.com/documentation/u...rd-solvers.php

Less exhaustive but easier to use, the simplified decision matrix :
https://www.openfoam.com/documentati...s-solvers.html

find 2 exctracts pictures of solvers list you should pick in.

I you're looking for a simple and robust User Interface (from CAD import to post-treament), think about ESI Visual-CFD ;)(https://openfoam.com/products/visualcfd.php)

bye
nderev

Hi, sorry to disturb you. I coupled chtMultiRegionFoam with the modified interCondensatingEvaporatingFoam, the solver compiles without problems (just judging from the compiled results), but once the temperature solves, there will be an error. I put up my TEqn.H and 0/T, please help me see where the problem is.

0/BASE/T

dimensions [0 0 0 1 0 0 0];
internalField uniform 0;
boundaryField
{
SJUERE
{
type zeroGradient;
}
SFRONTANDBACK
{
type empty;
}
HEAT
{
type fixedGradient;
gradient uniform 50;
}
BASE_FLUID
{
type compressible::turbulentTemperatureCoupledBaffleMix ed;
value $internalField;
Tnbr T;
kappaMethod solidThermo;
}
}

TEqn.H
{
Info<<"10"<<endl;
volScalarField& T = thermo->T();
Info<<"11"<<endl;
tmp<volScalarField> tcp(thermo->Cp());
const volScalarField& cp = tcp();
Info<<"12"<<endl;
const dimensionedScalar Cp1 = thermo->Cp1();
const dimensionedScalar Cp2 = thermo->Cp2();
Info<<"13"<<endl;
rhoCp = rho*cp;
Info<<"14"<<endl;
kappaEff = thermo->kappa() + rho*cp*turbulence.nut()/Prt;
Info<<"15"<<endl;
const surfaceScalarField rhoCpPhi
(
"rhoCpPhi",
rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2
);
Info<<"16"<<endl;
Pair<tmp<volScalarField>> vDotAlphal = mixture.mDot();
Info<<"17"<<endl;
fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T)
+ fvm::div(rhoCpPhi, T)
- fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoCpPhi), T)
- fvm::laplacian(kappaEff, T)
+ mixture.TSource()
);

Info<<"18"<<endl;
TEqn.relax();
TEqn.solve(mesh.solver(he.select(finalIter)));
thermo->correct();

Info<< "min/max(T) = " << min(T).value() << ", "
<< max(T).value() <<endl;
}

Error:
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3 Foam::mapDistributeBase::mapDistributeBase(Foam::U List<int> const&, Foam::UList<int> const&) at ??:?
#4 Foam::mapDistribute::mapDistribute(Foam::UList<int > const&, Foam::UList<int> const&) at ??:?
#5 Foam::mappedPatchBase::calcMapping() const at ??:?
#6 Foam::mappedPatchBase::map() const at ??:?
#7 void Foam::mappedPatchBase::distribute<double>(Foam::Li st<double>&) const at ??:?
#8 Foam::compressible::turbulentTemperatureCoupledBaf fleMixedFvPatchScalarField::updateCoeffs() at ??:?
#9 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in ~/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/chtConEvaFoam
#10 Foam::temperaturePhaseChangeTwoPhaseMixtures::cons tant::TSource() const at ??:?
#11 ? in ~/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/chtConEvaFoam
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? in ~/OpenFOAM/OpenFOAM-v2006/platforms/linux64GccDPInt32Opt/bin/chtConEvaFoam
Segmentation fault (core dumped)