Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffuse fraction”

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 325 records · Page 18

Understanding Viscoelasticity of an Entangled Silicone Copolymer via Coarse-Grained Molecular Dynamics Simulations

Entangled dynamics is important for understanding rheological properties of long-chain polymers. For entangled homopolymers, the classic tube-reptation model and its refinements have been successfully applied to quantify properties like diffusion coefficient and zero-rate viscosity. However, the application of such models to copolymers has been limited despite scientific and industrial importance. Here, we study the entangled melt dynamics of poly(dimethyl-co-diphenyl)siloxane random copolymer for a range of mean-composition-ratio ϕ of the diphenyl component via long-term molecular dynamics simulation with a recently developed coarse-grained model. We found that the segmental relaxation is heterogeneous at the monomeric level because of compositional fluctuations. However, at the chain-entanglement level and higher length scales, the viscoelastic response is homogeneous with compositional dependence only through the overall diphenyl fraction ϕ. The relaxation modulus of the entangled copolymer melt conforms to the Likhtman–McLeish model, and the viscosity predicted using our current coarse-grained parameters is in good quantitative agreement with experimental data.

Copolymers↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

On the high-temperature stability of the Al 8 Cu 3 Ce intermetallic in an additively manufactured Al-Cu-Ce-Zr alloy

High-temperature resistant eutectic Al alloys are crucial materials for lightweight and energy efficient design in the automotive and aviation industries. Additive manufacturing offers a pathway to refine eutectic microstructures and develop novel alloys with superior high-temperature strength. High-volume fraction intermetallic Al-Cu-Ce alloys have been developed to deliver high-temperature strength in combination with reduced hot-tearing susceptibility. Zr is added to provide additional strengthening via nanoscale Al 3 Zr precipitation, and to stabilize and avoid coarsening of the Al 8 Cu 3 Ce phase. However, the detailed interaction between Zr and Al 8 Cu 3 Ce remains unexplored. In this work, we show with synchrotron X-ray diffraction that laser powder bed fusion fabricated Al-Cu-Ce and Al-Cu-Ce-Zr alloys contain predominantly the Al 8 Cu 3 Ce intermetallic in the as-fabricated condition. Heat treatment of the Al-Cu-Ce alloy results in the Al 8 Cu 3 Ce → Al 8 Cu 4 Ce phase transformation. In the Al-Cu-Ce-Zr alloy, minor fractions of (Al,Cu,Si) 4 Ce and Al 2 Cu-θ are found in the as-fabricated condition, while Al 8 Cu 3 Ce remains stable during heat treatment. Atom probe microscopy quantifies intermetallic stoichiometries and reveals how Zr is enriched at the Al-matrix/Al 8 Cu 3 Ce interface acting as a diffusion barrier against solute exchange. Calibrated thermodynamic modeling underpins this as a kinetic effect. Here, a qualitative microstructural model summarizes, how Zr stabilizes Al 8 Cu 3 Ce against phase transformations and coarsening.

36 MATERIALS SCIENCE↗

Mechanism and kinetic study of rare earth extraction from allanite by direct acid leaching

Direct acid leaching was conducted to recover rare earth elements (REEs) from an allanite ore assaying 9,861.4 ppm of REEs. The effects of operational variables such as acid type, sulfuric acid (H 2 SO 4 ) concentration, temperature, solid/liquid (S/L) ratio, and particle size on the total REE (TREE) recovery were systematically examined. The REE mineralogy and the involved leaching mechanisms were investigated by microscopic characterization and residue re-leaching tests. It was found that approximately 80% of REEs was extracted using 1 M H 2 SO 4 at 75 °C for 2 h. Increasing temperature significantly improved the REE extraction, while changing acid type, H 2 SO 4 concentration, and particle size did not. The TREE recovery also gradually decreased at higher S/L ratios. The low fraction of Si leached (<4.0%) suggested that the extraction did not significantly alter the crystal structure of the dominant mineral. The leaching kinetics of REEs were fast within the first 10 min, followed by a much slower stage (i.e., 10-120 min). The kinetic modeling results showed that each stage may be controlled by both chemical reaction and diffusion (i.e., mix control). The activation energies for REE leaching using 1 M H 2 SO 4 in the 0-10 min and 10-120 min were determined to be 20.3 and 10.8 kJ/mol, respectively. Finally, the residue re-leaching results suggested that the fast leaching kinetics at the beginning may be explained by the preferential decomposition of metamict allanite rather than the formation of silica gel, while the following slower leaching stage was due to the dissolution of difficult-to-leach REE-bearing minerals such as well-crystallized allanite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Constraints on the Emission of Gamma-Rays from M31 with HAWC

