Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userContinuous data assimilation applied to a velocity-vorticity formulation of the 2D Navier-Stokes equations
63
0
0.0
(
0
)
Added by
Duygu Vargun
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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 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.
In two dimensions, we propose and analyze an a posteriori error estimator for finite element approximations of the stationary Navier Stokes equations with singular sources on Lipschitz, but not necessarily convex, polygonal domains. Under a smallness assumption on the continuous and discrete solutions, we prove that the devised error estimator is reliable and locally efficient. We illustrate the theory with numerical examples.
We investigate the spatio-temporal structure of the most likely configurations realising extremely high vorticity or strain in the stochastically forced 3D incompressible Navier-Stokes equations. Most likely configurations are computed by numerically finding the highest probability velocity field realising an extreme constraint as solution of a large optimisation problem. High-vorticity configurations are identified as pinched vortex filaments with swirl, while high-strain configurations correspond to counter-rotating vortex rings. We additionally observe that the most likely configurations for vorticity and strain spontaneously break their rotational symmetry for extremely high observable values. Instanton calculus and large deviation theory allow us to show that these maximum likelihood realisations determine the tail probabilities of the observed quantities. In particular, we are able to demonstrate that artificially enforcing rotational symmetry for large strain configurations leads to a severe underestimate of their probability, as it is dominated in likelihood by an exponentially more likely symmetry broken vortex-sheet configuration.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
