A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics

We present a high-order accurate weighted essentially non-oscillatory (WENO) finite difference scheme for solving the equations of ideal magnetohydrodynamics (MHD). This scheme is a direct extension of a WENO scheme, which has been successfully applied to hydrodynamic problems. The WENO scheme follows the same idea of an essentially non-oscillatory (ENO) scheme with an advantage of achieving higher-order accuracy with fewer computations. Both ENO and WENO can be easily applied to two and three spatial dimensions by evaluating the fluxes dimension-by-dimension. Details of the WENO scheme as well as the construction of a suitable eigen-system, which can properly decompose various families of MHD waves and handle the degenerate situations, are presented. Numerical results are shown to perform well for the one-dimensional Brio–Wu Riemann problems, the two-dimensional Kelvin–Helmholtz instability problems, and the two-dimensional Orszag–Tang MHD vortex system. They also demonstrate the importance of maintaining the divergence free condition for the magnetic field in achieving numerical stability. The tests also show the advantages of using the higher-order scheme. The new 5th-order WENO MHD code can attain an accuracy comparable with that of the second-order schemes with many fewer grid points.     

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.     

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.     

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.     

Convection-Diffusion Lattice Boltzmann Scheme for Irregular Lattices

In this paper, a lattice Boltzmann (LB) scheme for convection diffusion on irregular lattices is presented, which is free of any interpolation or coarse graining step. The scheme is derived using the axioma that the velocity moments of the equilibrium distribution equal those of the Maxwell–Boltzmann distribution. The axioma holds for both Bravais and irregular lattices, implying a single framework for LB schemes for all lattice types. By solving benchmark problems we have shown that the scheme is indeed consistent with convection diffusion. Furthermore, we have compared the performance of the LB schemes with that of finite difference and finite element schemes. The comparison shows that the LB scheme has a similar performance as the one-step second-order Lax–Wendroff scheme: it has little numerical diffusion, but has a slight dispersion error. By changing the relaxation parameter ω the dispersion error can be balanced by a small increase of the numerical diffusion.     

A General-Purpose Finite-Volume Advection Scheme for Continuous and Discontinuous Fields on Unstructured Grids

We present a new general-purpose advection scheme for unstructured meshes based on the use of a variation of the interface-tracking flux formulation recently put forward by O. Ubbink and R. I. Issa (J. Comput. Phys.153, 26 (1999)), in combination with an extended version of the flux-limited advection scheme of J. Thuburn (J. Comput. Phys.123, 74 (1996)), for continuous fields. Thus, along with a high-order mode for continuous fields, the new scheme presented here includes optional integrated interface-tracking modes for discontinuous fields. In all modes, the method is conservative, monotonic, and compatible. It is also highly shape preserving. The scheme works on unstructured meshes composed of any kind of connectivity element, including triangular and quadrilateral elements in two dimensions and tetrahedral and hexahedral elements in three dimensions. The scheme is finite-volume based and is applicable to control-volume finite-element and edge-based node-centered computations. An explicit–implicit extension to the continuous-field scheme is provided only to allow for computations in which the local Courant number exceeds unity. The transition from the explicit mode to the implicit mode is performed locally and in a continuous fashion, providing a smooth hybrid explicit–implicit calculation. Results for a variety of test problems utilizing the continuous and discontinuous advection schemes are presented.     

A Residual-Based Compact Scheme for the Compressible Navier–Stokes Equations

A simple and efficient time-dependent method is presented for solving the steady compressible Euler and Navier–Stokes equations with third-order accuracy. Owing to its residual-based structure, the numerical scheme is compact without requiring any linear algebra, and it uses a simple numerical dissipation built on the residual. The method contains no tuning parameter. Accuracy and efficiency are demonstrated for 2-D inviscid and viscous model problems. Navier–Stokes calculations are presented for a shock/boundary layer interaction, a separated laminar flow, and a transonic turbulent flow over an airfoil.     

A Composite Scheme for Gas Dynamics in Lagrangian Coordinates

One cycle of a composite finite difference scheme is defined as several time steps of an oscillatory scheme such as Lax–Wendroff followed by one step of a diffusive scheme such as Lax–Friedrichs. We apply this idea to gas dynamics in Lagrangian coordinates. We show numerical results in two dimensions for Noh's infinite strength shock problem and the Sedov blast wave problem, and for several one-dimensional problems including a Riemann problem with a contact discontinuity. For Noh's problem the composite scheme produces a better result than that obtained with a more conventional Lagrangian code.     

A Direct Spectral Collocation Poisson Solver in Polar and Cylindrical Coordinates

In this paper, we present a direct spectral collocation method for the solution of the Poisson equation in polar and cylindrical coordinates. The solver is applied to the Poisson equations for several different domains including a part of a disk, an annulus, a unit disk, and a cylinder. Unlike other Poisson solvers for geometries such as unit disks and cylinders, no pole condition is involved for the present solver. The method is easy to implement, fast, and gives spectral accuracy. We also use the weighted interpolation technique and nonclassical collocation points to improve the convergence.     