Cosmic rays, along with stellar radiation and magnetic fields, are known to make up a significant fraction of the energy density of galaxies such as the Milky Way. When cosmic rays interact in the interstellar medium, they produce gamma-ray emission which provides an important indication of how the cosmic rays propagate. Gamma rays from the Andromeda Galaxy (M31), located 785 kpc away, provide a unique opportunity to study cosmic-ray acceleration and diffusion in a galaxy with a structure and evolution very similar to the Milky Way. Using 33 months of data from the High Altitude Water Cherenkov Observatory, we search for TeV gamma rays from the galactic plane of M31. We also investigate past and present evidence of galactic activity in M31 by searching for Fermi Bubble-like structures above and below the galactic nucleus. No significant gamma-ray emission is observed, so we use the null result to compute upper limits on the energy density of cosmic rays $>10$ TeV in M31. In conclusion, the computed upper limits are approximately ten times higher than expected from the extrapolation of the Fermi LAT results.

79 ASTRONOMY AND ASTROPHYSICS↗

Microstructure and mechanical properties of hypoeutectic Al–6Ce–3Ni-0.7Fe (wt.%) alloy

Here, a hypoeutectic, Fe-modified Al–Ce–Ni alloy (Al–6Ce–3Ni-0.7Fe, wt.%) is studied in terms of microstructure, thermal stability, ambient temperature strengthening, and creep resistance. The as-cast microstructure consists of primary Al dendrites and interdendritic binary eutectic regions (Al–Al11Ce3 and/or Al–Al9(Fe,Ni)2), with micron/submicron lamellar spacing, depending on the location along the height of the ingot. The cast alloy exhibits excellent coarsening resistance at 400 °C, with mostly unchanged microstructure and microhardness after 6 weeks of aging, indicating good thermal stability of Al11Ce3 and Al9(Fe,Ni)2. Orowan strengthening and load transfer are identified as strengthening mechanisms at ambient and elevated temperature. A high volume fraction of the intermetallic phases (providing load transfer) and relatively coarse eutectic spacing (for modest Orowan strengthening) result in a moderate as-cast microhardness of 566 ± 32 MPa. Creep resistance at 300 and 350 °C is similar to a binary Al-12.5Ce eutectic alloy (with twice the Ce content) because of two countering effects: Al–6Ce–3Ni-0.7Fe shows a higher volume fraction of strengthening intermetallic phases, but it also exhibits a large fraction of primary Al dendrites which weaken the alloy. By contrast, the alloy, when laser-remelted at the surface, has a fully eutectic microstructure without primary aluminum dendrites achieved by high undercooling on solidification, with a refined network of eutectic phases that doubles the microhardness as compared to the cast alloy. Whereas coarsening is faster due to the shorter diffusion distances between the eutectic phases, hardness remains ~30% higher than the as-cast alloy after ~6 weeks aging at 400 °C.

36 MATERIALS SCIENCE↗

Identification of combustion mode under MILD conditions using Chemical Explosive Mode Analysis

Direct Numerical Simulations (DNS) data of Moderate or Intense Low-oxygen Dilution (MILD) combustion are analysed in this paper to identify the contributions of the autoignition and flame modes. This is performed using an extended Chemical Explosive Mode Analysis (CEMA) which accounts for diffusion effects allowing it to discriminate between deflagration and autoignition. This analysis indicates that in premixed MILD combustion conditions, the main combustion mode is ignition for all dilution and turbulence levels and for the two reactant temperature conditions considered. In non-premixed conditions, the preponderance of the ignition mode was observed to depend on the axial location and mixture fraction stratification. With a large mixture fraction lengthscale, ignition is more preponderant in the early part of the domain while the deflagrative mode increases further downstream. On the other hand, when the mixture fraction lengthscale is small, sequential autoignition is observed. Finally, the various combustion modes are observed to correlate strongly with mixture fraction where lean mixtures are more likely to autoignite while stoichiometric and rich mixtures are more likely to react as deflagrative structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalar Power Spectra and Scalar Structure Function Evolution in the Richtmyer–Meshkov Instability Upon Reshock

