This email from daniel_turci@id.uff.br originates from outside Imperial. Do not click on links and attachments unless you recognise the sender. If you trust the sender, add them to your safe senders list to disable email stamping for this address.

 

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. Turci
D. 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.