An Error Indicator Monitor Function for an r-Adaptive Finite-Element Method

An r-adaptive finite-element method based on moving-mesh partial differential equations (PDEs) and an error indicator is presented. The error indicator is obtained by applying a technique developed by Bank and Weiser to elliptic equations which result in this case from temporal discretization of the underlying physical PDEs on moving meshes. The construction of the monitor function based on the error indicator is discussed. Numerical results obtained with the current method and the commonly used method based on solution gradients are presented and analyzed for several examples.     

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.     

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.     

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.     

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.     

Coherent Raman and Infrared Studies of Sulfur Trioxide

High-resolution (0.001 cm⁻¹) coherent anti-Stokes Raman scattering (CARS) was used to observe the Q-branch structure of the IR-inactive ν₁ symmetric stretching mode of ³²S¹⁶O₃ and its various ¹⁸O isotopomers. The ν₁ spectrum of ³²S¹⁶O₃ reveals two intense Q-branches in the region 1065–1067 cm⁻¹, with surprisingly complex vibrational–rotational structure not resolved in earlier studies. Efforts to simulate this with a simple Fermi-resonance model involving ν₁ and 2ν₄ states do not reproduce the spectral detail, nor do they yield reasonable spectroscopic parameters. A more subtle combination of Fermi resonance and indirect Coriolis interactions with nearby states, 2ν₄(1=0, ±2), ν₂+ν₄(1=±1), 2ν₂(1=0), is suspected and a determination of the location of these coupled states by high-resolution infrared measurements is under way. At medium resolution (0.125 cm⁻¹), the infrared spectra reveal Q-branch features from which approximate band origins are estimated for the ν₂, ν₃, and ν₄ fundamental modes of ³²S¹⁸O₃, ³²S¹⁸O₂¹⁶O, and ³²S¹⁸O¹⁶O₂. These and literature data for ³²S¹⁶O₃ are used to calculate force constants for SO₃ and a comparison is made with similar values for SO₂ and SO. The frequencies and force constants are in excellent agreement with those obtained by Martin in a recent ab initio calculation.     

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.     

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.     

Separating Structure and Dynamics in CSA/DD Cross-Correlated Relaxation: A Case Study on Trehalose and Ubiquitin

We present two new sensitivity enhanced gradient NMR experiments for measuring interference effects between chemical shift anisotropy (CSA) and dipolar coupling interactions in a scalar coupled two-spin system in both the laboratory and rotating frames. We apply these methods for quantitative measurement of longitudinal and transverse cross-correlation rates involving interference of ¹³C CSA and ¹³C–¹H dipolar coupling in a disaccharide, α,α- -trehalose, at natural abundance of ¹³C as well as interference of amide ¹⁵N CSA and ¹⁵N–¹H dipolar coupling in uniformly ¹⁵N-labeled ubiquitin. We demonstrate that the standard heteronuclear T₁, T₂, and steady-state NOE autocorrelation experiments augmented by cross-correlation measurements provide sufficient experimental data to quantitatively separate the structural and dynamic contributions to these relaxation rates when the simplifying assumptions of isotropic overall tumbling and an axially symmetric chemical shift tensor are valid.     

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.     

Tagging of the Magnetization with the Transition Zones of 360° Rotations Generated by a Tandem of Two Adiabatic DANTE Inversion Sequences

A new technique is presented for generating myocardial tagging using the signal intensity minima of the transition zones between the bands of 0° and 360° rotations, induced by a tandem of two adiabatic delays alternating with nutations for tailored excitation (DANTE) inversion sequences. With this approach, the underlying matrix corresponds to magnetization that has experienced 0° or 360° rotations. The DANTE sequences were implemented from adiabatic parent pulses for insensitivity of the underlying matrix to B₁ inhomogeneity. The performance of the proposed tagging technique is demonstrated theoretically with computer simulations and experimentally on phantom and on the canine heart, using a surface coil for both RF transmission and signal reception. The simulations and the experimental data demonstrated uniform grid contrast and sharp tagging profiles over a twofold variation of the B₁ field magnitude.     

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.     

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.     

Millimeter-Wave Spectroscopy, High-Resolution Infrared Spectrum, ab Initio Calculations, and Molecular Geometry of FPS

The transient thiophosphenous fluoride FPS was produced by pyrolysis of 2.5% F₂PSPF₂ in Ar at 1300–1800°C. High-resolution (≥0.004 cm⁻¹) Fourier transform infrared spectra of the a-type ν₁ and b-type ν₂ bands, centered respectively at 803.249 and 726.268 cm⁻¹, were measured and fitted to rotational and quartic centrifugal distortion parameters. The millimeter-wave spectrum, essentially b-type, was measured between 300 and 370 GHz in the ground state and in the ν₃ excited state for FP³²S and in the ground state for FP³⁴S. The frequencies were fitted to a Watson-type A-reduced Hamiltonian up to sextic distortion terms. High level ab initio calculations with large basis sets were performed on FPS and supported the first identification of its infrared and millimeter wave spectra. The calculated anharmonic force field provided precise ab initio rovibrational α constants which were combined with the experimental molecular parameters to determine an accurate equilibrium structure of the molecule: r_e(PS)=188.86 pm, r_e(PF)=158.70 pm, θ(FPS)=109.28°. The collision-controlled 1/e lifetime measured in a 10-Pa (1 : 20) F₂PSPF₂/Ar mixture was 2 s, more than two orders of magnitude larger than that of FPO under the same experimental conditions.     

Probing Proteins in Solution by Xe NMR Spectroscopy

