A conservative SPH method for surfactant dynamics

In this paper, a Lagrangian particle method is proposed for the simulation of multiphase flows with surfactant. The model is based on the multiphase smoothed particle hydrodynamics (SPH) framework of Hu and Adams (2006) [1]. Surface-active agents (surfactants) are incorporated into our method by a scalar quantity describing the local concentration of molecules in the bulk phase and on the interface. The surfactant dynamics are written in conservative form, thus global mass of surfactant is conserved exactly. The transport model of the surfactant accounts for advection and diffusion. Within our method, we can simulate insoluble surfactant on an arbitrary interface geometry as well as interfacial transport such as adsorption or desorption. The flow-field dynamics and the surfactant dynamics are coupled through a constitutive equation, which relates the local surfactant concentration to the local surface-tension coefficient. Hence, the surface-tension model includes capillary and Marangoni-forces. The present numerical method is validated by comparison with analytic solutions for diffusion and for surfactant dynamics. More complex simulations of an oscillating bubble, the bubble deformation in a shear flow, and of a Marangoni-force driven bubble show the capabilities of our method to simulate interfacial flows with surfactants.     

An accurate adaptive solver for surface-tension-driven interfacial flows

A method combining an adaptive quad/octree spatial discretisation, geometrical Volume-Of-Fluid interface representation, balanced-force continuum-surface-force surface-tension formulation and height-function curvature estimation is presented. The extension of these methods to the quad/octree discretisation allows adaptive variable resolution along the interface and is described in detail. The method is shown to recover exact equilibrium (to machine accuracy) between surface-tension and pressure gradient in the case of a stationary droplet, irrespective of viscosity and spatial resolution. Accurate solutions are obtained for the classical test case of capillary wave oscillations. An application to the capillary breakup of a jet of water in air further illustrates the accuracy and efficiency of the method. The source code of the implementation is freely available as part of the Gerris flow solver.     

The Discontinuous Galerkin Time-Domain method for Maxwell’s equations with anisotropic materials

The Discontinuous Galerkin method is an accurate and efficient way to numerically solve the time-dependent Maxwell equations. In this paper, we extend the basic, two-dimensional formulation for isotropic materials to allow anisotropic permittivity tensors. Using a reference system with an analytical solution, we demonstrate that our extensions do not alter the superior convergence characteristics of the fundamental algorithm. We further apply our method to cylindrical invisibility cloaks to investigate the performance which can be achieved in experiments.     

A stable high-order Spectral Difference method for hyperbolic conservation laws on triangular elements

Numerical schemes using piecewise polynomial approximation are very popular for high order discretization of conservation laws. While the most widely used numerical scheme under this paradigm appears to be the Discontinuous Galerkin method, the Spectral Difference scheme has often been found attractive as well, because of its simplicity of formulation and implementation. However, recently it has been shown that the scheme is not linearly stable on triangles. In this paper we present an alternate formulation of the scheme, featuring a new flux interpolation technique using Raviart–Thomas spaces, which proves stable under a similar linear analysis in which the standard scheme failed. We demonstrate viability of the concept by showing linear stability both in the semi-discrete sense and for time stepping schemes of the SSP Runge–Kutta type. Furthermore, we present convergence studies, as well as case studies in compressible flow simulation using the Euler equations.     

Background field method for nonlinear σ-model in stochastic quantization

We formulate the background field method for the nonlinear σ-model in stochastic quantization. We demonstrate a one-loop calculation for a two-dimensional non-linear σ-model on a general riemannian manifold based on our formulation. The formulation is consistent with the known results in ordinary quantization. As a simple application, we also analyse the multiplicative renormalization of the O(N) nonlinear σ-model.     

三维声学多层快速多极子边界元及其应用

快速多极子边界元算法可以加速矩阵和向量乘法运算, 将传统边界元算法的计算量和内存占用量分别降为O(N log²N)和O(N), 适用于大型声学模型模拟计算. 本文发展了一种基于Burton-Miller方程的三维多层声学快速多极子边界元算法. 将新的自适应树状算法应用到对角形式的快速多极子边界元算法, 并使用最新提出的解析式源点矩计算公式, 进一步提高了快速多极子边界元的计算效率. 绝对软球体在内部共振频率处的散射声场计算, 验证了所发展算法在共振频率处求解的正确性. 与Bapat所提供的程序在多脉动球体辐射声场计算精度的比较, 验证了算法及程序在大型模型声学计算中的准确性, 同时显示了其求解的高效性. 最后, 将该算法用于车内声场及水下声学探测的分析计算.     

J-matrix method of scattering in one dimension: The nonrelativistic theory

We formulate a theory of nonrelativistic scattering in one dimension based on the J-matrix method. The scattering potential is assumed to have a finite range such that it is well represented by its matrix elements in a finite subset of a basis that supports a tridiagonal matrix representation for the reference wave operator. Contrary to our expectation, the 1D formulation reveals a rich and highly nontrivial structure compared to the 3D formulation. Examples are given to demonstrate the utility and accuracy of the method. It is hoped that this formulation constitutes a viable alternative to the classical treatment of 1D scattering problem and that it will help unveil new and interesting applications.     

