![]() | ||
The viscous vortex domains (VVD) method is a mesh-free method of computational fluid dynamics for directly numerically solving 2D Navier-Stokes equations in Lagrange coordinates It doesn't implement any turbulence model and free of arbitrary parameters. The main idea of this method is to present vorticity field with discrete regions (domains), which travel with diffusive velocity relatively to fluid and conserve their circulation. The same approach was used in Diffusion Velocity method of Ogami and Akamatsu, but VVD uses other discrete formulas
Contents
Features
The VVD method deals with viscous incompressible fluid. The viscosity and density of fluid is considered to be constant. Method can be extended for simulation of heat conductive fluid flows (viscous vortex-heat domains method)
The main features are:
Governing equations
The VVD method is based on a theorem, that circulation in viscous fluid is conserved on contours travelling with speed
where V is fluid velocity, Vd — diffusion velocity, ν — kinematic viscosity. This theorem shows resemblance with Kelvin's circulation theorem, but it works for viscid flows.
Basing on this theorem, flow region with non-zero circulation is presented with number of domains (small regions with finite volumes), which move with velocity u and thus their circulation
Discrete formulas
Fluid velocity V in point r can be calculated with help of Biot-savart law
where i indexes domains in flow, ri — tracking point of domain and γi — his circulation. δ is a so-called "radius of discreteness" — small value that smooths the vortex and helps to get rid of singularity in the domain tracking point. It equals to mean distance between domains.
Calculation of diffusion velocity is more difficult
First fraction produces vortex-vortex interaction (i — vortex index).
And second fraction represents vortex-boundary repulsion. It helps to calculate ∇Ω near body surface and properly describe boundary layer.
Here k indexes boundary segments, rk — its center, dSk — its normal.