Improved Broadband Inversion Performance for NMR in Liquids

New NMR broadband inversion pulses that compensate both for resonance offset and radiofrequency (RF) inhomogeneity are described. The approach described is a straightforward computer optimization of an initial digitized waveform generated from either a constant-amplitude frequency sweep or from an existing composite inversion pulse. Problems with convergence to local minima are alleviated by the way the optimization is carried out. For a given duration and maximum allowable RF field strength B₁ (but not necessarily given RMS power deposition), the resultant broadband inversion pulse (BIP) shows superior inversion compared to inversion pulses obtained from previous methods, including adiabatic inversion pulses. Any existing BIP can be systematically elaborated to build up longer inversion pulses that perform over larger and larger bandwidths. The resulting pulse need not be adiabatic throughout its duration or across the entire operational bandwidth.     

Nonlinear Landau Damping in Spherically Symmetric Vlasov Poisson Systems

The Vlasov Poisson system is a partial differential equation widely used to describe collisionless plasma. It is formulated in a six-dimensional phase space, this prohibits a numerical solution on a complete phase space grid. In some applications, however, spherical symmetry is given, which introduces singularities into the Vlasov Poisson equation. We focus on such problems and propose a stable algorithm using accommodating boundaries. At first, the method is tested in the linear regime, where analytical solutions are available. Thereafter it is applied to large disturbances from equilibrium.     

Half-Moment Closure for Radiative Transfer Equations

In this paper a moment method for radiative transfer equations is considered which has been developed and investigated using different approaches. Problems appearing for this moment system for boundary value problems using Maxwell-type boundary conditions are described. A new method based on the consideration of positive and negative half fluxes is developed and shown to overcome the above problems. Moreover, a numerical scheme and numerical results for the new moment system are presented.     

A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows

We propose a new model and a solution method for two-phase compressible flows. The model involves six equations obtained from conservation principles applied to each phase, completed by a seventh equation for the evolution of the volume fraction. This equation is necessary to close the overall system. The model is valid for fluid mixtures, as well as for pure fluids. The system of partial differential equations is hyperbolic. Hyperbolicity is obtained because each phase is considered to be compressible. Two difficulties arise for the solution: one of the equations is written in non-conservative form; non-conservative terms exist in the momentum and energy equations. We propose robust and accurate discretisation of these terms. The method solves the same system at each mesh point with the same algorithm. It allows the simulation of interface problems between pure fluids as well as multiphase mixtures. Several test cases where fluids have compressible behavior are shown as well as some other test problems where one of the phases is incompressible. The method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state.     

Design and implementation of a protection system for NMR spectrometers

We have implemented a scheme, SPECMON, for monitoring various parameters of a spectrometer, such as nitrogen pressure and sample temperature, and taking corrective action. The scheme is based on considerations of protection management which are of general application. Evaluation of the spectrometer state is incorporated in macros of the application software (VNMR) and is therefore very flexible. In contrast, corrective action is limited to the single one which is deemed fully safe: complete shutdown of the spectrometer and logging. Shutdown is implemented by a minor hardware modification of the spectrometer: the introduction of a second input to a relay already present for protection of the spectrometer power supply. Monitoring is handled by the host computer, and the shutdown command is transmitted via control lines of its series port, independent of the standard connection between the host computer and the NMR system console. The monitoring system (software and hardware) is unobtrusive in normal conditions, and it can be tested without affecting the operation of the spectrometer.     

Modeling Three-Dimensional Multiphase Flow Using a Level Contour Reconstruction Method for Front Tracking without Connectivity

Three-dimensional multiphase flow and flow with phase change are simulated using a simplified method of tracking and reconstructing the phase interface. The new level contour reconstruction technique presented here enables front tracking methods to naturally, automatically, and robustly model the merging and breakup of interfaces in three-dimensional flows. The method is designed so that the phase surface is treated as a collection of physically linked but not logically connected surface elements. Eliminating the need to bookkeep logical connections between neighboring surface elements greatly simplifies the Lagrangian tracking of interfaces, particularly for 3D flows exhibiting topology change. The motivation for this new method is the modeling of complex three-dimensional boiling flows where repeated merging and breakup are inherent features of the interface dynamics. Results of 3D film boiling simulations with multiple interacting bubbles are presented. The capabilities of the new interface reconstruction method are also tested in a variety of two-phase flows without phase change. Three-dimensional simulations of bubble merging and droplet collision, coalescence, and breakup demonstrate the new method's ability to easily handle topology change by film rupture or filamentary breakup. Validation tests are conducted for drop oscillation and bubble rise. The susceptibility of the numerical method to parasitic currents is also thoroughly assessed.     

