Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid water path”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

Liquid propulsion turbomachinery model testing

For the past few years an extensive experimental program to understand the fluid dynamics of the Space Shuttle Main Engine hot gas manifold has been in progress. This program includes models of the Phase II and II+ manifolds for each of the air and water flow facilities, as well as two different turbine flow paths and two simulated power levels for each manifold. All models are full-scale (geometric). The water models are constructed partially of acrylic to allow flow visualization. The intent of this paper is to discuss the concept, including the test objectives, facilities, and models, and to summarize the data for an example configuration, including static pressure data, flow visualization, and the solution of a specific flow problem.

Mcdaniels, David M.↗

Dimension-free path-integral molecular dynamics without preconditioning

Convergence with respect to imaginary-time discretization (i.e., the number of ring-polymer beads) is an essential part of any path-integral-based molecular dynamics (MD) calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular dynamics—including essentially all existing ring-polymer molecular dynamics (RPMD) and thermostatted RPMD (T-RPMD) methods—is that for a given MD time step, the overlap between the exact ring-polymer Boltzmann–Gibbs distribution and that sampled using MD becomes zero in the infinite-bead limit. This has clear implications for hybrid Metropolis Monte Carlo/MD sampling schemes, and it also causes the divergence with bead number of the primitive path-integral kinetic-energy expectation value when using standard RPMD or T-RPMD. We show that these and other problems can be avoided through the introduction of “dimension-free” numerical integration schemes for which the sampled ring-polymer position distribution has non-zero overlap with the exact distribution in the infinite-bead limit for the case of a harmonic potential. Most notably, we introduce the BCOCB integration scheme, which achieves dimension freedom via a particular symmetric splitting of the integration time step and a novel implementation of the Cayley modification [R. Korol et al., J. Chem. Phys. 151, 124103 (2019)] for the free ring-polymer half-steps. More generally, we show that dimension freedom can be achieved via mollification of the forces from the external physical potential. The dimension-free path-integral numerical integration schemes introduced here yield finite error bounds for a given MD time step, even as the number of beads is taken to infinity; these conclusions are proven for the case of a harmonic potential and borne out numerically for anharmonic systems that include liquid water. The numerical results for BCOCB are particularly striking, allowing for nearly three-fold increases in the stable time step for liquid water with respect to the Bussi–Parrinello (OBABO) and Leimkuhler (BAOAB) integrators, while introducing negligible errors in the calculated statistical properties and absorption spectrum. Importantly, the dimension-free, non-preconditioned integration schemes introduced here preserve ergodicity and global second-order accuracy, and they remain simple, black-box methods that avoid additional computational costs, tunable parameters, or system-specific implementations.

Korol, Roman (ORCID:0000000193076351)↗

Modeling Liquid Water by Climbing up Jacob’s Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials

Within the framework of Kohn–Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob’s ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric–isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Operational MODIS Cloud Optical and Microphysical Property Product: Overview of the Collection 6 Algorithm and Preliminary Results

Operational Moderate Resolution Imaging Spectroradiometer (MODIS) retrievals of cloud optical and microphysical properties (part of the archived products MOD06 and MYD06, for MODIS Terra and Aqua, respectively) are currently being reprocessed along with other MODIS Atmosphere Team products. The latest "Collection 6" processing stream, which is expected to begin production by summer 2012, includes updates to the previous cloud retrieval algorithm along with new capabilities. The 1 km retrievals, based on well-known solar reflectance techniques, include cloud optical thickness, effective particle radius, and water path, as well as thermodynamic phase derived from a combination of solar and infrared tests. Being both global and of high spatial resolution requires an algorithm that is computationally efficient and can perform over all surface types. Collection 6 additions and enhancements include: (i) absolute effective particle radius retrievals derived separately from the 1.6 and 3.7 !-lm bands (instead of differences relative to the standard 2.1 !-lm retrieval), (ii) comprehensive look-up tables for cloud reflectance and emissivity (no asymptotic theory) with a wind-speed interpolated Cox-Munk BRDF for ocean surfaces, (iii) retrievals for both liquid water and ice phases for each pixel, and a subsequent determination of the phase based, in part, on effective radius retrieval outcomes for the two phases, (iv) new ice cloud radiative models using roughened particles with a specified habit, (v) updated spatially-complete global spectral surface albedo maps derived from MODIS Collection 5, (vi) enhanced pixel-level uncertainty calculations incorporating additional radiative error sources including the MODIS L1 B uncertainty index for assessing band and scene-dependent radiometric uncertainties, (v) and use of a new 1 km cloud top pressure/temperature algorithm (also part of MOD06) for atmospheric corrections and low cloud non-unity emissivity temperature adjustments.

Platnick, Steven↗

