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.     

A Volume of Fluid Based Method for Fluid Flows with Phase Change

This paper presents a numerical method directed towards the simulation of flows with mass transfer due to changes of phase. We use a volume of fluid (VOF) based interface tracking method in conjunction with a mass transfer model and a model for surface tension. The bulk fluids are viscous, conducting, and incompressible. A one-dimensional test problem is developed with the feature that a thin thermal layer propagates with the moving phase interface. This test problem isolates the ability of a method to accurately calculate the thermal layers responsible for driving the mass transfer in boiling flows. The numerical method is tested on this problem and then is used in simulations of horizontal film boiling.     

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.     

An Incompressible Three-Dimensional Multiphase Particle-in-Cell Model for Dense Particle Flows

A three-dimensional, incompressible, multiphase particle-in-cell method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to an Eulerian grid and then mapping back computed stress tensors to particle positions. A subgrid particle, normal stress model for discrete particles which is robust and eliminates the need for an implicit calculation of the particle normal stress on the grid is presented. Interpolation operators and their properties are defined which provide compact support, are conservative, and provide fast solution for a large particle population. The solution scheme allows for distributions of types, sizes, and density of particles, with no numerical diffusion from the Lagrangian particle calculations. Particles are implicitly coupled to the fluid phase, and the fluid momentum and pressure equations are implicitly solved, which gives a robust solution.     

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.     

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.     

Multigrid Solution of the Incompressible Navier–Stokes Equations for Density-Stratified Flow past Three-Dimensional Obstacles

The steady incompressible Navier–Stokes equations in three dimensions are solved for neutral and stably stratified flow past three-dimensional obstacles of increasing spanwise width. The continuous equations are approximated using a finite volume discretisation on staggered grids with a flux-limited monotonic scheme for the advective terms. The discrete equations which arise are solved using a nonlinear multigrid algorithm with up to four grid levels using the SIMPLE pressure correction method as smoother. When at its most effective the multigrid algorithm is demonstrated to yield convergence rates which are independent of the grid density. However, it is found that the asymptotic convergence rate depends on the choice of the limiter used for the advective terms of the density equation, and some commonly used schemes are investigated. The variation with obstacle width of the influence of the stratification on the flow field is described and the results of the three-dimensional computations are compared with those of the corresponding computation of flow over a two-dimensional obstacle (of effectively infinite width). Also given are the results of time-dependent computations for three-dimensional flows under conditions of strong static stability when lee-wave propagation is present and the multigrid algorithm is used to compute the flow at each time step.     

MPDATA: A Finite-Difference Solver for Geophysical Flows

This article is a review of MPDATA, a class of methods for the numerical simulation of fluid flows based on the sign-preserving properties of upstream differencing. MPDATA was designed originally as an inexpensive alternative to flux-limited schemes for evaluating the advection of nonnegative thermodynamic variables (such as liquid water or water vapour) in atmospheric models. During the last decade, MPDATA has evolved from a simple advection scheme to a general approach for integrating the conservation laws of geophysical fluids on micro-to-planetary scales. The purpose of this paper is to summarize the basic concepts leading to a family of MPDATA schemes, to review existing MPDATA options, and to demonstrate the use of MPDATA to effectively construct two distinct types of models (elastic and anelastic) for complex geophysical flows.     

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.     

A Particle-Mesh Method for the Shallow Water Equations Near Geostrophic Balance

In this paper we outline a new particle-mesh method for rapidly rotating shallow water flows based on a set of regularized equations of motion. The time-stepping method uses an operator splitting of the equations into an Eulerian gravity wave part and a Lagrangian advection part. An essential ingredient is the advection of absolute vorticity by means of translated radial basis functions. We show that this implies exact conservation of enstrophy. The method is tested on two model problems based on the qualitative features of the solutions obtained (i.e., dispersion or smoothness of potential vorticity contours) as well as on the increase in mean divergence level.     

Unstructured Adaptive Grid Flow Simulations of Inert and Reactive Gas Mixtures

Unstructured adaptive grid flow simulation is applied to the calculation of high-speed compressible flows of inert and reactive gas mixtures. In the present case, the flowfield is simulated using the 2-D Euler equations, which are discretized in a cell-centered finite volume procedure on unstructured triangular meshes. Interface fluxes are calculated by a Liou flux vector splitting scheme which has been adapted to an unstructured grid context by the authors. Physicochemical properties are functions of the local mixture composition, temperature, and pressure, which are computed using the CHEMKIN-II subroutines. Computational results are presented for the case of premixed hydrogen–air supersonic flow over a 2-D wedge. In such a configuration, combustion may be triggered behind the oblique shock wave and transition to an oblique detonation wave is eventually obtained. It is shown that the solution adaptive procedure implemented is able to correctly define the important wave fronts. A parametric analysis of the influence of the adaptation parameters on the computed solution is performed.     

Comparison of Several Spatial Discretizations for the Navier–Stokes Equations

