Hi everyone:
i want to modify
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:http://www.cfd-online.com/OpenFOAM_D...part/proud.gifstream&) 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 |
BTW, i have add the definion o
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 |
Hi,
Could you find a solution for that error?! I see the same error when modifying a turbulence model. Thanks, |
All times are GMT -4. The time now is 18:24. |