459 research outputs found
New discretization and solution techniques for incompressible viscous flow problems
Several topics arising in the finite element solution of the incompressible Navier-Stokes equations are considered. Specifically, the question of choosing finite element velocity/pressure spaces is addressed, particularly from the viewpoint of achieving stable discretizations leading to convergent pressure approximations. The role of artificial viscosity in viscous flow calculations is studied, emphasizing work by several researchers for the anisotropic case. The last section treats the problem of solving the nonlinear systems of equations which arise from the discretization. Time marching methods and classical iterative techniques, as well as some modifications are mentioned
Proceedings for the ICASE Workshop on Heterogeneous Boundary Conditions
Domain Decomposition is a complex problem with many interesting aspects. The choice of decomposition can be made based on many different criteria, and the choice of interface of internal boundary conditions are numerous. The various regions under study may have different dynamical balances, indicating that different physical processes are dominating the flow in these regions. This conference was called in recognition of the need to more clearly define the nature of these complex problems. This proceedings is a collection of the presentations and the discussion groups
A split finite element algorithm for the compressible Navier-Stokes equations
An accurate and efficient numerical solution algorithm is established for solution of the high Reynolds number limit of the Navier-Stokes equations governing the multidimensional flow of a compressible essentially inviscid fluid. Finite element interpolation theory is used within a dissipative formulation established using Galerkin criteria within the Method of Weighted Residuals. An implicit iterative solution algorithm is developed, employing tensor product bases within a fractional steps integration procedure, that significantly enhances solution economy concurrent with sharply reduced computer hardware demands. The algorithm is evaluated for resolution of steep field gradients and coarse grid accuracy using both linear and quadratic tensor product interpolation bases. Numerical solutions for linear and nonlinear, one, two and three dimensional examples confirm and extend the linearized theoretical analyses, and results are compared to competitive finite difference derived algorithms
On the evolution of flow topology in turbulent Rayleigh-Bénard convection
Copyright 2016 AIP Publishing. This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing.Small-scale dynamics is the spirit of turbulence physics. It implicates many attributes of flow topology evolution, coherent structures, hairpin vorticity dynamics, and mechanism of the kinetic energy cascade. In this work, several dynamical aspects of the small-scale motions have been numerically studied in a framework of Rayleigh-Benard convection (RBC). To do so, direct numerical simulations have been carried out at two Rayleigh numbers Ra = 10(8) and 10(10), inside an air-filled rectangular cell of aspect ratio unity and pi span-wise open-ended distance. As a main feature, the average rate of the invariants of the velocity gradient tensor (Q(G), R-G) has displayed the so-calledPeer ReviewedPostprint (author's final draft
Adaptive rational Krylov methods for exponential Runge--Kutta integrators
We consider the solution of large stiff systems of ordinary differential
equations with explicit exponential Runge--Kutta integrators. These problems
arise from semi-discretized semi-linear parabolic partial differential
equations on continuous domains or on inherently discrete graph domains. A
series of results reduces the requirement of computing linear combinations of
-functions in exponential integrators to the approximation of the
action of a smaller number of matrix exponentials on certain vectors.
State-of-the-art computational methods use polynomial Krylov subspaces of
adaptive size for this task. They have the drawback that the required Krylov
subspace iteration numbers to obtain a desired tolerance increase drastically
with the spectral radius of the discrete linear differential operator, e.g.,
the problem size. We present an approach that leverages rational Krylov
subspace methods promising superior approximation qualities. We prove a novel
a-posteriori error estimate of rational Krylov approximations to the action of
the matrix exponential on vectors for single time points, which allows for an
adaptive approach similar to existing polynomial Krylov techniques. We discuss
pole selection and the efficient solution of the arising sequences of shifted
linear systems by direct and preconditioned iterative solvers. Numerical
experiments show that our method outperforms the state of the art for
sufficiently large spectral radii of the discrete linear differential
operators. The key to this are approximately constant rational Krylov iteration
numbers, which enable a near-linear scaling of the runtime with respect to the
problem size
A simple method for detecting chaos in nature
Chaos, or exponential sensitivity to small perturbations, appears everywhere
in nature. Moreover, chaos is predicted to play diverse functional roles in
living systems. A method for detecting chaos from empirical measurements should
therefore be a key component of the biologist's toolkit. But, classic
chaos-detection tools are highly sensitive to measurement noise and break down
for common edge cases, making it difficult to detect chaos in domains, like
biology, where measurements are noisy. However, newer tools promise to overcome
these limitations. Here, we combine several such tools into an automated
processing pipeline, and show that our pipeline can detect the presence (or
absence) of chaos in noisy recordings, even for difficult edge cases. As a
first-pass application of our pipeline, we show that heart rate variability is
not chaotic as some have proposed, and instead reflects a stochastic process in
both health and disease. Our tool is easy-to-use and freely available
A dynamically adaptive multigrid algorithm for the incompressible Navier-Stokes equations: Validation and model problems
An algorithm is described for the solution of the laminar, incompressible Navier-Stokes equations. The basic algorithm is a multigrid based on a robust, box-based smoothing step. Its most important feature is the incorporation of automatic, dynamic mesh refinement. This algorithm supports generalized simple domains. The program is based on a standard staggered-grid formulation of the Navier-Stokes equations for robustness and efficiency. Special grid transfer operators were introduced at grid interfaces in the multigrid algorithm to ensure discrete mass conservation. Results are presented for three models: the driven-cavity, a backward-facing step, and a sudden expansion/contraction
Efficient "black-box" multigrid solvers for convection-dominated problems
The main objective of this project is to develop a "black-box" multigrid preconditioner for the iterative solution of finite element discretisations of the convection-diffusion equation with dominant convection. This equation can be considered a stand alone scalar problem or as part of a more complex system of partial differential equations, such as the Navier-Stokes equations. The project will focus on the stand alone scalar problem. Multigrid is considered an optimal preconditioner for scalar elliptic problems. This strategy can also be used for convection-diffusion problems, however an appropriate robust smoother needs to be developed to achieve mesh-independent convergence. The focus of the thesis is on the development of such a smoother. In this context a novel smoother is developed referred to as truncated incomplete factorisation (tILU) smoother. In terms of computational complexity and memory requirements, the smoother is considerably less expensive than the standard ILU(0) smoother. At the same time, it exhibits the same robustness as ILU(0) with respect to the problem and discretisation parameters. The new smoother significantly outperforms the standard damped Jacobi smoother and is a competitor to the Gauss-Seidel smoother (and in a number of important cases tILU outperforms the Gauss-Seidel smoother). The new smoother depends on a single parameter (the truncation ratio). The project obtains a default value for this parameter and demonstrated the robust performance of the smoother on a broad range of problems. Therefore, the new smoothing method can be regarded as "black-box". Furthermore, the new smoother does not require any particular ordering of the nodes, which is a prerequisite for many robust smoothers developed for convection-dominated convection-diffusion problems. To test the effectiveness of the preconditioning methodology, we consider a number of model problems (in both 2D and 3D) including uniform and complex (recirculating) convection fields discretised by uniform, stretched and adaptively refined grids. The new multigrid preconditioner within block preconditioning of the Navier-Stokes equations was also tested. The numerical results gained during the investigation confirm that tILU is a scalable, robust smoother for both geometric and algebraic multigrid. Also, comprehensive tests show that the tILU smoother is a competitive method.EThOS - Electronic Theses Online ServiceEPSRCGBUnited Kingdo
- …