Explore the vast range of titles available.

Transport of viscous fluid through porous media is a direct consequence of the pore structure. Here we investigate transport through a specific class of two-dimensional porous geometries, namely those formed by fluid-mechanical erosion. We investigate the tortuosity and dispersion by analyzing the first two statistical moments of tracer trajectories. For most initial configurations, tortuosity decreases in time as a result of erosion increasing the porosity. However, we find that tortuosity can also increase transiently in certain cases. The porosity-tortuosity relationships that result from our simulations are compared with models available in the literature. Asymptotic dispersion rates are also strongly affected by the erosion process, as well as by the number and distribution of the eroding bodies. Finally, we analyze the pore size distribution of an eroding geometry. The simulations are performed by combining a boundary integral equation solver for the fluid equations, a second-order stable time-stepping method to simulate erosion, and high-order numerical methods to stably and accurately resolve nearly touching eroded bodies and particle trajectories near the eroding bodies.

An informative look at the theory, computer implementation, and application of the scaled boundary finite element method   This reliable resource, complete with MATLAB, is an easy-to-understand introduction to the fundamental principles of the scaled boundary finite element method.

The flow in a Hele-Shaw cell with a time-increasing gap poses a unique shrinking interface problem. When the upper plate of the cell is lifted perpendicularly at a prescribed speed, the exterior less viscous fluid penetrates the interior more viscous fluid, which generates complex, time-dependent interfacial patterns through the Saffman–Taylor instability. The pattern formation process sensitively depends on the lifting speed and is still not fully understood. For some lifting speeds, such as linear or exponential speed, the instability is transient and the interface eventually shrinks as a circle. However, linear stability analysis suggests there exist shape invariant shrinking patterns if the gap $b(t)$ is increased more rapidly: $b(t)=\\left (1-({7}/{2})\\tau \\mathcal {C} t\\right )^{-{2}/{7}}$, where $\\tau$ is the surface tension and $\\mathcal {C}$ is a function of the interface perturbation mode $k$. Here, we use a spectrally accurate boundary integral method together with an efficient time adaptive rescaling scheme, which for the first time makes it possible to explore the nonlinear limiting dynamical behaviour of a vanishing interface. When the gap is increased at a constant rate, our numerical results quantitatively agree with experimental observations (Nase et al., Phys. Fluids, vol. 23, 2011, 123101). When we use the shape invariant gap $b(t)$, our nonlinear results reveal the existence of $k$-fold dominant, one-dimensional, web-like networks, where the fractal dimension is reduced to almost unity at late times. We conclude by constructing a morphology diagram for pattern selection that relates the dominant mode $k$ of the vanishing interface and the control parameter $\\mathcal {C}$.

In this paper, we investigate a sink-driven three-layer flow in a radial Hele-Shaw cell. The three fluids are of different viscosities, with one fluid occupying an annulus-like domain, forming two interfaces with the other two fluids. Using a boundary integral method and a semi-implicit time stepping scheme, we alleviate the numerical stiffness in updating the interfaces and achieve spectral accuracy in space. The interaction between the two interfaces introduces novel dynamics leading to rich pattern formation phenomena, manifested by two typical events: either one of the two interfaces reaches the sink faster than the other (forming cusp-like morphology), or they come very close to each other (suggesting a possibility of interface merging). In particular, the inner interface can be wrapped by the other to have both scenarios. We find that multiple parameters contribute to the dynamics, including the width of the annular region, the location of the sink, and the mobilities of the fluids.

Based on the kernel-free boundary integral method proposed by Ying and Henriquez (J Comput Phys 227(2):1046–1074, 2007 ), which is a second-order accurate method for general elliptic partial differential equations, this work develops it to be a fourth-order accurate version for the modified Helmholtz equation. The updated method is in line with the original one. Unlike the traditional boundary integral method, it does not need to know any analytical expression of the fundamental solution or Green’s function in evaluation of boundary or volume integrals. Boundary value problems under consideration are reformulated into Fredholm boundary integral equations of the second kind, whose corresponding discrete forms are solved with the simplest Krylov subspace iterative method, the Richardson iteration. During each iteration, a Cartesian grid based nine-point compact difference scheme is used to discretize the simple interface problem whose solution is the boundary or volume integral in the BIEs. The resulting linear system is solved by a fast Fourier transform based solver, whose computational work is roughly proportional to the number of grid nodes in the Cartesian grid used. As the discrete boundary integral equations are well-conditioned, the iteration converges within an essentially fixed number of steps, independent of the mesh parameter. Numerical results are presented to verify the solution accuracy and demonstrate the algorithm efficiency.

