Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

Accurate solution of the boundary problem for collisionless gas in the field of a plane gravitational wave

An accurate solution of the Cauchy problem is found for a general-relativity collisionless kinetic equation against the background of the metric of a nonlinear plane gravitational field with an arbitrary law of gas-particle reflection from a boundary of any specified form. It is shown that in the field of a gravitational wave, interaction of the gas with the boundary necessarily leads to the appearance of shock waves.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 12, pp. 41–45, December, 1985.It remains to thank the participants of the gravitation seminar at Kazan' State University for discussions of the work.     

Open and traction boundary conditions for the incompressible Navier–Stokes equations

We present numerical schemes for the incompressible Navier–Stokes equations (NSE) with open and traction boundary conditions. We use pressure Poisson equation (PPE) formulation and propose new boundary conditions for the pressure on the open or traction boundaries. After replacing the divergence free constraint by this pressure Poisson equation, we obtain an unconstrained NSE. For Stokes equation with open boundary condition on a simple domain, we prove unconditional stability of a first order semi-implicit scheme where the pressure is treated explicitly and hence is decoupled from the computation of velocity. Using either boundary condition, the schemes for the full NSE that treat both convection and pressure terms explicitly work well with various spatial discretizations including spectral collocation and C⁰ finite elements. Moreover, when Reynolds number is of O(1) and when the first order semi-implicit time stepping is used, time step size of O(1) is allowed in benchmark computations for the full NSE. Besides standard stability and accuracy check, various numerical results including flow over a backward facing step, flow past a cylinder and flow in a bifurcated tube are reported. Numerically we have observed that using PPE formulation enables us to use the velocity/pressure pairs that do not satisfy the standard inf–sup compatibility condition. Our results extend that of Johnston and Liu [H. Johnston, J.-G. Liu, Accurate, stable and efficient Navier–Stokes solvers based on explicit treatment of the pressure term. J. Comp. Phys. 199 (1) (2004) 221–259] which deals with no-slip boundary conditions only.     

移动边界问题的时空域正则区域迭代配点法

考虑热传导方程的移动边界问题,其定解区域随着时间而变化.构造一种时空域上的高精度数值算法求解1+1维移动边界问题.在时空域上假设一个初始移动边界位置,构成移动边界问题的不规则计算区域,选择一个适当的正则区域(矩形区域)完全覆盖所计算的不规则区域,在正则区域上利用移动边界约束条件和固定边界条件,采用时空域重心插值配点法求...     

Half-Range Expansions for an Astrophysical Problem

Problems of spectral analysis are studied for an indefinite singular boundary value problem coming from astrophysical theory of particle acceleration around shocks. This leads to a nonclassical initial-boundary value problem for a partial differential equation that can bereduced by separation of variables to an indefinite Sturm–Liouville problem for which we establish Riesz basis properties of the eigen- and associated functions and formulate completeness and expansion theorems.     

Method of solving diffusion boundary-value problems for a region with a boundary moving in accordance with an arbitrary law

A general method is proposed for solving the boundary-value problem of the diffusion equation in a limited region with a boundary that moves in accordance with an arbitrary law. The method is used to solve the first linear diffusion problem. Other boundary-value problems can be solved in similar fashion.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 12, pp. 97–101, December, 1970.     

Wave Functions for Arbitrary Operator Ordering in the de Sitter Minisuperspace Approximation

We derive exact series solutions for the Wheeler–DeWitt equation corresponding to a spatially closed Friedmann–Robertson–Walker universe with cosmological constant for arbitrary operator ordering of the scale factor of the universe. The resulting wave functions are those relevant to the approximation which has been widely used in two-dimensional minisuperspace models with an inflationary scalar field for the purpose of predicting the period of inflation which results from competing boundary condition proposals for the wave function of the universe. The problem that Vilenkin's tunneling wave function is not normalizable for general operator orderings, is shown to persist for other values of the spatial curvature, and when additional matter degrees of freedom such as radiation are included.     

Computation of discrete transparent boundary conditions for the 2D Helmholtz equation

