skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Finite Element Method for Two-Phase Flow with Material Viscous Interface
Abstract This paper studies a model of two-phase flow with an immersed material viscous interface and a finite element method for the numerical solution of the resulting system of PDEs.The interaction between the bulk and surface media is characterized by no-penetration and slip with friction interface conditions.The system is shown to be dissipative, and a model stationary problem is proved to be well-posed.The finite element method applied in this paper belongs to a family of unfitted discretizations.The performance of the method when model and discretization parameters vary is assessed.Moreover, an iterative procedure based on the splitting of the system into bulk and surface problems is introduced and studied numerically.  more » « less
Award ID(s):
2011444 1953535
PAR ID:
10320567
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Computational Methods in Applied Mathematics
Volume:
22
Issue:
2
ISSN:
1609-4840
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The paper studies a higher order unfitted finite element method for the Stokes system posed on a surface in ℝ 3 . The method employs parametric P k - P k −1 finite element pairs on tetrahedral bulk mesh to discretize the Stokes system on embedded surface. Stability and optimal order convergence results are proved. The proofs include a complete quantification of geometric errors stemming from approximate parametric representation of the surface. Numerical experiments include formal convergence studies and an example of the Kelvin--Helmholtz instability problem on the unit sphere. 
    more » « less
  2. Summary In this paper, we propose and analyze two stabilized mixed finite element methods for the dual‐porosity‐Stokes model, which couples the free flow region and microfracture‐matrix system through four interface conditions on an interface. The first stabilized mixed finite element method is a coupled method in the traditional format. Based on the idea of partitioned time stepping, the four interface conditions, and the mass exchange terms in the dual‐porosity model, the second stabilized mixed finite element method is decoupled in two levels and allows a noniterative splitting of the coupled problem into three subproblems. Due to their superior conservation properties and convenience of the computation of flux, mixed finite element methods have been widely developed for different types of subsurface flow problems in porous media. For the mixed finite element methods developed in this article, no Lagrange multiplier is used, but an interface stabilization term with a penalty parameter is added in the temporal discretization. This stabilization term ensures the numerical stability of both the coupled and decoupled schemes. The stability and the convergence analysis are carried out for both the coupled and decoupled schemes. Three numerical experiments are provided to demonstrate the accuracy, efficiency, and applicability of the proposed methods. 
    more » « less
  3. null (Ed.)
    Abstract Shock waves from underwater and air explosions are significant threats to surface and underwater vehicles and structures. Recent studies on the mechanical and thermal properties of various phase-separated elastomers indicate the possibility of applying these materials as a coating to mitigate shock-induced structural failures. To demonstrate this approach and investigate its efficacy, this paper presents a fluid-structure coupled computational model capable of predicting the dynamic response of air-backed bilayer (i.e. elastomer coating – metal substrate) structures submerged in water to hydrostatic and underwater explosion loads. The model couples a three-dimensional multiphase finite volume computational fluid dynamics model with a nonlinear finite element computational solid dynamics model using the FIVER (FInite Volume method with Exact multi-material Riemann solvers) method. The kinematic boundary condition at the fluid-structure interface is enforced using an embedded boundary method that is capable of handling large structural deformation and topological changes. The dynamic interface condition is enforced by formulating and solving local, one-dimensional fluid-solid Riemann problems, which is well-suited for transferring shock and impulsive loads. The capability of this computational model is demonstrated through a numerical investigation of hydrostatic and shock-induced collapse of aluminum tubes with polyurea coating on its inner surface. The thickness of the structure is resolved explicitly by the finite element mesh. The nonlinear material behavior of polyurea is accounted for using a hyper-viscoelastic constitutive model featuring a modified Mooney-Rivlin equation and a stress relaxation function in the form of prony series. Three numerical experiments are conducted to simulate and compare the collapse of the structure in different loading conditions, including a constant pressure, a fluid environment initially in hydrostatic equilibrium, and a two-phase fluid flow created by a near-field underwater explosion. 
    more » « less
  4. Interface problems have wide applications in modern scientific research. Obtaining accurate numerical solutions of multi-domain problems involving triple junction conditions remains a significant challenge. In this paper, we develop an efficient finite element method based on non-body-fitting meshes for solving multi-domain elliptic interface problems. We follow the idea of immersed finite element by modifying local basis functions to accommodate interface conditions. We enrich the local finite element space by adding new basis functions for handling non-homogeneous flux jump. The numerical scheme is symmetric and positive definite. Numerical experiments are provided to demonstrate the features of our method. 
    more » « less
  5. null (Ed.)
    Numerical modelling of deformation in hydromechanical systems can be time-consuming using fully coupled classical numerical methods for large representative porous media samples. In this paper, we present a new two-way coupled partitioned fluid–solid system. The coupled system is applied for simulating geomechanical processes at the pore-scale. We track the deformation of the solid resulting from the drainage of resident fluids in the pores, as well as the evolution of fluid properties from dynamic loading. The finite element method is responsible for capturing the structural deformation in the coupled system while the dynamic pore network is used for modelling multiphase flow in the fluid subsystem. A fictitious fluid–solid interface is created at each pore network-finite element node junction via convex hulling, followed by data exchange using linear interpolation. The results show good agreement with a pre-existing coupled finite volume model and the computations are completed in much less time. 
    more » « less