Academic literature on the topic 'Raviart-Thomas vector space'

AbstractIdentities that relate projections of Raviart–Thomas finite element vector fields to discrete gradients of Crouzeix–Raviart finite element functions are derived under general conditions. Various implications such as discrete convex duality results and a characterization of the image of the projection of the Crouzeix–Ravaiart space onto elementwise constant functions are deduced.

AbstractA general approach was proposed in this article to develop high-order exponentially fitted basis functions for finite element approximations of multi-dimensional drift-diffusion equations for modeling biomolecular electrodiffusion processes. Such methods are highly desirable for achieving numerical stability and efficiency. We found that by utilizing the one-to-one correspondence between the continuous piecewise polynomial space of degree k + 1 and the divergencefree vector space of degree k, one can construct high-order two-dimensional exponentially fitted basis functions that are strictly interpolative at a selected node set but are discontinuous on edges in general, spanning nonconforming finite element spaces. First order convergence was proved for the methods constructed from divergence-free Raviart-Thomas space RT 00 at two different node sets.

AbstractWe combine theoretical results from polytope domain meshing, generalized barycentric coordinates, and finite element exterior calculus to construct scalar- and vector-valued basis functions for conforming finite element methods on generic convex polytope meshes in dimensions 2 and 3. Our construction recovers well-known bases for the lowest order Nédélec, Raviart–Thomas, and Brezzi–Douglas–Marini elements on simplicial meshes and generalizes the notion of Whitney forms to non-simplicial convex polygons and polyhedra. We show that our basis functions lie in the correct function space with regards to global continuity and that they reproduce the requisite polynomial differential forms described by finite element exterior calculus. We present a method to count the number of basis functions required to ensure these two key properties.

We introduce and analyse an augmented mixed variational formulation for the coupling of the Stokes and heat equations. More precisely, the underlying model consists of the Stokes equation suggested by the Oldroyd model for viscoelastic flow, coupled with the heat equation through a temperature-dependent viscosity of the fluid and a convective term. The original unknowns are the polymeric part of the extra-stress tensor, the velocity, the pressure, and the temperature of the fluid. In turn, for convenience of the analysis, the strain tensor, the vorticity, and an auxiliary symmetric tensor are introduced as further unknowns. This allows to join the polymeric and solvent viscosities in an adimensional viscosity, and to eliminate the polymeric part of the extra-stress tensor and the pressure from the system, which, together with the solvent part of the extra-stress tensor, are easily recovered later on through suitable postprocessing formulae. In this way, a fully mixed approach is applied, in which the heat flux vector is incorporated as an additional unknown as well. Furthermore, since the convective term in the heat equation forces both the velocity and the temperature to live in a smaller space than usual, we augment the variational formulation by using the constitutive and equilibrium equations, the relation defining the strain and vorticity tensors, and the Dirichlet boundary condition on the temperature. The resulting augmented scheme is then written equivalently as a fixed-point equation, so that the well-known Schauder and Banach theorems, combined with the Lax-Milgram theorem and certain regularity assumptions, are applied to prove the unique solvability of the continuous system. As for the associated Galerkin scheme, whose solvability is established similarly to the continuous case by using the Brouwer fixed-point and Lax–Milgram theorems, we employ Raviart–Thomas approximations of order k for the stress tensor and the heat flux vector, continuous piecewise polynomials of order ≤ k + 1 for velocity and temperature, and piecewise polynomials of order ≤ k for the strain tensor and the vorticity. Finally, we derive optimal a priori error estimates and provide several numerical results illustrating the good performance of the scheme and confirming the theoretical rates of convergence.

AbstractIn this paper, we introduce and analyze a high-order, fully-mixed finite element method for the free convection of n-dimensional fluids, {n\in\{2,3\}}, with temperature-dependent viscosity and thermal conductivity. The mathematical model is given by the coupling of the equations of continuity, momentum (Navier–Stokes) and energy by means of the Boussinesq approximation, as well as mixed thermal boundary conditions and a Dirichlet condition on the velocity. Because of the dependence on the temperature of the fluid properties, several additional variables are defined, thus resulting in an augmented formulation that seeks the rate of strain, pseudostress and vorticity tensors, velocity, temperature gradient and pseudoheat vectors, and temperature of the fluid. Using a fixed-point approach, smallness-of-data assumptions and a slight higher-regularity assumption for the exact solution provide the necessary well-posedness results at both continuous and discrete levels. In addition, and as a result of the augmentation, no discrete inf-sup conditions are needed for the well-posedness of the Galerkin scheme, which provides freedom of choice with respect to the finite element spaces. In particular, we suggest a combination based on Raviart–Thomas, Lagrange and discontinuous elements for which we derive optimal a priori error estimates. Finally, several numerical examples illustrating the performance of the method and confirming the theoretical rates of convergence are reported.

