Convergence of a non-stiff boundary integral method for interfacial flows with surface tension

Boundary integral methods to simulate interfacial flows are very sensitive to numerical instabilities. In addition, surface tension introduces nonlinear terms with high order spatial derivatives into the interface dynamics. This makes the spatial discretization even more difficult and, at the same time, imposes a severe time step constraint for stable explicit time integration methods.

A proof of the convergence of a reformulated boundary integral method for two-density fluid interfaces with surface tension is presented. The method is based on a scheme introduced by Hou, Lowengrub and Shelley [ J. Comp. Phys. 114 (1994), pp. 312-338] to remove the high order stability constraint or stiffness. Some numerical filtering is applied carefully at certain places in the discretization to guarantee stability. The key of the proof is to identify the most singular terms of the method and to show, through energy estimates, that these terms balance one another.

The analysis is at a time continuous-space discrete level but a fully discrete case for a simple Hele-Shaw interface is also studied. The time discrete analysis shows that the high order stiffness is removed and also provides an estimate of how the CFL constraint depends on the curvature and regularity of the solution.

The robustness of the method is illustrated with several numerical examples. A numerical simulation of an unstably stratified two-density interfacial flow shows the roll-up of the interface; the computations proceed up to a time where the interface is about to pinch off and trapped bubbles of fluid are formed. The method remains stable even in the full nonlinear regime of motion. Another application of the method shows the process of drop formation in a falling single fluid.

     

The Neumann problem for the Laplace equation on general domains

The solution of the weak Neumann problem for the Laplace equation with a distribution as a boundary condition is studied on a general open set G in the Euclidean space. It is shown that the solution of the problem is the sum of a constant and the Newtonian potential corresponding to a distribution with finite energy supported on ∂G. If we look for a solution of the problem in this form we get a bounded linear operator. Under mild assumptions on G a necessary and sufficient condition for the solvability of the problem is given and the solution is constructed. The research was supported by the Academy of Sciences of the Czech Republic, Institutional Research Plan No. AV0Z10190503.     

弹性力学问题解唯一的边界积分方程

从积分方程式出发,应用基本解的特性分析,说明在力边值问题中,位移边界积分方程和面力边界积分方程的位移解不唯一.提出了位移解唯一的条件,建立了唯一解的位移边界积分方程和面力边界积分方程.实例计算结果表明唯一解的边界积分方程是有效的.     

A stable boundary elements method for magnetohydrodynamic channel flows at high Hartmann numbers

The article is devoted to extension of boundary element method (BEM) for solving coupled equations in velocity and induced magnetic field for time dependent magnetohydrodynamic (MHD) flows through a rectangular pipe. The BEM is equipped with finite difference approach to solve MHD problem at high Hartmann numbers up to 10⁶. In fact, the finite difference approach is used to approximate partial derivatives of unknown functions at boundary points respect to outward normal vector. It yields a numerical method with no singular boundary integrals. Besides, a new approach is suggested in this article where transforms 2D singular BEM's integrals to 1D nonsingular ones. The new approach reduces computational cost, significantly. Note that the stability of the numerical scheme is proved mathematically when computational domain is discretized uniformly and Hartmann number is 40 times bigger than length of boundary elements. Numerical examples show behavior of velocity and induced magnetic field across the sections.     

Electromagnetic scattering by a homogeneous chiral obstacle in a chiral environment

Time-harmonic electromagnetic waves are scattered by a homogeneouschiral obstacle embedded in a chiral environment. The correspondingtransmission problem is reduced, via Bohren's decomposition,to an integral equation over the interface between the obstacleand the surrounding medium. This integral equation is shownto be uniquely solvable except for a discrete set of electromagneticparameters of the obstacle.     

Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R~2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.     

The thermostat problem with a nonlocal nonlinear boundary condition

The thermostat controller for an air-conditioning system isusually placed in a position at some distance from the unitand this can lead to large swings in temperature. This paperaddresses this question by studying a paradigm—a one-dimensionalheat conduction equation with and without heat loss, and wherethe flux of heat extracted or input by the unit is consideredto be a function of the temperature at the other end. The essential results are that the system can be unstable andthat this is exacerbated both by a more powerful air-conditioningunit and by more efficient insulation.     

An Artificial Boundary Condition for the Vortex Movements in Two Dimensions

1 Motion of vortices and cloud in cell method The motion of incompressible inviscid ?ow in two dimensions can be described by the equations ?u ?t (u ?) ρ1 ?P = f (1) ?u = 0 , (2) where u = (u, v), ρ, P , and f = (f1, f2) denote ?uid velocity, densit…     

Asymptotic expansions for currents caused by small interface changes of an electromagnetic inclusion

