We prove that for nonlinear elastic energies with strong enough energetic control of the outer distortion of admissible deformations, almost everywhere global invertibility as constraint can be obtained in the $\Gamma$-limit of the elastic energy with an added nonlocal self-repulsion term with asymptocially vanishing coefficient. The self-repulsion term considered here formally coincides with a Sobolev–Slobodeckiĭ seminorm of the inverse deformation. Variants near the boundary or on the surface of the domain are also studied.

This article surveys research on the application of compatible finite element methods to large-scale atmosphere and ocean simulation. Compatible finite element methods extend Arakawa’s C-grid finite difference scheme to the finite element world. They are constructed from a discrete de Rham complex, which is a sequence of finite element spaces linked by the operators of differential calculus. The use of discrete de Rham complexes to solve partial differential equations is well established, but in this article we focus on the specifics of dynamical cores for simulating weather, oceans and climate. The most important consequence of the discrete de Rham complex is the Hodge–Helmholtz decomposition, which has been used to exclude the possibility of several types of spurious oscillations from linear equations of geophysical flow. This means that compatible finite element spaces provide a useful framework for building dynamical cores. In this article we introduce the main concepts of compatible finite element spaces, and discuss their wave propagation properties. We survey some methods for discretizing the transport terms that arise in dynamical core equation systems, and provide some example discretizations, briefly discussing their iterative solution. Then we focus on the recent use of compatible finite element spaces in designing structure preserving methods, surveying variational discretizations, Poisson bracket discretizations and consistent vorticity transport.

The aim of this paper is to provide and numerically test in the presence of measurement noise a procedure for target classification in wave imaging based on comparing frequency-dependent distribution descriptors with precomputed ones in a dictionary of learned distributions. Distribution descriptors for inhomogeneous objects are obtained from the scattering coefficients. First, we extract the scattering coefficients of the (inhomogeneous) target from the perturbation of the reflected waves. Then, for a collection of inhomogeneous targets, we build a frequency-dependent dictionary of distribution descriptors and use a matching algorithm in order to identify a target from the dictionary up to some translation, rotation and scaling.

We present a JASMIN-based two-dimensional parallel implementation of an adaptive combined preconditioner for the solution of linear problems arising in the finite volume discretisation of one-group and multi-group radiation diffusion equations. We first propose the attribute of patch-correlation for cells of a two-dimensional monolayer piecewise rectangular structured grid without any suspensions based on the patch hierarchy of JASMIN, classify and reorder these cells via their attributes, and derive the conversion of cell-permutations. Using two cell-permutations, we then construct some parallel incomplete LU factorisation and substitution algorithms, to provide our parallel -GMRES solver with the help of the default BoomerAMG in the HYPRE library. Numerical results demonstrate that our proposed parallel incomplete LU preconditioner (ILU) is of higher efficiency than the counterpart in the Euclid library, and that the proposed parallel -GMRES solver is more robust and more efficient than the default BoomerAMG-GMRES solver.

A second order Ghost Fluid method is proposed for the treatment of interface problems of elliptic equations with discontinuous coefficients. By appropriate use of auxiliary virtual points, physical jump conditions are enforced at the interface. The signed distance function is used for the implicit description of irregular domain. With the additional unknowns, high order approximation considering the discontinuity can be built. To avoid the ill-conditioned matrix, the interpolation stencils are selected adaptively to balance the accuracy and the numerical stability. Additional equations containing the jump restrictions are assembled with the original discretized algebraic equations to form a new sparse linear system. Several Krylov iterative solvers are tested for the newly derived linear system. The results of a series of 1-D, 2-D tests show that the proposed method possesses second order accuracy in L∞ norm. Besides, the method can be extended to the 3-D problems straightforwardly. Numerical results reveal the present method is highly efficient and robust in dealing with the interface problems of elliptic equations.

The computation efficiency of high dimensional (3D and 4D) B-spline interpolation, constructed by classical tensor product method, is improved greatly by precomputing the B-spline function. This is due to the character of NLT code, i.e. only the linearised characteristics are needed so that the unperturbed orbit as well as values of the B-spline function at interpolation points can be precomputed at the beginning of the simulation. By integrating this fixed point interpolation algorithm into NLT code, the high dimensional gyro-kinetic Vlasov equation can be solved directly without operator splitting method which is applied in conventional semi-Lagrangian codes. In the Rosenbluth-Hinton test, NLT runs a few times faster for Vlasov solver part and converges at about one order larger time step than conventional splitting code.

In this article, we study theoretically and numerically the interaction of a vortex induced by a rotating cylinder with a perpendicular plane. We show the existence of weak solutions to the swirling vortex models by using the Hopf extension method, and by an elegant contradiction argument, respectively. We demonstrate numerically that the model could produce phenomena of swirling vortex including boundary layer pumping and two-celled vortex that are observed in potential line vortex interacting with a plane and in a tornado.