浸润接触线的摩擦性质与固体表面张力的Wenzel行为

依据Adam和Jessop关于固-液-气三相接触线具有静摩擦性的观点，能从浸润接触角数据推算 出固体的表面张力系数，但结果显然会与摩擦条件的具体形式有关.以报道的不锈钢和聚丙 烯实验数据为例，通过对比Mises和Amonton两种摩擦定律给出的表面张力系数，后者给出的 张力才具有总是随表面粗糙性的增大而增加的Wenzel效应，以及当固-液，固-气界面的张力 系数之差等于液-气界面的张力系数时，退后角将变为零的全浸润条件.这似乎表明用Amonto n定律描写接触线的静摩擦要更为合理. 关键词： 浸润角 静摩擦 固体表面张力     

Understanding quantum entanglement by thermo field dynamics

We propose a new method to understand quantum entanglement using the thermo field dynamics (TFD) described by a double Hilbert space. The entanglement states show a quantum-mechanically complicated behavior. Our new method using TFD makes it easy to understand the entanglement states, because the states in the tilde space in TFD play a role of tracer of the initial states. For our new treatment, we define an extended density matrix on the double Hilbert space. From this study, we make a general formulation of this extended density matrix and examine some simple cases using this formulation. Consequently, we have found that we can distinguish intrinsic quantum entanglement from the thermal fluctuations included in the definition of the ordinary quantum entanglement at finite temperatures. Through the above examination, our method using TFD can be applied not only to equilibrium states but also to non-equilibrium states. This is shown using some simple finite systems in the present paper.     

Efficient numerical methods for computing ground states and dynamics of dipolar Bose–Einstein condensates

New efficient and accurate numerical methods are proposed to compute ground states and dynamics of dipolar Bose–Einstein condensates (BECs) described by a three-dimensional (3D) Gross–Pitaevskii equation (GPE) with a dipolar interaction potential. Due to the high singularity in the dipolar interaction potential, it brings significant difficulties in mathematical analysis and numerical simulations of dipolar BECs. In this paper, by decoupling the two-body dipolar interaction potential into short-range (or local) and long-range interactions (or repulsive and attractive interactions), the GPE for dipolar BECs is reformulated as a Gross–Pitaevskii–Poisson type system. Based on this new mathematical formulation, we prove rigorously existence and uniqueness as well as nonexistence of the ground states, and discuss the existence of global weak solution and finite time blow-up of the dynamics in different parameter regimes of dipolar BECs. In addition, a backward Euler sine pseudospectral method is presented for computing the ground states and a time-splitting sine pseudospectral method is proposed for computing the dynamics of dipolar BECs. Due to the adoption of new mathematical formulation, our new numerical methods avoid evaluating integrals with high singularity and thus they are more efficient and accurate than those numerical methods currently used in the literatures for solving the problem. Extensive numerical examples in 3D are reported to demonstrate the efficiency and accuracy of our new numerical methods for computing the ground states and dynamics of dipolar BECs.     

Terahertz Thermometry: Combining Hyperspectral Imaging and Temperature Mapping at Terahertz Frequencies

The accurate and non‐invasive determination of multiple physical parameters, with well‐defined spatial resolution, is crucial for applications in manufacturing, chemistry, medicine and biology. Specifically, the ability to simultaneously measure both temperature and spectral signatures is still experimentally unavailable. To this end, we propose a mapping technique for biological systems, which exploits a linear correlation between terahertz wave reflectivity and temperature, and allows to spatially and spectrally resolve thermal distributions. This method is applied to a model biological system in two relevant cases where in one example, nanoplasmonic‐induced photothermal effects are imaged gaining new insights into collective heating phenomena. In the second example, we demonstrate a joint thermal‐hyperspectral imaging approach to chemically map the presence of a model drug formulation and simultaneously investigate its thermal stability in our biological system. This concept can be easily extended and widely applied to all materials that demonstrate a measurable change in their dielectric properties.     

Generalized Foldy-Lax formulation

We present a method for numerical wave propagation in a heterogeneous medium. The medium is defined in terms of an extended scatterer or target which is surrounded by many small scatterers. By extending the classic Foldy-Lax formulation we developed an efficient algorithm for numerical wave propagation in two dimension. In the method that we set forth multiple scattering among the point scatterers and the extended target is fully taken into account via a boundary integral formulation coupled with the Foldy-Lax formulation. This formulation forms the basis for our numerical procedure.     

Localization of diffusion sources in complex networks with sparse observations

Locating sources in a large network is of paramount importance to reduce the spreading of disruptive behavior. Based on the backward diffusion-based method and integer programming, we propose an efficient approach to locate sources in complex networks with limited observers. The results on model networks and empirical networks demonstrate that, for a certain fraction of observers, the accuracy of our method for source localization will improve as the increase of network size. Besides, compared with the previous method (the maximum–minimum method), the performance of our method is much better with a small fraction of observers, especially in heterogeneous networks. Furthermore, our method is more robust against noise environments and strategies of choosing observers.     

Wavelet frame based scene reconstruction from range data