We consider solutions to the Helmholtz equation in two dimensions. The aim of this article is to advance the development of high-order asymptotic expansions for boundary perturbations of currents caused by small perturbations of the shape of an inhomogeneity with 𝒞²-boundary. The work represents a natural completion of Ammari et al. [H. Ammari, H. Kang, M. Lim, and H. Zribi, Conductivity interface problems. Part I: Small perturbations of an interface, Trans. Am. Math. Soc. 363 (2010), pp. 2901–2922], where the solution for the Helmholtz equation is represented by a system and the proof of our asymptotic expansion is radically different from Ammari et al. (2010). Our derivation is rigorous and is based on the field expansion method. Its proof relies on layer potential techniques. It plays a key role in developing effective algorithms to determine certain properties of the shape of an inhomogeneity based on boundary measurements.     

关于薄板的无网格局部边界积分方程方法中的友解

无网格局部边界积分方程方法是最近发展起来的一种新的数值方法,这种方法综合了伽辽金有限元、边界元和无单元伽辽金法的优点,是一种具有广阔应用前景的、真正的无网格方法．把无网格局部边界积分方程方法应用于求解薄板问题,给出了薄板无网格局部边界积分方程方法所需要的友解及其全部公式．     

Treatment of singularities in 2-D domains using BIEM

The singularities which arise when there is a sudden change of boundary conditions are modelled using spectral shape interpolation functions. The procedure can be used for elasticity as well as potential theory and to any degree of accuracy with respect to the smooth part of the curve.     

Boundary integral method for scattering problems with cracks buried in a piecewise homogeneous medium

We consider the scattering of time‐harmonic acoustic plane waves by a crack buried in a piecewise homogeneous medium. The integral representation for a solution is obtained in the form of potentials by using Green's formula. The density in potentials satisfies the uniquely solvable Fredholm integral equation. Then we obtain the existence and uniqueness of the solution. Copyright © 2011 John Wiley & Sons, Ltd.     

Boundary integral method for Stokes flow past a porous body

In this paper we obtain an indirect boundary integral method in order to prove existence and uniqueness of the classical solution to a boundary value problem for the Stokes–Brinkman-coupled system, which describes an unbounded Stokes flow past a porous body in terms of Brinkman's model. Therefore, one assumes that the flow inside the body is governed by the continuity and Brinkman equations. Some asymptotic results in both cases of large and, respectively, of low permeability are also obtained. Copyright © 2007 John Wiley & Sons, Ltd.     

The method of integral equation for scattering problem with a mixed crack

We consider the scattering of an electromagnetic time‐harmonic plane wave by an infinite cylinder having a mixed open crack (or arc) in R² as the cross section. The crack is made up of two parts, and one of the two parts is (possibly) coated by a material with surface impedance λ. We transform the scattering problem into a system of boundary integral equations by adopting a potential approach, and establish the existence and uniqueness of a weak solution to the system by the Fredholm theory. Copyright © 2011 John Wiley & Sons, Ltd.     

Transonic potential flows in a convergent-divergent approximate nozzle

In this paper we prove existence, uniqueness and regularity of certain perturbed (subsonic-supersonic) transonic potential flows in a two-dimensional Riemannian manifold with “convergent-divergent” metric, which is an approximate model of the de Laval nozzle in aerodynamics. The result indicates that transonic flows obtained by quasi-one-dimensional flow model in fluid dynamics are stable with respect to the perturbation of the velocity potential function at the entry (i.e., tangential velocity along the entry) of the nozzle. The proof is based upon linear theory of elliptic-hyperbolic mixed type equations in physical space and a nonlinear iteration method.     

Stochastic flows of mappings

In this paper,the stochastic flow of mappings generated by a Feller convolution semigroup on a compact metric space is studied.This kind of flow is the generalization of superprocesses of stochastic flows and stochastic diffeomorphism induced by the strong solutions of stochastic differential equations.     

泛函积分方程的正解与应用

运用锥拉伸与压缩不动点定理,研究一类泛函积分方程的正解．把所得结果应用于研究常微分方程积分边值问题正解的存在性,所得结果推广和本质上改进了马如云等作者最近的一些结果．     

Computation of axisymmetric MHD flows in a channel with an external longitudinal magnetic field

A numerical study of two-dimensional plasma flows in coaxial channels of plasma accelerators is presented. Two new results are obtained. First, for the computation of MHD problems belonging to the class under consideration, Zalesak’s method is used. It is based on an explicit finite difference scheme with flux correction. This method is free of space splitting, and, therefore, is well suited for parallel computations on multiprocessors. Second, the statement of the problem is extended so that the acceleration of the plasma by the azimuth magnetic self-field goes on in the presence of an external longitudinal field. The results of test computations demonstrate the efficiency of the method and made it possible to investigate the influence of the longitudinal field on the properties of the plasma flows.     

Numerical simulation of pressure-driven startup laminar flows through a planar T-junction channel

A start-up flow of a viscous incompressible fluid in a T-junction channel is studied numerically. The flow starting from rest is driven by a constant pressure drops suddenly applied between the entries and exits of a planar T-junction channel. The Navier-Stokes equations in primitive variables are solved numerically using finite-volume techniques. Predicted variations with time of the volume flow rates and the flow patterns are presented for several values of pressure drops. It has been shown that a start-up flow can pass through different regimes (or different flow direction) before asymptotically reaching steady state distribution.