In this paper, a novel penalty method based on the immersed boundary formulation is proposed for simulating the transient Stokes flow with an inextensible interface enclosing a suspended solid particle. The main idea of this approach relies on the penalty techniques by modifying the constitutive equation of Stokes flow to weaken the incompressibility condition, relating the surface divergence to the elastic tension $\sigma$ to relax the interface's inextensibility, and connecting the particle surface-velocity with the particle surface force ${\mbf F}$ to regularize the particle's rigid motion. The advantage of these regularized governing equations is that when they are discretized by the standard centered difference scheme on a staggered grid, the resulting linear system can easily be reduced by eliminating the unknowns $p_h$, $\sigma_h$ and ${\mbf F}_h$ directly, so that we just need to solve a smaller linear system of the velocity approximation ${\mbf u}_h$. This advantage is preserved and even enhanced when such approach is applied to the transient Stokes flow with multiple compound vesicles. Moreover, this smaller linear system is symmetric and negative-definite, which enables us to use efficient linear solvers. Another important feature of the proposed method is that the discretization scheme is unconditionally stable in the sense that an appropriately defined energy functional associated with the discrete system is decreasing and hence bounded in time. We numerically test the accuracy and stability of the immersed boundary discretization scheme. The tank-treading and tumbling motions of inextensible interface with a suspended solid particle in the simple shear flow will be studied extensively. The simulation of the motion of multiple compound vesicles will be performed as well. Numerical results illustrate the superior performance of the proposed penalty method.

In this paper, we combine deep learning concepts and some proper orthogonal decomposition (POD) model reduction methods for predicting flow in heterogeneous porous media. Nonlinear flow dynamics is studied, where the dynamics is regarded as a multi-layer network. The solution at the current time step is regarded as a multi-layer network of the solution at the initial time and input parameters. As for input, we consider various sources, which include source terms (well rates), permeability fields, and initial conditions. We consider the flow dynamics, where the solution is known at some locations and the data is integrated to the flow dynamics by modifying the reduced-order model. This approach allows modifying the reduced-order formulation of the problem. Because of the small problem size, limited observed data can be handled. We consider enriching the observed data using the computational data in deep learning networks. The basis functions of the global reduced order model are selected such that the degrees of freedom represent the solution at observation points. This way, we can avoid learning basis functions, which can also be done using neural networks. We present numerical results, where we consider channelized permeability fields, where the network is constructed for various channel configurations. Our numerical results show that one can achieve a good approximation using forward feed maps based on multi-layer networks.

Using the vorticity and stream function variables is an effective way to compute 2-D incompressible flow due to the facts that the incompressibility constraint for the velocity is automatically satisfied, the pressure variable is eliminated, and high order schemes can be efficiently implemented. However, a difficulty arises in a multi-connected computational domain in determining the constants for the stream function on the boundary of the “holes”. This is an especially challenging task for the calculation of unsteady flows, since these constants vary with time to reflect the total fluxes of the flow in each sub-channel. In this paper, we propose an efficient method in a finite difference setting to solve this problem and present some numerical experiments, including an accuracy check of a Taylor vortex-type flow, flow past a non-symmetric square, and flow in a heat exchanger.

In this paper, we discuss a general multiscale model reduction framework based on multiscale finite element methods. We give a brief overview of related multiscale methods. Due to page limitations, the overview focuses on a few related methods and is not intended to be comprehensive. We present a general adaptive multiscale model reduction framework, the Generalized Multiscale Finite Element Method. Besides the method's basic outline, we discuss some important ingredients needed for the method's success. We also discuss several applications. The proposed method allows performing local model reduction in the presence of high contrast and no scale separation.

We study the generalized eigenvalue problems (GEPs) derived from modeling the surface acoustic wave in piezoelectric materials with periodic inhomogeneity. The eigenvalues appear in the reciprocal pairs due to periodic boundary conditions in the modeling. By transforming the GEP into a T-palindromic quadratic eigenvalue problem (TPQEP), the reciprocal relationship of the eigenvalues can be maintained. In this paper, we outline four recently developed structure-preserving algorithms, SA, SDA, TSHIRA and GTSHIRA, for solving the TPQEP. Numerical comparisons on the accuracy and the computational costs of these algorithm are presented. The eigenvalues close to unit circle on the complex plane are of interest in the area of filter and sensor designs. Our numerical results show that the Arnoldi-type structure-preserving algorithms TSHIRA and GTSHIRA with re-symplectic and re-bi-isotropic processes