AbstractA general setting for a new approach to constructing vector splitting schemes in mixed FEM for heat transfer problem is considered. For space approximation Raviart–Thomas finite elements of lowest order are implemented on rectangular (parallelepiped) mesh. The main question discussed is how to implement time discretization so as to obtain efficient numerical algorithms.The key idea of the proposed approach is to use scalar splitting schemes for heat flux divergence. This allows one to carry out accuracy and stability analysis on the basis of the well-known results for underlying scalar splitting schemes.

AbstractIn this paper, we present a novel arbitrary-order discrete de Rham (DDR) complex on general polyhedral meshes based on the decomposition of polynomial spaces into ranges of vector calculus operators and complements linked to the spaces in the Koszul complex. The DDR complex is fully discrete, meaning that both the spaces and discrete calculus operators are replaced by discrete counterparts, and satisfies suitable exactness properties depending on the topology of the domain. In conjunction with bespoke discrete counterparts of $$\text {L}^2$$ L 2 -products, it can be used to design schemes for partial differential equations that benefit from the exactness of the sequence but, unlike classical (e.g., Raviart–Thomas–Nédélec) finite elements, are nonconforming. We prove a complete panel of results for the analysis of such schemes: exactness properties, uniform Poincaré inequalities, as well as primal and adjoint consistency. We also show how this DDR complex enables the design of a numerical scheme for a magnetostatics problem, and use the aforementioned results to prove stability and optimal error estimates for this scheme.

Dans cette thèse, une nouvelle approche est développée pour les problèmes de viscoplasticité et de dynamique non linéaire. En particulier, les équations variationnelles sont élaborées selon le principe de Helligner-Reissner, de sorte que les champs de contrainte et de déplacement apparaissent comme des champs inconnus sous la forme faible. Trois nouveaux éléments finis sont développés. Le premier élément fini est formulé pour le problème axisymétrique, dans lequel le champ de contraintes est approximé par des polynômes d’ordre inférieur tels que des fonctions linéaires. Cette approche donne des solutions précises spécifiquement dans les problèmes incompressibles et rigides. De plus, un élément fini de flexion de membrane et de plaque est nouvellement conçu en discrétisant le champ de contraintes en utilisant l’espace vectoriel de Raviart-Thomas d’ordre le plus bas RT0. Cette approche garantit la continuité du champ de contraintes sur tout un domaine discret, ce qui est un avantage significatif dans la méthode numérique, notamment pour les problèmes de propagation des ondes. Les développements sont effectués pour le comportement constitutif visco-plastique des matériaux, où les équations d’évolution correspondantes sont obtenues en faisant appel au principe de dissipation maximale. Pour résoudre les équations d’équilibre dynamique, des schémas de conservation et de décroissance de l’énergie sont formulés en conséquence. Le schéma de conservation de l’énergie est inconditionnellement stable, car il peut préserver l’énergie totale d’un système donné sous une vibration libre, tandis que le schéma décroissant peut dissiper des modes de vibration à plus haute fréquence. La dernière partie de cette thèse présente les procédures d’upscaling du comportement des matériaux visco-plastiques. Plus précisément, la mise à l’échelle est effectuée par une méthode d’identification stochastique via une mise à jour baysienne en utilisant le filtre de Gauss-Markov-Kalman pour l’assimilation des propriétés importantes des matériaux dans les régimes élastique et inélastique
In this thesis, a novel approach is developed for visco-plasticity and nonlinear dynamics problems. In particular, variational equations are elaborated following the Helligner-Reissner principle, so that both stress and displacement fields appear as unknown fields in the weak form. Three novel finite elements are developed. The first finite element is formulated for the axisymmetric problem, in which the stress field is approximated by low-order polynomials such as linear functions. This approach yields accurate solutions specifically in incompressible and stiff problems. In addition, a membrane and plate bending finite element are newly designed by discretizing the stress field using the lowest order Raviart-Thomas vector space RT0. This approach guarantees the continuity of the stress field over an entire discrete domain, which is a significant advantage in the numerical method, especially for the wave propagation problems. The developments are carried out for the viscoplastic constitutive behavior of materials, where the corresponding evolution equations are obtained by appealing to the principle of maximum dissipation. To solve the dynamic equilibrium equations, energy conserving and decaying schemes are formulated correspondingly. The energy conserving scheme is unconditional stable, since it can preserve the total energy of a given system under a free vibration, while the decaying scheme can dissipate higher frequency vibration modes. The last part of this thesis presents procedures for upscaling of the visco-plastic material behavior. Specifically, the upscaling is performed by stochastic identification method via Baysian updating using the Gauss-Markov-Kalman filter for assimilation of important material properties in the elastic and inelastic regimes