We present a family of non-local transparent boundary conditions for the 2D Helmholtz equation. The whole domain, on which the Helmholtz equation is defined, is decomposed into an interior and an exterior domain. The corresponding interior Helmholtz problem is formulated as a variational problem in a standard manner, representing a boundary value problem, whereas the exterior problem is posed as an initial value problem in the radial variable. This problem is then solved approximately by means of the Laplace transformation. The derived boundary conditions are asymptotically correct, model inhomogeneous exterior domains and are simple to implement.     

On a (3 + 1)-dimensional solution of the Landau-Lifshitz equation

We obtain a new nonstationary solution of the nonlinear Landau-Liftshitz equation in three spatial dimensions for easy-axis ferromagnets and antiferromagnets in an external magnetic field. The solution has discrete axial symmetry and depending on the values of the free parameters, it can have the form of a multidomain periodic exchange structure or a domain boundary separating two semi-infinite regions of space. The singular behavior of the solution near the symmetry axis is pointed out.V. V. Kuibyshev State University, Tomsk. Institute of Applied Mathematics, Far-Eastern Branch, Russian Academy of Sciences. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 6, pp. 68–71, June, 1994.     

The moment method for boundary layer problems in Brownian motion theory

We apply Grad's moment method, with Hermite moments and Marshak-type boundary conditions, to several boundary layer problems for the Klein-Kramers equation, the kinetic equation for noninteracting Brownian particles, and study its convergence properties as the number of moments is increased. The errors in various quantities of physical interest decrease asymptotically as inverse powers of this number; the exponent is roughly three times as large as in an earlier variational method, based on an expansion in the exact boundary layer eigenfunctions. For the case of a fully absorbing wall (the Milne problem) we obtain full agreement with the recent exact solution of Marshall and Watson; the relevant slip coefficient, the Milne length, is reproduced with an accuracy better than 10^–6. We also consider partially absorbing walls, with specular or diffuse reflection of nonabsorbed particles. In the latter case we allow for a temperature difference between the wall and the medium in which the particles move. There is noa priori reason why our method should work only for Brownian dynamics; one may hope to extend it to a broad class of linear transport equations. As a first test, we looked at the Milne problem for the BGK equation. In spite of the completely different analytic structure of the boundary layer eigenfunctions, the agreement with the exact solution is almost as good as for the Klein-Kramers equation.     

Radiation Reaction of the Classical Off-Shell Relativistic Charged Particle

It has been shown by Gupta and Padmanabhan that the radiation reaction force of the Abraham–Lorentz–Dirac equation can be obtained by a coordinate transformation from the inertial frame of an accelerating charged particle to that of the laboratory. We show that the problem may be formulated in a flat space of five dimensions, with five corresponding gauge fields in the framework of the classical version of a fully gauge covariant form of the Stueckelberg–Feynman–Schwinger covariant mechanics (the zero mode fields of the 0, 1, 2, 3 components correspond to the Maxwell fields). Without additional constraints, the particles and fields are not confined to their mass shells. We show that in the mass-shell limit, the generalized Lorentz force obtained by means of the retarded Green's functions for the five dimensional field equations provides the classical Abraham–Lorentz–Dirac radiation reaction terms (with renormalized mass and charge). We also obtain general coupled equations for the orbit and the off-shell dynamical mass during the evolution as well as an autonomous non-linear equation of third order for the off-shell mass. The theory does not admit radiation if the particle does not move off-shell. The structure of the equations implies that mass-shell deviation is bounded when the external field is removed.     

Numerical Simulation of Grain-Boundary Grooving by Level Set Method

A numerical investigation of grain-boundary grooving by means of a level set method is carried out. An idealized polycrystalline interconnect which consists of grains separated by parallel grain boundaries aligned normal to the average orientation of the surface is considered. Initially, the surface diffusion is the only physical mechanism assumed. The surface diffusion is driven by surface-curvature gradients, while a fixed surface slope and zero atomic flux are assumed at the groove root. The corresponding mathematical system is an initial boundary value problem for a two-dimensional equation of Hamilton–Jacobi type. The results obtained are in good agreement with both Mullins analytical “small-slope” solution of the linearized problem (W. W. Mullins, 1957, j. Appl. Phys. 28, 333) (for the case of an isolated grain boundary) and with the solution for a periodic array of grain boundaries (S. A. Hackney, 1988, Scripta Metall. 22, 1731). Incorporation of an electric field changes the problem to one of electromigration. Preliminary results of electromigration drift velocity simulations in copper lines are presented and discussed.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

