A New Characteristic Nonconforming Mixed Finite Element Scheme for Convection-Dominated Diffusion Problem

A characteristic nonconforming mixed finite element method (MFEM) is proposed for the convection-dominated diffusion problem based on a new mixed variational formulation. The optimal order error estimates for both the original variable u and the auxiliary variable σ with respect to the space are obtained by employing some typical characters of the interpolation operator instead of the mixed (or expanded mixed) elliptic projection which is an indispensable tool in the traditional MFEM analysis. At last, we give some numerical results to confirm the theoretical analysis

Unified convergence analysis of numerical schemes for a miscible displacement problem

This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the framework of the gradient discretisation method for diffusion operators on generic grids. We use it to establish a novel convergence result in $L^\infty(0,T; L^2(\Omega))$ of the approximate concentration using minimal regularity assumptions on the solution to the continuous problem. The convection term in the concentration equation is discretised using a centred scheme. We present a variety of numerical tests from the literature, as well as a novel analytical test case. The performance of two schemes are compared on these tests; both are poor in the case of variable viscosity, small diffusion and medium to small time steps. We show that upstreaming is not a good option to recover stable and accurate solutions, and we propose a correction to recover stable and accurate schemes for all time steps and all ranges of diffusion

Numerical simulation of incompressible fluid flow by the spectral element method

Tato diplomová práce prezentuje metodu spektrálních prvků. Tato metoda je použita k řešení stacionárního 2-D laminárního proudění Newtonovské nestlačitelné tekutiny. Proudění je popsáno stacionarní Navier-Stokesovou rovnicí. Dohromady s okrajovou pod- mínkou tvoří Navier-Stokesův problém. Na slabou formulaci této úlohy je aplikována metoda spektrálních prvků. Touto discretizací se získá soustava nelineárních rovnic. K obrdžení lineární soustavy je použita Newtonova iterační metoda. Podorobný algorit- mus tvoří jádro Navier-Stokeseva solveru, který je naprogramován v Matlabu. Na závěr jsou pomocí tohoto solveru řešeny dva příklady: proudění v kavitě a obtékání válce. Přík- lady jsou řešeny pro různé Reynoldsovy čísla. První od 1 do 1000 a druhý od 1 do 100.The thesis presents the spectral element method and its application to a steady 2-D laminar flow of an incompressible Newtonian fluid. Main features of this method are presented in the thesis. The flow is governed by the steady Navier-Stokes equation. Together with boundary data they form the steady Navier-Stokes problem. Its weak form is a starting point for the method. A space discretization is applied and it results into a nonlinear system of equations. Due to this, the nonlinearity has to be treated. To obtain a linear system of equations is the Newton iteration method used. This algorithm forms the kernel of a Navier-Stokes solver that is implemented in Matlab. Finally, there are presented two examples: the lid driven cavity flow and the flow over a cylinder. The first one is solved for Reynolds numbers from 1 to 1000 and the second one for Reynolds numbers from 1 to 100.

On the convergence order of the finite element error in the kinetic energy for high Reynolds number incompressible flows

The kinetic energy of a flow is proportional to the square of the norm of the velocity. Given a sufficient regular velocity field and a velocity finite element space with polynomials of degree , then the best approximation error in is of order . In this survey, the available finite element error analysis for the velocity error in is reviewed, where is a final time. Since in practice the case of small viscosity coefficients or dominant convection is of particular interest, which may result in turbulent flows, robust error estimates are considered, i.e., estimates where the constant in the error bound does not depend on inverse powers of the viscosity coefficient. Methods for which robust estimates can be derived enable stable flow simulations for small viscosity coefficients on comparatively coarse grids, which is often the situation encountered in practice. To introduce stabilization techniques for the convection-dominated regime and tools used in the error analysis, evolutionary linear convection–diffusion equations are studied at the beginning. The main part of this survey considers robust finite element methods for the incompressible Navier–Stokes equations of order , , and for the velocity error in . All these methods are discussed in detail. In particular, a sketch of the proof for the error bound is given that explains the estimate of important terms which determine finally the order of convergence. Among them, there are methods for inf–sup stable pairs of finite element spaces as well as for pressure-stabilized discretizations. Numerical studies support the analytic results for several of these methods. In addition, methods are surveyed that behave in a robust way but for which only a non-robust error analysis is available. The conclusion of this survey is that the problem of whether or not there is a robust method with optimal convergence order for the kinetic energy is still open

Theoretical and numerical studies of chaotic mixing

Theoretical and numerical studies of chaotic mixing are performed to circumvent the difficulties of efficient mixing, which come from the lack of turbulence in microfluidic devices. In order to carry out efficient and accurate parametric studies and to identify a fully chaotic state, a spectral element algorithm for solution of the incompressible Navier-Stokes and species transport equations is developed. Using Taylor series expansions in time marching, the new algorithm employs an algebraic factorization scheme on multi-dimensional staggered spectral element grids, and extends classical conforming Galerkin formulations to nonconforming spectral elements. Lagrangian particle tracking methods are utilized to study particle dispersion in the mixing device using spectral element and fourth order Runge-Kutta discretizations in space and time, respectively. Comparative studies of five different techniques commonly employed to identify the chaotic strength and mixing efficiency in microfluidic systems are presented to demonstrate the competitive advantages and shortcomings of each method. These are the stirring index based on the box counting method, Poincare sections, finite time Lyapunov exponents, the probability density function of the stretching field, and mixing index inverse, based on the standard deviation of scalar species distribution. Series of numerical simulations are performed by varying the Peclet number (Pe) at fixed kinematic conditions. The mixing length (lm) is characterized as function of the Pe number, and lm ∝ ln(Pe) scaling is demonstrated for fully chaotic cases. Employing the aforementioned techniques, optimum kinematic conditions and the actuation frequency of the stirrer that result in the highest mixing/stirring efficiency are identified in a zeta potential patterned straight micro channel, where a continuous flow is generated by superposition of a steady pressure driven flow and time periodic electroosmotic flow induced by a stream-wise AC electric field. Finally, it is shown that the invariant manifold of hyperbolic periodic point determines the geometry of fast mixing zones in oscillatory flows in two-dimensional cavity

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

Schnelle Löser für partielle Differentialgleichungen

[no abstract available