The Richtmyer–Meshkov instability of a twice-shocked gas interface is studied using high-speed planar laser-induced fluorescence in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. Furthermore, this initial condition is accelerated by a shock of nominal strength M = 1.9, and then accelerated again by the transmitted shock that reflects off the end wall of the tube. Three individual experiments are analyzed, the energy spectrum and the structure functions of the light gas mole fraction field are calculated and compared.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robustness of the Galactic Center excess morphology against masking

The Galactic Center excess (GCE) remains an enduring mystery, with leading explanations being annihilating dark matter or an unresolved population of millisecond pulsars. Analyzing the morphology of the GCE provides critical clues to identify its exact origin. We investigate the robustness of the inferred GCE morphology against the effects of masking, an important step in the analysis where the gamma-ray emission from point sources and the galactic disk are excluded. Using different masks constructed from Fermi point source catalogs and a wavelet method, we find that the GCE morphology, particularly its ellipticity and cuspiness, is relatively independent of the choice of mask for energies above 2–3 GeV. The GCE morphology systematically favors an approximately spherical shape, as expected for dark matter annihilation. Compared to various stellar bulge profiles, a spherical dark matter annihilation profile better fits the data across different masks and galactic diffuse emission backgrounds, except for the stellar bulge profile which provides a similar fit to the data. Modeling the GCE with two components, one from dark matter annihilation and one tracing the Coleman bulge, we find this two-component model outperforms any single component or combinations of dark matter annihilation and other stellar bulge profiles. Uncertainty remains about the exact fraction contributed by each component across different background models and masks. Furthermore, when the Coleman bulge dominates, its corresponding spectrum lacks characteristics typically associated with millisecond pulsars, suggesting that it mostly models the emission from other sources instead of the GCE that is still present and spherically symmetric.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

(Towards) DNS of a Laboratory Lean CH4/H2 Low-Swirl Flame Impinging on an Inclined Wall

Due to downsizing trends, flame-wall interaction (FWI) is increasingly prominent in gas turbines (GTs). FWI has direct consequences on flame stabilization and pollutant emissions, but it is not well understood in turbulent flows representative of GTs. We present results from a direct numerical simulation (DNS) of a turbulent CH4/H2 model GT low-swirl laboratory-scale flame interacting with an inclined wall. The results from the laboratory flame include simultaneous measurements of velocity using stereo particle imaging velocimetry and OHxCH2O planar laser induced fluorescence. The adaptive-mesh refinement solver PeleLMeX is used, with 24-species, 105-reaction reduced Aramco chemical kinetics mechanism. The premixed fuel-air mixture consists of hydrogen-enriched methane with 70% hydrogen volume fraction and 0.4 equivalence ratio. The inflow is prescribed to match experimental measurements at the burner exit. Karlovitz and turbulent Reynolds numbers are 300 and 400, respectively. The simulation and experimental results show excellent agreement. The flame features a bowl-shape stabilization, with a corrugated, continuous flame front at the leading edge, followed by fragmented reaction zones downstream. A large diffuse cloud of CH2O is formed downstream of the quenching point. The simulation results indicate that the cloud of CH2O is the result of incomplete methane combustion, with CH2O "leaking" from the locally quenched reaction zones.The DNS provides fine-grain resolution of turbulence-flame-wall interaction that cannot be captured with experimental measurements. With access to the entire solution vector at each cell of the computational domain, the local quenching.

flame-wall interactions↗

Optical Investigation of a Partial Fuel Stratification Strategy to Stabilize Overall Lean Operation of a DISI Engine Fueled with Gasoline and E30

This paper offers new insights into a partial fuel stratification (PFS) combustion strategy that has proven to be effective at stabilizing overall lean combustion in direct injection spark ignition engines. To this aim, high spatial and temporal resolution optical diagnostics were applied in an optically accessible engine working in PFS mode for two fuels and two different durations of pilot injection at the time of spark: 210 µs and 330 µs for E30 (gasoline blended with ethanol by 30% volume fraction) and gasoline, respectively. In both conditions, early injections during the intake stroke were used to generate a well-mixed lean background. The results were compared to rich, stoichiometric and lean well-mixed combustion with different spark timings. In the PFS combustion process, it was possible to detect a non-spherical and highly wrinkled blue flame, coupled with yellow diffusive flames due to the combustion of rich zones near the spark plug. The initial flame spread for both PFS cases was faster compared to any of the well-mixed cases (lean, stoichiometric and rich), suggesting that the flame propagation for PFS is enhanced by both enrichment and enhanced local turbulence caused by the pilot injection. Different spray evolutions for the two pilot injection durations were found to strongly influence the flame kernel inception and propagation. PFS with pilot durations of 210 µs and 330 µs showed some differences in terms of shapes of the flame front and in terms of extension of diffusive flames. Yet, both cases were highly repeatable.

