Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum compression”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulating a pulsed-power-driven plasma with ideal MHD

We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a recipe for coupling pulsed power generators to the magnetohydrodynamics (MHD) domain by approximating the perfectly insulating vacuum as a light, perfectly conducting, inviscid MHD fluid and discuss the applicability of this counter-intuitive technique. This much more affordable ideal MHD representation is particularly useful in situations where a plasma exhibits interesting three-dimensional phenomena, either due to the design of the experiment or due to developing instabilities. We verified that this coupling recipe works by modeling an exactly solvable flux compression generator as well as a self-similar Noh-like solution and demonstrated convergence to the theoretical solution. We also showed examples of simulating complex three-dimensional pulsed power devices with this technique. We release our code implementation to the public.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamical tuning of the chemical potential to achieve a target particle number in grand canonical Monte Carlo simulations

We present a method to facilitate Monte Carlo simulations in the grand canonical ensemble given a target mean particle number. The method imposes a fictitious dynamics on the chemical potential, to be run concurrently with the Monte Carlo sampling of the physical system. Corrections to the chemical potential are made according to time-averaged estimates of the mean and variance of the particle number, with the latter being proportional to thermodynamic compressibility. Here, we perform a variety of tests, and in all cases find rapid convergence of the chemical potential—inexactness of the tuning algorithm contributes only a minor part of the total measurement error for realistic simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Equations of State in Computational Physics: What, Why, and How

Summary: 1. Equations of state furnish thermodynamic properties of a material. 2. Usually involve P, p (or V), and T, or P, r, and E 3. Necessary ingredient to solve compressible flow problems. 4. Varieties: i. Analytical expressions (ideal gas, Grüneisen) – fast, convenient for verification and test problems ii. Graphs and charts – good for hand calculations iii. Tables and databases – generally most accurate, require interpolation i. Steam tables ii. NIST database iii. SESAME EOS Library 5. Constructed by using semi-empirical models from Statistical Mechanics and fitting to available data

36 MATERIALS SCIENCE↗
DOE OSTI2020

Automated Fitting of a Semi-empirical Multiphase Equation of State for Carbon

The equation of state (EOS) of carbon is important in high explosive, geophysical, and inertial confinement fusion applications. Within the semi-empirical Sesame framework, the EOS of each phase is represented by a sum of cold, vibrational, and thermal electronic Helmholtz free energy contributions. Each phase has ~5-10 independent parameters that are adjusted to reproduce single-phase data (e.g., thermal expansion, isothermal compression) as well as experimentally- and computationally- derived phase boundaries. Manual calibration of the full multiphase EOS is arduous. We present our progress in development of automated EOS parameterization based on minimization of an objective function. Here, this function encodes deviation of model EOS results from experimental/computational benchmarks. Optimization is implemented as a combination of global (particle swarm) and local optimization techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Chapman–Jouguet deflagration criteria and compressibility dynamics of turbulent fast flames for turbulence-induced deflagration-to-detonation transition

This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Understanding how chemical structure affects ignition-delay-time $\phi$-sensitivity

$\phi$-sensitivity is the change in ignition delay time (IDT) with respect to the fuel-to-air equivalence ratio ($\phi$). High $\phi$-sensitivity is a desirable fuel property for applications in advanced compression ignition and multi-mode engine designs. Understanding how $\phi$-sensitivity depends on chemical structure is essential for selecting promising biofuels from the ever-growing list of proposed candidates. Here, we investigate the effect of chemical structure on $\phi$-sensitivity with experiment, simulation, and theory. Experimental Advanced Fuel Ignition Delay Analyzer (AFIDA) measurements for 2,4-dimethylpentane and diisopropyl ether provide evidence that branching and functional groups strongly impact $\phi$-sensitivity. Further insights into this dependence are obtained with 0-D kinetic simulations with existing mechanisms for n-pentane, diethyl ether, 3-pentanone, n-heptane, 2-methylhexane, 2,4-dimethylpentane, and 2,2,3-trimethylbutane. Quantum mechanical (QM) G4 calculations of low-temperature reactions help explain the observed experimental and simulation trends. Specifically, these QM calculations provide theoretical estimates of the ketohydroperoxide (KHP) dissociation rates, the HO 2 formation rates from peroxy radical (ROO), and the “cross-over” temperatures, i.e., the temperature at which ROO dissociation is favored compared to hydroperoxyl radical (QOOH) formation. Each of these reaction rates is compared to the n-alkane reference point to determine the impact of branching and different functional groups. Although kinetic mechanisms typically assume that KHP dissociation rates are invariant of chemical environment, our QM results suggest that this rate can span a range of roughly two orders of magnitude. We also discuss the importance of including the peroxy-hydroperoxy (OO-OOH) hydrogen transfer reaction for branched ethers. Finally, the insights gained assist in proposing a highly $\phi$-sensitive compound, namely, isopropyl propyl ether.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Structure and pressure dependence of the Fermi surface of lithium

