A Nyström method for weakly singular integral operators on surfaces

We describe a modified Nyström method for the discretization of the weakly singular boundary integral operators which arise from the formulation of linear elliptic boundary value problems as integral equations. Standard Nyström and collocation schemes proceed by representing functions via their values at a collection of quadrature nodes. Our method uses appropriately scaled function values in lieu of such representations. This results in a scheme which is mathematically equivalent to Galerkin discretization in that the resulting matrices are related to those obtained by Galerkin methods via conjugation with well-conditioned matrices, but which avoids the evaluation of double integrals. Moreover, we incorporate a new mechanism for approximating the singular integrals which arise from the discretization of weakly singular integral operators which is considerably more efficient than standard methods. We illustrate the performance of our method with numerical experiments.     

A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov–Poisson equations

The Vlasov–Poisson equations describe the evolution of a collisionless plasma, represented through a probability density function (PDF) that self-interacts via an electrostatic force. One of the main difficulties in numerically solving this system is the severe time-step restriction that arises from parts of the PDF associated with moderate-to-large velocities. The dominant approach in the plasma physics community for removing these time-step restrictions is the so-called particle-in-cell (PIC) method, which discretizes the distribution function into a set of macro-particles, while the electric field is represented on a mesh. Several alternatives to this approach exist, including fully Lagrangian, fully Eulerian, and so-called semi-Lagrangian methods. The focus of this work is the semi-Lagrangian approach, which begins with a grid-based Eulerian representation of both the PDF and the electric field, then evolves the PDF via Lagrangian dynamics, and finally projects this evolved field back onto the original Eulerian mesh. In particular, we develop in this work a method that discretizes the 1 + 1 Vlasov–Poisson system via a high-order discontinuous Galerkin (DG) method in phase space, and an operator split, semi-Lagrangian method in time. Second-order accuracy in time is relatively easy to achieve via Strang operator splitting. With additional work, using higher-order splitting and a higher-order method of characteristics, we also demonstrate how to push this scheme to fourth-order accuracy in time. We show how to resolve all of the Lagrangian dynamics in such a way that mass is exactly conserved, positivity is maintained, and high-order accuracy is achieved. The Poisson equation is solved to high-order via the smallest stencil local discontinuous Galerkin (LDG) approach. We test the proposed scheme on several standard test cases.     

Nyström Method for the Coulomb and Screened Coulomb Potentials

In this paper, we show that the simple Nyström method can yield very accurate eigenvalues and eigenfunctions not only for large principal quantum number but also for large angular momentum quantum number. We demonstrate that the furcation phenomenon emerging in the calculated eigenfunctions can be regarded as an indicator for the bad behavior of the integral equation and the unreliability of the obtained results.     

A local energy-preserving scheme for Klein Gordon Schrdinger equations

A local energy conservation law is proposed for the Klein–Gordon–Schr ¨odinger equations, which is held in any local time–space region. The local property is independent of the boundary condition and more essential than the global energy conservation law. To develop a numerical method preserving the intrinsic properties as much as possible, we propose a local energy-preserving(LEP) scheme for the equations. The merit of the proposed scheme is that the local energy conservation law can hold exactly in any time–space region. With the periodic boundary conditions, the scheme also possesses the discrete change and global energy conservation laws. A nonlinear analysis shows that the LEP scheme converges to the exact solutions with order O(τ2+ h2). The theoretical properties are verified by numerical experiments.     

Darboux transformations in integral and differential forms for generalized Schrödinger equations

The Darboux transformation operator technique is applied to the generalized Schrödinger equation with a position-dependent effective mass and with linearly energy-dependent potentials. Intertwining operators are obtained in an explicit form and used for constructing generalized Darboux transformations. An interrelation is established between the differential and integral transformation operators. It is shown how to construct the quantum well potentials in nanoelectronic with a given spectrum.     

A note on the integral formulation of Einstein’s equations induced on a braneworld

We revisit the integral formulation (or Greens function approach) of Einsteins equations in the context of braneworlds. The integral formulation has been proposed independently by several authors in the past, based on the assumption that it possible to give a reinterpretation of the local metric field in curved spacetimes as an integral expression involving sources and boundary conditions. This allows one to separate source-generated and source-free contributions to the metric field. As a consequence, an exact meaning to Machs Principle can be achieved in the sense that only source-generated (matter fields) contributions to the metric are allowed for; universes which do not obey this condition would be non-Machian. In this paper, we revisit this idea concentrating on a Randall–Sundrum-type model with a non-trivial cosmology on the brane. We argue that the role of the surface term (the source-free contribution) in the braneworld scenario may be quite subtler than in the 4D formulation. This may pose, for instance, an interesting issue to the cosmological constant problem.     

A mass,energy, vorticity,and potential enstrophy conserving lateral fluid–land boundary scheme for the shallow water equations