How to reconstruct the scene (a visible surface) from a set of scattered, noisy and possibly sparse range data is a challenging problem in robotic navigation and computer graphics. As most real scenes can be modeled by piecewise smooth surfaces, traditional surface fitting techniques (e.g. smoothing spline) generally can not preserve sharp discontinuities of surfaces. Based on sparse approximation of piecewise smooth functions in frame domain, we propose a new tight frame based formulation for reconstructing a piecewise smooth surface from a sparse range data set, which is robust to both additive noise and outliers. Furthermore, the resulting minimization problem from our formulation can be efficiently solved by the split Bregman method [1], [2]. The numerical experiments show that the proposed approach is capable of reconstructing a piecewise smooth surface with sharp edges from sparse range data corrupted with noise and outliers.     

Uniformly stable numerical schemes for the Boltzmann equation preserving the compressible Navier–Stokes asymptotics

In this paper, we develop a numerical method to solve Boltzmann like equations of kinetic theory which is able to capture the compressible Navier–Stokes dynamics at small Knudsen numbers. Our approach is based on the micro/macro decomposition technique, which applies to general collision operators. This decomposition is performed in all the phase space and leads to an equivalent formulation of the Boltzmann (or BGK) equation that couples a kinetic equation with macroscopic ones. This new formulation is then discretized with a semi-implicit time scheme combined with a staggered grid space discretization. Finally, several numerical tests are presented in order to illustrate the efficiency of our approach. Incidentally, we also introduce in this paper a modification of a standard splitting method that allows to preserve the compressible Navier–Stokes asymptotics in the case of the simplified BGK model. Up to our knowledge, this property is not known for general collision operators.     

A hybrid variational front tracking-level set mesh generator for problems exhibiting large deformations and topological changes

We present a method for generating 2-D unstructured triangular meshes that undergo large deformations and topological changes in an automatic way. We employ a method for detecting when topological changes are imminent via distance functions and shape skeletons. When a change occurs, we use a level set method to guide the change of topology of the domain mesh. This is followed by an optimization procedure, using a variational formulation of active contours, that seeks to improve boundary mesh conformity to the zero level contour of the level set function. Our method is advantageous for Arbitrary-Lagrangian–Eulerian (ALE) type methods and directly allows for using a variational formulation of the physics being modeled and simulated, including the ability to account for important geometric information in the model (such as for surface tension driven flow). Furthermore, the meshing procedure is not required at every time-step and the level set update is only needed during a topological change. Hence, our method does not significantly affect computational cost.     

Provably unconditionally stable,second-order time-accurate,mixed variational methods for phase-field models

We introduce provably unconditionally stable mixed variational methods for phase-field models. Our formulation is based on a mixed finite element method for space discretization and a new second-order accurate time integration algorithm. The fully-discrete formulation inherits the main characteristics of conserved phase dynamics, namely, mass conservation and nonlinear stability with respect to the free energy. We illustrate the theory with the Cahn–Hilliard equation, but our method may be applied to other phase-field models. We also propose an adaptive time-stepping version of the new time integration method. We present some numerical examples that show the accuracy, stability and robustness of the new method.     

A novel video-based microsphere localization algorithm for low contrast silica particles under white light illumination

On the basis of a brief review of four common image recognition algorithms for microspheres made of polystyrene or melamine resin, we present a new microsphere localization method for low-contrast silica beads under white light illumination. We compare both the polystyrene and silica procedures with respect to accuracy and precision by means of an optical tweezers setup providing CMOS video microscopy capability. By that we demonstrate that our new silica algorithm achieves a relative position uncertainty of less than ±1 nm for micron-sized microspheres, significantly exceeding the precision of the other silica approaches studied. Second, we present an advancement of our single microsphere tracking method to scenarios where two polystyrene, melamine resin or silica microspheres are in close-to-contact proximity. While the majority of the analysis algorithms studied generate artefacts due to interference effects under these conditions, we show that our new approach yields accurate and precise results.     

On the problem of reality conditions for the Ashtekar formulation

In this paper we discuss the reality conditions for Lorentzian and Euclidean gravity in the Ashtekar formulation by introducing a double conformal transformation.We generalize Marugan‘s results and demonstrate that the values of the double conformal factor have to be either real or double complex numbers.Either Lorentzian or Enclidean gravitational theory is up to the different values of the double conformal factor.Furthermore,the reality conditions of Lorentzian and Euclidean gravitational theory can be expressed in a unified way be use of the double complex function method.     

A fourth-order auxiliary variable projection method for zero-Mach number gas dynamics

A fourth-order numerical method for the zero-Mach-number limit of the equations for compressible flow is presented. The method is formed by discretizing a new auxiliary variable formulation of the conservation equations, which is a variable density analog to the impulse or gauge formulation of the incompressible Euler equations. An auxiliary variable projection method is applied to this formulation, and accuracy is achieved by combining a fourth-order finite-volume spatial discretization with a fourth-order temporal scheme based on spectral deferred corrections. Numerical results are included which demonstrate fourth-order spatial and temporal accuracy for non-trivial flows in simple geometries.