Here, we report studies of the Fermi surface (FS) of isotopically pure polycrystalline 7 Li from ambient pressure to 4.7 GPa. Shubnikov–de Haas (SdH) oscillations at 300 mK in external magnetic fields up to 35 T are measured to map the spherical parts of the FS of lithium. Our ambient pressure data show that the principal SdH frequencies consist of three distinct peaks at 41.25 ± 0.05, 41.65 ± 0.05, and 42.05 ± 0.05 kT, which are theoretically consistent with the presence of two crystal structure domains exhibiting nearly spherical FS. The size of the spherical parts of the FS is compatible with bcc and fcc crystal structures. The measured frequencies at 41.25 and 42.65 kT present direct quantitative evidence for the spherical deformation of the FS in fcc Li. Our high-pressure data show that while the FS of Li deforms under compression, it remains mostly spherical up to 4.7 GPa and the pressure dependence of the SdH frequency is consistent with the theoretically calculated pressure dependence in the fcc structure. Finally, we find that the electron effective mass does not deviate under pressure significantly from its ambient pressure value.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

High-performance cementitious composites containing nanostructured carbon additives made from charred coal fines

Carbon-based nanomaterials, such as carbon nanoplatelets, graphene oxide, and carbon quantum dots, have many possible end-use applications due to their ability to impart unique mechanical, electrical, thermal, and optical properties to cement composites. Despite this potential, these materials are rarely used in the construction industry due to high material costs and limited data on performance and durability. In this study, domestic coal is used to fabricate low-cost carbon nanomaterials that can be used economically in cement formulations. A range of chemical and physical processing approaches are employed to control the size, morphology, and chemical functionalization of the carbon nanomaterial, which improves its miscibility with cement formulations and its impact on mechanical properties and durability. At loadings of 0.01 to 0.07 wt.% of coal-derived carbon nanomaterial, the compressive and flexural strength of cement samples are enhanced by 24% and 23%, respectively, in comparison to neat cement. At loadings of 0.02 to 0.06 wt.%, the compressive and flexural strength of concrete composites increases by 28% and 21%, respectively, in comparison to neat samples. Additionally, the carbon nanomaterial additives studied in this work reduce cement porosity by 36%, permeability by 86%, and chloride penetration depth by 60%. These results illustrate that low-loadings of coal-derived carbon nanomaterial additives can improve the mechanical properties, durability, and corrosion resistance of cement composites.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of the surface termination and the oxygen vacancy position on LaNiO 3 ultra-thin films: First-principles study

In this study, while ultra-thin layers of the LaNiO 3 film exhibit a remarkable metal-insulator transition as the film thickness becomes smaller than a few unit cell (u.c.), the formation of possible oxygen vacancies and their effects on the correlated electronic structure have been rarely studied using first principles. Here, we investigate the effects of the surface termination and the oxygen vacancy position on the electronic properties and vacancy energetics of LaNiO 3 ultra-thin films under the compressive strain using density functional theory plus U (DFT + U). We find that oxygen vacancies can be easily formed in the Ni layers with the NiO 2 terminated surface (0.5 u.c. and 1.5 u.c. thickness) compared to the structures with the LaO terminated surface and the in-plane vacancy is energetically favored than the out-of-plane vacancy. When two vacancy sites are allowed, the Ni square plane geometry is energetically more stable in most cases as two oxygen vacancies tend to stay near a Ni ion. Strong anisotropy between the in-plane and out-of-plane vacancy formation as well as the layer and orbital dependent electronic structure occur due to strain, surface termination, charge reconstruction, and quantum confinement effects. The in-plane vacancy of the NiO 2 terminated structure is favored since the released charge due to the oxygen vacancy can be easily accommodated in the d x 2 -y 2 orbital, which is less occupied than the d z 2 orbital. Remarkably, the oxygen vacancy structure containing the Ni square-plane geometry becomes an insulating state in DFT + U with a sizable band gap of 1.2 eV because the large crystal field splitting between d z 2 and d x 2 -y 2 orbitals in the square-plane favors an insulating state and the Mott insulating state is induced in other Ni sites due to strong electronic correlations. In the thin-film structure without oxygen vacancies, the strong correlation effect in DFT + U drives a pseudogap ground state at the Fermi energy, similarly as the experimental photo-emission spectra; however, the variation of the DFT + U electronic structure depending on the surface termination becomes weaker compared to those obtained in DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Advanced Ab Initio Methods for Nuclear Structure (Final Report)

Over the past decade, there has been enormous progress in the description of nuclear structure from first principles, using interactions from Chiral Effective Field Theory that are rooted in Quantum Chromodynamics, the fundamental theory of strong interactions, and many-body methods that solve the Schrödinger equation with systematically improvable approximations. Amongst those are the family of In-Medium Similarity Renormalization Group (IMSRG) framework developed by the PI and his co-workers. While these methods scale polynomially in the size N of the single-particle basis, computational efforts still grows dramatically as we increase the degrees of freedom for the nucleons by relaxing symmetries or introducing continuum couplings (see below), or as we push to improved truncations to provide precise inputs for experimental efforts, in particular in fundamental symmetry searches. In order to address this growing computational cost, one focus area of this award was the exploration of compression and factorization methods. The key to success or failure is the presence of low-rank structures within the matrix elements of NN and 3N interactions or the IMSRG evolution operator, and a means to reformulate the method that will let us exploit them.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Electronic structure of Co 3 d states in the Kitaev material candidate honeycomb cobaltate Na 3 Co 2 SbO 6 probed with x-ray dichroism