Metastable precipitation and ion–extractant transport in liquid–liquid separations of trivalent elements

The extractant-assisted transport of metal ions from aqueous to organic environments by liquid–liquid extraction has been widely used to separate and recover critical elements on an industrial scale. While current efforts focus on designing better extractants and optimizing process conditions, the mechanism that underlies ionic transport remains poorly understood. Here, we report a nonequilibrium process in the bulk aqueous phase that influences interfacial ion transport: the formation of metastable ion–extractant precipitates away from the liquid–liquid interface, separated from it by a depletion region without precipitates. Although the precipitate is soluble in the organic phase, the depletion region separates the two and ions are sequestered in a long-lived metastable state. Since precipitation removes extractants from the aqueous phase, even extractants that are sparingly soluble in water will continue to be withdrawn from the organic phase to feed the aqueous precipitation process. Solute concentrations in both phases and the aqueous pH influence the temporal evolution of the process and ionic partitioning between the precipitate and organic phase. Aqueous ion–extractant precipitation during liquid–liquid extraction provides a reaction path that can influence the extraction kinetics, which plays an important role in designing advanced processes to separate rare earths and other minerals.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Observed and Simulated Variability of Droplet Spectral Dispersion in Convective Clouds Over the Amazon

Abstract In this study, the variability of the spectral dispersion of droplet size distributions (DSDs) in convective clouds is investigated. Analyses are based on aircraft measurements of growing cumuli near the Amazon basin, and on numerical simulations of an idealized ice‐free cumulus. In cleaner clouds, the relative dispersion , defined as the ratio of the standard deviation to the mean value of the droplet diameter, is negatively correlated with the ratio of the cloud water content ( ) to the adiabatic liquid water content ( ), while no strong correlation between and is seen in polluted clouds. Bin microphysics numerical simulations suggest that these contrasting behaviors are associated with the effect of collision‐coalescence in cleaner clouds, and secondary droplet activation in polluted clouds, in addition to the turbulent mixing of parcels that experienced different paths within the cloud. Collision‐coalescence simultaneously broadens the DSDs and decreases , explaining the inverse relationship between and in cleaner clouds. Secondary droplet activation broadens the DSDs but has little direct impact on . The combination of a rather modest DSD broadening due to weak collision‐coalescence with enhanced droplet activation in both diluted and highly undiluted cloud regions may contribute to maintain a relatively uniform within polluted clouds. These findings can be useful for parameterizing the shape parameter ( ) of gamma DSDs in bulk microphysics cloud‐resolving models. It is shown that emulating the observed relationship improves the estimation of the collision‐coalescence rate in bulk microphysics simulations compared to the bin simulations.

54 ENVIRONMENTAL SCIENCES↗

Flow Straightener for a Rotating-Drum Liquid Separator

A flow straightener has been incorporated into a rotary liquid separator that originally comprised an inlet tube, a shroud plate, an impeller, an inner drum, an outer drum, a housing, a pitot tube, and a hollow shaft motor. As a consequence of the original geometry of the impeller, shroud, inner drum, and hollow shaft, swirl was created in the airflow inside the hollow shaft during operation. The swirl speed was large enough to cause a significant pressure drop. The flow straightener consists of vanes on the back side of the shroud plate. These vanes compartmentalize the inside of the inner drum in such a way as to break up the flow path and thereby stop the air from swirling; as a result, the air enters the hollow shaft with a predominantly axial velocity instead of a swirl. Tests of the rotary liquid separator at an airflow rate of 10 cu ft/min (0.0047 cu m/s) revealed that the dynamic pressure drop was 8 in. of water (approx.=2 kPa) in the absence of the flow straightener and was reduced to 1 in. of water (approx.=0.25 kPa) in the presence of the flow straightener.

O'Coin, James R.↗

A Machine‐Learning‐Assisted Stochastic Cloud Population Model as a Parameterization of Cumulus Convection

Abstract A machine‐learning‐assisted stochastic cloud population model is coupled with the Advanced Research Weather Research and Forecasting (WRF) model to represent fluctuations in the cloud‐base mass flux associated with the life cycles and interactions among cumulus convection cells. In this cloud population model, the size distribution and the associated cloud‐base mass flux of the convective cells are related to their previous state and to the change in the total convective area via a transition function. The convective area tendency in turn is assumed to depend on the cloud‐base mass flux that is resolved by the host WRF model. The transition function is represented by a single hidden‐layer neural network trained by the evolution of convective cell size distributions in a 1‐km grid‐spacing WRF simulation run over the Australian Monsoon region. At every grid point of the host model, the cloud population model predicts the cell size and cloud‐base mass flux distributions from which a random sample of cells is fed to an entraining parcel model that calculates precipitation as well as the associated liquid water potential temperature and total moisture tendencies. These tendencies are averaged over the cells and provided to the host model. Several regional simulations are performed over tropical and midlatitude domains to test this as a potential approach to scale‐aware parameterization. It is shown that such an approach could be a new promising path to simulating realistic precipitation statistics and propagation of precipitation associated with the Madden‐Julian Oscillation while maintaining realistic depictions of the diurnal cycle over both land and ocean.

