A General Deterministic Treatment of Derivatives in Particle Methods

A unified approach to approximating spatial derivatives in particle methods using integral operators is presented. The approach is an extension of particle strength exchange, originally developed for treating the Laplacian in advection–diffusion problems. Kernels of high order of accuracy are constructed that can be used to approximate derivatives of any degree. A new treatment for computing derivatives near the edge of particle coverage is introduced, using “one-sided” integrals that only look for information where it is available. The use of these integral approximations in wave propagation applications is considered and their error is analyzed in this context using Fourier methods. Finally, simple tests are performed to demonstrate the characteristics of the treatment, including an assessment of the effects of particle dispersion, and their results are discussed.     

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.     

Imaging the long-range dipolar field in structured liquid state samples

We describe imaging experiments in which the pattern of the dipolar field generated by spatially modulated nuclear magnetization is directly visualized in simply structured phantoms. Two types of experiment have been carried out at 11.7 T using (1)H NMR signals. In the first, the field from a single spin species is imaged via its own NMR signal. In the second, the NMR signal from one spin species is used to image the field generated by a second species. The field patterns measured in these experiments correspond well with those calculated using simple theoretical expressions for the dipolar field. The results also directly demonstrate the spatial sensitivity of the signal generated using dipolar field effects, indicating that the range of the field depends upon the inverse of the spatial frequency with which the magnetization is modulated.     

Wave Simulation in Frozen Porous Media

We propose a numerical algorithm for simulation of wave propagation in frozen porous media, where the pore space is filled with ice and water. The model, based on a Biot-type three-phase theory, predicts three compressional waves and two shear waves and models the attenuation level observed in rocks. Attenuation is modeled with exponential relaxation functions which allow a differential formulation based on memory variables. The wavefield is obtained using a grid method based on the Fourier differential operator and a Runge–Kutta time-integration algorithm. Since the presence of slow quasistatic modes makes the differential equations stiff, a time-splitting integration algorithm is used to solve the stiff part analytically. The modeling is second-order accurate in the time discretization and has spectral accuracy in the calculation of the spatial derivatives.     

A New Version of the Fast Multipole Method for Screened Coulomb Interactions in Three Dimensions

We present a new version of the fast multipole method (FMM) for screened Coulomb interactions in three dimensions. Existing schemes can compute such interactions in O(N) time, where N denotes the number of particles. The constant implicit in the O(N) notation, however, is dominated by the expense of translating far-field spherical harmonic expansions to local ones. For each box in the FMM data structure, this requires 189p⁴ operations per box, where p is the order of the expansions used. The new formulation relies on an expansion in evanescent plane waves, with which the amount of work can be reduced to 40p²+6p³ operations per box.     

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.     

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.     

Relaxation Behavior of Laser-Polarized Xe Gas: Size Dependency and Wall Effect of the T1 Relaxation Time in Glass and Gelatin Bulbs

Size dependency of the relaxation time T₁ was measured for laser-polarized ¹²⁹Xe gas encapsulated in different sized cavities made by glass bulbs or gelatin capsules. The use of laser-polarized gas enhances the sensitivity a great deal, making it possible to measure the longer ¹²⁹Xe relaxation time in quite a short time. The size dependency is analyzed on the basis of the kinetic theory of gases and a relationship is derived in which the relaxation rate is connected with the square inverse of the diameter of the cavity. Such an analysis provides a novel parameter which denotes the wall effect on the relaxation rate when a gas molecule collides with the surface once in a second. The relaxation time of ¹²⁹Xe gas is also dependent on the material which forms the cavity. This dependency is large and the relaxation study using polarized ¹²⁹Xe gas is expected to offer important information about the state of the matter of the cavity wall.     

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.     

Lattice dynamics,third order elastic constants and low temperature lattice thermal expansion of zirconium

The lattice dynamics and the anharmonic properties of the hexagonal Zirconium are worked out using Keating's approach. The dispersion curves are fitted using twelve second order parameters and the six second order elastic constants are evaluated. The ten third order elastic constants are calculated using five third order parameters. The experimental measurements on the pressure derivatives of the second order elastic constants in Zirconium are in good agreement with the calculated values. The low-temperature limit of the lattice thermal expansion is calculated which agrees well with the value obtained from thermal expansion data. The variation of the generalised GPs of the elastic modes with the direction of propagation is illustrated by polar diagram.     

Generation and Application of Ultra-High-Intensity Magnetic Field Gradient Pulses for NMR Spectroscopy

Two different concepts of gradient current power supplies are introduced, which are suitable for the generation of ultra-high intensity pulsed magnetic field gradients of alternating polarity. The first system consists of a directly binary coded current source (DBCCS). It yields current pulses of up to ±120 A and a maximum voltage across the gradient coil of ±400 V. The second system consists of two TECHRON 8606 power supplies in push–pull configuration (PSPPC). It yields current pulses of up to ±100 A and a maximum voltage across the gradient coil of ±300 V. In combination with actively shielded anti-Helmholtz gradient coils, both systems are used routinely in NMR diffusion studies with unipolar pulsed field gradients of up to 35 T/m. Until now, alternating pulsed field gradient experiments were successfully performed with gradient intensities of up to ±25 T/m (DBCCS) and ±35 T/m (PSPPC), respectively. Based on the observation of the NMR spin echo in the presence of a small read gradient, procedures to test the stability and the matching of such ultra-high pulsed field gradient intensities as well as an automated routine for the compensation of possible mismatches are introduced. The results of these procedures are reported for the PSPPC system.     

The δ-function expansion of the modified two-particle Ursell function of a hard-sphere fluid

