Cut finite element methods for coupled bulk–surface problems

We develop a cut finite element method for a second order elliptic coupled bulk-surface model problem. We prove a priori estimates for the energy and L2L2 norms of the error. Using stabilization terms we show that the resulting algebraic system of equations has a similar condition number as a standard fitted finite element method. Finally, we present a numerical example illustrating the accuracy and the robustness of our approach

A Stabilized Cut Streamline Diffusion Finite Element Method for Convection-Diffusion Problems on Surfaces

We develop a stabilized cut finite element method for the stationary convection diffusion problem on a surface embedded in R d . The cut finite element method is based on using an embedding of the surface into a three dimensional mesh consisting of tetrahedra and then using the restriction of the standard piecewise linear continuous elements to a piecewise linear approximation of the surface. The stabilization consists of a standard streamline diffusion stabilization term on the discrete surface and a so called normal gradient stabilization term on the full tetrahedral elements in the active mesh. We prove optimal order a priori error estimates in the standard norm associated with the streamline diffusion method and bounds for the condition number of the resulting stiffness matrix. The condition number is of optimal order for a specific choice of method parameters. Numerical examples supporting our theoretical results are also included

Stabilized CutFEM for the convection problem on surfaces

We develop a stabilized cut finite element method for the convection problem on a surface based on continuous piecewise linear approximation and gradient jump stabilization terms. The discrete piecewise linear surface cuts through a background mesh consisting of tetrahedra in an arbitrary way and the finite element space consists of piecewise linear continuous functions defined on the background mesh. The variational form involves integrals on the surface and the gradient jump stabilization term is defined on the full faces of the tetrahedra. The stabilization term serves two purposes: first the method is stabilized and secondly the resulting linear system of equations is algebraically stable. We establish stability results that are analogous to the standard meshed flat case and prove h3 / 2order convergence in the natural norm associated with the method and that the full gradient enjoys h3 / 4order of convergence in L2. We also show that the condition number of the stiffness matrix is bounded by h- 2. Finally, our results are verified by numerical examples

Effect of anisotropy and destructuration on behavior of Haarajoki test embankment

This paper investigates the influence of anisotropy and destructuration on the behavior of Haarajoki test embankment, which was built by the Finnish National Road Administration as a noise barrier in 1997 on a soft clay deposit. Half of the embankment is constructed on an area improved with prefabricated vertical drains, while the other half is constructed on the natural deposit without any ground improvement. The construction and consolidation of the embankment is analyzed with the finite-element method using three different constitutive models to represent the soft clay. Two recently proposed constitutive models, namely S-CLAY1 which accounts for initial and plastic strain induced anisotropy, and its extension, called S-CLAY1S which accounts, additionally, for interparticle bonding and degradation of bonds, were used in the analysis. For comparison, the problem is also analyzed with the isotropic modified cam clay model. The results of the numerical analyses are compared with the field measurements. The simulations reveal the influence that anisotropy and destructuration have on the behavior of an embankment on soft clay

A large-strain radial consolidation theory for soft clays improved by vertical drains

A system of vertical drains with combined vacuum and surcharge preloading is an effective solution for promoting radial flow, accelerating consolidation. However, when a mixture of soil and water is deposited at a low initial density, a significant amount of deformation or surface settlement occurs. Therefore, it is necessary to introduce large-strain theory, which has been widely used to manage dredged disposal sites in one-dimensional theory, into radial consolidation theory. A governing equation based on Gibson's large-strain theory and Barron's free-strain theory incorporating the radial and vertical flows, the weight of the soil, variable hydraulic conductivity and compressibility during the consolidation process is therefore presented