Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNumerical analysis of a new formulation for the Oseen equations in terms of vorticity and Bernoulli pressure
56
0
0.0
(
0
)
Added by
Ricardo Ruiz Baier
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
A variational formulation is introduced for the Oseen equations written in terms of vor-ti-city and Bernoulli pressure. The velocity is fully decoupled using the momentum balance equation, and it is later recovered by a post-process. A finite element method is also proposed, consisting in equal-order Nedelec finite elements and piecewise continuous polynomials for the vorticity and the Bernoulli pressure, respectively. The {it a priori} error analysis is carried out in the $mathrm{L}^2$-norm for vorticity, pressure, and velocity; under a smallness assumption either on the convecting velocity, or on the mesh parameter. Furthermore, an {it a posteriori} error estimator is designed and its robustness and efficiency are studied using weighted norms. Finally, a set of numerical examples in 2D and 3D is given, where the error indicator serves to guide adaptive mesh refinement. These tests illustrate the behaviour of the new formulation in typical flow conditions, and they also confirm the theoretical findings.
rate research
Read More
We introduce a family of mixed methods and discontinuous Galerkin discretisations designed to numerically solve the Oseen equations written in terms of velocity, vorticity, and Bernoulli pressure. The unique solvability of the continuous problem is addressed by invoking a global inf-sup property in an adequate abstract setting for non-symmetric systems. The proposed finite element schemes, which produce exactly divergence-free discrete velocities, are shown to be well-defined and optimal convergence rates are derived in suitable norms. In addition, we establish optimal rates of convergence for a class of discontinuous Galerkin schemes, which employ stabilisation. A set of numerical examples serves to illustrate salient features of these methods.
We propose and analyse an augmented mixed finite element method for the Oseen equations written in terms of velocity, vorticity, and pressure with non-constant viscosity and homogeneous Dirichlet boundary condition for the velocity. The weak formulation includes least-squares terms arising from the constitutive equation and from the incompressibility condition, and we show that it satisfies the hypotheses of the Babuvska-Brezzi theory. Repeating the arguments of the continuous analysis, the stability and solvability of the discrete problem are established. The method is suited for any Stokes inf-sup stable finite element pair for velocity and pressure, while for vorticity any generic discrete space (of arbitrary order) can be used. A priori and a posteriori error estimates are derived using two specific families of discrete subspaces. Finally, we provide a set of numerical tests illustrating the behaviour of the scheme, verifying the theoretical convergence rates, and showing the performance of the adaptive algorithm guided by residual a posteriori error estimation.
We study a continuous data assimilation (CDA) algorithm for a velocity-vorticity formulation of the 2D Navier-Stokes equations in two cases: nudging applied to the velocity and vorticity, and nudging applied to the velocity only. We prove that under a typical finite element spatial discretization and backward Euler temporal discretization, application of CDA preserves the unconditional long-time stability property of the velocity-vorticity method and provides optimal long-time accuracy. These properties hold if nudging is applied only to the velocity, and if nudging is also applied to the vorticity then the optimal long-time accuracy is achieved more rapidly in time. Numerical tests illustrate the theory, and show its effectiveness on an application problem of channel flow past a flat plate.
Infinite-dimensional Newton methods can be effectively used to derive numerical proofs of the existence of solutions to partial differential equations (PDEs). In computer-assisted proofs of PDEs, the original problem is transformed into the infinite Newton-type fixed point equation $w = - {mathcal L}^{-1} {mathcal F}(hat{u}) + {mathcal L}^{-1} {mathcal G}(w)$, where ${mathcal L}$ is a linearized operator, ${mathcal F}(hat{u})$ is a residual, and ${mathcal G}(w)$ is a local Lipschitz term. Therefore, the estimations of $| {mathcal L}^{-1} {mathcal F}(hat{u}) |$ and $| {mathcal L}^{-1}{mathcal G}(w) |$ play major roles in the verification procedures. In this paper, using a similar concept as the `Schur complement for matrix problems, we represent the inverse operator ${mathcal L}^{-1}$ as an infinite-dimensional operator matrix that can be decomposed into two parts, one finite dimensional and one infinite dimensional. This operator matrix yields a new effective realization of the infinite-dimensional Newton method, enabling a more efficient verification procedure compared with existing methods for the solution of elliptic PDEs. We present some numerical examples that confirm the usefulness of the proposed method. Related results obtained from the representation of the operator matrix as ${mathcal L}^{-1}$ are presented in the appendix.
We introduce a high-order numerical scheme for fractional ordinary differential equations with the Caputo derivative. The method is developed by dividing the domain into a number of subintervals, and applying the quadratic interpolation on each subinterval. The method is shown to be unconditionally stable, and for general nonlinear equations, the uniform sharp numerical order $3- u$ can be rigorously proven for sufficiently smooth solutions at all time steps. The proof provides a general guide for proving the sharp order for higher-order schemes in the nonlinear case. Some numerical examples are given to validate our theoretical results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
