Publications

Coupling the mimetic finite difference method with the tensor-train format results in a very effective method for low-rank numerical approximations of the solutions of the time-dependent Maxwell wave propagation equations in three dimensions. To this end, we discretize the curl operators on the primal/dual tensor product grid complex and we couple the space discretization with a staggered-in-time second-order accurate time-marching scheme. The resulting solver is accurate to the second order in time and space, and still compatible, so that the approximation of the magnetic flux field has zero discrete divergence with a discrepancy close to the machine precision level. Our approach is not limited to the second-order of accuracy. We can devise higher-order formulations in space through suitable extensions of the tensor-train stencil to compute the derivatives of the mimetic differential operators. Employing the tensor-train format improves the solver performance by orders of magnitude in terms of CPU time and memory storage. A final set of numerical experiments confirms this expectation.

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO₂, predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO2 preferentially into their interlayers and therefore potentially act as a material for CO2 separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO2 into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO2 intercalation. Our free energy calculations indicate that CO2 favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO2/H2O ratio in clay interlayers as compared to that in bulk water. CO2 molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO2 intercalation. The structural charge distribution significantly affects the amount of CO2 intercalated in the interlayers. These results provide a mechanistic understanding of CO2 intercalation in clays for CO2 separation, capture, and storage.

Alkyl-substituted oxetanes are cyclic ethers formed via unimolecular reactions of QOOH radicals produced via a six-membered transition state in the preceding isomerization step of organic peroxy radicals, ROO. Owing to radical isomer-specific formation pathways, cyclic ethers are unambiguous proxies for inferring QOOH reaction rates. Therefore, accounting for subsequent oxidation of cyclic ethers is important in order to accurately determine rates for QOOH → products. Cyclic ethers can react via unimolecular reaction (ring-opening) or via bimolecular reaction with O2 to form cyclic ether-peroxy adducts. The computations herein provide reaction mechanisms and theoretical rate coefficients for the former type in order to determine competing pathways for the cyclic ether radicals. Rate coefficients of unimolecular reactions of 2,4-dimethyloxetanyl radicals were computed using master equation modeling from 0.01 to 100 atm and from 300 to 1000 K. Coupled-cluster methods were utilized for stationary-point energy calculations, and uncertainties in the computed rate coefficients were accounted for using variation in barrier heights and in well depths. The potential energy surfaces reveal accessible channels to several species via crossover reactions, such as 2-methyltetrahydrofuran-5-yl and pentanonyl isomers. For the range of temperature over which 2,4-dimethyloxetane forms during n-pentane oxidation, the following are the major channels: 2,4-dimethyloxetan-1-yl → acetaldehyde + allyl, 2,4-dimethyloxetan-2-yl → propene + acetyl, and 2,4-dimethyloxetan-3-yl → 3-butenal + methyl, or, 1-penten-3-yl-4-ol. Well-skipping reactions were significant in a number of channels and also exhibited a markedly different pressure dependence. The calculations show that rate coefficients for ring-opening are approximately an order of magnitude lower for the tertiary 2,4-dimethyloxetanyl radicals than for the primary and secondary 2,4-dimethyloxetanyl radicals. Unlike for reactions of the corresponding ROO radicals, however, unimolecular rate coefficients are independent of the stereochemistry. Moreover, rate coefficients of cyclic ether radical ring-opening are of the same order of magnitude as O2 addition, underscoring the point that a competing network of reactions is necessary to include for accurate chemical kinetics modeling of species profiles for cyclic ethers.

Here we present a classical molecular-spin dynamics (MSD) methodology that enables accurate computations of the temperature dependence of the magnetocrystalline anisotropy as well as magnetoelastic properties of magnetic materials. The nonmagnetic interactions are accounted for by a spectral neighbor analysis potential (SNAP) machine-learned interatomic potential, whereas the magnetoelastic contributions are accounted for using a combination of an extended Heisenberg Hamiltonian and a Néel pair interaction model, representing both the exchange interaction and spin-orbit-coupling effects, respectively. All magnetoelastic potential components are parameterized using a combination of first-principles and experimental data. Our framework is applied to the α phase of iron. Initial testing of our MSD model is done using a 0 K parametrization of the Néel interaction model. After this, we examine how individual Néel parameters impact the $B$₁ and $B$₂ magnetostrictive coefficients using a moment-independent δ sensitivity analysis. The results from this study are then used to initialize a genetic algorithm optimization which explores the Néel parameter phase space and tries to minimize the error in the B₁ and B₂ magnetostrictive coefficients in the range of 0–1200 K. Our results show that while both the 0 K and genetic algorithm optimized parametrization provide good experimental agreement for $B$₁ and $B$₂, only the genetic algorithm optimized results can capture the second peak in the $B$₁ magnetostrictive coefficient which occurs near approximately 800 K.