The expansion of the modified two-particle Ursell functionU(r) of a hardsphere quantal fluid is obtained in terms of a series of derivatives of δ-function. This expansion has been used to expand the second virial co-efficientB ₂ of the fluid. The expansion is correct up to the fourth power in thermal wavelength and the terms of the order of λ⁸ and λ⁴ in the first expa nsion are new.     

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.     

Expansions of the solutions of the biconfluent Heun equation in terms of incomplete Beta and Gamma functions

Considering the equations for some functions involving the first or the second derivatives of the biconfluent Heun function, we construct two expansions of the solutions of the biconfluent Heun equation in terms of incomplete Beta functions. The first series applies single Beta functions as expansion functions, while the second one involves a combination of two Beta functions. The coefficients of expansions obey four- and five-term recurrence relations, respectively. It is shown that the proposed technique is potent to produce series solutions in terms of other special functions. Two examples of such expansions in terms of the incomplete Gamma functions are presented.     

The effects of dissolved oxygen upon amide proton relaxation and chemical shift in a perdeuterated protein

The effects of dissolved molecular oxygen upon amide proton ((1)H(N)) longitudinal and transverse relaxation rates and chemical shifts were studied for a small protein domain, the second type 2 module of fibronectin ((2)F2)-isotopically enriched to 99% (2)H, 98% (15)N. Longitudinal relaxation rate enhancements, R(O(2))((1)H(N)), of individual backbone (1)H(N) nuclei varied up to 14 fold between a degassed and oxygenated (1 bar) solution, indicating that the oxygen distribution within the protein is inhomogeneous. On average, smaller relaxation rate enhancements were observed for (1)H(N) nuclei associated with the core of the protein compared to (1)H(N) nuclei closer to the surface, suggesting restricted oxygen accessibility to some regions. In agreement with an O(2)-(1)H(N) hyperfine interaction in the extreme narrowing limit, the (1)H(N) transverse relaxation rates showed no significant change, up to an oxygen pressure of 9.5 bar (the maximum pressure used in this study). For most (1)H(N) resonances, small deltadelta(O(2))((1)H(N)) hyperfine chemical shifts could be detected between oxygen pressures of 1 bar and 9.5 bar.     

Dynamical Study of the Empty Bianchi Type I Model in Generalised Scalar-tensor Theory

A dynamical study of the generalised scalar-tensor theory in the empty Bianchi type I model is made. We use a method from which we derive the sign of the first and second derivatives of the metric functions and examine three different theories that can all tend towards relativistic behaviours at late time. We determine conditions so that the dynamic is in expansion and decelerated at late time.     

Lattice dynamics,third order elastic constants and thermal expansion of rhenium

A systematic investigation of the lattice dynamics, second and third order elastic constants and thermal expansion has been carried out in rhenium using Keating's approach. The ten third-order elastic constants of rhenium are calculated using four anharmonic parameters. The present model reproduces the measured pressure derivatives of the second order elastic constants well. The low and high temperature limits γ̄L and γ̄H of the lattice thermal expansion are calculated. The agreement between the calculated γ̄H and that obtained by Gschneidner from the thermal expansion and specific heat data of rhenium is satisfactory.     

Dispersion relations and lattice thermal expansion of dysprosium

The lattice dynamics, second and third order elastic constants and the lattice thermal expansion of dysprosium have been calculated using Keating's approach. The ten third order elastic constants are calculated using four anharmonic parameters. The present model reproduces the measured pressure derivatives of the second order elastic constants of dysprosium well. The low and high temperature limits γ?_I and γ?_H of the lattice thermal expansion are evaluated and the agreement between the calculated γ?_H and that obtained from the thermal expansion and specific heat data is satisfactory.     

Approximation and Reconstruction of the Electrostatic Field in Wire–Plate Precipitators by a Low-Order Model

The numerical computation of the ionic space charge and electric field produced by corona discharge in a wire–plate electrostatic precipitator (ESP) is considered. The electrostatic problem is defined by a reduced set of the Maxwell equations. Since self-consistent conditions at the wire and at the plate cannot be specified a priori, a time-consuming iterative numerical procedure is required. The efficiency of all numerical solvers of the reduced Maxwell equations depends in particular on the accuracy of the initial guess solution. The objectives of this work are two: first, we propose a semianalytical technique based on the Karhunen–Loève (KL) decomposition of the current density field J, which can significantly improve the performance of a numerical solver; second, we devise a procedure to reconstruct the complete electric field from a given J. The approximate solution of the current density field is based on the derivation of an analytical approximation , which, added to a linear combination of few KL basis functions, constitutes an accurate approximation of J. In the first place, this result is useful for optimization procedures of the current density field, which involve the computation of many different configurations. Second, we show that from the current density field we can obtain an accurate estimate for the complete electrostatic field which can be used to speed up the convergence of the iterative procedure of standard numerical solvers.     

Exact renormalization group study of fermionic theories

The exact renormalization group approach (ERG) is developed for the case of pure fermionic theories by deriving a Grassmann version of the ERG equation and applying it to the study of fixed point solutions and critical exponents of the two-dimensional chiral Gross-Neveu model. An approximation based on the derivative expansion and a further truncation in the number of fields is used. Two solutions are obtained analytically in the limit N → ∞, with N being the number of fermionic species. For finite N some fixed point solutions, with their anomalous dimensions and critical exponents, are computed numerically. The issue of separation of physical results from the numerous spurious ones is discussed. We argue that one of the solutions we find can be identified with that of Dashen and Frishman, whereas the others seem to be new ones.