13,372 research outputs found
A new discrete velocity method for Navier-Stokes equations
The relation between Latttice Boltzmann Method, which has recently become
popular, and the Kinetic Schemes, which are routinely used in Computational
Fluid Dynamics, is explored. A new discrete velocity model for the numerical
solution of the Navier-Stokes equations for incompressible fluid flow is
presented by combining both the approaches. The new scheme can be interpreted
as a pseudo-compressibility method and, for a particular choice of parameters,
this interpretation carries over to the Lattice Boltzmann Method.Comment: 28 pages, 8 figure
Sobolev gradients and image interpolation
We present here a new image inpainting algorithm based on the Sobolev
gradient method in conjunction with the Navier-Stokes model. The original model
of Bertalmio et al is reformulated as a variational principle based on the
minimization of a well chosen functional by a steepest descent method. This
provides an alternative of the direct solving of a high-order partial
differential equation and, consequently, allows to avoid complicated numerical
schemes (min-mod limiters or anisotropic diffusion). We theoretically analyze
our algorithm in an infinite dimensional setting using an evolution equation
and obtain global existence and uniqueness results as well as the existence of
an -limit. Using a finite difference implementation, we demonstrate
using various examples that the Sobolev gradient flow, due to its smoothing and
preconditioning properties, is an effective tool for use in the image
inpainting problem
Continuous, Semi-discrete, and Fully Discretized Navier-Stokes Equations
The Navier--Stokes equations are commonly used to model and to simulate flow
phenomena. We introduce the basic equations and discuss the standard methods
for the spatial and temporal discretization. We analyse the semi-discrete
equations -- a semi-explicit nonlinear DAE -- in terms of the strangeness index
and quantify the numerical difficulties in the fully discrete schemes, that are
induced by the strangeness of the system. By analyzing the Kronecker index of
the difference-algebraic equations, that represent commonly and successfully
used time stepping schemes for the Navier--Stokes equations, we show that those
time-integration schemes factually remove the strangeness. The theoretical
considerations are backed and illustrated by numerical examples.Comment: 28 pages, 2 figure, code available under DOI: 10.5281/zenodo.998909,
https://doi.org/10.5281/zenodo.99890
The Navier-Stokes-alpha model of fluid turbulence
We review the properties of the nonlinearly dispersive Navier-Stokes-alpha
(NS-alpha) model of incompressible fluid turbulence -- also called the viscous
Camassa-Holm equations and the LANS equations in the literature. We first
re-derive the NS-alpha model by filtering the velocity of the fluid loop in
Kelvin's circulation theorem for the Navier-Stokes equations. Then we show that
this filtering causes the wavenumber spectrum of the translational kinetic
energy for the NS-alpha model to roll off as k^{-3} for k \alpha > 1 in three
dimensions, instead of continuing along the slower Kolmogorov scaling law,
k^{-5/3}, that it follows for k \alpha < 1. This rolloff at higher wavenumbers
shortens the inertial range for the NS-alpha model and thereby makes it more
computable. We also explain how the NS-alpha model is related to large eddy
simulation (LES) turbulence modeling and to the stress tensor for second-grade
fluids. We close by surveying recent results in the literature for the NS-alpha
model and its inviscid limit (the Euler-alpha model).Comment: 22 pages, 1 figure. Dedicated to V. E. Zakharov on the occasion of
his 60th birthday. To appear in Physica
Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations
A general solution adaptive scheme based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required
- …