AeroForces definition

29 May 2024

      *******************
This email 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 https://spam.ic.ac.uk/SpamConsole/Senders.aspx 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.*

Daniel Rubano Barretto Turci

Alexandra Liosi

tags

participants (2)