Hi all,
I am running an 2D incompressible Navier Stokes simulation with a dimensionless setup and using the AerForces filter.
I am trying to reproduce results from the literature that use the drag coefficient defined as Cd = 2*Fd/(rho*U*U*D), where Fd is the drag force per unit of length and D is the diameter of a cylinder. I wish to know whether the value exported to the file DragLift.fce is half the Cd as defined previously or the Cd itself. I think it is the former, but I would like a confirmation.
For reference, I am sending a snippet of the session xml file.
<SOLVERINFO>
<I PROPERTY="SolverType" VALUE="VelocityCorrectionScheme" />
<I PROPERTY="EQTYPE" VALUE="UnsteadyNavierStokes" />
<I PROPERTY="EvolutionOperator" VALUE="Nonlinear" />
<I PROPERTY="Projection" VALUE="Galerkin" />
<I PROPERTY="GlobalSysSoln" VALUE="DirectStaticCond" />
<I PROPERTY="TimeIntegrationMethod" VALUE="IMEXOrder2" />
<I PROPERTY="Driver" VALUE="Standard" />
<I PROPERTY="SPECTRALHPDEALIASING" VALUE="True" />
</SOLVERINFO>
<PARAMETERS>
<P> TimeStep = 0.01 </P>
<P> FinTime = 20 </P>
<P> NumSteps = FinTime/TimeStep </P>
<P> IO_CheckSteps = NumSteps/20 </P>
<P> IO_InfoSteps = 5 </P>
<P> IO_FiltersInfoSteps = 100 </P>
<P> IO_CFLSteps = 10 </P>
<P> Re = 66 </P>
<P> Kinvis = 1/Re </P>
</PARAMETERS>
<FILTERS>
<FILTER TYPE="AeroForces">
<PARAM NAME="OutputFile"> DragLift </PARAM>
<PARAM NAME="OutputFrequency"> 5 </PARAM>
<PARAM NAME="Boundary"> B[3] </PARAM>
</FILTER>
-- Daniel Rubano B. TurciD. Sc. Student at Universidade Federal Fluminense (UFF).
Laboratory of Thermal Sciences - LATERMO.
Rua Passo da Pátria, 156, Bloco D, sala 302, Niterói. RJ - 24210-240, Brazil.