33 ADVANCED PROPULSION SYSTEMS↗

Spatially Distributed Ramp Reversal Memory in VO 2

Ramp‐reversal memory has recently been discovered in several insulator‐to‐metal transition materials where a non‐volatile resistance change can be set by repeatedly driving the material partway through the transition. This study uses optical microscopy to track the location and internal structure of accumulated memory as a thin film of VO 2 is temperature cycled through multiple training subloops. These measurements reveal that the gain of insulator phase fraction between consecutive subloops occurs primarily through front propagation at the insulator‐metal boundaries. By analyzing transition temperature maps, it is found, surprisingly, that the memory is also stored deep inside both insulating and metallic clusters throughout the entire sample, making the metal‐insulator coexistence landscape more rugged. This non‐volatile memory is reset after heating the sample to higher temperatures, as expected. Diffusion of point defects is proposed to account for the observed memory writing and subsequent erasing over the entire sample surface. By spatially mapping the location and character of non‐volatile memory encoding in VO 2 , this study results enable the targeting of specific local regions in the film where the full insulator‐to‐metal resistivity change can be harnessed in order to maximize the working range of memory elements for conventional and neuromorphic computing applications.

36 MATERIALS SCIENCE↗

Role of ring-enlargement reactions in the formation of aromatic hydrocarbons

Ring-enlargement reactions can provide a fast route towards the formation of six-membered single-ring or polycyclic aromatic hydrocarbons (PAHs). To investigate the participation of the cyclopentadienyl (C 5 H 5 ) radical in ring-enlargement reactions in high-temperature environments, a mass-spectrometric study was conducted. Experimental access to the C 5 H 5 high-temperature chemistry was provided by two counterflow diffusion flames. Cyclopentene was chosen as a primary fuel given the large amount of resonantly stabilized cyclopentadienyl radicals produced by its decomposition and its high tendency to form PAHs. In a second experiment, methane was added to the fuel stream to promote methyl addition pathways and to assess the importance of ring-enlargement reactions for PAH growth. The experimental dataset includes mole fraction profiles of small intermediate hydrocarbons and of several larger species featuring up to four condensed aromatic rings. Results show that, while the addition of methane enhances the production of methylcyclopentadiene and benzene, the concentration of larger polycyclic hydrocarbons is reduced. The increase of benzene is probably attributable to the interaction between the methyl and the cyclopentadienyl radicals. However, the formation of larger aromatics seems to be dominated only by the cyclopentadienyl driven molecular-growth routes which are hampered by the addition of methane. Furthermore, in addition to the experimental work, two chemical mechanisms were tested and newly calculated reaction rates for cyclopentadiene reactions were included. In an attempt to assess the impact of cyclopentadienyl ring-enlargement chemistry on the mechanisms’ predictivity, pathways to form benzene, toluene, and ethylbenzene were investigated. Results show that the updated mechanism provides an improved agreement between the computed and measured aromatics concentrations. Nevertheless, a detailed study of the single reaction steps leading to toluene, styrene, and ethylbenzene would be certainly beneficial.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cosmogenic 3 He exposure dating in mafic rocks by “Virtual mineral separation” of pyroxene

We describe a “virtual mineral separation” method for measuring the cosmogenic 3 He concentration in pyroxene in mafic rocks that consist mainly of plagioclase and pyroxene, without physically separating the minerals. This approach is significantly faster and more cost-effective than the conventional method, which requires physical separation and purification of pyroxene grains by time-consuming and labor-intensive crushing, acid cleaning, magnetic separation, HF etching, and handpicking under a microscope. The premise of the method is that helium diffusivity is much higher in plagioclase than in pyroxene, so controlled preheating of a mixed whole-rock sample can degas 3 He from plagioclase while retaining all 3 He in pyroxene. A second heating step releases all 3 He from pyroxene for measurement. To then obtain a 3 He concentration in pyroxene rather than the whole rock, we determine the pyroxene weight fraction in the sample using X-ray computed tomography (CT). A comparison of 3 He concentrations in pyroxene measured using virtual mineral separation with those measured in the same samples by physical mineral separation in previous work shows no evidence of systematic bias between the methods. Virtual mineral separation greatly simplifies the workflow for 3 He exposure-dating of mafic rocks, reduces time, effort, and cost, and permits measurements on very small samples. This enables new emerging applications of exposure dating, such as quantifying stochastic surface processes, ecosystem studies, and potential subglacial bedrock exposure dating.