The recent prediction that honeycomb lattices of Co 2+ (3d 7 ) ions could host dominant Kitaev interactions provides an exciting direction for exploration of new routes to stabilizing Kitaev’s quantum spin liquid in real materials. Na 3 Co 2 SbO 6 has been singled out as a potential material candidate provided that spin and orbital moments couple into a J eff = $\frac{1}{2}$ ground state, and that the relative strength of trigonal crystal field and spin-orbit coupling acting on Co ions can be tailored. Using x-ray linear dichroism (XLD) and x-ray magnetic circular dichroism (XMCD) experiments, alongside configuration interaction calculations, we confirm the counterintuitive positive sign of the trigonal crystal field acting on Co 2+ ions and test the validity of the J eff = $\frac{1}{2}$ description of the electronic ground state. Furthermore, the results lend experimental support to recent theoretical predictions that a compression (elongation) of CoO 6 octahedra along (perpendicular to) the trigonal axis would drive this cobaltate toward the Kitaev limit, assuming the J eff = $\frac{1}{2}$ character of the electronic ground state is preserved.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Structure-Related Electronic and Magnetic Properties in Ultrathin Epitaxial Ni x Fe 3–x O 4 Films on MgO(001)

Off-stoichiometric Ni x Fe 3–x O 4 ultrathin films (x < 2.1) with varying Ni content x and thickness 16 (±2) nm were grown on MgO(001) by reactive molecular beam epitaxy. Synchrotron-based high-resolution X-ray diffraction measurements reveal vertical compressive strain for all films, resulting from a lateral pseudomorphic adaption of the film to the substrate lattice without any strain relaxation. Complete crystallinity with smooth interfaces and surfaces is obtained independent of the Ni content x. For x < 1 an expected successive conversion from Fe 3 O 4 to NiFe 2 O 4 is observed, whereas local transformation into NiO structures is observed for films with Ni content x > 1. However, angle-resolved hard X-ray photoelectron spectroscopy measurements indicate homogeneous cationic distributions without strictly separated phases independent of the Ni content, while X-ray absorption spectroscopy shows that also for x > 1, not all Fe 2+ cations are substituted by Ni 2+ cations. The ferrimagnetic behavior, as observed by superconducting quantum interference device magnetometry, is characterized by decreasing saturation magnetization due to the formation of antiferromagnetic NiO parts.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multi-millijoule hollow-core fiber compression of short-wave infrared pulses to a single cycle

Approaches for efficient pulse compression can enable dramatic increases in the available peak power, as well as enable the generation of isolated attosecond X-ray pulses. Achieving high compression ratios for longer wavelength drivers has, however, been challenging. We present the compression of few-cycle 2.1 µm central wavelength short-wave infrared laser pulses to 6.9 fs with 2.35 mJ pulse energy at a 10 kHz repetition rate. Electric field resolved measurements reveal a single cycle light field oscillation. With a carrier-envelope-phase stability of 131 mrad and average power fluctuations below 1 %, the system constitutes an excellent light source for strong-field experiments and attosecond physics.

Blöchl, Johannes [Ludwig Maximilian Univ. of Munic↗
DOE OSTI2020

Impact of Laboratory-Accelerated Aging Methods to Study Alkali–Silica Reaction and Reinforcement Corrosion on the Properties of Concrete

This study focuses on two separate investigations of the main aging mechanisms: alkali–silica reactivity (ASR) and the corrosion of reinforcing steel (rebar) concrete, both of which may result in a premature failure to meet the serviceability or strength requirements of a concrete structure. However, these processes occur very slowly, spanning decades. The impact of direct chemical additives to fresh concrete to accelerate ASR and the corrosion of reinforcing steel on the fresh and hardened properties of the ensuing material are investigated to inform the potential use of chemicals in large-scale studies. The deterioration of reinforced concrete (RC) is determined by means of expansion, cracking, bulk diffusivity and surface resistivity measurements, and compressive, split tensile and flexural strength tests. The results indicate that the addition of sodium hydroxide and calcium chloride can effectively accelerate the crack formation and propagation in concrete due to ASR and the corrosion of rebar, respectively. The ASR-induced cracks maintained a constant crack width from 0.05 mm to 0.1 mm over the measurement period regardless of the intensity of aging acceleration. Adding 4% chloride by weight of cement for accelerating rebar corrosion resulted in an average crack that was 82% larger than in the case of ASR accelerated with the addition of sodium hydroxide. The addition of alkali resulted in an increase in early-age (7-day) strength. At a total alkali loading of 2.98 kg/m3, 3.84 kg/m3 and 5.57 kg/m3, the 28-day compressive strength of concrete decreased by 3%, 10% and 24%, respectively. Similarly, a higher early-age strength and a lower later-age strength was observed for the concrete in the presence of corrosive calcium chloride. The results from this research are expected to inform future studies on the long-term performance of RC structures under accelerated ASR and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