SEARCH · Engineering Papers
Results for “ring model”
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.
A control-oriented combustion model framework for compression ignition engines operating on low-reactivity fuel
This work focuses on zero-dimensional modeling of the heat release rate in a compression ignition engine operating on gasoline-like fuels. Due to the properties of gasoline, such as high volatility and longer ignition delay than diesel, the injection strategies can vary significantly from the operation with conventional diesel fuel. Different injection strategies are commonly used to achieve varying degrees of in-cylinder stratification in order to shape the combustion event and maximize efficiency. The proposed zero-dimensional combustion model was developed to account for the different stages in combustion caused by the fuel stratification. As the ignition delay model is an integral part of the entire combustion process and significantly affects the prediction accuracy, special attention has been paid to local phenomena influencing ignition delay. A one-dimensional spray model by Musculus and Kattke was employed in conjunction with a Lagrangian tracking approach in order to estimate the local air–fuel ratio within the spray tip, as a proxy for reactivity. The local air–fuel ratio, in-cylinder temperature and pressure were used in an integral fashion to estimate the ignition delay. Heat release rates were modeled using first-order non-linear differential equations. The proposed combustion model was validated against experimental data of a heavy-duty compression ignition engine with up to three injection events at mostly 1038 r/min and 14 bar brake mean effective pressure. Further validation of the model was carried out at other engine loads and speeds. Model prediction errors in CA50 of less than 1 °CA across all conditions were found. Modeling results of other combustion metrics such as combustion duration and indicated mean effective pressure are also highly satisfactory. In addition, the model has been shown to be capable of estimating the ringing intensity for most conditions.
Simulating Isotropic Rayleigh Wave Horizontal to Vertical Ratios at Rock Valley
Elastic anisotropy in the subsurface can manifest from pre-existing structure or stress-induced cracking – both phenomena are of interest at the Rock Valley site in Nevada. The azimuthal dependence of Rayleigh wave ellipticity, also known as H/V, can reveal elastic anisotropy, but the effect of isotropic variability in H/V must first be accounted for. To that end, we simulate synthetic seismograms from a ring of sources in a fully 3-D isotropic velocity model that captures the geology at Rock Valley. We measure H/V from these waveforms and compare them with H/V predictions obtained using 1D models extracted at specific locations from the 3D model. We find that an average of the azimuthal H/V measurements from the 3-D simulations and 1-D predictions agree, but that the azimuthal variability of H/V is a large fraction of the mean value. These results estimate the noise character in azimuthal H/V so it can be factored out of future anisotropic simulations.
Ringing Artifact Reduction Method for Ultrasound Reconstruction Using Multi-Agent Consensus Equilibrium
Non-destructive characterization of multi-layered structures that can be accessed from only a single side is important for applications such as well-bore integrity inspection. Existing methods related to Synthetic Aperture Focusing Technique (SAFT) rapidly produce acceptable results but with significant artifacts. Recently, ultrasound model-based iterative reconstruction (UMBIR) approaches have shown significant improvements over SAFT. However, even these methods produce ringing artifacts due to the high fractional-bandwidth of the excitation signal.In this paper, we propose a ringing artifact reduction method for ultrasound image reconstruction that uses a multi-agent consensus equilibrium (RARE-MACE) framework. Our approach integrates a physics-based forward model that accounts for the propagation of a collimated ultrasonic beam in multi-layered media, a spatially varying image prior, and a denoiser designed to suppress the ringing artifacts that are characteristic of reconstructions from high-fractional bandwidth ultrasound sensor data. We test our method on simulated and experimental measurements and show substantial improvements in image quality compared to SAFT and UMBIR.
Accurate Thermochemistry of Complex Lignin Structures via Density Functional Theory, Group Additivity, and Machine Learning
A molecular-level understanding of lignin structures and bond dissociation energies could facilitate depolymerization technologies. Still, this information is currently limited due to the lack of databases and the simplification of surrogate models. Here, substitution effects on seven common linkages in lignin polymers are systematically investigated. An automated reaction network generator is employed to create a database of structures. A new group additivity (GA) model based on principal component analysis (PCA) descriptors is introduced and trained on gas-phase density functional theory data of 4100 species at the M06-2X/6-311++G(d,p) level. Hydrogen bonds, local steric, and nonaromatic ring contributions are also incorporated. Lastly, we improve the accuracy of the group additivity model to reach the G4 theory by computing a data set of 770 species at this level and using a data fusion approach.
Nonlinear Shear Rheology of Entangled Polymer Rings
Steady-state shear viscosity η($\bar{γ}$) of unconcatenated ring polymer melts as a function of the shear rate $\bar{γ}$ is studied by a combination of experiments, simulations, and theory. Experiments using polystyrenes with Z ≈ 5 and Z ≈ 11 entanglements indicate weaker shear thinning for rings compared to linear polymers exhibiting power law scaling of shear viscosity η ~ $\bar{γ}$ –0.56 ± 0.02 , independent of chain length, for Weissenberg numbers up to about 10 2 . Nonequilibrium molecular dynamics simulations using the bead-spring model reveal a similar behavior with η ~ $\bar{γ}$ -0.57 ± 0.08 for 4 ≤ Z ≤ 57. Viscosity decreases with chain length for high $\bar{γ}$. In our experiments, we see the onset of this regime, and in simulations, which we extended to Wi ~ 10 4 , the nonuniversality is fully developed. In addition to a naive scaling theory yielding for the universal regime η ~ $\bar{γ}$ –0.57 , we developed a novel shear slit model explaining many details of observed conformations and dynamics as well as the chain length-dependent behavior of viscosity at large $\bar{γ}$. The signature feature of the model is the presence of two distinct length scales: the size of tension blobs and much larger thickness of a shear slit in which rings are self-consistently confined in the velocity gradient direction and which is dictated by the size of a chain section with relaxation time 1/$\bar{γ}$. These two length scales control the two normal stress differences. In this model, the chain length-dependent onset of nonuniversal behavior is set by tension blobs becoming as small as about one Kuhn segment. This model explains the approximate applicability of the Cox–Merz rule for ring polymers.
Mimicking Extradiol Dioxygenase Reactivity on Iridium
Extradiol dioxygenases catalyze the cleavage of benzenediol (catechol) or vicinal aminophenols via oxygen atom insertion into the 2,3- C–C bond. These reactions are most often proposed to proceed through the migratory rearrangement of a d6 alkylperoxide, generating the corresponding d 6 ring expanded product. However, regiospecific insertion remains a rare outcome among synthetic model complexes. Here, a dioxygenated Ir complex (2) converts to (a) the paramagnetic metallatrioxolane (3) and (b) oxygen atom inserted products (4) and (5); all three complexes are third-row metal analogues of enzymatic intermediates. The conversion of 2 to 3 was triggered by an H• abstraction, generating a third-row metallatrioxolane that is one electron reduced from the canonical d 6 alkylperoxide. Alternatively, photolysis of 2 (467 nm) results in ring-expanded product 4, from which an H• can be abstracted to generate 5. This latter complex is also one electron reduced relative to the canonical ring expansion product. Because extradiol mechanisms were largely defined using Fe(II) metallocofactors, 3 and 5 may be especially relevant to the known Co(II) accepting variants. We propose these states became synthetically accessible via the incorporation of a catechol-like substrate into the larger, multidentate ligand L1. This perturbation enhances the affinity of L1 (and related intermediates) to the metal. The same perturbation may have also been key in characterizing the first κ 2 -bound, dianionic ortho ester ligand – the observed binding mode in the x-ray structure of 5. Broadly, we propose through this case study that connecting abiological dioxygen complexes to non-heme oxygenase-like intermediates might be a “roadmap” towards regiospecific aerobic oxygenations.
Mean-Field Ring Polymer Rates Using a Population Dividing Surface
Mean-field Ring Polymer Molecular Dynamics (MF-RPMD) offers a computationally efficient method for the simulation of reaction rates in multi-level systems. Previous work has established that, to model a nonadiabatic state-to-state reaction accurately, the dividing surface must be chosen to explicitly sample kinked ring polymer configurations where at least one bead is in a different electronic state than the others. Building on this, we introduce a population difference coordinate and a kink-constrained dividing surface, and we test the accuracy of the resulting mean-field rate theory on a series of bilinear coupling model systems. Here, we demonstrate that this new MF-RPMD rate approach is efficient to implement and quantitatively accurate for models over a wide range of driving forces, coupling strengths, and temperatures.
First-Principles Study of n -Butane Monomolecular Cracking and Dehydrogenation on Two-Dimensional-Zeolite Model Systems: Reaction Mechanisms and Effects of Spatial Confinement
Two-dimensional (2D) ultrathin (~0.5 nm) aluminosilicate bilayer films, consisting of hexagonal prisms (a.k.a. double 6-membered rings D6R) with acidic bridging hydroxyl groups exposed on the surface, have been previously synthesized on a Ru(0001) surface as a zeolite model system. These structures are helpful for mimicking zeolite catalysts with D6R building blocks, such as chabazite. We performed density functional theory calculations to investigate the monomolecular cracking and dehydrogenation of n-butane molecules over the acidic hydroxyl groups of the 2D model system and compared the reaction energetics with that in bulk chabazite. The intrinsic activation energy barrier is the highest for dehydrogenation and lowest for central C–C bond cracking in bulk chabazite. The trend of intrinsic energy barriers for dehydrogenation and terminal and central C–C bond cracking is reproduced on the 2D aluminosilicate film. Overall, the activation barriers are higher on the 2D film than in bulk chabazite due to the lack of confinement in the former. We further explored the effects of the zeolite channel size on the n-butane adsorption and monomolecular cracking using different bulk nanoporous zeolite frameworks (TON, MEL, MEI, and VFI). We found that as the confinement of channels decreases, n-butane adsorption becomes weaker, and the intrinsic energy barrier of terminal C–C cracking increases. The activation energy barriers (dehydrogenation and terminal and central C–C cracking) on the 2D bilayer film surface, which may be considered as zeolite cages at the infinite cage size limit, are close to that in VFI with a relatively large channel size. Comparing the reaction pathway of n-butane terminal C–C cracking in 3D nanocages and on the surface of the 2D aluminosilicate film revealed that stabilizing the transition states in the 3D nanocages is responsible for the decrease in the intrinsic energy barriers for bulk zeolites.
How Aqueous Solvation Impacts the Frequencies and Intensities of Infrared Absorption Bands in Flavin: The Quest for a Suitable Solvent Model
FTIR spectroscopy accompanied by quantum chemical simulations can reveal important information about molecular structure and intermolecular interactions in the condensed phase. Simulations typically account for the solvent either through cluster quantum mechanical (QM) models, polarizable continuum models (PCM), or hybrid quantum mechanical/molecular mechanical (QM/MM) models. Recently, we studied the effect of aqueous solvent interactions on the vibrational frequencies of lumiflavin, a minimal flavin model, using cluster QM and PCM models. Those models successfully reproduced the relative frequencies of four prominent stretching modes of flavin’s isoalloxazine ring in the diagnostic 1450–1750 cm−1 range but poorly reproduced the relative band intensities. Here, we extend our studies on this system and account for solvation through a series of increasingly sophisticated models. Only by combining elements of QM clusters, QM/MM, and PCM approaches do we obtain an improved agreement with the experiment. The study sheds light more generally on factors that can impact the computed frequencies and intensities of IR bands in solution.
Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry
The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.
Effects of Magnetic Fields on Gas Dynamics and Star Formation in Nuclear Rings
Nuclear rings at the centers of barred galaxies are known to be strongly magnetized. To explore the effects of magnetic fields on star formation in these rings and nuclear gas flows, we run magnetohydrodynamic simulations in which there is a temporally constant magnetized inflow to the ring, representing a bar-driven inflow. The mass inflow rate is 1M ⊙ yr –1 , and we explore models with a range of field strength in the inflow. We adopt the TIGRESS framework developed by Kim & Ostriker to handle radiative heating and cooling, star formation, and resulting supernova (SN) feedback. We find that magnetic fields are efficiently amplified in the ring due to rotational shear and SN feedback. Within a few 100 Myr, the turbulent component Btrb in the ring saturates at ~35 μG (in rough equipartition with the turbulent kinetic energy density), while the regular component Breg exceeds 50 μG. Expanding superbubbles created by clustered SN explosions vertically drag predominantly toroidal fields from near the midplane to produce poloidal fields in high-altitude regions. The growth of magnetic fields greatly suppresses star formation at late times. Simultaneously, strong magnetic tension in the ring drives radially inward accretion flows from the ring to form a circumnuclear disk in the central region; this feature is absent in the unmagnetized model.
S 8 ring content in liquid sulfur determined from x-ray diffraction
The S 8 ring content in liquid sulfur has been determined from the 4.47 Å peak in the x-ray pair distribution function over the temperature range −55 °C –330 °C. This peak represents the characteristic 3rd and 4th nearest atom-atom distances across an S 8 -ring. Three procedures were used to calculate the S 8 -ring fraction; scaling to crystalline β-S, a linear combination fit and a direct integration method. All methods gave similar values, indicating that approximately one third of the S 8 -rings break up upon melting into a distribution of short chains with an average length of 5 ± 2 atoms. Our results contradict established models of low temperature sulfur that indicate below the λ-transition the liquid comprises almost entirely of S8-rings. In addition, we find the S 8 -ring fraction increases to 76 ± 10% at −30 °C in the supercooled liquid, and decreases to 23% ± 3% at 330 °C in the high temperature melt. First principle molecular dynamics calculations confirm the onset of dynamic ring distortion via broad bond angle distributions, associated with the interstitial region between the first and second peaks in the pair distribution function. Our results are consistent with the presence of a propensity of short chains in the low temperature liquid, that are required for the bond switching mechanism responsible for the polymerization process to occur at the λ-transition temperature, T λ .
A unique, ring-like radio source with quadrilateral structure detected with machine learning
ABSTRACT We report the discovery of a unique object in the MeerKAT Galaxy Cluster Legacy Survey (MGCLS) using the machine learning anomaly detection framework astronomaly. This strange, ring-like source is 30′ from the MGCLS field centred on Abell 209, and is not readily explained by simple physical models. With an assumed host galaxy at redshift 0.55, the luminosity (1025 W Hz−1) is comparable to powerful radio galaxies. The source consists of a ring of emission 175 kpc across, quadrilateral enhanced brightness regions bearing resemblance to radio jets, two ‘ears’ separated by 368 kpc, and a diffuse envelope. All of the structures appear spectrally steep, ranging from −1.0 to −1.5. The ring has high polarization (25 per cent) except on the bright patches (<10 per cent). We compare this source to the Odd Radio Circles recently discovered in ASKAP data and discuss several possible physical models, including a termination shock from starburst activity, an end-on radio galaxy, and a supermassive black hole merger event. No simple model can easily explain the observed structure of the source. This work, as well as other recent discoveries, demonstrates the power of unsupervised machine learning in mining large data sets for scientifically interesting sources.
First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way
We present the first Event Horizon Telescope (EHT) observations of Sagittarius A* (Sgr A*), the Galactic center source associated with a supermassive black hole. These observations were conducted in 2017 using a global interferometric array of eight telescopes operating at a wavelength of λ = 1.3 mm. The EHT data resolve a compact emission region with intrahour variability. A variety of imaging and modeling analyses all support an image that is dominated by a bright, thick ring with a diameter of 51.8 ± 2.3 μas (68% credible interval). The ring has modest azimuthal brightness asymmetry and a comparatively dim interior. Using a large suite of numerical simulations, we demonstrate that the EHT images of Sgr A* are consistent with the expected appearance of a Kerr black hole with mass ∼4 × 10$^{6}$ M $_{⊙}$, which is inferred to exist at this location based on previous infrared observations of individual stellar orbits, as well as maser proper-motion studies. Our model comparisons disfavor scenarios where the black hole is viewed at high inclination (i > 50°), as well as nonspinning black holes and those with retrograde accretion disks. Our results provide direct evidence for the presence of a supermassive black hole at the center of the Milky Way, and for the first time we connect the predictions from dynamical measurements of stellar orbits on scales of 10$^{3}$–10$^{5}$ gravitational radii to event-horizon-scale images and variability. Furthermore, a comparison with the EHT results for the supermassive black hole M87* shows consistency with the predictions of general relativity spanning over three orders of magnitude in central mass.
A molecular mediator for reductive concerted proton-electron transfers via electrocatalysis
Electrocatalytic approaches to the activation of unsaturated substrates via reductive concerted proton-electron transfer (CPET) must overcome competing, often kinetically dominant hydrogen evolution. We introduce the design of a molecular mediator for electrochemically triggered reductive CPET through the synthetic integration of a Brønsted acid and a redox mediator. Cathodic reduction at the cobaltocenium redox mediator substantially weakens the homolytic nitrogen-hydrogen bond strength of a Brønsted acidic anilinium tethered to one of the cyclopentadienyl rings. The electrochemically generated molecular mediator is demonstrated to transform a model substrate, acetophenone, to its corresponding neutral α-radical via a rate-determining CPET.
The Muon g-2 Experiment at Fermilab: Precision Measurements and Nonlinear Beam Dynamics
The Muon g-2 Experiment has achieved a significant milestone by measuring the positive muon anomalous magnetic moment to an unprecedented accuracy of 0.20 ppm, paving the way to a better understanding of particle physics, including beyond-the-Standard-Model possibilities. In the experiment, muons are circulated within a storage ring, and their anomalous precession frequency — the spin precession relative to momentum — is determined from decay positron time and energy data, captured using calorimeters. In support of accurate storage ring simulations, we performed high-order calculations of the field of the experiment's high-voltage quadrupoles. We also carried out accurate, transfer map–based calculations of the storage ring's chromaticity. These calculations are essential for understanding the beam dynamics effects in the Muon g-2 storage ring.
Data-driven chaos indicator for nonlinear dynamics and applications on storage ring lattice design
A data-driven chaos indicator concept is introduced to characterize the degree of chaos for nonlinear dynamical systems. The indicator is represented by the prediction accuracy of surrogate models established purely from data. It provides a metric for the predictability of nonlinear motions in a given system. When using the indicator to implement a tune-scan for a quadratic Hénon map, the main resonances and their asymmetric stop-band widths can be identified. When applied to particle transportation in a storage ring, as particle motion becomes more chaotic, its surrogate model prediction accuracy decreases correspondingly. So, the prediction accuracy, acting as a chaos indicator, can be used directly as the objective for nonlinear beam dynamics optimization. This method provides a different perspective on nonlinear beam dynamics and an efficient method for nonlinear lattice optimization. Applications in dynamic aperture optimization are demonstrated as real world examples.