45 research outputs found
Analysis of the Brinkman-Forchheimer equations with slip boundary conditions
In this work, we study the Brinkman-Forchheimer equations driven under slip
boundary conditions of friction type. We prove the existence and uniqueness of
weak solutions by means of regularization combined with the Faedo-Galerkin
approach. Next we discuss the continuity of the solution with respect to
Brinkman's and Forchheimer's coefficients. Finally, we show that the weak
solution of the corresponding stationary problem is stable
Error estimates for the discretization of the velocity tracking problem
In this paper we are continuing our work (Casas and Chrysafinos, SIAM J Numer Anal 50(5):2281–2306, 2012), concerning a priori error estimates for the velocity tracking of two-dimensional evolutionary Navier–Stokes flows. The controls are of distributed type, and subject to point-wise control constraints. The discretization scheme of the state and adjoint equations is based on a discontinuous time-stepping scheme (in time) combined with conforming finite elements (in space) for the velocity and pressure. Provided that the time and space discretization parameters, t and h respectively, satisfy t = Ch2, error estimates of order O(h2) and O(h 3/2 – 2/p ) with p > 3 depending on the regularity of the target and the initial velocity, are proved for the difference between the locally optimal controls and their discrete approximations, when the controls are discretized by the variational discretization approach and by using piecewise-linear functions in space respectively. Both results are based on new duality arguments for the evolutionary Navier–Stokes equations
HDG-NEFEM with Degree Adaptivity for Stokes Flows
This paper presents the first degree adaptive procedure able to directly use the geometry given by a CAD model. The technique uses a hybridisable discontinuous Galerkin discretisation combined with a NURBS-enhanced rationale, completely removing the uncertainty induced by a polynomial approximation of curved boundaries that is common within an isoparametric approach. The technique is compared against two strategies to perform degree adaptivity currently in use. This paper demonstrates, for the first time, that the most extended technique for degree adaptivity can easily lead to a non-reliable error estimator if no communication with CAD software is introduced whereas if the communication with the CAD is done, it results in a substantial computing time. The proposed technique encapsulates the CAD model in the simulation and is able to produce reliable error estimators irrespectively of the initial mesh used to start the adaptive process. Several numerical examples confirm the findings and demonstrate the superiority of the proposed technique. The paper also proposes a novel idea to test the implementation of high-order solvers where different degrees of approximation are used in different elements
Systems of Differential Algebraic Equations in Computational Electromagnetics
Starting from space-discretisation of Maxwell's equations, various classical
formulations are proposed for the simulation of electromagnetic fields. They
differ in the phenomena considered as well as in the variables chosen for
discretisation. This contribution presents a literature survey of the most
common approximations and formulations with a focus on their structural
properties. The differential-algebraic character is discussed and quantified by
the differential index concept