Search

We derive a new kind of boundary conditions for the vorticity equation with solid wall boundaries for fluid flow problems. The formulation uses a Dirichlet condition for the normal component of vorticity and Neumann type conditions for the tangential components. In a Galerkin (integral) formulation the tangential condition is natural, i.e., it is enforced by a right-hand side functional and does not impose a boundary constraint on trial and test spaces. The functional involves the pressure variable, and we discuss several velocity–vorticity formulations where the proposed condition is appropriate. Several numerical experiments are given that illustrate the validity of the new approach. [ABSTRACT FROM AUTHOR]

In this paper, we consider phase-field based fracture propagation in elastic media. The main purpose is the development of a robust and efficient numerical scheme. To enforce crack irreversibility as a constraint, we use a primal-dual active set strategy, which can be identified as a semi-smooth Newton method. The active set iteration is merged with the Newton iteration for solving the fully-coupled nonlinear partial differential equation discretized using finite elements, resulting in a single, rapidly converging nonlinear scheme. It is well known that phase-field models require fine meshes to accurately capture the propagation dynamics of the crack. Because traditional estimators based on adaptive mesh refinement schemes are not appropriate, we develop a predictor-corrector scheme for local mesh adaptivity to reduce the computational cost. This method is both robust and efficient and allows us to treat temporal and spatial refinements and to study the influence of model regularization parameters. Finally, our proposed approach is substantiated with different numerical tests for crack propagation in elastic media and pressurized fracture propagation in homogeneous and heterogeneous media. [ABSTRACT FROM AUTHOR]

Abstract: In this paper the grad–div stabilization for the incompressible Navier–Stokes finite element approximations is considered from two different viewpoints: (i) as a variational multiscale approach for the pressure subgrid modeling and (ii) as a stabilization procedure of least-square type. Some new error estimates for the linearized problem with the grad–div stabilization are proved with the help of norms induced by the pressure Schur complement operator. We discuss the stabilization parameter choice arising in the frameworks of least-square and multiscale methods and consider assumptions which allow to relate both approaches. [Copyright &y& Elsevier]