[waynezw0618]

Hi everyone:
i want to modify the SST model with AIAA paper A98-32817.so i change the definion of member function F1 like this
tmp<volscalarfield> kOmegaSSTRC::F1(const volScalarField& CDkOmega) const
{
//**modification** define the limitation of the CDkw by the maximum value of cross section term in the flow field (AIAA98-2554)
double CDkOmegalim = 1.0e-8*gMax(CDkOmega);

volScalarField CDkOmegaPlus = max
(
CDkOmega,
dimensionedScalar("CDkOmegalim", dimless/sqr(dimTime), CDkOmegalim)
);//**modification** use the CDkOmegalim instead of 1.0e-10,

volScalarField arg1 = min
(
min
(
max
(
(scalar(1)/betaStar)*sqrt(k_)/(omega_*y_),
scalar(500)*nu()/(sqr(y_)*omega_)
),
(4*alphaOmega2)*k_/(CDkOmegaPlus*sqr(y_))
),
scalar(10)
);

return tanh(pow4(arg1));
}

but it was run with error that

-> FOAM FATAL ERROR : Arguments of max have different dimensions
dimensions : [0 0 -1 0 0 0 0] and [0 0 -0.5 0 0 0 0]
#0 Foam::error::printStack(Foam:stream&) in "/home/waynezw0618/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/waynezw0618/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::max(Foam::dimensionSet const&, Foam::dimensionSet const&) in "/home/waynezw0618/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#3 Foam::tmp<foam::geometricfield<double,> > Foam::max<double,>(Foam::tmp<foam::geometricfield< double,> > const&, Foam::tmp<foam::geometricfield<double,> > const&) in "/home/waynezw0618/OpenFOAM/waynezw0618-1.4.1/lib/linux64GccDPOpt/libincompressi bleTurbulenceModels.so"
#4 Foam::turbulenceModels::kOmegaSSTRC::correct() in "/home/waynezw0618/OpenFOAM/waynezw0618-1.4.1/lib/linux64GccDPOpt/libincompressi bleTurbulenceModels.so"
#5 main in "/home/waynezw0618/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/simp leFoam"
#6 __libc_start_main in "/lib64/libc.so.6"
#7 __gxx_personality_v0 in "/home/waynezw0618/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/simp leFoam"


From function max(const dimensionSet& ds1, const dimensionSet& ds2)
in file dimensionSet/dimensionSet.C at line 213.

FOAM aborting

could you please tell me what is wrong with my modification?

thanks !
yours wayne

[waynezw0618]

BTW, i have add the definion of private date F1_ declaration and initialize it in the constructor of the model like this :
omega_
(
IOobject
(
"omega",
runTime_.timeName(),
mesh_,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh_
),


F1_
(
IOobject
(
"F1",
runTime_.timeName(),
mesh_,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh_
),
nut_(a1*k_/max(a1*(omega_ + omegaSmall_), F2()*mag(symm(fvc::grad(U_)))))

and i give a file named F1 in "0" document :
class volScalarField;
object F1;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{

BOTTOM
{
type fixedValue;
value uniform 0;
}

INLET1
{
type fixedValue;
value uniform 0;
}

OUTLET
{
type fixedValue;
value uniform 0;
}
SIDE
{
type fixedValue;
value uniform 0;
}

TOP
{
type fixedValue;
value uniform 0;
}

WING
{
type fixedValue;
value uniform 0;
}


}


// ************************************************** *********************** //

can anyone help me!!
thanks!

wayne

[sam1364]

Hi,

Could you find a solution for that error?!
I see the same error when modifying a turbulence model.

Thanks,