CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Dimentsions error in my modification of SST model

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 19, 2008, 05:50
Default Hi everyone: i want to modify
  #1
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 300
Rep Power: 8
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to 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 is offline   Reply With Quote

Old   June 19, 2008, 06:03
Default BTW, i have add the definion o
  #2
Senior Member
 
wayne.zhang
Join Date: Mar 2009
Location: Shanghai, Shanghai, P.R.China
Posts: 300
Rep Power: 8
waynezw0618 is on a distinguished road
Send a message via MSN to waynezw0618 Send a message via Skype™ to 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
waynezw0618 is offline   Reply With Quote

Old   April 19, 2012, 12:36
Default
  #3
Member
 
pooyan
Join Date: Nov 2011
Posts: 61
Rep Power: 4
sam1364 is on a distinguished road
Hi,

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

Thanks,
sam1364 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Turbulence Model Modification and Implement qtian OpenFOAM Running, Solving & CFD 1 January 25, 2009 06:53
design modification alireza CFX 5 February 18, 2008 07:46
modification of UDF ammi FLUENT 2 January 18, 2007 21:35
modification for CFX? caty CFX 5 May 19, 2005 18:48
UDF modification merac FLUENT 4 March 1, 2001 06:42


All times are GMT -4. The time now is 16:13.