A numerical scheme for treating fluid–land boundaries in inviscid shallow water flows is derived that conserves the domain-summed mass, energy, vorticity, and potential enstrophy in domains with arbitrarily shaped boundaries. The boundary scheme is derived from a previous scheme that conserves all four domain-summed quantities only in periodic domains without boundaries. It consists of a method for including land in the model along with evolution equations for the vorticity and extrapolation formulas for the depth at fluid–land boundaries. Proofs of mass, energy, vorticity, and potential enstrophy conservation are given. Numerical simulations are carried out demonstrating the conservation properties and accuracy of the boundary scheme for inviscid flows and comparing its performance with that of four alternative boundary schemes. The first of these alternatives extrapolates or finite-differences the velocity to obtain the vorticity at boundaries; the second enforces the free-slip boundary condition; the third enforces the super-slip condition; and the fourth enforces the no-slip condition. Comparisons of the conservation properties demonstrate that the new scheme is the only one of the five that conserves all four domain-summed quantities, and it is the only one that both prevents a spurious energy cascade to the smallest resolved scales and maintains the correct flow orientation with respect to an external forcing. Comparisons of the accuracy demonstrate that the new scheme generates vorticity fields that have smaller errors than those generated by any of the alternative schemes, and it generates depth and velocity fields that have errors about equal to those in the fields generated by the most accurate alternative scheme.     

A new proposal for the discretization of the Griffin—Wheeler—Hartree—Fock equations

A polynomial expansion is proposed as a new way to discretize the Griffin-Wheeler-Hartree-Fock equations of the Generator Coordinate Hartree-Fock method. The implementation of the polynomial expansion in the Generator Coordinate Hartree-Fock method discretizes the Griffin-Wheeler-Hartree-Fock equations through a numerical mesh which is not equally spaced. This procedure makes the optimization of Gaussian exponents in the Generator Coordinate Hartree-Fock method more flexible and more efficient. The results obtained with the polynomial expansion for atomic Hartree-Fock energies show this technique is very powerful when employed in the design of compact and high accurate Gaussian basis sets used in ab initio non-relativistic (Hartree-Fock) and relativistic (Dirac-Fock) calculations.     

A positivity-preserving ALE finite element scheme for convection–diffusion equations in moving domains

A new high-resolution scheme is developed for convection–diffusion problems in domains with moving boundaries. A finite element approximation of the governing equation is designed within the framework of a conservative Arbitrary Lagrangian Eulerian (ALE) formulation. An implicit flux-corrected transport (FCT) algorithm is implemented to suppress spurious undershoots and overshoots appearing in convection-dominated problems. A detailed numerical study is performed for P₁ finite element discretizations on fixed and moving meshes. Simulation results for a Taylor dispersion problem (moderate Peclet numbers) and for a convection-dominated problem (large Peclet numbers) are presented to give a flavor of practical applications.     

A semi-discrete central scheme for magnetohydrodynamics on orthogonal–curvilinear grids

This work describes a novel scheme for the equations of magnetohydrodynamics on orthogonal–curvilinear grids within a finite-volume framework. The scheme is based on a combination of central-upwind techniques for hyperbolic conservation laws and projection–evolution methods originally developed for Hamilton–Jacobi equations. The scheme is derived in semi-discrete form, and a full-fledged version is obtained by applying any stable and accurate solver for integration in time. The divergence-free condition of the magnetic field is a built-in property of the scheme by virtue of a constrained-transport ansatz for the induction equation. From the general formulation second-order accurate schemes for cylindrical grids and spherical grids are introduced in some more detail pointing out their potential importance in many applications. Special emphasis in this context is put to a treatment of the geometric axis implying severe complications because of the presence of coordinate singularities and associated grid degeneracy. An attempt to tackle these problems is presented. Numerical experiments illustrate the overall robustness and performance of the scheme for a small suite of tests.     

Implicit LU-SGS algorithm for high-order methods on unstructured grid with p-multigrid strategy for solving the steady Navier–Stokes equations

The fluid dynamic equations are discretized by a high-order spectral volume (SV) method on unstructured tetrahedral grids. We solve the steady state equations by advancing in time using a backward Euler (BE) scheme. To avoid the inversion of a large matrix we approximate BE by an implicit lower–upper symmetric Gauss–Seidel (LU-SGS) algorithm. The implicit method addresses the stiffness in the discrete Navier–Stokes equations associated with stretched meshes. The LU-SGS algorithm is then used as a smoother for a p-multigrid approach. A Von Neumann stability analysis is applied to the two-dimensional linear advection equation to determine its damping properties. The implicit LU-SGS scheme is used to solve the two-dimensional (2D) compressible laminar Navier–Stokes equations. We compute the solution of a laminar external flow over a cylinder and around an airfoil at low Mach number. We compare the convergence rates with explicit Runge–Kutta (E-RK) schemes employed as a smoother. The effects of the cell aspect ratio and the low Mach number on the convergence are investigated. With the p-multigrid method and the implicit smoother the computational time can be reduced by a factor of up to 5–10 compared with a well tuned E-RK scheme.     