We adopt a boundary integral method to study the dynamics of a translating droplet confined in a Hele-Shaw cell in the Stokes regime. The droplet is driven by the motion of the ambient fluid with the same viscosity. We characterize the three-dimensional (3D) nature of the droplet interface and of the flow field. The interface develops an arc-shaped ridge near the rear-half rim with a protrusion in the rear and a laterally symmetric pair of higher peaks; this pair of protrusions has been identified by recent experiments (Huerre et al., Phys. Rev. Lett., vol. 115 (6), 2015, 064501) and predicted asymptotically (Burgess & Foster, Phys. Fluids A, vol. 2 (7), 1990, pp. 1105–1117). The mean film thickness is well predicted by the extended Bretherton model (Klaseboer et al., Phys. Fluids, vol. 26 (3), 2014, 032107) with fitting parameters. The flow in the streamwise wall-normal middle plane is featured with recirculating zones, which are partitioned by stagnation points closely resembling those of a two-dimensional droplet in a channel. Recirculation is absent in the wall-parallel, unconfined planes, in sharp contrast to the interior flow inside a moving droplet in free space. The preferred orientation of the recirculation results from the anisotropic confinement of the Hele-Shaw cell. On these planes, we identify a dipolar disturbance flow field induced by the travelling droplet and its $1/r^{2}$ spatial decay is confirmed numerically. We pinpoint counter-rotating streamwise vortex structures near the lateral interface of the droplet, further highlighting the complex 3D flow pattern.

In this paper, the finite element boundary integral (FEBI) method, for the first time, is applied to solve the 3-D arbitrary shaped eddy current nondestructive testing (ECNDT) problems. FEBI alleviates the extra computational costs of truncation region in finite element method (FEM) and the difficulty in derivation of the Green’s function in boundary element method (BEM). The boundary integral equation (BIE) selected is the combined field integral equation (CFIE) in the TENH (tangential testing of electric field and normal testing of magnetic field) form, which shows better convergence compared with other forms. In BEM, the equivalent electric and magnetic surface currents are expanded by Rao-Wilton-Glisson (RWG) vector basis functions. While in FEM, the electric field and electric surface current are expanded by tetrahedron-based edge elements and RWG vector basis functions, respectively. The discretized matrix achieved by BEM and FEM is coupled by the field continuity conditions. For ECNDT problems, inhomogeneous meshes are required due to the small size of cracks or slots than the whole solution domain. It makes the convergence for solving the coupled matrix formed by the sparse matrix generated by FEM and the dense matrix produced by BEM worse, thus, precondition is required for FEBI solution in iterative method, which complicates the solving procedure. To alleviate the cumbersome solving process, the inward-looking formulation method is studied to work as precondition by solving the inverse of FEM matrix directly, and then the coupled discretized matrix is solved iteratively. By evaluating several ECNDT benchmark cases involving the cylindrical flaws and surface slots, the predicted impedance changes achieved by FEBI method are compared with those by semi-analytical method, FEM, and experiment which demonstrates that the proposed FEBI method based forward solver can simulate the ECNDT problems both accurately and efficiently.

We present a kernel-free boundary integral method (KFBIM) for solving variable coefficients partial differential equations (PDEs) in a doubly connected domain. We focus our study on boundary value problems (BVP) and interface problems. A unique feature of the KFBIM is that the method does not require an analytical form of the Green’s function for designing quadratures but rather computes boundary or volume integrals by solving an equivalent interface problem on Cartesian mesh. We first decompose the problem defined in a doubly connected into two separate interface problems. The system of boundary integral equations is solved using the Krylov method. The method is second-order accurate in space, and its complexity is linearly proportional to the number of mesh points. Numerical examples demonstrate that the method is robust for variable coefficients PDEs, even for cases with large diffusion coefficients ratio and complex geometries where two interfaces are close.

We study the pairwise interactions of drops in an applied uniform DC electric field within the framework of the leaky dielectric model. We develop three-dimensional numerical simulations using the boundary integral method and an analytical theory assuming small drop deformations. We apply the simulations and the theory to explore the electrohydrodynamic interactions between two identical drops with arbitrary orientation of their line of centres relative to the applied field direction. Our results show a complex dynamics depending on the conductivities and permittivities of the drops and suspending fluids, and the initial drop pair alignment with the applied electric field.

A generalized boundary integral equation method for boundary value problems of two-dimensional isotropic lattice Laplacian is proposed in this paper. The proposed method is an extension of the classical boundary integral equation method with notable advantage. By utilizing the asymptotic expression of the fundamental solution at infinity, this method effectively addresses the challenge of numerical integration involving singular integral kernels. The introduction of Green’s formulas, Dirichlet and Neumann traces, and other tools which are parallel to the traditional integral equation method, form a solid foundation for the development of the generalized boundary integral equation method. The solvability of boundary integral equations and the solvability of lattice interface problem are important guarantees for the feasibility of this method, and these are emphasized in this paper. Subsequently, the generalized boundary integral equation method is applied to boundary value problems equipped with either Dirichlet or Neumann boundary conditions. Simple numerical examples demonstrate the accuracy and effectiveness of the generalized boundary integral equation method.