Geosciences↗

A simple, validated approach for design of two-dimensional periodic particle patterns via acoustophoresis

Two-dimensional patterning of microparticles enables a wide range of functional materials, including patterned energy storage electrodes, flexible electronics, and sensor arrays. Particle patterning via acoustics offers an attractive path to generate a wide variety of 2D periodic patterns that introduce tailorable hierarchical porosity, useful for controlling surface area, transport distances, and other properties. This method is most effective with micron scale particles and patterns of tens to hundreds of microns. To enable systematic exploration of the broad design space for such patterns, this work develops a model of 2D and 3D assembly of particles at high loadings and validates the obtained patterns against both experiments and more computationally intensive modeling techniques. Using this simple model, connections are mapped between input parameters (like actuation conditions, particle volume fraction, material properties) and output geometrical features (like void size and shape, pattern connectivity, and surface area) so that they can be tailored to given applications. The utility of this simple model is illustrated by predicting and then experimentally demonstrating new hierarchical patterns resulting from multiple waves of different frequencies interacting. These multiscale patterns offer the potential to lift the limits on surface area, diffusion distances, and other features.

42 ENGINEERING↗

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗

Enhancement and maximum in the isobaric specific-heat capacity measurements of deeply supercooled water using ultrafast calorimetry

Knowledge of the temperature dependence of the isobaric specific heat (C p ) upon deep supercooling can give insights regarding the anomalous properties of water. If a maximum in C p exists at a specific temperature, as in the isothermal compressibility, it would further validate the liquid–liquid critical point model that can explain the anomalous increase in thermodynamic response functions. The challenge is that the relevant temperature range falls in the region where ice crystallization becomes rapid, which has previously excluded experiments. Here, we have utilized a methodology of ultrafast calorimetry by determining the temperature jump from femtosecond X-ray pulses after heating with an infrared laser pulse and with a sufficiently long time delay between the pulses to allow measurements at constant pressure. Evaporative cooling of ~15-µm diameter droplets in vacuum enabled us to reach a temperature down to ~228 K with a small fraction of the droplets remaining unfrozen. We observed a sharp increase in C p , from 88 J/mol/K at 244 K to about 218 J/mol/K at 229 K where a maximum is seen. The C p maximum is at a similar temperature as the maxima of the isothermal compressibility and correlation length. From the C p measurement, we estimated the excess entropy and self-diffusion coefficient of water and these properties decrease rapidly below 235 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Control-Oriented Modeling of Lithium-Ion Batteries

Battery management systems (BMSs), which monitor and optimize performance while ensuring safety, require control-oriented models, i.e., models tailored to the design and implementation of estimation and control algorithms. Physics-based electrochemical models describe detailed battery phenomena, but are too computationally intensive for use in estimation and control applications. Single particle models (SPMs), which retain some of the physics of electrochemical models, are often used for control-oriented battery modeling since they are computationally efficient; however, they are only valid over very low frequency ranges and C-rates. Empirical equivalent circuit models (ECMs) are also used for control-oriented battery modeling since they are computationally efficient and can describe battery behavior over wide frequency ranges; however, they provide no physical understanding of the battery and, therefore, have limited applicability. Further, fractional order terms (e.g., Warburg impedances) are often employed, making the models unwieldy for use in the time domain. This work provides a control-oriented battery model that combines the benefits of SPM and ECM models, while overcoming their limitations. The proposed model incorporates some of the battery physics found in electrochemical models, can easily be used in both the time and frequency domains, and describes battery behavior over its entire frequency range. A linearized single particle model, which incorporates key electrochemical parameters, is used for modeling battery physics at very low frequencies. For low frequencies, integer-order linear systems are used to approximate diffusion physics described by Warburg impedances, and high frequency behavior is modeled by the double layer capacitance effect. The proposed battery model is more computationally efficient than full electrochemical models since it does not require the solution of PDEs, is accurate for a wider frequency range than the SPM considered in this paper, and does not suffer from the unwieldiness and limited applicability of empirical ECMs. Finally, the model is validated in the time and frequency domains via a comparison to pseudo-two-dimensional (P2D) model simulations and experimental data.

25 ENERGY STORAGE↗