BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jun Lai

Jun Lai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysisphysics.flu-dynphysics.optics

Catalog footprint

What is connected

11works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

A spectral boundary integral method for the elastic obstacle scattering problem in three dimensions

In this paper, we consider the scattering of a plane wave by a rigid obstacle embedded in a homogeneous and isotropic elastic medium in three dimensions. Based on the Helmholtz decomposition, the elastic scattering problem is reduced to a coupled boundary value problem for the Helmholtz and Maxwell equations. A novel system of boundary integral equations is formulated and a spectral boundary integral method is developed for the coupled boundary value problem. Numerical experiments are presented to demonstrate the superior performance of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Study of a droplet breakup process in decaying homogeneous isotropic turbulence based on the phase-field DUGKS approach

The breakup of a spherical droplet in a decaying homogeneous isotropic turbulence is studied by solving the Cahn-Hilliard-Navier-Stokes equations, using the discrete unified gas kinetic scheme combined with the free-energy-based phase-field model. We focus on the combined effects of turbulence and surface tension on the breakup process by assuming that the two fluid phases have the same density and same viscosity. The key physical parameters of the system include the volume fraction (6.54%), the initial Weber number (21.7), and the initial Taylor microscale Reynolds number (58). Three distinct stages of droplet evolution are identified, namely, the deformation stage when the initially spherical droplet evolves into an irregular geometric shape with complex structures, the breakup stage when many daughter droplets are formed, and the restoration stage when the droplets relax towards spherical shape. These three stages are analyzed systematically from several perspectives: (1) a geometric perspective concerning the maximum equivalent diameter, the total number of droplets, total interface area, and probability distribution of droplet diameters, (2) a dynamic perspective concerning the evolution of local velocity and vorticity at the fluid-fluid interface, (3) a global perspective concerning the evolution of average kinetic energy / dissipation rate and their Fourier spectra, (4) spherical harmonics based energetics concerning simultaneous transfer of kinetic energy across different length scales and different radii relative to initial droplet center, and (5) the time evolution of global kinetic energy and free energy of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A highly accurate boundary integral method for the elastic obstacle scattering problem

Consider the scattering of a time-harmonic plane wave by a rigid obstacle embedded in a homogeneous and isotropic elastic medium in two dimensions. In this paper, a novel boundary integral formulation is proposed and its highly accurate numerical method is developed for the elastic obstacle scattering problem. More specifically, based on the Helmholtz decomposition, the model problem is reduced to a coupled boundary integral equation with singular kernels. A regularized system is constructed in order to handle the degenerated integral operators. The semi-discrete and full-discrete schemes are studied for the boundary integral system by using the trigonometric collocation method. Convergence is established for the numerical schemes in some appropriate Sobolev spaces. Numerical experiments are presented for both smooth and nonsmooth obstacles to demonstrate the superior performance of the proposed method.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An inverse acoustic-elastic interaction problem with phased or phaseless far-field data

Consider the scattering of a time-harmonic acoustic plane wave by a bounded elastic obstacle which is immersed in a homogeneous acoustic medium. This paper concerns an inverse acoustic-elastic interaction problem, which is to determine the location and shape of the elastic obstacle by using either the phased or phaseless far-field data. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The jump relations are studied for the second derivatives of the single-layer potential in order to establish the corresponding boundary integral equations. The well-posedness is discussed for the solution of the coupled boundary integral equations. An efficient and high order Nyström-type discretization method is proposed for the integral system. A numerical method of nonlinear integral equations is developed for the inverse problem. For the case of phaseless data, we show that the modulus of the far-field pattern is invariant under a translation of the obstacle. To break the translation invariance, an elastic reference ball technique is introduced. We prove that the inverse problem with phaseless far-field pattern has a unique solution under certain conditions. In addition, a numerical method of the reference ball technique based nonlinear integral equations is also proposed for the phaseless inverse problem. Numerical experiments are provided to demonstrate the effectiveness and robustness of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A Fast Solver for the Elastic Scattering of Multiple Particles

Consider the elastic scattering of a time-harmonic wave by multiple well separated rigid particles in two dimensions. To avoid using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert the boundary value problem of the elastic wave equation into a coupled boundary value problem of Helmholtz equations. Based on single, double, and combined layer potentials with the simpler Green's function of the Helmholtz equation, we present three different boundary integral equations for the coupled boundary value problem. The well-posedness of the new integral equations are established. Computationally, a scattering matrix based method is proposed to evaluate the elastic wave for arbitrarily shaped particles. The method uses the local expansion for the incident wave and the multipole expansion for the scattered wave. The linear system of algebraic equations is solved by GMRES with fast multipole method (FMM) acceleration. Numerical results show that the method is fast and highly accurate for solving the elastic scattering problem with multiple particles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Quantitative Study of the Effect of Cladding Thickness on Modal Confinement Loss in Photonic Waveguides

