4,078 research outputs found
Spectral element discretization of the vorticity, velocity and pressure formulation of the Navier-Stokes problem
The two-dimensional Navier–Stokes equations, when provided with non standard boundary conditions which involve the normal component of the velocity and the vorticity, admit a variational formulation with three independent unknowns, the vorticity, the velocity and the pressure. We propose a discretization of this problem by spectral element methods. A detailed numerical analysis leads to optimal error estimates for the three unknowns and numerical experiments confirm the interest of the discretization
Effects of discrete energy and helicity conservation in numerical simulations of helical turbulence
Helicity is the scalar product between velocity and vorticity and, just like
energy, its integral is an in-viscid invariant of the three-dimensional
incompressible Navier-Stokes equations. However, space-and time-discretization
methods typically corrupt this property, leading to violation of the inviscid
conservation principles. This work investigates the discrete helicity
conservation properties of spectral and finite-differencing methods, in
relation to the form employed for the convective term. Effects due to
Runge-Kutta time-advancement schemes are also taken into consideration in the
analysis. The theoretical results are proved against inviscid numerical
simulations, while a scale-dependent analysis of energy, helicity and their
non-linear transfers is performed to further characterize the discretization
errors of the different forms in forced helical turbulence simulations
Effects of discrete energy and helicity conservation in numerical simulations of helical turbulence
Helicity is the scalar product between velocity and vorticity and, just like
energy, its integral is an in-viscid invariant of the three-dimensional
incompressible Navier-Stokes equations. However, space-and time-discretization
methods typically corrupt this property, leading to violation of the inviscid
conservation principles. This work investigates the discrete helicity
conservation properties of spectral and finite-differencing methods, in
relation to the form employed for the convective term. Effects due to
Runge-Kutta time-advancement schemes are also taken into consideration in the
analysis. The theoretical results are proved against inviscid numerical
simulations, while a scale-dependent analysis of energy, helicity and their
non-linear transfers is performed to further characterize the discretization
errors of the different forms in forced helical turbulence simulations
FluSI: A novel parallel simulation tool for flapping insect flight using a Fourier method with volume penalization
FluSI, a fully parallel open source software for pseudo-spectral simulations
of three-dimensional flapping flight in viscous flows, is presented. It is
freely available for non-commercial use under
[https://github.com/pseudospectators/FLUSI]. The computational framework runs
on high performance computers with distributed memory architectures. The
discretization of the three-dimensional incompressible Navier--Stokes equations
is based on a Fourier pseudo-spectral method with adaptive time stepping. The
complex time varying geometry of insects with rigid flapping wings is handled
with the volume penalization method. The modules characterizing the insect
geometry, the flight mechanics and the wing kinematics are described.
Validation tests for different benchmarks illustrate the efficiency and
precision of the approach. Finally, computations of a model insect in the
turbulent regime demonstrate the versatility of the software
A fully semi-Lagrangian discretization for the 2D Navier--Stokes equations in the vorticity--streamfunction formulation
A numerical method for the two-dimensional, incompressible Navier--Stokes
equations in vorticity--streamfunction form is proposed, which employs
semi-Lagrangian discretizations for both the advection and diffusion terms,
thus achieving unconditional stability without the need to solve linear systems
beyond that required by the Poisson solver for the reconstruction of the
streamfunction. A description of the discretization of Dirichlet boundary
conditions for the semi-Lagrangian approach to diffusion terms is also
presented. Numerical experiments on classical benchmarks for incompressible
flow in simple geometries validate the proposed method
Spectral method for the unsteady incompressible Navier-Stokes equations in gauge formulation
A spectral method which uses a gauge method, as opposed to a projection method, to decouple the computation of velocity and pressure in the unsteady incompressible Navier-Stokes equations, is presented. Gauge methods decompose velocity into the sum of an auxilary field and the gradient of a gauge variable, which may, in principle, be assigned arbitrary boundary conditions, thus overcoming the issue of artificial pressure boundary conditions in projection methods. A lid-driven cavity flow is used as a test problem. A subtraction method is used to reduce the pollution effect of singularities at the top corners of the cavity. A Chebyshev spectral collocation method is used to discretize spatially. An exponential time differencing method is used to discretize temporally. Matrix diagonalization procedures are used to compute solutions directly and efficiently. Numerical results for the flow at Reynolds number Re = 1000 are presented, and compared to benchmark results. It is shown that the method, called the spectral gauge method, is straightforward to implement, and yields accurate solutions if Neumann boundary conditions are imposed on the gauge variable, but suffers from reduced convergence rates if Dirichlet boundary conditions are imposed on the gauge variable
Computation of Steady Incompressible Flows in Unbounded Domains
In this study we revisit the problem of computing steady Navier-Stokes flows
in two-dimensional unbounded domains. Precise quantitative characterization of
such flows in the high-Reynolds number limit remains an open problem of
theoretical fluid dynamics. Following a review of key mathematical properties
of such solutions related to the slow decay of the velocity field at large
distances from the obstacle, we develop and carefully validate a
spectrally-accurate computational approach which ensures the correct behavior
of the solution at infinity. In the proposed method the numerical solution is
defined on the entire unbounded domain without the need to truncate this domain
to a finite box with some artificial boundary conditions prescribed at its
boundaries. Since our approach relies on the streamfunction-vorticity
formulation, the main complication is the presence of a discontinuity in the
streamfunction field at infinity which is related to the slow decay of this
field. We demonstrate how this difficulty can be overcome by reformulating the
problem using a suitable background "skeleton" field expressed in terms of the
corresponding Oseen flow combined with spectral filtering. The method is
thoroughly validated for Reynolds numbers spanning two orders of magnitude with
the results comparing favourably against known theoretical predictions and the
data available in the literature.Comment: 39 pages, 12 figures, accepted for publication in "Computers and
Fluids
- …