An experimental investigation and the optical modeling of the structural coloration produced from total internal reflection interference within 3D microstructures are described. Ray-tracing simulations coupled with color visualization and spectral analysis techniques are used to model, examine, and rationalize the iridescence generated for a range of microgeometries, including hemicylinders and truncated hemispheres, under varying illumination conditions. An approach to deconstruct the observed iridescence and complex far-field spectral features into its elementary components and systematically link them to ray trajectories that emanate from the illuminated microstructures is demonstrated. The results are compared with experiments, wherein microstructures are fabricated with methods such as chemical etching, multiphoton lithography, and grayscale lithography. Microstructure arrays patterned on surfaces with varying orientation and size lead to unique color-traveling optical effects and highlight opportunities for how total internal reflection interference can be used to create customizable reflective iridescence. The findings herein provide a robust conceptual framework for rationalizing this multibounce interference mechanism and establish approaches for characterizing and tailoring the optical and iridescent properties of microstructured surfaces.

Tungsten (W) is a material of choice for the divertor material due to its high melting temperature, thermal conductivity, and sputtering threshold. However, W has a very high brittle-to-ductile transition temperature, and at fusion reactor temperatures (≥1000 K), it may undergo recrystallization and grain growth. Dispersion-strengthening W with zirconium carbide (ZrC) can improve ductility and limit grain growth, but much of the effects of the dispersoids on microstructural evolution and thermomechanical properties at high temperatures are still unknown. We present a machine learned Spectral Neighbor Analysis Potential for W-ZrC that can now be used to study these materials. In order to construct a potential suitable for large-scale atomistic simulations at fusion reactor temperatures, it is necessary to train on ab initio data generated for a diverse set of structures, chemical environments, and temperatures. Further accuracy and stability tests of the potential were achieved using objective functions for both material properties and high temperature stability. Validation of lattice parameters, surface energies, bulk moduli, and thermal expansion is confirmed on the optimized potential. Tensile tests of W/ZrC bicrystals show that although the W(110)-ZrC(111) C-terminated bicrystal has the highest ultimate tensile strength (UTS) at room temperature, observed strength decreases with increasing temperature. At 2500 K, the terminating C layer diffuses into the W, resulting in a weaker W-Zr interface. Meanwhile, the W(110)-ZrC(111) Zr-terminated bicrystal has the highest UTS at 2500 K.

Here in this work, we investigate the applicability of the harmonic balance method (HBM) to predict periodic solutions of a single degree-of-freedom forced Duffing oscillator with freeplay nonlinearity. By studying the route to impact, which refers to a parametric study as the contact stiffness increases from soft to hard, the convergence behavior of the HBM can be understood in terms of the strength of the non-smooth forcing term. HBM results are compared to time-integration results to facilitate an evaluation of the accuracy of nonlinear periodic responses. An additional contribution of this study is to perform convergence and stability analysis specifically for isolas generated by the non-smooth nonlinearity. Residual error analysis is used to determine the approximate number of harmonics required to get results accurate to a given error tolerance. Hill’s method and Floquet theory are employed to compute the stability of periodic solutions and identify the types of bifurcations in the system.

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. Successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. The solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

The Z machine is a current driver producing up to 30 MA in 100 ns that utilizes a wide range of diagnostics to assess accelerator performance and target behavior conduct experiments that use the Z target as a source of radiation or high pressures. Here, we review the existing suite of diagnostic systems, including their locations and primary configurations. The diagnostics are grouped in the following categories: pulsed power diagnostics, x-ray power and energy, x-ray spectroscopy, x-ray imaging (including backlighting, power flow, and velocimetry), and nuclear detectors (including neutron activation). We will also briefly summarize the primary imaging detectors we use at Z: image plates, x-ray and visible film, microchannel plates, and the ultrafast x-ray imager. The Z shot produces a harsh environment that interferes with diagnostic operation and data retrieval. We term these detrimental processes “threats” of which only partial quantifications and precise sources are known. Finally, we summarize the threats and describe techniques utilized in many of the systems to reduce noise and backgrounds.

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