Grid convergence studies for subsonic and transonic flows over airfoils are presented in order to compare the accuracy of several spatial discretizations for the compressible Navier–Stokes equations. The discretizations include the following schemes for the inviscid fluxes: (1) second-order-accurate centered differences with third-order matrix numerical dissipation, (2) the second-order convective upstream split pressure scheme (CUSP), (3) third-order upwind-biased differencing with Roe's flux-difference splitting, and (4) fourth-order centered differences with third-order matrix numerical dissipation. The first three are combined with second-order differencing for the grid metrics and viscous terms. The fourth discretization uses fourth-order differencing for the grid metrics and viscous terms, as well as higher-order approximations near boundaries and for the numerical integration used to calculate forces and moments. The results indicate that the discretization using higher-order approximations for all terms is substantially more accurate than the others, producing less than two percent numerical error in lift and drag components on grids with less than 13,000 nodes for subsonic cases and less than 18,000 nodes for transonic cases. Since the cost per grid node of all of the discretizations studied is comparable, the higher-order discretization produces solutions of a given accuracy much more efficiently than the others.     

Mass Flux Schemes and Connection to Shock Instability

We analyze numerical mass fluxes with an emphasis on their capability for accurately capturing shock and contact discontinuities. The study of mass flux is useful because it is the term common to all conservation equations and the numerical diffusivity introduced in it bears a direct consequence to the prediction of contact (stationary and moving) discontinuities, which are considered to be the limiting case of the boundary layer. We examine several prominent numerical flux schemes and analyze the structure of numerical diffusivity. This leads to a detailed investigation into the cause of certain catastrophic breakdowns by some numerical flux schemes. In particular, we identify the dissipative terms that are responsible for shock instabilities, such as the odd–even decoupling and the so-called “carbuncle phenomenon”. As a result, we propose a conjecture stating the connection of the pressure difference term to these multidimensional shock instabilities and hence a cure to those difficulties. The validity of this conjecture has been confirmed by examining a wide class of upwind schemes. The conjecture is useful to the flux function development, for it indicates whether the flux scheme under consideration will be afflicted with these kinds of failings. Thus, a class of shock-stable schemes can be identified. Interestingly, a shock-stable scheme's self-correcting capability is demonstrated with respect to carbuncle-contaminated profiles for flows at both low supersonic and high Mach numbers.     

A Higher-Order Boundary Treatment for Cartesian-Grid Methods

The Euler equations describe the flow phenomena of compressible inviscid gas dynamics. We simulate such flows using a higher-order Cartesian-grid method, together with a special treatment for the cells cut by the boundary of an object. A new method for the treatment of the boundary is described where these cut boundary cells are maintained as whole cells rather than as cut cells, thus avoiding stability problems. The method is second-order accurate in one dimension and higher-order accurate in two dimensions but not strictly conservative; however, we show that this error in the conservation does not lead to spurious phenomena on some representative test calculations. The advantages of the new boundary treatment are that it is higher-order accurate, that it is independent of the applied method, and that it is simple.     

An open magnet utilizing ferro-refraction current magnification

Ferro-refraction is the field magnification that is obtained when a current segment is near a high magnetic permeable boundary. It is shown that ferro-refraction may be used in the design of magnets for NMR or MRI to increase the efficiency of these magnets. The field may be modeled analytically with the Biot--Savart law and the inclusion of mirror image currents. Ferro-refraction is particularly useful in the design of monohedral magnets, magnets producing a remote homogeneous region which have the magnetic sources arranged to one side. These magnets have also been called planar magnets. Two designs for a monohedral magnet which produce good agreement between experimental and analytic results are presented.     

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.     

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.     

Improved Spin-Echo-Edited NMR Diffusion Measurements

The need for simple and robust schemes for the analysis of ligand-protein binding has resulted in the development of diffusion-based NMR techniques that can be used to assay binding in protein solutions containing a mixture of several ligands. As a means of gaining spectral selectivity in NMR diffusion measurements, a simple experiment, the gradient modified spin-echo (GOSE), has been developed to reject the resonances of coupled spins and detect only the singlets in the (1)H NMR spectrum. This is accomplished by first using a spin echo to null the resonances of the coupled spins. Following the spin echo, the singlet magnetization is flipped out of the transverse plane and a dephasing gradient is applied to reduce the spectral artifacts resulting from incomplete cancellation of the J-coupled resonances. The resulting modular sequence is combined here with the BPPSTE pulse sequence; however, it could be easily incorporated into any pulse sequence where additional spectral selectivity is desired. Results obtained with the GOSE-BPPSTE pulse sequence are compared with those obtained with the BPPSTE and CPMG-BPPSTE experiments for a mixture containing the ligands resorcinol and tryptophan in a solution of human serum albumin.     

GAMMA simulations of stray field responses: slice thickness and pulse calibration

A stray field (STRAFI) module has been added to the GAMMA magnetic resonance simulation platform in order to facilitate computational investigations of NMR experiments in large static field gradients that are on the order of 50 T/m. The package has been used to examine system response during echo trains generated by the application of shaped pulses. The associated echo amplitude maxima and effective slice thickness are presented. A new accurate method for STRAFI pulse calibration based on relative echo amplitudes is proposed.     

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.