Correspondence-principle condition on the Reissner-Nordström solution for charged leptons

Current experimental evidence supports the picture of three charged leptons (e ^±, ^±, ^±), a triplet with mass splitting that bears a resemblance to the Gell-Mann-Okubo form. To elucidate the overall mass scale, a charged lepton is viewed as a mass point that engenders a local Reissner-Nordström spacetime geometry, and the Einstein-Maxwell action is evaluated through an invariant space-time region associated with the particle's radiation reaction interval 2e ²/3m. What emerge for the electron, muon, and tau are values of the Einstein-Maxwell action in the neighborhood of. The mean value of the three action integrals is =(140.93)e ²=(1.0284), and thus the apparent quantum condition sets the mass scale for the three charged leptons.     

Integrable discretisations for a class of nonlinear Schrödinger equations on Grassmann algebras

Integrable discretisations for a class of coupled (super) nonlinear Schrödinger (NLS) type of equations are presented. The class corresponds to a Lax operator with entries in a Grassmann algebra. Elementary Darboux transformations are constructed. As a result, Grassmann generalisations of the Toda lattice and the NLS dressing chain are obtained. The compatibility (Bianchi commutativity) of these Darboux transformations leads to integrable Grassmann generalisations of the difference Toda and NLS equations. The resulting systems will have discrete Lax representations provided by the set of two consistent elementary Darboux transformations. For the two discrete systems obtained, initial value and initial-boundary problems are formulated.     

A matched interface and boundary method for solving multi-flow Navier–Stokes equations with applications to geodynamics

We have developed a second-order numerical method, based on the matched interface and boundary (MIB) approach, to solve the Navier–Stokes equations with discontinuous viscosity and density on non-staggered Cartesian grids. We have derived for the first time the interface conditions for the intermediate velocity field and the pressure potential function that are introduced in the projection method. Differentiation of the velocity components on stencils across the interface is aided by the coupled fictitious velocity values, whose representations are solved by using the coupled velocity interface conditions. These fictitious values and the non-staggered grid allow a convenient and accurate approximation of the pressure and potential jump conditions. A compact finite difference method was adopted to explicitly compute the pressure derivatives at regular nodes to avoid the pressure–velocity decoupling. Numerical experiments verified the desired accuracy of the numerical method. Applications to geophysical problems demonstrated that the sharp pressure jumps on the clast-Newtonian matrix are accurately captured for various shear conditions, moderate viscosity contrasts and a wide range of density contrasts. We showed that large transfer errors will be introduced to the jumps of the pressure and the potential function in case of a large absolute difference of the viscosity across the interface; these errors will cause simulations to become unstable.     

Absorbing boundary conditions for nonlinear Euler and Navier–Stokes equations based on the perfectly matched layer technique

Absorbing boundary conditions for the nonlinear Euler and Navier–Stokes equations in three space dimensions are presented based on the perfectly matched layer (PML) technique. The derivation of equations follows a three-step method recently developed for the PML of linearized Euler equations. To increase the efficiency of the PML, a pseudo mean flow is introduced in the formulation of absorption equations. The proposed PML equations will absorb exponentially the difference between the nonlinear fluctuation and the prescribed pseudo mean flow. With the nonlinearity in flux vectors, the proposed nonlinear absorbing equations are not formally perfectly matched to the governing equations as their linear counter-parts are. However, numerical examples show satisfactory results. Furthermore, the nonlinear PML reduces automatically to the linear PML upon linearization about the pseudo mean flow. The validity and efficiency of proposed equations as absorbing boundary conditions for nonlinear Euler and Navier–Stokes equations are demonstrated by numerical examples.     

An efficient method for the incompressible Navier–Stokes equations on irregular domains with no-slip boundary conditions,high order up to the boundary

Common efficient schemes for the incompressible Navier–Stokes equations, such as projection or fractional step methods, have limited temporal accuracy as a result of matrix splitting errors, or introduce errors near the domain boundaries (which destroy uniform convergence to the solution). In this paper we recast the incompressible (constant density) Navier–Stokes equations (with the velocity prescribed at the boundary) as an equivalent system, for the primary variables velocity and pressure. equation for the pressure. The key difference from the usual approaches occurs at the boundaries, where we use boundary conditions that unequivocally allow the pressure to be recovered from knowledge of the velocity at any fixed time. This avoids the common difficulty of an, apparently, over-determined Poisson problem. Since in this alternative formulation the pressure can be accurately and efficiently recovered from the velocity, the recast equations are ideal for numerical marching methods. The new system can be discretized using a variety of methods, including semi-implicit treatments of viscosity, and in principle to any desired order of accuracy. In this work we illustrate the approach with a 2-D second order finite difference scheme on a Cartesian grid, and devise an algorithm to solve the equations on domains with curved (non-conforming) boundaries, including a case with a non-trivial topology (a circular obstruction inside the domain). This algorithm achieves second order accuracy in the L^∞ norm, for both the velocity and the pressure. The scheme has a natural extension to 3-D.