Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Constants”

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 343 records · Page 19

Strong coupling constant from 1-loop improved static energy

The static energy is an excellent observable for extracting the strong coupling $α_s$ on the lattice. For short distances, the static energy can be calculated both on the lattice using Wilson line correlators, and with perturbation theory up to three loop accuracy with leading ultrasoft log resummation. Comparing the perturbative expression and lattice data allows for precise determination of $α_s$. We present early results for 1-loop lattice perturbation theory improvement of the Wilson loop and show how it improves the $α_s$ extraction. We present a preliminary reanalysis of the TUMQCD (2+1)-flavor QCD data.

Leino, Viljami [U. Southern Denmark, QTC; Southern↗

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗

Cosmological Constraints from Gas Mass Fractions of Massive, Relaxed Galaxy Clusters

We present updated cosmological constraints from measurements of the gas mass fractions (f(gas)) of massive, dynamically relaxed galaxy clusters. Our new data set has greater leverage on models of dark energy, thanks to the addition of the Perseus cluster at low redshifts, two new clusters at redshifts z greater than or similar to 1, and significantly longer observations of four clusters at 0.6 < z < 0.9. Our low-redshift (z < 0.16) f(gas) data, combined with the cosmic baryon fraction measured from the cosmic microwave background (CMB), imply a Hubble constant of h = 0.722 +/- 0.067. Combining the full f(gas) data set with priors on the cosmic baryon density and the Hubble constant, we constrain the dark energy density to be Omega(Lambda) = 0.865 +/- 0.119 in non-flat Lambda cold dark matter (cosmological constant) models, and its equation of state to be in flat, constant-w models, respectively 41percent and 29percent tighter than our previous work, and comparable to the best constraints available from other probes. Combining f(gas), CMB, supernova, and baryon acoustic oscillation data, we also constrain models with global curvature and evolving dark energy. For the massive, relaxed clusters employed here, we find the scaling of f(gas) with mass to be consistent with a constant, with an intrinsic scatter that corresponds to just similar to 3percent in distance.

clusters: general – X-rays:↗

Particle dynamics and geometric optics in Chern–Simons black holes

Highlights: • Black holes in Chern–Simons modified gravity have been studied. • The effects of the coupling constant of the Chern–Simons gravity on some physical quantities have been investigated. • The dependence of Hawking mass and Penrose process on the coupling constant is studied. • Rotation of the polarization vector is also studied for dependence on the Chern–Simons coupling constant. In this paper we study the effects of the coupling constant of the Chern–Simons modified gravity on some physical properties of black holes. The Hawking mass is one of the proposed definitions of quasilocal mass. We find that, for slowly rotating Chern–Simons black holes, the Hawking mass is independent of the coupling constant. Next, we show the dependence on the centre of mass energy, for two neutral colliding particles, of coupling constant and the rotation parameter. We also investigate energy extraction through Penrose process and find that the energy gain and efficiency of the Penrose process are independent of this coupling constant. Rotation of the polarization vector is also studied for dependence on the Chern–Simons coupling constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Study of the Radiation Tolerance of CVD Diamond to 70 MeV Protons, Fast Neutrons and 200 MeV Pions

We measured the radiation tolerance of commercially available diamonds grown by the Chemical Vapor Deposition process by measuring the charge created by a 120 GeV hadron beam in a 50 μm pitch strip detector fabricated on each diamond sample before and after irradiation. We irradiated one group of samples with 70 MeV protons, a second group of samples with fast reactor neutrons (defined as energy greater than 0.1 MeV), and a third group of samples with 200 MeV pions, in steps, to (8.8±0.9) × 1015 protons/cm2, (1.43±0.14) × 1016 neutrons/cm2, and (6.5±1.4) × 1014 pions/cm2, respectively. By observing the charge induced due to the separation of electron–hole pairs created by the passage of the hadron beam through each sample, on an event-by-event basis, as a function of irradiation fluence, we conclude all datasets can be described by a first-order damage equation and independently calculate the damage constant for 70 MeV protons, fast reactor neutrons, and 200 MeV pions. We find the damage constant for diamond irradiated with 70 MeV protons to be 1.62±0.07(stat)±0.16(syst)× 10−18 cm2/(p μm), the damage constant for diamond irradiated with fast reactor neutrons to be 2.65±0.13(stat)±0.18(syst)× 10−18 cm2/(n μm), and the damage constant for diamond irradiated with 200 MeV pions to be 2.0±0.2(stat)±0.5(syst)× 10−18 cm2/(π μm). The damage constants from this measurement were analyzed together with our previously published 24 GeV proton irradiation and 800 MeV proton irradiation damage constant data to derive the first comprehensive set of relative damage constants for Chemical Vapor Deposition diamond. We find 70 MeV protons are 2.60 ± 0.29 times more damaging than 24 GeV protons, fast reactor neutrons are 4.3 ± 0.4 times more damaging than 24 GeV protons, and 200 MeV pions are 3.2 ± 0.8 more damaging than 24 GeV protons. We also observe the measured data can be described by a universal damage curve for all proton, neutron, and pion irradiations we performed of Chemical Vapor Deposition diamond. Finally, we confirm the spatial uniformity of the collected charge increases with fluence for polycrystalline Chemical Vapor Deposition diamond, and this effect can also be described by a universal curve.

36 MATERIALS SCIENCE↗

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗

General Method for Determining Redox Potentials without Electrolyte

A novel method to determine redox potentials without electrolyte is presented here. The method is based on a new ability to determine the dissociation constant, $K°_d$, for ion pairs formed between any radical anion and any inert electrolyte counterion. These dissociation constants can be used to determine relative shifts of redox potential as a function of electrolyte concentration, connecting referenced potentials determined with electrochemistry (with 0.1 M electrolyte) to electrolyte-free values. Pulse radiolysis created radical anions enabling determination of equilibrium constants for electron transfer between anions of donor and acceptor molecules as a function of electrolyte concentration in THF. The measurements determined “composite equilibrium constants”, $K_{eqC}$, which contain information about the dissociation constant for the electrolyte cations, X + , with the radical anions of both the donor, $K°_d$(D –• ,X + ) and the acceptor, $K°_d$(A –• ,X + ). Dissociation constants were obtained for a selection of radical anions with tetrabutylammonium (TBA + ). The electrolyte was found to shift the reduction potentials of small molecules 1-methylpyrene and trans-stilbene by close to +130 mV whereas oligo-fluorenes and polyfluorenes experienced shifts of only (+25 ± 6) mV due to charge delocalization weakening the ion pair. These shifts for reduction of aromatic hydrocarbon molecules are smaller than shifts of +232 and +451 mV seen previously for benzophenone radical anion with TBA + and Na + respectively where the charge on the radical anion is localized largely on one C=O bond, thus forming a more tightly bound ion pair.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface Dynamics of Charge Transport in LaAlO 3 /SrTiO 3 with Time-Resolved Kelvin Probe Force Microscopy

Time-resolved Kelvin probe force microscopy (tr-KPFM) was used to map surface potential dynamics in LAO/STO under its critical thickness, where it acts as an insulator, as polarizing and depolarizing electric fields were applied to planar electrodes in ambient conditions. Two time constants, for polarization and relaxation (i.e., depolarization) each, were derived and were shown to be independent of the direction of electrical field, revealing the reproducibility and recovery of the surfaces. Ambient aging or recent electrical history did not change the surface kinetics as determined by tr-KPFM. During polarization, the faster time constant (~3–6 s) is likely associated with the transport of H+ ions through the full length of the channel. The slower time constant (~40–50 s), on the other hand, is observed across only half the channel width near the higher-potential–more positive–electrode, which is likely associated with an electrochemical reaction such as the dissociation of water. This slower time constant drops to ~10–20 s toward the middle of the channel which could represent the transport of OH– ions. Furthermore, during relaxation, the fast time constant is consistent with polarization and is correlated to the same mechanism while the slower time constant is higher than it was during polarization and may represent the reabsorption of moisture from the atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical Understanding of Effects of Operating Modes on the Performance Durability of Solid Oxide Cells: A Comparison between Potentiostatic and Galvanostatic Operations

In solid oxide cells (SOCs), the choice between galvanostatic (constant current) and potentiostatic (constant voltage) modes does not significantly affect the performance of SOCs as long as the cell components remain unchanged and intact. However, the degradation of cell components, which leads to changes in electrochemical and physical properties of the cell, elevates the importance of the selected operating mode. This paper aims to investigate the effects of galvanostatic and potentiostatic operating modes on the evolving properties of SOCs and their subsequent influence on performance durability. Employing non-equilibrium thermodynamic analysis, a crucial approach for understanding the degradation phenomena within an active electrochemical system, this study aims to provide in-depth insights into how these operating modes affect the longevity and efficacy of SOCs. Key findings include: In cases where oxygen electrode (OE) degradation is accelerated by higher partial pressure of oxygen ( p O 2 ), operating under constant voltage electrolysis can mitigate the high p O 2 at the OE|electrolyte (OE|EL) interface. Conversely, if OE degradation occurs more rapidly under a lower p O 2 , constant current electrolysis is more effective in suppressing degradation by achieving a high p O 2 at the OE|EL interface. For degradation of the fuel electrode (FE) due to higher p O 2 , constant current electrolysis is beneficial for more stable performance, which helps maintain low p O 2 at the FE|EL interface. When FE degradation is accelerated by lower p O 2 , constant voltage electrolysis can avert low p O 2 at the FE|EL interface. In practical scenarios, more complex degradation mechanisms come into play, especially when p O 2 significantly deviates from initial conditions. Degradation in one electrode can influence p O 2 in the other electrode, a phenomenon more pronounced in potentiostatic than in galvanostatic electrolysis.

Electrochemistry↗

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of Critical Compression Ratio on Rating Gasoline Knock Propensity

It is common practice in the automotive industry to explore the knock limits of fuels on an engine by a comparison of the knock limited spark advance (KLSA) at threshold knock intensity. However, the knock propensity of gasolines can be rated by changing one of three metrics on a variable compression ratio Cooperative Fuels Research (CFR) octane rating engine while holding the other two variables constant: knock intensity, spark timing, and critical compression ratio. The operational differences between the standard research octane number (RON) rating and modern engine operation have been explored in three parts. The first part focused on the effects of lambda and knock characterization. The second part studied the effects of spark timing. This third part explores the knock ratings of several gasolines by comparing the critical compression ratios at constant combustion phasing and knock intensity. The threshold knock intensity was based on the standard octane rating D1 pickup or by maximum amplitude of pressure oscillations (MAPO) measured by a piezoelectric cylinder pressure transducer. Several Fuels for Advanced Combustion Engines (FACE) gasolines, primary reference fuels (PRFs), and toluene standardization fuels (TSFs) were tested on a CFR octane rating engine with advanced data acquisition equipment and a piezoelectric cylinder pressure transducer. These tests deviated from the ASTM D2699 standard octane rating procedure. For each test fuel, the CFR engine was operated at stoichiometry at a constant combustion phasing (CA50) and the compression ratio was modified until a threshold knock intensity was realized. It was found that the chemical composition of the fuels affected the relationship of critical compression ratios between the D1 knockmeter and piezoelectric pressure transducer knock intensity thresholds, as well as the measured combustion maximum pressure rise rate and spark timing setting for constant CA50. For highly aromatic fuels tested at a constant MAPO knock intensity threshold, it was found that the maximum pressure rise rate was two to three times higher than that of highly paraffinic fuels with similar RON and the spark advance was several crank angle degrees less for constant combustion phasing.

Kolodziej, Christopher P↗

Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models

Abstract. Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics–dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.

58 GEOSCIENCES↗