The interaction of xenon with different proteins in aqueous solution is investigated by ¹²⁹Xe NMR spectroscopy. Chemical shifts are measured in horse metmyoglobin, hen egg white lysozyme, and horse cytochrome c solutions as a function of xenon concentration. In these systems, xenon is in fast exchange between all possible environments. The results suggest that nonspecific interactions exist between xenon and the protein exteriors and the data are analyzed in term of parameters which characterize the protein surfaces. The experimental data for horse metmyoglobin are interpreted using a model in which xenon forms a 1:1 complex with the protein and the chemical shift of the complexed xenon is reported (Locci et al., Keystone Symposia “Frontiers of NMR in Molecular Biology VI”, Jan. 9–15, 1999, Breckenridge, CO, Abstract E216, p. 53; Locci et al., XeMAT 2000 “Optical Polarization and Xenon NMR of Materials”, June 28–30, 2000, Sestri Levante, Italy, p. 46).     

The Adjoint Method for an Inverse Design Problem in the Directional Solidification of Binary Alloys

This paper presents a systematic procedure based on the adjoint method for solving a class of inverse directional alloy solidification design problems in which a desired growth velocityv_fis achieved under stable growth conditions. To the best of our knowledge, this is the first time that a continuum adjoint formulation is proposed for the solution of an inverse problem with simultaneous heat and mass transfer, thermo-solutal convection, and phase change. In this paper, the interfacial stability is considered to imply a sharp solid–liquid freezing interface. This condition is enforced using the constitutional undercooling criterion in the form of an inequality constraint between the thermal and solute concentration gradients,GandG_c, respectively, at the freezing front. The main unknowns of the design problem are the heating and/or cooling boundary conditions on the mold walls. The inverse design problem is formulated as a functional optimization problem. The cost functional is defined by the square of theL₂norm of the deviation of the freezing interface temperature from the temperature corresponding to thermodynamic equilibrium. A continuum adjoint system is derived to calculate the adjoint temperature, concentration, and velocity fields such that the gradient of the cost functional can be expressed analytically. The cost functional minimization process is realized by the conjugate gradient method via the finite element method solutions of the continuum direct, sensitivity, and adjoint problems. The developed formulation is demonstrated with an example of designing the directional solidification of a binary aqueous solution in a rectangular mold such that a stable vertical interface advances from left to right with a desired growth velocity.     

H NMR Characterization of the Product from Single Solid-Phase Resin Beads Using Capillary NMR Flow Probes

A capillary NMR flow probe was designed to generate high-resolution ¹H NMR spectra at 600 MHz from the cleaved product of individual 160-μm Tentagel combinatorial chemistry beads. By injecting a dissolved sample sandwiched between an immiscible, perfluorinated organic liquid directly into the probe, NMR spectra of the product cleaved from single beads were acquired in just 1 h of spectrometer time without diffusional dilution. Sample handling efficiency on the single bead scale was comparable to that obtained with a bulk sample. Using the relative intensity of the DMSO-d₅H versus the analyte signals in a fully relaxed CPMG spectrum, the amount of product cleaved from a single bead was determined to be 540±170 pmol in one of the samples. Following the NMR data collection, the samples were examined with electrospray ionization mass spectrometry to provide additional structural information. By coupling with microliter-volume fluidic capabilities, the capillary flow probe described here will enable multidimensional characterization of single solid-phase resin products in an online manner.     

H Spectroscopy without Solvent Suppression: Characterization of Signal Modulations at Short Echo Times

While most proton (¹H) spectra acquired in vivo utilize selective suppression of the solvent signal for more sensitive detection of signals from the dilute metabolites, recent reports have demonstrated the feasibility and advantages of collecting in vivo data without solvent attenuation. When these acquisitions are performed at short echo times, the presence of frequency modulations of the water resonance may become an obstacle to the identification and quantitation of metabolite resonances. The present report addresses the characteristics, origin, and elimination of these sidebands. Sideband amplitudes were measured as a function of delay time between gradient pulse and data collection, as a function of gradient pulse amplitude, and as a function of spatial location of the sample for each of the three orthogonal gradient sets. Acoustic acquisitions were performed to demonstrate the correlation between mechanical vibration resonances and the frequencies of MR sidebands. A mathematical framework is developed and compared with the experimental results. This derivation is based on the theory that these frequency modulations are induced by magnetic field fluctuations generated by the transient oscillations of gradient coils.     

Linewidth-Resolved N HSQC, a Simple 3D Method to Measure N Relaxation Times from T1 and T2 Linewidths

A three-dimensional approach for measuring ¹⁵N relaxation times is described. Instead of selecting particular values for the relaxation period, in the proposed method the relaxation period is incremented periodically in order to create a 3D spectrum. This additional frequency domain of the transformed spectrum contains the relaxation time information in the T₁ and T₂ linewidths, and thus the longitudinal and transverse ¹⁵N relaxation times can be measured without determination of 2D cross peak volumes/intensities and subsequent curve fitting procedures.     

Angular Flow in Toroid Cavity Probes

NMR signals from samples that rotate uniformly about the central conductor of a TCD (toroid cavity detector) exhibit frequency shifts that are directly proportional to the sample's angular velocity. This newly observed effect is based on the unique radiofrequency field inside TCDs, which is variable in direction. If a liquid sample is pumped through a capillary tube wound about the central conductor, the frequency shift is proportional to the flow rate. A mathematical relationship between a volumetric flow rate and the frequency shift is established and experimentally verified to high precision. Additionally, two-dimensional flow-resolved NMR spectroscopy for discrimination between components with different flow velocities yet retaining chemical shift information for structural analysis is presented. The application of the two-dimensional method in chromatographic NMR is suggested. Furthermore, utilization of the frequency-shift effect for rheologic studies if combined with toroid-cavity rotating-frame imaging is proposed.