[uyelken]

Hi guys, im using chtmultiregionfoam mixed with ddtfoam. When i define the BCs for conjugate heat transfer, im taking error. On the other hand, when i give zerogradient for BCs, it is working well

This is from xieqn.h

#include "compute_sL.H"
  

    alphaUnburned = oneKGMS * (-5.47712e-01*pow(fH,8)+1.44848e+00*pow(fH,7)-1.59433e+00*pow(fH,6) +9.49494e-01*pow(fH,5)-3.33735e-01*pow(fH,4)+7.10672e-02*pow(fH,3)-9.01563e-03*pow(fH,2)+5.74169e-04*fH+2.59170e-05);
    alphaUnburned *= pow(p/pref,0.200);

    const volScalarField& rhoUnburned = rhou;
    
    aUnburned = alphaUnburned / rhoUnburned;
    aUnburned.max(1e-10); // avoid division by zero in calculation of gCr

    dimensionedScalar eps0("eps0",dimensionSet(0,2,-3,0,0),1e-9);
    volScalarField epsilon = max(turbulence.epsilon(),eps0);

    dimensionedScalar k0("k0",dimensionSet(0,2,-2,0,0),1e-9);
    volScalarField k = max(turbulence.k(),k0);
    volScalarField up = uPrimeCoef*sqrt(2.0/3.0*k);
    up.max(1e-9);

    
    volScalarField lInt_Peters = 0.37*pow(up,3)/epsilon;	// Peters (CFX)

volScalarField b = scalar(1.0) - c;

muUnburned = oneKGMS * (-1.08702e-04*pow(fH,2)+1.70856e-04*fH+1.59871e-05);
p.max(1.0); // avoid p=0 in the following expression;
muUnburned *= pow(p/pref,0.205);
muUnburned.max(1e-8);	// bounding


// quenching:
    nu = thermo.mu()/rho;
    gCr = pow(sL,2)/aUnburned;
    epsilonCr = 15*nu*pow(gCr,2);
    epsilonCr.max(0.1); // can become zero on BCs
    volScalarField& lInt = lInt_Peters;
    volScalarField lKolmogorov = pow(pow(nu,3)/epsilon,0.25); 
    volScalarField sigma = max(0.28*Foam::log(lInt/lKolmogorov), scalar(1e-15));
    G = 0.5 * erfc(-pow(2*sigma,-0.5)*(log(epsilonCr/epsilon) + sigma/2.0));

The error:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2   in "/lib64/libc.so.6"
#3   in "/lib64/libm.so.6"
#4  Foam::pow(Foam::Field<double>&, Foam::UList<double> const&, double const&) at ??:?
#5  void Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensioned<double> const&) at ??:?
#6  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::pow<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::dimensioned<double> const&) at ??:?
#7
 at ??:?
#8
 at ??:?
#9  __libc_start_main in "/lib64/libc.so.6"
#10
 at ??:?