On Backtracking Failure in Newton–GMRES Methods with a Demonstration for the Navier–Stokes Equations

In an earlier study of inexact Newton methods, we pointed out that certain counterintuitive behavior may occur when applying residual backtracking to the Navier–Stokes equations with heat and mass transport. Specifically, it was observed that a Newton–GMRES method globalized by backtracking (linesearch, damping) may be less robust when high accuracy is required of each linear solve in the Newton sequence than when less accuracy is required. In this brief discussion, we offer a possible explanation for this phenomenon, together with an illustrative numerical experiment involving the Navier–Stokes equations.     

A subspace time-domain algorithm for automated NMR spectral normalization

Recently, two methods have been proposed for quantitatively comparing NMR spectra of control and treated samples, in order to examine the possible occurring variations in cell metabolism and/or structure in response to numerous physical, chemical, and biological agents. These methods are the maximum superposition normalization algorithm (MaSNAl) and the minimum rank normalization algorithm (MiRaNAl). In this paper a new subspace-based time-domain normalization algorithm, denoted by SuTdNAl (subspace time-domain normalization algorithm), is presented. By the determination of the intersection of the column spaces of two Hankel matrices, the common signal poles and further on the components having proportionally varying amplitudes are detected. The method has the advantage that it is computationally less intensive than the MaSNAl and the MiRaNAl. Furthermore, no approximate estimate of the normalization factor is required. The algorithm was tested by Monte Carlo simulations on a set of simulation signals. It was shown that the SuTdNAl has a statistical performance similar to that of the MiRaNAl, which itself is an improvement over the MaSNAl. Furthermore, two samples of known contents are compared with the MiRaNAl, the SuTdNAl, and an older method using a standard. Finally, the SuTdNAl is tested on a realistic simulation example derived from an in vitro measurement on cells.     

Tailored HCCH-TOCSY experiment for resonance assignment in the proximity of a paramagnetic center

The presence of a paramagnetic center may disturb both coherent and incoherent communication between nuclear spins that are affected, to some extent, by the hyperfine interaction. This is a limiting factor to an extensive use of paramagnetic probes in NMR spectroscopy to enhance partial alignment and to exploit cross correlation effects and pseudocontact shifts. We propose here an HCCH-TOCSY experiment tailored to identify spin systems involving resonances that are partly or completely affected by hyperfine interaction. The efficiency of polarization transfer steps when fast relaxing nuclei are involved is discussed. The sequence is tested for the protein Calbindin D(9k), in which one of the two native Ca2+ ions is replaced by the paramagnetic Ce3+ ion as well as for the oxidized form of cytochrome b(562).     

Multifrequency Probe for Pulsed EPR and ENDOR Spectroscopy

The design, construction, and performance of a multifrequency pulsed EPR and ENDOR probe for use at cryogenic temperatures are described. Interchangeable resonators based on a folded strip line design allow variation of the resonance frequency over a range of 5-11 GHz. Variable coupling to the resonator is achieved capacitively via a simple mechanical adjustment which is thermally and mechanically stable. The entire assembly is robust and easily fabricated. Common methods of analyzing the resonator parameters such as the Q-factor and coupling coefficient are discussed quantitatively. Probe performance data and multifrequency pulsed ENDOR spectra are presented.     

Improved narrowband dipolar recoupling for homonuclear distance measurements in rotating solids

Recovery of the magnetic dipolar interaction between nuclei bearing the same gyromagnetic ratio in rotating solids can be promoted by synchronous rf irradiation. Determination of the dipolar interaction strength can serve as a tool for structural elucidation in polycrystalline powders. Spinning frequency dependent narrow-band (nb) RFDR and SEDRA experiments are utilized as simple techniques for the determination of dipolar interactions between the nuclei in coupled homonuclear spin pairs. The magnetization exchange and coherence dephasing due to a fixed number of rotor-synchronously applied pi-pulses is monitored at spinning frequencies in the vicinity of the rotational resonance (R(2)) conditions. The powder nbRFDR and nbSEDRA decay curves of spin magnetizations and coherences, respectively, as a function of the spinning frequency can be measured and analyzed using simple rate equations providing a quantitative measure of the dipolar coupling. The effects of the phenomenological relaxation parameters in these rate equations are discussed and an improved methodology is suggested for analyzing nbRFDR data for small dipolar couplings. The distance between the labeled nuclei in the 1,3-(13)C(2)-hydroxybutyric acid molecule is rederived using existing nbRFDR results and the new simulation procedure. A nbSEDRA experiment has been performed successfully on a powder sample of singly labeled 1-(13)C-L-leucine measuring the dipolar interaction between the labeled carboxyl carbon and the natural abundant beta-carbon. Both narrowband techniques are employed for the determination of the nuclear distances between the side-chain carbons of leucine and its carbonyl carbon in a tripeptide Leu-Gly-Phe that is singly (13)C-labeled at the leucine carbonyl carbon position.