BaCe0.25Mn0.75O3−δ (BCM), a non-stoichiometric oxide with a layered perovskite-like crystal structure, has recently emerged as a prospective contender for application in renewable energy harvesting by solar thermochemical hydrogen generation. Using solar-thermal energy and a reducing environment, oxygen vacancies can be created in high-temperature BCM, and the reduced crystal so obtained can, in turn, produce H2 by stripping oxygen from H2O. Therefore, a first step toward understanding the working mechanism and optimizing the performance of BCM is a thorough and comparative analysis of the electronic structure of the pristine and the reduced material. In this paper, we probe the electronic structure of BCM using the combined effort of first-principles calculations and experimental O K-edge X-ray absorption spectroscopy (XAS). The computed projected density of states (PDOS) and orbital plots are used to propose a simplified model for orbital mixing between the oxygen and metal atoms. With the help of state-of-the-art simulations, we are able to find the origins of the XAS peaks and categorize them on the basis of contribution from Ce and Mn. For the reduced crystal, the calculations show that the change in electron density resulting from the reduction is strongly localized around the oxygen vacancy. Experimental measurements reveal a marked lowering of the first O K-edge peak in the reduced crystal. Using theoretical analysis, this is shown to result from lifting of spin degeneracy in the absorption peaks as well as from a diminished O 2p contribution to the frontier unoccupied orbitals, in accordance with the tight binding scheme. The simulated results serve as a reference for the extent of spectral change as a function of the percentage of oxygen vacancies in the reduced crystal. Our study paves the way for the investigation of the working mechanism of BCM and for computational and experimental efforts aimed at design and discovery of efficient water-splitting oxides.

The present contribution emphasizes the formation of oligomeric products in various depolymerization approaches of lignin, namely reductive catalytic fractionation, oxidative catalytic fractionation, and pyrolysis. Three possible routes to form such oligomers in these depolymerization processes are summarized and compared from various studies conducted on model compounds. Next, the main identification techniques for characterizing oligomeric products are highlighted. Particular focus is given to 2D-HSQC-NMR, GPC, Maldi-TOF-MS and FT-ICR-MS, which represent the state-of-art characterization of lignin. Special attention was paid to the transferability of these techniques for depolymerized oligomeric lignin. Finally, both the existing and expected potential lignin valorization routes are discussed for these oligomers, and technical hurdles and recommendations are provided in an attempt to catalyze the development of new discoveries and enabling technologies.

Herein, the International Commission on Illumination (CIE) designed its color space to be perceptually uniform so that a given numerical change in the color code corresponds to perceived change in color. This color encoding is demonstrated to be advantageous in scientific visualization and analysis of vector fields. The specific application is analysis of ice motion in the Arctic where patterns in smooth monthly-averaged ice motion are seen. Furthermore, fractures occurring in the ice cover result in discontinuities in the ice motion. This vector jump in displacement can also be visualized. We then analyze modeled and observed fractures through the use of a metric on the color space, and image amplitude and phase metrics. Amplitude and phase metrics arise from image registration that is accomplished by sampling images using space filling curves, thus reducing the image registration problem to the more reliable functional alignment problem. We demonstrate this through an exploration of the metrics to compare model runs to an observed ice crack.

Here we used aerosol mass spectrometry coupled with tunable synchrotron photoionization to measure radical and closed-shell species associated with particle formation in premixed flames and during pyrolysis of butane, ethylene, and methane. We analyzed photoionization (PI) spectra for the C₇H₇ radical to identify the isomers present during particle formation. For the combustion and pyrolysis of all three fuels, the PI spectra can be fit reasonably well with contributions from four radical isomers: benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl. Although there are significant experimental uncertainties in the isomeric speciation of C7H₇, the results clearly demonstrate that the isomeric composition of C₇H₇ strongly depends on the combustion or pyrolysis conditions and the fuel or precursors. Fits to the PI spectra using reference curves for these isomers suggest that all of these isomers may contribute to m/z 91 in butane and methane flames, but only benzyl and vinylcyclopentadienyl contribute to the C₇H₇ isomer signal in the ethylene flame. Only tropyl and benzyl appear to play a role during pyrolytic particle formation from ethylene, and only tropyl, vinylcyclopentadienyl, and o-tolyl appear to participate during particle formation from butane pyrolysis. There also seems to be a contribution from an isomer with an ionization energy below 7.5 eV for the flames but not for the pyrolysis conditions. Kinetic models with updated and new reactions and rate coefficients for the C₇H₇ reaction network predict benzyl, tropyl, vinylcyclopentadienyl, and o-tolyl to be the primary C₇H₇ isomers and predict negligible contributions from other C₇H₇ isomers. These updated models provide better agreement with the measurements than the original versions of the models but, nonetheless, underpredict the relative concentrations of tropyl, vinylcyclopentadienyl, and o-tolyl in both flames and pyrolysis and overpredict benzyl in pyrolysis. Our results suggest that there are additional important formation pathways for the vinylcyclopentadienyl, tropyl, and o-tolyl radicals and/or loss pathways for the benzyl radical that are currently unaccounted for in the present models.

Tritium population thermodynamics and transport kinetics critically define the tritium storage performance of zirconium tritides that can be used for a variety of nuclear applications including tritium-producing burnable absorber rods. Both thermodynamic and kinetic properties can be sensitive to grain sizes of materials and can be significantly altered by irradiated defects during operation under the reactor environments. A thorough experimental characterization of how these properties evolve under different reactor conditions and different initial grain structures is extremely challenging. Here molecular dynamics simulations are used to investigate tritium population and diffusion in zirconium with and without different planar symmetric and asymmetric tilt grain boundaries and irradiated defects. Here, we found that in addition to trapping tritium, the most significant effect of planar grain boundaries is to increase tritium diffusivity on the boundary plane. Furthermore, fine grain structures are found to mitigate the change of tritium diffusivity due to irradiated point defects as these point defects are likely to migrate to and sink at grain boundaries.