Integral equation methods for elliptic problems with boundary conditions of mixed type

Laplace’s equation with mixed boundary conditions, that is, Dirichlet conditions on parts of the boundary and Neumann conditions on the remaining contiguous parts, is solved on an interior planar domain using an integral equation method. Rapid execution and high accuracy is obtained by combining equations which are of Fredholm’s second kind with compact operators on almost the entire boundary with a recursive compressed inverse preconditioning technique. Then an elastic problem with mixed boundary conditions is formulated and solved in an analogous manner and with similar results. This opens up for the rapid and accurate solution of several elliptic problems of mixed type.     

A free-surface hydrodynamic model for density-stratified flow in the weakly to strongly non-hydrostatic regime

A non-hydrostatic density-stratified hydrodynamic model with a free surface has been developed from the vorticity equations rather than the usual momentum equations. This approach has enabled the model to be obtained in two different forms, weakly non-hydrostatic and fully non-hydrostatic, with the computationally efficient weakly non-hydrostatic form applicable to motions having horizontal scales greater than the local water depth. The hydrodynamic model in both its weakly and fully non-hydrostatic forms is validated numerically using exact nonlinear non-hydrostatic solutions given by the Dubriel–Jacotin–Long equation for periodic internal gravity waves, internal solitary waves, and flow over a ridge. The numerical code is developed based on a semi-Lagrangian scheme and higher order finite-difference spatial differentiation and interpolation. To demonstrate the applicability of the model to coastal ocean situations, the problem of tidal generation of internal solitary waves at a shelf-break is considered. Simulations carried out with the model obtain the evolution of solitary wave generation and propagation consistent with past results. Moreover, the weakly non-hydrostatic simulation is shown to compare favorably with the fully non-hydrostatic simulation. The capability of the present model to simulate efficiently relatively large scale non-hydrostatic motions suggests that the weakly non-hydrostatic form of the model may be suitable for application in a large-area domain while the computationally intensive fully non-hydrostatic form of the model may be used in an embedded sub-domain where higher resolution is needed.     

Adaptive artificial boundary condition for the two-level Schrödinger equation with conical crossings

In this paper, we present an adaptive approach to design the artificial boundary conditions for the two-level Schrödinger equation with conical crossings on the unbounded domain. We use the windowed Fourier transform to obtain the local wave number information in the vicinity of artificial boundaries, and adopt the operator splitting method to obtain an adaptive local artificial boundary condition. Then reduce the original problem into an initial boundary value problem on the bounded computational domain, which can be solved by the finite difference method. By this numerical method, we observe the surface hopping phenomena of the two-level Schrödinger equation with conical crossings. Several numerical examples are provided to show the accuracy and convergence of the proposed method.     

Approximate method for calculating the temperature field in a growing crystal with account of thermal radiation

An engineering method is described for solving the heat-conduction problem in a growing single crystal for the case in which heat is lost at the surface through radiation. An integral substitution is used to transfer the nonlinearity in the boundary condition into the basic equation for the process; the nonlinear complex in the transformed equation can be evaluated. The approximate solution is compared with the results of a computer numerical calculation.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii Fizika, No. 2, pp. 89–95, February, 1971.     

Diffraction of a Plane Electromagnetic Wave by an Infinite Double-Ridge Wedge with a Dielectric Cylinder at the Crest

We reduce the considered problem to solving a matrix equation of the second kind for unknown coefficients of expansion of a diffracted field into a Fourier–Bessel series. This expansion was obtained by imposing boundary conditions on the diffracted field with the subsequent re-expansion of the field function over basis functions in a given interval. The expansion coefficients were determined analytically in the case where the electric diameter of the cylinder is less than unity as well as numerically with a high accuracy by solving the obtained matrix equation using the reduction method. We derived expressions for the pattern of the far-zone field scattered by the studied structure and the backscattering cross section and give exact numerical results for the case of an E-polarized incident wave.