Search

Abstract In this paper, we propose and analyze the stability and the dissipative structure of a new dynamic term-by-term stabilized finite element formulation for the Navier–Stokes problem that can be viewed as a variational multiscale (VMS) method under some assumptions. The essential point of the formulation is the time dependent nature of the subscales and, contrary to residual-based formulations, the introduction of two velocity subscale components. They represent the components of the convective and the pressure gradient terms, respectively, of the momentum equation that cannot be captured by the finite element mesh. A key idea of the proposed method is that the convective subscale is close to a solenoidal field and the pressure gradient subscale is close to a potential field. The method ensures stability in anisotropic space–time discretizations, which is proved using numerical analysis for a linearized problem and demonstrated in classical numerical tests. The work includes a detailed description of the proposed formulation and several numerical examples that serve to justify our claims. Highlights • A new dynamic term-by-term stabilized finite element formulation is proposed. • The proposed method ensures stability in anisotropic space–time discretizations. • The dissipative structure of the method is analyzed and some stability results are presented. • Several numerical examples are presented to justify our claims. [ABSTRACT FROM AUTHOR]

The theories for thick plates and beams, namely Reissner–Mindlin's and Timoshenko's theories, are well known to suffer numerical locking when approximated using the standard Galerkin finite element method for small thicknesses. This occurs when the same interpolations are used for displacement and rotations, reason for which stabilization becomes necessary. To overcome this problem, a Variational Multiscale stabilization method is analyzed in this paper. In this framework, two different approaches are presented: the Algebraic Sub-Grid Scale formulation and the Orthogonal Sub-Grid Scale formulation. Stability and convergence is proved for both approaches, explaining why the latter performs much better. Although the numerical examples show that the Algebraic Sub-Grid Scale approach is in some cases able to overcome the numerical locking, it is highly sensitive to stabilization parameters and presents difficulties to converge optimally with respect to the element size in the L 2 norm. In this regard, the Orthogonal Sub-Grid Scale approach, which considers the space of the sub-grid scales to be orthogonal to the finite element space, is shown to be stable and optimally convergent independently of the thickness of the solid. The final formulation is similar to approaches developed previously, thus justifying them in the frame of the Variational Multiscale concept. • Variational Multiscale Stabilization for plates and beams Finite Elements. • Models designed for thick theories Reissner–Mindlin for plates Timoshenko for beams. • The OSGS formulation does not present shear locking when the thickness goes to zero. • A complete numerical analysis is performed for the OSGS formulation. • Numerical results confirm the theoretical findings. [ABSTRACT FROM AUTHOR]

In this article, we analyse several discontinuous Galerkin (DG) methods for the Stokes problem under minimal regularity on the solution. We assume that the velocity u belongs to and the pressure . First, we analyse standard DG methods assuming that the right-hand side f belongs to . A DG method that is well defined for f belonging to is then investigated. The methods under study include stabilized DG methods using equal-order spaces and inf–sup stable ones where the pressure space is one polynomial degree less than the velocity space. [ABSTRACT FROM AUTHOR]

In this paper we provide an error analysis of a fractional-step method for the numerical solution of the incompressible Navier–Stokes equations. Under mild regularity assumptions on the continuous solution, we obtain first order error estimates in the time step size, both for the intermediate and the end-of-step velocities of the method; we also give some error estimates for the pressure solution. [Copyright &y& Elsevier]

In this work we compare two apparently different stabilization procedures for the finite element approximation of the incompressible Navier–Stokes equations. The first is the characteristic-based split (CBS). It combines the characteristic Galerkin method to deal with convection dominated flows with a classical splitting technique, which in some cases allows us to use equal velocity–pressure interpolations. The second approach is the Galerkin-least-squares (GLS) method, in which a least-squares form of the element residual is added to the basic Galerkin equations. It is shown that both formulations display similar stabilization mechanisms, provided the stabilization parameter of the GLS method is identified with the time step of the CBS approach. This identification can be understood from a formal Fourier analysis of the linearized problem. Copyright © 2001 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

This work presents a conservative scheme for iteration-by-subdomain domain decomposition (DD) strategies applied to the finite element solution of flow problems. The DD algorithm is based on the iterative update of the boundary conditions on the interfaces between the subregions, the so-called transmission conditions. The transmission conditions involve the essential and natural boundary conditions of the weak form of the problem, and should ensure strong continuity of the velocity and weak continuity of the traction. As a first approach, the transmission conditions are interpolated using the classical Lagrange interpolation functions. Conservation problems might arise when two adjacent subdomains have a sensibly different mesh spacing. In order to conserve any desired quantity of interest, an interface constraining is introduced: continuity of the transmission conditions are constrained under a scalar conservation equation. An example of mass conservation illustrates the algorithm. Copyright © 2001 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

Abstract: This paper deals with the question of strain localization associated with materials which exhibit softening due to tensile straining. A standard local isotropic Rankine damage model with strain-softening is used as exemplary constitutive model. Both the irreducible and mixed forms of the problem are examined and stability and solvability conditions are discussed. Lack of uniqueness and convergence difficulties related to the strong material nonlinearities involved are also treated. From this analysis, the issue of local discretization error in the pre-localization regime is deemed as the main difficulty to be overcome in the discrete problem. Focus is placed on low order finite elements with continuous strain and displacement fields (triangular P1P1 and quadrilateral Q1Q1), although the presented approach is very general. Numerical examples show that the resulting procedure is remarkably robust: it does not require the use of auxiliary tracking techniques and the results obtained do not suffer from spurious mesh-bias dependence. [Copyright &y& Elsevier]

Abstract: This paper exploits the concept of stabilized finite element methods to formulate stable mixed stress/displacement and strain/displacement finite elements for the solution of nonlinear solid mechanics problems. The different assumptions and approximations used to derive the methods are exposed. The proposed procedure is very general, applicable to 2D and 3D problems. Implementation and computational aspects are also discussed, showing that a robust application of the proposed formulation is feasible. Numerical examples show that the results obtained compare favorably with those obtained with the corresponding irreducible formulation. [Copyright &y& Elsevier]