The main scope of this work is to carry out a mathematical framework and its corresponding finite element (FE) discretization for the partially saturated soil consolidation modelling in presence of an immiscible pollutant. A multiphase system with the interstitial voids in the grain matrix...

In this paper we are concerned with the elliptic PDEs with highly oscillating coefficients which model the behaviour of composite linear elastic materials. Analytical expressions for the effective coefficients are obtained for the case of layered materials using the theory of homogenisation....

We consider a system of Maxwell’s and Landau-Lifshitz-Gilbert equations describing magnetization dynamics in micromagnetism. The problem is discretized by a convergent, unconditionally stable finite element method. A multigrid preconditioned Uzawa type method for the solution of the algebraic...

This paper deals with the development of efficient numerical solvers for the primitive equations of the ocean in turbulent regime. We derive the numerical approximation of a reduced model by the Smagorinsky turbulence model that includes stabilization of the pressure discretization by a penalty...

A method for constructing Whitney forms is proposed, which applies to tetrahedra, hexahedra, triangular prisms, and pyramids in a similar way, and proceeds from a unique generating principle, thus unifying their presentation. The principle automatically enforces conformity (i.e.,...

A parallel implementation of the Balancing Domain Decomposition based on Constraints (BDDC) method using the frontal solver is employed to solve systems of linear equations from finite element analysis, and incorporated into a standard finite element system for engineering analysis by linear...

A dyadic grid is a d-dimensional hierarchical mesh where a cell at level k is partitioned into two equal children at level k+1 by a hyperplane perpendicular to coordinate axis (kmodd). We consider here the finite element approach on adaptive grids, static and dynamic, for various functional...

The mechanical behavior of magnetic reading devices is mainly governed by compressible Reynolds equations when the air bearing modeling approximation is considered. First, the convection dominated feature motivates the use of a characteristics scheme adapted to steady state problems. Secondly, a...

In this paper the problem of blending parametric surfaces is discussed. We present a constructing method of blending surfaces by parametric discrete interpolation PDE splines. These functions are obtained from some boundary conditions and a given interpolation data point set, by minimizing a...

Two classes of nodal methods, weakly and strongly discontinuous, are introduced and applied to the numerical solution of the neutron transport equation in two-dimensional Cartesian geometry and discrete ordinates. These methods are then applied for the approximation of the solution of a...