There has been increasing interest in making the photonic devices more and more compact in the integrated photonics industry, and one of the important questions for manufacturers and design engineers is how to quantify the effect of the finite cladding thickness on the modal confinement loss of photonic waveguides. This requires at least six to seven digits accuracy for the computation of propagation constant $β$ since the modal confinement loss is proportional to the imaginary part of $β$ that is six to seven orders of magnitude smaller than its real part by the industrial standard. In this paper, we present an accurate and efficient method to compute the propagation constant of electromagnetic modes of photonic waveguides with arbitrary number of (nonsmooth) inclusions in a layered media. The method combines a well-conditioned boundary integral equation formulation for photonic waveguides which requires the discretization of the material interface only, and efficient Sommerfeld integral representations to treat the effect of the layered medium. Our scheme is capable of calculating the propagation loss of the electromagnetic modes with high fidelity, even for waveguides with corners imbedded in a cladding material of finite thickness. The numerical results, with more than $10$-digit accuracy, show quantitatively that the modal confinement loss of the rectangular waveguide increases exponentially fast as the cladding thickness decreases.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A new hybrid integral representation for frequency domain scattering in layered media

A variety of problems in acoustic and electromagnetic scattering require the evaluation of impedance or layered media Green's functions. Given a point source located in an unbounded half-space or an infinitely extended layer, Sommerfeld and others showed that Fourier analysis combined with contour integration provides a systematic and broadly effective approach, leading to what is generally referred to as the Sommerfeld integral representation. When either the source or target is at some distance from an infinite boundary, the number of degrees of freedom needed to resolve the scattering response is very modest. When both are near an interface, however, the Sommerfeld integral involves a very large range of integration and its direct application becomes unwieldy. Historically, three schemes have been employed to overcome this difficulty: the method of images, contour deformation, and asymptotic methods of various kinds. None of these methods make use of classical layer potentials in physical space, despite their advantages in terms of adaptive resolution and high-order accuracy. The reason for this is simple: layer potentials are impractical in layered media or half-space geometries since they require the discretization of an infinite boundary. In this paper, we propose a hybrid method which combines layer potentials (physical-space) on a finite portion of the interface together with a Sommerfeld-type (Fourier) correction. We prove that our method is efficient and rapidly convergent for arbitrarily located sources and targets, and show that the scheme is particularly effective when solving scattering problems for objects which are close to the half-space boundary or even embedded across a layered media interface.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Second Kind Integral Equation Formulation for the Mode Calculation of Optical Waveguides

We present a second kind integral equation (SKIE) formulation for calculating the electromagnetic modes of optical waveguides, where the unknowns are only on material interfaces. The resulting numerical algorithm can handle optical waveguides with a large number of inclusions of arbitrary irregular cross section. It is capable of finding the bound, leaky, and complex modes for optical fibers and waveguides including photonic crystal fibers (PCF), dielectric fibers and waveguides. Most importantly, the formulation is well conditioned even in the case of nonsmooth geometries. Our method is highly accurate and thus can be used to calculate the propagation loss of the electromagnetic modes accurately, which provides the photonics industry a reliable tool for the design of more compact and efficient photonic devices. We illustrate and validate the performance of our method through extensive numerical studies and by comparison with semi-analytical results and previously published results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fast and robust solver for the scattering from a layered periodic structure containing multi-particle inclusions

We present a solver for plane wave scattering from a periodic dielectric grating with a large number $M$ of inclusions lying in each period of its middle layer.Such composite material geometries have a growing role in modern photonic devices and solar cells. The high-order scheme is based on boundary integral equations, and achieves many digits of accuracy with ease. The usual way to periodize the integral equation---via the quasi-periodic Green's function---fails at Wood's anomalies. We instead use the free-space Green's kernel for the near field, add auxiliary basis functions for the far field, and enforce periodicity in an expanded linear system; this is robust for all parameters. Inverting the periodic and layer unknowns, we are left with a square linear system involving only the inclusion scattering coefficients. Preconditioning by the single-inclusion scattering matrix, this is solved iteratively in $O(M)$ time using a fast matrix-vector product. Numerical experiments show that a diffraction grating containing $M=1000$ inclusions per period can be solved to 9-digit accuracy in under 5 minutes on a laptop.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fast direct solver for high frequency scattering from a large cavity in two dimensions

We present a fast direct solver for the simulation of electromagnetic scattering from an arbitrarily-shaped, large, empty cavity embedded in an infinite perfectly conducting half space. The governing Maxwell equations are reformulated as a well-conditioned second kind integral equation and the resulting linear system is solved in nearly linear time using a hierarchical matrix factorization technique. We illustrate the performance of the scheme with several numerical examples for complex cavity shapes over a wide range of frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A fast solver for multi-particle scattering in a layered medium

In this paper, we consider acoustic or electromagnetic scattering in two dimensions from an infinite three-layer medium with thousands of wavelength-size dielectric particles embedded in the middle layer. Such geometries are typical of microstructured composite materials, and the evaluation of the scattered field requires a suitable fast solver for either a single configuration or for a sequence of configurations as part of a design or optimization process. We have developed an algorithm for problems of this type by combining the Sommerfeld integral representation, high order integral equation discretization, the fast multipole method and classical multiple scattering theory. The efficiency of the solver is illustrated with several numerical experiments.

0 citations0 reviewsNo code linkedOpen access