High-Order Nonreflecting Boundary Conditions without High-Order Derivatives

A wave problem in an unbounded domain is often treated numerically by truncating the infinite domain via an artificial boundary , imposing a so-called nonreflecting boundary condition (NRBC) on , and then solving the problem numerically in the finite domain bounded by . A general approach is devised here to construct high-order local NRBCs with a symmetric structure and with only low (first- or second-) order spatial and/or temporal derivatives. This enables the practical use of NRBCs of arbitrarily high order. In the case of time-harmonic waves with finite element discretization, the approach yields a symmetric C⁰ finite element formulation in which standard elements can be employed. The general methodology is presented for both the time-harmonic case (Helmholtz equation) and the time-dependent case (the wave equation) and is demonstrated numerically in the former case.     

The Three Dimensional Non-conforming Finite Element Solution of the Chapman–Ferraro Problem

We demonstrate the feasibility of using a non-conforming, piecewise harmonic finite element method on an unstructured grid in solving a magnetospheric physics problem. We use this approach to construct a global discrete model of the magnetic field of the magnetosphere that includes the effects of shielding currents at the outer boundary (the magnetopause). As in the approach of F. R. Toffolettoet al.(1994,Geophys. Res. Lett.21, 7) the internal magnetospheric field model is that of R. V. Hilmer and G.-H. Voigt (1995,J. Geophys. Res.) while the magnetopause shape is based on an empirically determined approximation (1997, J. Shueet al.,J. Geophys. Res.102, 9497). The results is a magnetic field model whose field lines are completely confined within the magnetosphere. The presented numerical results indicate that the discrete non-conforming finite element model is well-suited for magnetospheric field modeling.     

Physics-Based Preconditioning and the Newton–Krylov Method for Non-equilibrium Radiation Diffusion

An algorithm is presented for the solution of the time dependent reaction-diffusion systems which arise in non-equilibrium radiation diffusion applications. This system of nonlinear equations is solved by coupling three numerical methods, Jacobian-free Newton–Krylov, operator splitting, and multigrid linear solvers. An inexact Newton's method is used to solve the system of nonlinear equations. Since building the Jacobian matrix for problems of interest can be challenging, we employ a Jacobian–free implementation of Newton's method, where the action of the Jacobian matrix on a vector is approximated by a first order Taylor series expansion. Preconditioned generalized minimal residual (PGMRES) is the Krylov method used to solve the linear systems that come from the iterations of Newton's method. The preconditioner in this solution method is constructed using a physics-based divide and conquer approach, often referred to as operator splitting. This solution procedure inverts the scalar elliptic systems that make up the preconditioner using simple multigrid methods. The preconditioner also addresses the strong coupling between equations with local 2×2 block solves. The intra-cell coupling is applied after the inter-cell coupling has already been addressed by the elliptic solves. Results are presented using this solution procedure that demonstrate its efficiency while incurring minimal memory requirements.     

New High-Resolution Semi-discrete Central Schemes for Hamilton–Jacobi Equations

We introduce a new high-resolution central scheme for multidimensional Hamilton–Jacobi equations. The scheme retains the simplicity of the non-oscillatory central schemes developed by C.-T. Lin and E. Tadmor (in press, SIAM J. Sci. Comput.), yet it enjoys a smaller amount of numerical viscosity, independent of 1/Δt. By letting Δt↓0 we obtain a new second-order central scheme in the particularly simple semi-discrete form, along the lines of the new semi-discrete central schemes recently introduced by the authors in the context of hyperbolic conservation laws. Fully discrete versions are obtained with appropriate Runge–Kutta solvers. The smaller amount of dissipation enables efficient integration of convection-diffusion equations, where the accumulated error is independent of a small time step dictated by the CFL limitation. The scheme is non-oscillatory thanks to the use of nonlinear limiters. Here we advocate the use of such limiters on second discrete derivatives, which is shown to yield an improved high resolution when compared to the usual limitation of first derivatives. Numerical experiments demonstrate the remarkable resolution obtained by the proposed new central scheme.     

A Roe Scheme for Ideal MHD Equations on 2D Adaptively Refined Triangular Grids