54 ENVIRONMENTAL SCIENCES↗

Porous mesh manifold for enhanced boiling performance

High-performance electronics are continuously demanding cooling of higher heat fluxes. Phase-change cooling, including pool boiling, is a useful approach to address this challenge; however, competition between liquid and vapor flows generally limit the heat fluxes that can be dissipated. A range of strategies to control these flows have been investigated previously, including capillary guides. Here a manifold structure formed from a metallic mesh is investigated to control the disposition of liquid and vapor phases above a pool fed boiling surface enhanced with porous structures. Copper mesh forms defined liquid flow paths, using capillary action to guide and distribute liquid evenly over the heated surface, along with open channels to facilitate vapor escape. The mesh provides a novel structure for liquid guidance that imposes low resistance to liquid flow while occluding a minimal area of heated surface underneath. The manifold performance is characterized in boiling fed by a pool of water above a laser-textured aluminum nitride heat dissipation surface with pin–fin structures having heights of 110 µm and spacing of 30 µm with a heated area of 5 mm x 5 mm. A maximum heat flux of 490 W/cm 2 is reached with the manifold in the pool fed configuration, representing an increase of more than 65% over the porous pin fin surface alone. The maximum stable superheat observed for the manifold of 36K is 14K higher than that for the porous surface without the manifold. The factors limiting performance of the manifold are analyzed. High superheat is attributed to partial flooding of the boiling surface as suggested by the reduction in superheat using external suction. Similar systems and structures for enhanced two-phase cooling are compared.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗

Nanofluidic Size-Exclusion Chromatograph

Efforts are under way to develop a nanofluidic size-exclusion chromatograph (SEC), which would be a compact, robust, lightweight instrument for separating molecules of interest according to their sizes and measuring their relative abundances in small samples. About as large as a deck of playing cards, the nanofluidic SEC would serve, in effect, as a laboratory on a chip that would perform the functions of a much larger, conventional, bench-top SEC and ancillary equipment, while consuming much less power and much smaller quantities of reagent and sample materials. Its compactness and low power demand would render it attractive for field applications in which, typically, it would be used to identify and quantitate a broad range of polar and nonpolar organic compounds in soil, ice, and water samples. Size-exclusion chromatography is a special case of high-performance liquid chromatography. In a conventional SEC, a sample plug is driven by pressure along a column packed with silica or polymer beads that contain uniform nanopores. The interstices between, and the pores in, the beads collectively constitute a size-exclusion network. Molecules follow different paths through the size-exclusion network, such that characteristic elution times can be related to sizes of molecules: basically, smaller molecules reach the downstream end of the column after the larger ones do because the smaller ones enter minor pores and stay there for a while, whereas the larger ones do not enter the pores. The volume accessible to molecules gradually diminishes as their size increases. All molecules bigger than a pore size elute together. For most substances, the elution times and sizes of molecules can be correlated directly with molecular weights. Hence, by measuring the flux of molecules arriving at the downstream end as a function of time, one can obtain a liquid mass spectrum for the molecules present in a sample over a broad range of molecular weights.

Feldman, Sabrina↗

Characterization of Liquid Fuel Evaporation of a Lifted Methanol Spray Flame in a Vitiated Coflow Burner

An experimental investigation of lifted spray flames in a coflow of hot, vitiated gases is presented. The vitiated coflow burner is a spray flame that issues into a coaxial flow of hot combustion products from a lean, premixed H2/Air flame. The spray flame in a vitiated coflow emulates the combustion that occurs in many advanced combustors without the detailed fluid mechanics. Two commercially available laser diagnostic systems are used to characterize the spray flame and to demonstrate the vitiated coflow burner's amenability to optical investigation. The Ensemble Particle Concentration and Size (EPCS) system is used to measure the path-average droplet size distribution and liquid volume fraction at several axial locations while an extractive probe instrument named the Real-time Fuel-air Analyzer (RFA) is used to measure the air to fuel ratio downstream of the spray nozzle with high temporal and spatial resolution. The effect of coflow conditions (stoichiometry) and dilution of the fuel with water was studied with the EPCS optical system. As expected, results show that water retards the evaporation and combustion of fuels. Measurements obtained by the RFA extractive probe show that while the Delavan manufactured nozzle does distribute the fuel over the manufacturer specified spray angle, it unfortunately does not distribute the fuel uniformly, providing conditions that may result in the production of unwanted NOx. Despite some limitations due to the inherent nature of the experimental techniques, the two diagnostics can be readily applied to spray flames in the vitiated coflow environment.

