måndag 4 juli 2022

Euler CFD as Parameter Free CFD as ToE

The Euler equations in velocity-pressure $(u,p)$ and $(x,t)$-coordinates are invariant under a rescaling of velocity $u$ into $\bar u =\frac{u}{U}$ with $U$ a reference speed such as free stream speed in bluff body flow with corresponding rescaling of pressure $p$ into $\bar p=\frac{p}{U^2}$ and time $t$ into $\bar t =Ut$ without rescaling of space with thus $\bar x = x$. The scaling of pressure with $U^2$ conforms with Bernoulli's Law and the scaling of drag force $\sim C_DU^2$ from a drag coefficient $C_D$. The propulsion power to balance drag thus scales with $U^3$. The Euler equations are thus formally invariant under change of velocity scale as an expression of formally zero viscosity or infinite Reynolds number.

The basic energy estimate of Euler CFD expresses a balance between rate of loss of kinetic energy and computational residual-based turbulent dissipation of the basic simplified form $C\frac{h}{\vert u\vert}|\vert u\cdot\nabla u\vert^2$ both scaling with $u^3$. The propulsion power is balanced by the rate of loss of kinetic energy and so by turbulent dissipation. The drag coefficient can thus alternatively be computed from total turbulent dissipation.


The (remarkable) fact that the drag coefficient $C_D$ does not include dependence of the Reynolds number $Re$, expresses observations that drag depends little on $Re$ beyond drag crisis, which connects to Kolmogorov's conjecture of finite limit of turbulent dissipation as well as mesh and stabilisation independence in computation. The functionality of the drag coefficient supports Euler's Dream that Euler CFD offers a Theory of Everything ToE for slightly viscous incompressible flow with independence of $Re$ beyond drag crisis. Since total drag shows little dependence on $Re$ while in principle it has a contribution from skin friction with a skin friction coefficient (scaling with $U^2$) decreasing with $Re$, the skin friction contribution appears to be small, in contradiction to a common conception of major contribution: If major drag indeed would come from skin friction, then drag would decrease with increasing $Re$, but it does not beyond drag crisis.


Notice that the Navier-Stokes equations with constant viscosity $\nu$ with turbulent dissipation intensity $\nu\vert\nabla u\vert^2$ scaling with $u^2$, are not velocity scale invariant and thus carry a dependence on $Re$ possibly making computational solution impossible for large $Re$. 


Recall that the definition of $Re =\frac{UL}{\nu}$ with $U$ a reference speed and $L$ a reference length and $\nu$ a viscosity is not well determined and so independence of mean value quantities such as drag, lift and pitch moment is a necessary requirement to make CFD predictable.


Here is experimental evidence that $C_D$ for NACA0012 at zero angle of attack does not depend on $Re$ beyond drag crisis:



Notice the reduction of $C_D$ by a factor 2 from $Re =100.00$ to $Re > 500.00$ as an expression of drag crisis.


Inga kommentarer:

Skicka en kommentar