In this paper we present a second order finite volume method for the resolution of the bidimensional ideal MHD equations on adaptively refined triangular meshes. Our numerical flux function is based on a multidimensional extension of the Roe scheme proposed by Cargo and Gallice for the 1D MHD system. If the mesh is only composed of triangles, our scheme is proved to be weakly consistent with the condition …B=0. This property fails on a cartesian grid. The efficiency of our refinement procedure is shown on 2D MHD shock capturing simulations. Numerical results are compared in case of the interaction of a supersonic plasma with a cylinder on the adapted grid and several non-refined grids. We also present a mass loading simulation which corresponds to a 2D version of the interaction between the solar wind and a comet.     

Laser-polarized Xe NMR and MRI at Ultralow Magnetic Fields

Laser-polarized ¹²⁹Xe and a high-T_csuperconducting quantum interference device (SQUID) are used to obtain magnetic resonance images in porous materials at a magnetic field of 2.3 mT, corresponding to a Larmor frequency of 27 kHz. Image resolution of 1 mm is obtained with gradients of only 1 mT/m. The resolution of xenon chemical shifts in different physicochemical environments at ultralow fields is also demonstrated. Details of the circulating flow optical pumping apparatus and the SQUID spectrometer are presented.     

Accurate ω-ψ Spectral Solution of the Singular Driven Cavity Problem

This article provides accurate spectral solutions of the driven cavity problem, calculated in the vorticity–stream function representation without smoothing the corner singularities—a prima facie impossible task. As in a recent benchmark spectral calculation by primitive variables of Botella and Peyret, closed-form contributions of the singular solution for both zero and finite Reynolds numbers are subtracted from the unknown of the problem tackled here numerically in biharmonic form. The method employed is based on a split approach to the vorticity and stream function equations, a Galerkin–Legendre approximation of the problem for the perturbation, and an evaluation of the nonlinear terms by Gauss–Legendre numerical integration. Results computed for Re=0, 100, and 1000 compare well with the benchmark steady solutions provided by the aforementioned collocation–Chebyshev projection method. The validity of the proposed singularity subtraction scheme for computing time-dependent solutions is also established.     

High-Order Scheme for a Nonlinear Maxwell System Modelling Kerr Effect

This paper is devoted to the derivation of an efficient numerical scheme for the Kerr–Maxwell system. We begin by studying the 1-D Riemann problem. We obtain a result of existence and uniqueness for large data. Then we develop a high-order Roe solver and exhibit solutions in 1-D and 2-D simulations.     

Measurement of Small One-Bond Proton–Carbon Residual Dipolar Coupling Constants in Partially Oriented C Natural Abundance Oligosaccharide Samples: Analysis of Heteronuclear JCH-Modulated Spectra with the BIRD Inversion Pulse

Two 2D J-modulated HSQC-based experiments were designed for precise determination of small residual dipolar one-bond carbon–proton coupling constants in ¹³C natural abundance carbohydrates. Crucial to the precision of a few hundredths of Hz achieved by these methods was the use of long modulation intervals and BIRD pulses, which acted as semiselective inversion pulses. The BIRD pulses eliminated effective evolution of all but ¹J_CH couplings, resulting in signal modulation that can be described by simple modulation functions. A thorough analysis of such modulation functions for a typical four-spin carbohydrate spin system was performed for both experiments. The results showed that the evolution of the ¹H–¹H and long-range ¹H–¹³C couplings during the BIRD pulses did not necessitate the introduction of more complicated modulation functions. The effects of pulse imperfections were also inspected. While weakly coupled spin systems can be analyzed by simple fitting of cross peak intensities, in strongly coupled spin systems the evolution of the density matrix needs to be considered in order to analyse data accurately. However, if strong coupling effects are modest the errors in coupling constants determined by the “weak coupling” analysis are of similar magnitudes in oriented and isotropic samples and are partially cancelled during dipolar coupling calculation. Simple criteria have been established as to when the strong coupling treatment needs to be invoked.     

Moving Mesh Methods in Multiple Dimensions Based on Harmonic Maps

In practice, there are three types of adaptive methods using the finite element approach, namely the h-method, p-method, and r-method. In the h-method, the overall method contains two parts, a solution algorithm and a mesh selection algorithm. These two parts are independent of each other in the sense that the change of the PDEs will affect the first part only. However, in some of the existing versions of the r-method (also known as the moving mesh method), these two parts are strongly associated with each other and as a result any change of the PDEs will result in the rewriting of the whole code. In this work, we will propose a moving mesh method which also contains two parts, a solution algorithm and a mesh-redistribution algorithm. Our efforts are to keep the advantages of the r-method (e.g., keep the number of nodes unchanged) and of the h-method (e.g., the two parts in the code are independent). A framework for adaptive meshes based on the Hamilton–Schoen–Yau theory was proposed by Dvinsky. In this work, we will extend Dvinsky's method to provide an efficient solver for the mesh-redistribution algorithm. The key idea is to construct the harmonic map between the physical space and a parameter space by an iteration procedure. Each iteration step is to move the mesh closer to the harmonic map. This procedure is simple and easy to program and also enables us to keep the map harmonic even after long times of numerical integration. The numerical schemes are applied to a number of test problems in two dimensions. It is observed that the mesh-redistribution strategy based on the harmonic maps adapts the mesh extremely well to the solution without producing skew elements for multi-dimensional computations.     