Cabra, Ricardo↗

Surface and Atmospheric Contributions to Passive Microwave Brightness Temperatures

Physically-based passive microwave precipitation retrieval algorithms require a set of relationships between satellite observed brightness temperatures (TB) and the physical state of the underlying atmosphere and surface. These relationships are typically non-linear, such that inversions are ill-posed especially over variable land surfaces. In order to better understand these relationships, this work presents a theoretical analysis using brightness temperature weighting functions to quantify the percentage of the TB resulting from absorption/emission/reflection from the surface, absorption/emission/scattering by liquid and frozen hydrometeors in the cloud, the emission from atmospheric water vapor, and other contributors. The results are presented for frequencies from 10 to 874 GHz and for several individual precipitation profiles as well as for three cloud resolving model simulations of falling snow. As expected, low frequency channels (<89 GHz) respond to liquid hydrometeors and the surface, while the higher frequency channels become increasingly sensitive to ice hydrometeors and the water vapor sounding channels react to water vapor in the atmosphere. Low emissivity surfaces (water and snow-covered land) permit energy downwelling from clouds to be reflected at the surface thereby increasing the percentage of the TB resulting from the hydrometeors. The slant path at a 53deg viewing angle increases the hydrometeor contributions relative to nadir viewing channels and show sensitivity to surface polarization effects. The TB percentage information presented in this paper answers questions about the relative contributions to the brightness temperatures and provides a key piece of information required to develop and improve precipitation retrievals over land surfaces.

Gail Skofronick-Jackson↗

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗

Monitoring of Ionic Silver Biocide with Capacitively-Coupled Contactless Conductivity Detection (C4D)

Monitoring of biocide dosing in spacecraft potable water systems is desirable to ensure dosing system operation and resultant microbial control. Electrolytic conductivity is linearly related to the concentration ionic silver (Ag+, the baseline biocide for future NASA missions) chemically dosed into highly purified water, and could serve as a proxy measurement. Traditional conductivity cells rely on immersed metal or graphite electrodes, which may potentially be subject to significant current-induced or auto-galvanic deposition of Ag+. Furthermore, oxide fouling may potentially induce measurement errors. Capacitively-coupled contactless conductivity detection (C4D) is a promising alternative. C4D uses electrodes placed outside an inert, non-conductive flow path, potentially mitigating such problems. We identified and purchased a commercial C4D system with claimed performance (conductivity range: 20 nS/cm to 4 µS/cm) appropriate for this application. We characterized the instrument’s performance with conductivity standards and Ag+ solutions. We began investigations considering practical and long-term detector operation, including potential current-induced deposition of conductive deposits at the tubing wall. In addition, we designed a side-stream flow system to provide liquid flow through the narrow tubing (1/16” outer diameter) required by the detector and characterized the performance of the combined C4D detector/flow system. Finally, we completed theoretical modeling of expected conductivity change with multiple dosing techniques and water chemistries.

Lance Dean Delzeit↗

Simplified Method for the In Situ Collection and Laboratory Analysis of Cosmogenic Tracers (Sulfur-35 and Sodium-22) to Determine Residence Time Distributions and Water Ages

The use of cosmogenically produced sulfur-35 (T 1/2 = 87 days) and sodium-22 (T 1/2 = 2.6 years) as intrinsic tracers can provide valuable information on catchment hydrology, flow paths, and subsurface storage. A new and straightforward method was created to determine the activities of both 35 S and 22 Na in various water sources by pumping large volumes (up to 1000 L) of water through cation- and anion-exchange resin columns in the field to collect sodium and sulfate ions and simple chemistry in the lab. Samples are counted for 35 S using liquid scintillation counting (LSC) and for 22 Na via γ spectroscopy. Our novel in situ method provides faster sample throughput as well as better counting statistics and lower detection limits. Both methods were successfully applied at the Southern Sierra Critical Zone Observatory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spiral Flow Separator

Proposed liquid-separating device relies on centrifugal force in liquid/liquid or liquid/solid mixture in spiral path. Operates in continuous flow at relatively high rates. Spiral tubes joined in sequence, with outlet tubes connected to joints. Cross-sectional areas of successive spiral tubes decreases by cross-sectional areas of outlet tubes. Centrifugal force pushes denser particles or liquids to outer edge of spiral, where removed from flow. Principle exploited to separate solids from wastewater, oil from fresh or salt water, or contaminants from salt water before evaporation. Also used to extract such valuable materials as precious metals from slurries.

Robertson, Glen A.↗

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