5,618 research outputs found
Spectral discretization of the time-dependent Navier-Stokes problem coupled with the heat equation
29 pages avec calculsThe aim of this work is to present the unsteady Navier{Stokes equations coupled with the heat equation where the viscosity depends on the temperature. We propose a discretization of theses equations that combines Euler's implicit scheme in time and spectral methods in space. We prove optimal error estimates between the continuous and discrete solutions. Some numerical experiments con rm the interest of this approach
Volume penalization for inhomogeneous Neumann boundary conditions modeling scalar flux in complicated geometry
We develop a volume penalization method for inhomogeneous Neumann boundary
conditions, generalizing the flux-based volume penalization method for
homogeneous Neumann boundary condition proposed by Kadoch et al. [J. Comput.
Phys. 231 (2012) 4365]. The generalized method allows us to model scalar flux
through walls in geometries of complex shape using simple, e.g. Cartesian,
domains for solving the governing equations. We examine the properties of the
method, by considering a one-dimensional Poisson equation with different
Neumann boundary conditions. The penalized Laplace operator is discretized by
second order central finite-differences and interpolation. The discretization
and penalization errors are thus assessed for several test problems.
Convergence properties of the discretized operator and the solution of the
penalized equation are analyzed. The generalized method is then applied to an
advection-diffusion equation coupled with the Navier-Stokes equations in an
annular domain which is immersed in a square domain. The application is
verified by numerical simulation of steady free convection in a concentric
annulus heated through the inner cylinder surface using an extended square
domain.Comment: 32 pages, 19 figure
Spectral methods for CFD
One of the objectives of these notes is to provide a basic introduction to spectral methods with a particular emphasis on applications to computational fluid dynamics. Another objective is to summarize some of the most important developments in spectral methods in the last two years. The fundamentals of spectral methods for simple problems will be covered in depth, and the essential elements of several fluid dynamical applications will be sketched
A matrix-free high-order discontinuous Galerkin compressible Navier-Stokes solver: A performance comparison of compressible and incompressible formulations for turbulent incompressible flows
Both compressible and incompressible Navier-Stokes solvers can be used and
are used to solve incompressible turbulent flow problems. In the compressible
case, the Mach number is then considered as a solver parameter that is set to a
small value, , in order to mimic incompressible flows.
This strategy is widely used for high-order discontinuous Galerkin
discretizations of the compressible Navier-Stokes equations. The present work
raises the question regarding the computational efficiency of compressible DG
solvers as compared to a genuinely incompressible formulation. Our
contributions to the state-of-the-art are twofold: Firstly, we present a
high-performance discontinuous Galerkin solver for the compressible
Navier-Stokes equations based on a highly efficient matrix-free implementation
that targets modern cache-based multicore architectures. The performance
results presented in this work focus on the node-level performance and our
results suggest that there is great potential for further performance
improvements for current state-of-the-art discontinuous Galerkin
implementations of the compressible Navier-Stokes equations. Secondly, this
compressible Navier-Stokes solver is put into perspective by comparing it to an
incompressible DG solver that uses the same matrix-free implementation. We
discuss algorithmic differences between both solution strategies and present an
in-depth numerical investigation of the performance. The considered benchmark
test cases are the three-dimensional Taylor-Green vortex problem as a
representative of transitional flows and the turbulent channel flow problem as
a representative of wall-bounded turbulent flows
Immersed Boundary Smooth Extension: A high-order method for solving PDE on arbitrary smooth domains using Fourier spectral methods
The Immersed Boundary method is a simple, efficient, and robust numerical
scheme for solving PDE in general domains, yet it only achieves first-order
spatial accuracy near embedded boundaries. In this paper, we introduce a new
high-order numerical method which we call the Immersed Boundary Smooth
Extension (IBSE) method. The IBSE method achieves high-order accuracy by
smoothly extending the unknown solution of the PDE from a given smooth domain
to a larger computational domain, enabling the use of simple Cartesian-grid
discretizations (e.g. Fourier spectral methods). The method preserves much of
the flexibility and robustness of the original IB method. In particular, it
requires minimal geometric information to describe the boundary and relies only
on convolution with regularized delta-functions to communicate information
between the computational grid and the boundary. We present a fast algorithm
for solving elliptic equations, which forms the basis for simple, high-order
implicit-time methods for parabolic PDE and implicit-explicit methods for
related nonlinear PDE. We apply the IBSE method to solve the Poisson, heat,
Burgers', and Fitzhugh-Nagumo equations, and demonstrate fourth-order pointwise
convergence for Dirichlet problems and third-order pointwise convergence for
Neumann problems
Spectral collocation methods
This review covers the theory and application of spectral collocation methods. Section 1 describes the fundamentals, and summarizes results pertaining to spectral approximations of functions. Some stability and convergence results are presented for simple elliptic, parabolic, and hyperbolic equations. Applications of these methods to fluid dynamics problems are discussed in Section 2
- …