In this paper we present a fast and accurate numerical algorithm for the computation of hyperspherical Bessel functions of large order and real arguments. For the hyperspherical Bessel functions of closed type, no stable algorithm existed so far due to the lack of a backwards recurrence. We solved this problem by establishing a relation to Gegenbauer polynomials. All our algorithms are written in C and are publicly available at Github [https://github.com/lesgourg/class_public]. A Python wrapper is available upon request.

We introduce PyCFTBoot, a wrapper designed to reduce the barrier to entry in conformal bootstrap calculations that require semidefinite programming. Symengine and SDPB are used for the most intensive symbolic and numerical steps respectively. After reviewing the built-in algorithms for conformal blocks, we explain how to use the code through a number of examples that verify past results. As an application, we show that the multi-correlator bootstrap still appears to single out the Wilson-Fisher fixed points as special theories in dimensions between 3 and 4 despite the recent proof that they violate unitarity.

In magnetically confined plasma research, the understandings of small and large perturbations at equilibrium are both critical for plasma controlling and steady state operation. Numerical simulations using original MHD model can hardly give clear picture for small perturbations, while non-conservative perturbed MHD model may break conservation law, and give unphysical results when perturbations grow large after long-time computation. In this paper, we present a nonlinear conservative perturbed MHD model by splitting primary variables in original MHD equations into equilibrium part and perturbed part, and apply an approach in the framework of discontinuous Galerkin (DG) spatial discretization for numerical solutions. This enables high resolution of very small perturbations, and also gives satisfactory non-smooth solutions for large perturbations, which are both broadly concerned in magnetically confined plasma research. Numerical examples demonstrate satisfactory performance of the proposed model clearly. For small perturbations, the results have higher resolution comparing with the original MHD model; for large perturbations, the non-smooth solutions match well with existing references, confirming reliability of the model for instability investigations in magnetically confined plasma numerical research.

How to effectively simulate the interaction between fluid and solid elements of different sizes remains to be challenging. The discrete element method (DEM) has been used to deal with the interactions between solid elements of various shapes and sizes, while the material point method (MPM) has been developed to handle the multiphase (solid-liquid-gas) interactions involving failure evolution. A combined MPM-DEM procedure is proposed to take advantage of both methods so that the interaction between solid elements and fluid particles in a container could be better simulated. In the proposed procedure, large solid elements are discretized by the DEM, while the fluid motion is computed using the MPM. The contact forces between solid elements and rigid walls are calculated using the DEM. The interaction between solid elements and fluid particles are calculated via an interfacial scheme within the MPM framework. With a focus on the boundary condition effect, the proposed procedure is illustrated by representative examples, which demonstrates its potential for a certain type of engineering problems.

A Lagrangian surface hopping algorithm is implemented to study the two dimensional massless Dirac equation for Graphene with an electrostatic potential, in the semiclassical regime. In this problem, the crossing of the energy levels of the system at Dirac points requires a particular treatment in the algorithm in order to describe the quantum transition—characterized by the Landau-Zener probability— between different energy levels. We first derive the Landau-Zener probability for the underlying problem, then incorporate it into the surface hopping algorithm. We also show that different asymptotic models for this problem derived in [O. Morandi, F. Schurrer, J. Phys. A:Math. Theor. 44 (2011) 265301]may give different transition probabilities. We conduct numerical experiments to compare the solutions to the Dirac equation, the surface hopping algorithm, and the asymptotic models of [O. Morandi, F. Schurrer, J. Phys. A: Math. Theor. 44 (2011) 265301].

In this paper, a hybrid dynamic mesh generation method for multi-block structured grid is presented based on inverse distance weighting (IDW) interpolation and transfinite interpolation (TFI). The major advantage of the algorithm is that it maintains the effectiveness of TFI, while possessing the ability to deal with multi-block structured grid from the IDW method. In this approach, dynamic mesh generation is made in two steps. At first, all domain vertexes with known deformation are selected as sample points and IDW interpolation is applied to get the grid deformation on domain edges. Then, an arc-length-based TFI is employed to efficiently calculate the grid deformation on block faces and inside each block. The present approach can be well applied to both two-dimensional (2D) and three-dimensional (3D) problems. The proposed method has been well-validated by several test cases. Numerical results show that dynamic meshes with high quality can be generated in an accurate and efficient manner.

In this paper, a double-passage shape correction (DPSC) method is presented for simulation of unsteady flows around vibrating blades and aeroelastic prediction. Based on the idea of phase-lagged boundary conditions, the shape correction method was proposed aimed at efficiently dealing with unsteady flow problems in turbomachinery. However, the original single-passage shape correction (SPSC) may show the disadvantage of slow convergence of unsteady solutions and even produce nonphysical oscillation. The reason is found to be related with the disturbances on the circumferential boundaries that can not be damped by numerical schemes. To overcome these difficulties, the DPSC method is adopted here, in which the Fourier coefficients are computed from flow variables at implicit boundaries instead of circumferential boundaries in the SPSC method. This treatment actually reduces the interaction between the calculation of Fourier coefficients and the update of flow variables. Therefore a faster convergence speed could be achieved and also the solution stability is improved. The present method is developed to be suitable for viscous and turbulent flows. And for real three-dimensional (3D) problems, the rotating effects are also considered. For validation, a 2D oscillating turbine cascade, a 3D oscillating flat plate cascade and a 3D practical transonic fan rotor are investigated. Comparisons with experimental data or other solutions and relevant discussions are presented in detail. Numerical results show that the solution accuracy of DPSC method is favorable and at least comparable to the SPSC method. However, fewer iteration cycles are needed to get a converged and stable unsteady solution, which greatly improves the computational efficiency.

In this paper, we propose a new partitioned approach to compute fluid-structure interaction (FSI) by extending the original direct-forcing technique and integrating it with the immersed boundary method. The fluid and structural equations are calculated separately via their respective disciplinary algorithms, with the fluid motion solved by the immersed boundary method on a uniform Cartesian mesh and the structural motion solved by a finite element method, and their solution data only communicate at the fluid-structure interface. This computational framework is capable of handling FSI problems with sophisticated structures described by detailed constitutive laws. The proposed methods are thoroughly tested through numerical simulations involving viscous fluid flow interacting with rigid, elastic solid, and elastic thin-walled structures.

In this paper, we present accurate and economic integration quadratures for hypersingular functions over three simple geometric shapes in ℝ 3 (spheres, cubes, and cylinders). The quadrature nodes are made of the tensor-product of 1-D Gauss nodes on [–1,1] for non-periodic variables or uniform nodes on [0,2π] or [0,π] for periodic ones. The quadrature weights are converted from a brute-force integration of the hypersingular function through interpolating the smooth component of the integrand. Numerical results are presented to validate the accuracy and efficiency of computing hypersingular integrals, as in the computations of Cauchy principal values, with a minimum number of quadrature nodes. The pre-calculated quadrature tables can be then readily used to implement Nyström collocation methods of hypersingular volume integral equations such as the one for Maxwell equations.

The level set method is one of the most successful methods for the simulation of multi-phase flows. To keep the level set function close the signed distance function, the level set function is constantly reinitialized by solving a Hamilton-Jacobi type of equation during the simulation. When the fluid interface intersects with a solid wall, a moving contact line forms and the reinitialization of the level set function requires a boundary condition in certain regions on the wall. In this work, we propose to use the dynamic contact angle, which is extended from the contact line, as the boundary condition for the reinitialization of the level set function. The reinitialization equation and the equation for the normal extension of the dynamic contact angle form a coupled system and are solved simultaneously. The extension equation is solved on the wall and it provides the boundary condition for the reinitialization equation; the level set function provides the directions along which the contact angle is extended from the contact line. The coupled system is solved using the 3rd order TVD Runge-Kutta method and the Godunov scheme. The Godunov scheme automatically identifies the regions where the angle condition needs to be imposed. The numerical method is illustrated by examples in three dimensions.

In this paper, we develop a novel energy-preserving wavelet collocation method for solving general multi-symplectic formulations of Hamiltonian PDEs. Based on the autocorrelation functions of Daubechies compactly supported scaling functions, the wavelet collocation method is conducted for spatial discretization. The obtained semi-discrete system is shown to be a finite-dimensional Hamiltonian system, which has an energy conservation law. Then, the average vector field method is used for time integration, which leads to an energy-preserving method for multi-symplectic Hamiltonian PDEs. The proposed method is illustrated by the nonlinear Schrödinger equation and the Camassa-Holm equation. Since differentiation matrix obtained by the wavelet collocation method is a cyclic matrix, we can apply Fast Fourier transform to solve equations in numerical calculation. Numerical experiments show the high accuracy, effectiveness and conservation properties of the proposed method.

Multiscale modelling is a powerful technique, which allows for computational efficiency while retaining small-scale details when they are essential for understanding a finer behaviour of the studied system. In the case of materials modelling, one of the effective multiscaling concepts is domain partitioning, which implies the existence of an explicit interface between various material descriptions, for instance atomistic and continuum regions. When dynamic material behaviour is considered, the major problem for this technique is dealing with reflections of high frequency waves from the interface separating two scales. In this article, a new method is suggested, which overcomes this problem for the case of magnetisation dynamics. The introduction of a damping band at the interface between scales, which absorbs high frequency waves, is suggested. The idea is verified using a number of one-dimensional examples with fine/coarse scale discretisation of a continuum problem of spin wave propagation. This work is the first step towards establishing a reliable atomistic/continuum multiscale transition for the description of the evolution of magnetic properties of ferromagnets.

The main purpose of this work is to contrast and analyze a large time-stepping numerical method for the Swift-Hohenberg (SH) equation. This model requires very large time simulation to reach steady state, so developing a large time step algorithm becomes necessary to improve the computational efficiency. In this paper, a semi-implicit Euler schemes in time is adopted. An extra artificial term is added to the discretized system in order to preserve the energy stability unconditionally. The stability property is proved rigorously based on an energy approach. Numerical experiments are used to demonstrate the effectiveness of the large time-stepping approaches by comparing with the classical scheme.