A High-Order Discontinuous Galerkin Method for 2D Incompressible Flows

In this paper we introduce a high-order discontinuous Galerkin method for two-dimensional incompressible flow in the vorticity stream-function formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method. The stream function is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability. The method is efficient for inviscid or high Reynolds number flows. Optimal error estimates are proved and verified by numerical experiments.     

Two-Level Hierarchical FEM Method for Modeling Passive Microwave Devices

In recent years multigrid methods have been proven to be very efficient for solving large systems of linear equations resulting from the discretization of positive definite differential equations by either the finite difference method or theh-version of the finite element method. In this paper an iterative method of the multiple level type is proposed for solving systems of algebraic equations which arise from thep-version of the finite element analysis applied to indefinite problems. A two-levelV-cycle algorithm has been implemented and studied with a Gauss–Seidel iterative scheme used as a smoother. The convergence of the method has been investigated, and numerical results for a number of numerical examples are presented.     

Solution of Time-Dependent Diffusion Equations with Variable Coefficients Using Multiwavelets

A new numerical algorithm is developed for the solution of time-dependent differential equations of diffusion type. It allows for an accurate and efficient treatment of multidimensional problems with variable coefficients, nonlinearities, and general boundary conditions. For space discretization we use the multiwavelet bases introduced by Alpert (1993,SIAM J. Math. Anal.24, 246–262), and then applied to the representation of differential operators and functions of operators presented by Alpert, Beylkin, and Vozovoi (Representation of operators in the multiwavelet basis, in preparation). An important advantage of multiwavelet basis functions is the fact that they are supported only on non-overlapping subdomains. Thus multiwavelet bases are attractive for solving problems in finite (non periodic) domains. Boundary conditions are imposed with a penalty technique of Hesthaven and Gottlieb (1996,SIAM J. Sci. Comput., 579–612) which can be used to impose rather general boundary conditions. The penalty approach was extended to a procedure for ensuring the continuity of the solution and its first derivative across interior boundaries between neighboring subdomains while time stepping the solution of a time dependent problem. This penalty procedure on the interfaces allows for a simplification and sparsification of the representation of differential operators by discarding the elements responsible for interactions between neighboring subdomains. Consequently the matrices representing the differential operators (on the finest scale) have block-diagonal structure. For a fixed order of multiwavelets (i.e., a fixed number of vanishing moments) the computational complexity of the present algorithm is proportional to the number of subdomains. The time discretization method of Beylkin, Keiser, and Vozovoi (1998, PAM Report 347) is used in view of its favorable stability properties. Numerical results are presented for evolution equations with variable coefficients in one and two dimensions.     

MUSIC and Aromatic Residues: Amino Acid Type-Selective H–N Correlations, III

Amino acid type-selective experiments can help to remove ambiguities in automated assignment procedures for ¹⁵N/¹³C-labeled proteins. Here we present five triple-resonance experiments that yield amino acid type-selective ¹H–¹⁵N correlations for aromatic amino acids. Four of the novel experiments are based on the MUSIC coherence transfer scheme that replaces the initial INEPT transfer and is selective for CH₂. The MUSIC sequence is combined with selective excitation pulses to create experiments for Trp (W-HSQC) as well as Phe, Tyr, and His (FYH-HSQC). In addition, an experiment selective for Trp H¹–N¹ is presented. The new experiments are recorded as two-dimensional experiments and their performance is demonstrated with the application to a protein domain of 115 amino acids.     

A Triangulated Vortex Method for the Axisymmetric Euler Equations