Here, the effectiveness of continuous vibro-impact forcing representations for the cantilevered pipe that conveys fluid is explored and analyzed. The previously accepted forcing model utilizing a smoothened trilinear spring is estimated using three continuous forcing representations, namely, polynomial, rational polynomial, and hyperbolic tangent. The accuracy of the estimated forcing functions is investigated and analyzed by calculating the root mean square error, and bifurcation diagrams are generated and compared to the nominal system. Additionally, the dynamic response of the system is further characterized using Poincare maps, power spectra, and basins of attraction. Once all continuous forcing representations are analyzed and compared to the nominal system, the computational cost of each method is examined, and further limitations of the hyperbolic tangent method are discovered. It is proved that the hyperbolic tangent forcing representation most accurately captures the dynamic response of the pipeline, and the least accurate representation is the rational polynomial representation. Additionally, considerable computational cost is saved when employing the hyperbolic tangent representation compared to the discontinuous representation.

Abstract not provided.

Mixtures of gas-phase hydrogen isotopologues (diatomic combinations of protium, deuterium, and tritium) can be separated using columns containing a solid such as palladium that reversibly absorbs hydrogen. A temperature-swing process can transport hydrogen into or out of a column by inducing temperature-dependent absorption or desorption reactions. We consider two designs: a thermal cycling absorption process, which moves hydrogen back and forth between two columns, and a simulated moving bed (SMB), where columns are in a circular arrangement. We present a numerical mass and heat transport model of absorption columns for hydrogen isotope separation. It includes a detailed treatment of the absorption-desorption reaction for palladium. By comparing the isotope concentrations within the columns as a function of position and time, we observe that SMB can lead to sharper separations for a given number of thermal cycles by avoiding the remixing of isotopes.

Vertical-axis wind turbines (VAWTs) have a long history, with a wide variety of turbine archetypes that have been designed and tested since the 1970s. While few utility-scale VAWTs currently exist, the placement of the generator near the turbine base could make VAWTs advantageous over tradition horizontal-axis wind turbines for floating offshore wind applications via reduced platform costs and improved scaling potential. However, there are currently few numerical design and analysis tools available for VAWTs. One existing engineering toolset for aero-hydro-servo-elastic simulation of VAWTs is the Offshore Wind ENergy Simulator (OWENS), but its current modeling capability for floating systems is non-standard and not ideal. This article describes how OWENS has been coupled to several OpenFAST modules to update and improve modeling of floating offshore VAWTs and discusses the verification of these new capabilities and features. The results of the coupled OWENS verification test agree well with a parallel OpenFAST simulation, validating the new modeling and simulation capabilities in OWENS for floating VAWT applications. These developments will enable the design and optimization of floating offshore VAWTs in the future.

Abstract not provided.

Titanium alloys are used in a large array of applications. In this work we focus our attention on the most used alloy, Ti-6Al-4V (Ti64), which has excellent mechanical and biocompatibility properties with applications in aerospace, defense, biomedical, and other fields. Here we present high-fidelity experimental shock compression data measured on Sandia's Z machine. We extend the principal shock Hugoniot for Ti64 to more than threefold compression, up to over 1.2 TPa. We use the data to validate our ab initio molecular dynamics simulations and to develop a highly reliable, multiphase equation of state (EOS) for Ti64, spanning a broad range of temperature and pressures. The first-principles simulations show very good agreement with Z data and with previous three-stage gas gun data from Sandia's STAR facility. The resulting principal Hugoniot and the broad-range EOS and phase diagram up to 10 TPa and 105 K are suitable for use in shock experiments and in hydrodynamic simulations. The high-precision experimental results and high-fidelity simulations demonstrate that the Hugoniot of the Ti64 alloy is stiffer than that of pure Ti and reveal that Ti64 melts on the Hugoniot at a significantly lower pressure and temperature than previously modeled.

Fireballs produced from the detonation of high explosives often contain particulates primarily composed of various phases of carbon soot. The transport and concentration of these particulates is of interest for model validation and emission characterization. This work proposes ultra-high-speed imaging techniques to observe a fireball's structure and optical depth. An extinction-based diagnostic applied at two wavelengths indicates that extinction scales inversely with wavelength, consistent with particles in the Rayleigh limit and dimensionless extinction coefficients which are independent of wavelength. Within current confidence bounds, the extinction-derived soot mass concentrations agree with expectations based upon literature reported soot yields. Results also identify areas of high uncertainty where additional work is recommended.