A new method is presented for the prediction of unsteady axisymmetric inviscid flows. By combining a triangulated vortex approach with a novel evaluation technique for the Biot–Savart integrals, a Lagrangian vortex method is developed which eliminates the singularities usually present in axisymmetric methods, without recourse to normalizations or other approximations. Furthermore, the computational effort scales as the number of control points N and, in the large N limit, depends only on the order of quadrature employed. The accuracy and computational effort are assessed by comparison with the velocity field of a Gaussian core vortex ring and the use of the technique is illustrated by computation of the motion of Norbury rings and of vortex ring pairing.     

Unitary representations of superconformal algebra

superconformal algebra is algebra with two Virasoro operators. The Kac determinant is calculated and the complete list of unitary representations is determined. Two types of extensions of algebra are discussed. A new approach to construction of algebras from rational conformal field theories is proposed.     

Unconditionally Stable Methods for Hamilton–Jacobi Equations

We present new numerical methods for constructing approximate solutions to the Cauchy problem for Hamilton–Jacobi equations of the form u_t+H(D_xu)=0. The methods are based on dimensional splitting and front tracking for solving the associated (non-strictly hyperbolic) system of conservation laws p_t+D_xH(p)=0, where p=D_xu. In particular, our methods depend heavily on a front tracking method for one-dimensional scalar conservation laws with discontinuous coefficients. The proposed methods are unconditionally stable in the sense that the time step is not limited by the space discretization and they can be viewed as “large-time-step” Godunov-type (or front tracking) methods. We present several numerical examples illustrating the main features of the proposed methods. We also compare our methods with several methods from the literature.     

Implementation Aspects of 3D Lattice-BGK: Boundaries, Accuracy, and a New Fast Relaxation Method

In many realistic fluid-dynamical simulations the specification of the boundary conditions, the error sources, and the number of time steps to reach a steady state are important practical considerations. In this paper we study these issues in the case of the lattice-BGK model. The objective is to present a comprehensive overview of some pitfalls and shortcomings of the lattice-BGK method and to introduce some new ideas useful in practical simulations. We begin with an evaluation of the widely used bounce-back boundary condition in staircase geometries by simulating flow in an inclined tube. It is shown that the bounce-back scheme is first-order accurate in space when the location of the non-slip wall is assumed to be at the boundary nodes. Moreover, for a specific inclination angle of 45 degrees, the scheme is found to be second-order accurate when the location of the non-slip velocity is fitted halfway between the last fluid nodes and the first solid nodes. The error as a function of the relaxation parameter is in that case qualitatively similar to that of flat walls. Next, a comparison of simulations of fluid flow by means of pressure boundaries and by means of body force is presented. A good agreement between these two boundary conditions has been found in the creeping-flow regime. For higher Reynolds numbers differences have been found that are probably caused by problems associated with the pressure boundaries. Furthermore, two widely used 3D models, namelyD₃Q₁₅andD₃Q₁₉, are analysed. It is shown that theD₃Q₁₅model may induce artificial checkerboard invariants due to the connectivity of the lattice. Finally, a new iterative method, which significantly reduces the saturation time, is presented and validated on different benchmark problems.     

Minimization of the Truncation Error by Grid Adaptation

A new grid adaptation strategy, which minimizes the truncation error of a pth-order finite difference approximation, is proposed. The main idea of the method is based on the observation that the global truncation error associated with discretization on nonuniform meshes can be minimized if the interior grid points are redistributed in an optimal sequence. The method does not explicitly require the truncation error estimate, and at the same time, it allows one to increase the design order of approximation globally by one, so that the same finite difference operator reveals superconvergence properties on the optimal grid. Another very important characteristic of the method is that if the differential operator and the metric coefficients are evaluated identically by some hybrid approximation, then the single optimal grid generator can be employed in the entire computational domain independently of points where the hybrid discretization switches from one approximation to another. Generalization of the present method to multiple dimensions is presented. Numerical calculations of several one-dimensional and one two-dimensional test examples demonstrate the performance of the method and corroborate the theoretical results.     

A Finite Volume Method for the Approximation of Diffusion Operators on Distorted Meshes

A new finite volume method is presented for discretizing general linear or nonlinear elliptic second-order partial-differential equations with mixed boundary conditions. The advantage of this method is that arbitrary distorted meshes can be used without the numerical results being altered. The resulting algorithm has more unknowns than standard methods like finite difference or finite element methods. However, the matrices that need to be inverted are positive definite, so the most powerful linear solvers can be applied. The method has been tested on a few elliptic and parabolic equations, either linear, as in the case of the standard heat diffusion equation, or nonlinear, as in the case of the radiation diffusion equation and the resistive diffusion equation with Hall term.