Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Identifying Bounds of Inorganic Content in Solventless Processing of Hybrid Solid Electrolytes

Solid-state lithium batteries require safe, robust electrolytes to enable higher energy densities and improved safety over conventional cells. Hybrid polymer–ceramic electrolytes are a promising solution, combining the processability of polymers with the high ionic conductivity and mechanical strength of inorganic fillers. In this work, we demonstrate a solventless, UV-curing method to produce hybrid solid electrolytes using a poly(ethylene glycol) dimethyl ether (PEGDME)-based photocurable matrix incorporating Li 1.5 Al 0.5 Ge 1.5 (PO 4 ) 3 (LAGP) or Li 7 La 3 Zr 2 O 12 (LLZO) ceramic electrolyte. Inorganic filler loadings up to ∼55 wt.% could be successfully incorporated via this process which was the highest inorganic content at which the slurry remains processable and cured into a uniform film. The resulting UV-cured composite electrolytes remain flexible and exhibit room-temperature ionic conductivities on the order of 10 −4 S·cm −1 , along with notably improved lithium-ion transference numbers compared to conventional polymer electrolytes. Similar performance and processing limits were observed for both LAGP and LLZO, indicating that ceramic filler chemistry does not significantly affect the UV-curing process or the electrolyte's ion transport properties in this regime. Eliminating solvents from fabrication not only simplifies processing and mitigates environmental concerns but also enables higher solid contents that enhance mechanical strength and help suppress lithium dendrite formation. In conclusion, this scalable approach thus paves the way for manufacturing robust composite solid electrolytes for next-generation solid-state batteries (SSBs).

Batteries↗
DOE OSTI2024

Separation of Lipoproteins for Quantitative Analysis of 14 C-Labeled Lipid-Soluble Compounds by Accelerator Mass Spectrometry

To date, 14 C tracer studies using accelerator mass spectrometry (AMS) have not yet resolved lipid-soluble analytes into individual lipoprotein density subclasses. The objective of this work was to develop a reliable method for lipoprotein separation and quantitative recovery for biokinetic modeling purposes. The novel method developed provides the means for use of small volumes (10–200 µL) of frozen plasma as a starting material for continuous isopycnic lipoprotein separation within a carbon- and pH-stable analyte matrix, which, following post-separation fraction clean up, created samples suitable for highly accurate 14 C/ 12 C isotope ratio determinations by AMS. Manual aspiration achieved 99.2 ± 0.41% recovery of [5- 14 CH 3 ]-(2R, 4'R, 8'R)-α-tocopherol contained within 25 µL plasma recovered in triacylglycerol rich lipoproteins (TRL = Chylomicrons + VLDL), LDL, HDL, and infranatant (INF) from each of 10 different sampling times for one male and one female subject, n = 20 total samples. Small sample volumes of previously frozen plasma and high analyte recoveries make this an attractive method for AMS studies using newer, smaller footprint AMS equipment to develop genuine tracer analyses of lipophilic nutrients or compounds in all human age ranges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Water vapor profile retrieval simulation studies for the HIRS/MSU and AMTS/MSU sounders

A combined physical statistical method for the retrieval of atmospheric water vapor column densities has been developed at the Goddard Laboratory for Atmospheric Sciences. In this method differences between observed and calculated clear column brightness temperatures are used in a regression formalism to obtain estimates of integrated column densities in given atmospheric layers. The physics of the radiative transfer problem is contained in the calculated clear column brightness temperatures and in accounting for the effects of clouds on the observed radiances. The regression matrix is equivalent to a set of empirical sensitivity factors.

Reuter, D.↗
NASA NTRS2005

E-Beam Processing of Polymer Matrix Composites for Multifunctional Radiation Shielding

Aliphatic polymers were identified as optimum radiation shielding polymeric materials for building multifunctional structural elements for in-space habitats. Conceptual damage tolerant configurations of polyolefins have been proposed, but many manufacturing issues relied on methods and materials which have sub-optimal radiation shielding characteristics (for example, epoxy matrix and adhesives). In the present approach, we shall investigate e-beam processing technologies for inclusion of high-strength aliphatic polymer reinforcement structures into a highly cross-linked polyolefin matrix. This paper reports the baseline thermo-mechanical properties of low density polyethylene and highly crystallized polyethylene.

Hou, Tan-Hung↗
NASA NTRS2001

Retrieval of Polar Stratospheric Cloud Microphysical Properties From Lidar Measurements: Dependence on Particle Shape Assumptions

Knowledge of particle sizes and number densities of polar stratospheric clouds (PSCs) is highly important, because they are critical parameters for the modeling of the ozone chemistry of the stratosphere. In situ measurements of PSC particles are rare. the main instrument for the accumulation of PSC data are lidar systems. Therefore the derivation of some microphysical properties of PSCS from the optical parameters measured by lidars would be highly beneficial for ozone research. Inversion of lidar data obtained in the presence of PSCs formed from crystalline particles type 11 and the various nitric acid tri Ydrrate (NAT) types cannot be easily accomplished, because a suitable scattering theory for small faceted crystals has not been readily available tip to now. As a consequence, the T-matrix method is commonly used for the interpretation of these PSC lidar data. Here the assumption is made that the optical properties of an ensemble of spheroids resemble those of crystalline PSCs, and microphysical properties of the PSC are inferred from the optical signatures of the PSC at two or more wavelengths. The problem with the T-matrix approach is that the assumption of spheroidal instead of faceted particles can lead to dramatically wrong results: Usually cloud particle properties are deduced from analysis of lidar profiles of backscatter ratio and depolarization ratio. The particle contribution to the backscatter ratio is given by the product of the particle number density and the backscattering cross section. The latter is proportional to the value of the particle's scattering phase function at 180 degrees scattering angle. At 180 degrees however, the phase functions of rough, faceted crystals and of spheroids with same maximum dimension differ by a factor of 6. From this it follows that for a PSC consisting of faceted crystals, the particle number density is underestimated by roughly the same factor if spheroidal particles are unrealistically assumed. We are currently developing a retrieval technique for determining the microphysical parameters of crystalline PSCs that takes into account the faceted shape of the PSC particles. This approach utilizes finite-difference time-domain (FDTD) calculations of particle optical properties. The accuracy and the free choice of the shape of the scattering particle make the FDTD technique a promising tool for the inversion of PSC lidar data. A first comparison of FDTD and T-matrix calculations will be presented.

Reichardt, Susanne↗
DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A Monte Carlo Model of Irradiation-Induced Recrystallization in Polycrystalline UMo Fuels

Experiments show that irradiation-induced recrystallization speeds up the swelling kinetics in U10wt%Mo fuels. However, recrystallization mechanisms and the effect of initial grain microstructures on recrystallization kinetics are still unclear. In this work a Monte Carlo model coupling the rate theory of defect evolution has been developed to study the irradiation-induced recrystallization. The rate theory is used to describe the spatial evolution of gas bubbles, interstitials and interstitial loops; First-Passage Kinetic Monte Carlo (FPKMC) approach is used to describe the fast and strongly anisotropic migration of interstitials, and a Cellular Automata method is used to model the formation of recrystallized grains. With the assumption that recrystallization may occur when the local interstitial loop density is larger than a given critical value, simulation results reveal that 1) recrystallized grains first nucleate on grain boundaries and the recrystallization zone front moves to the center of original coarse grains in the UMo matrix, and 2) recrystallization starts earlier in coarse polycrystalline structures, while the overall recrystallization kinetics decreases with increasing grain size. These results agree with experimental observations. The comparison of recrystallization kinetics obtained from experiments and modeling suggests that the interstitial loop accumulation leads to the recrystallization and the interstitial loop growth is suppressed inside coarse grains due to the over-pressured intra-granular gas bubbles.

Frazier, William E.↗
NASA NTRS1976

Radiance, polarization, and ellipticity of the radiation in the earth's atmosphere

The complete radiation field including polarization is calculated for a model of the real atmosphere by the matrix operator method. The radiance, direction and amount of polarization, and ellipticity are obtained at the top and bottom of the atmosphere for three values of the surface albedo (0; 0.15 0.90) and five solar zenith angles. Scattering and absorption by molecules (including ozone) and by aerosols are taken into account together with the variation of the number density of these substances with height. All results are calculated for both a normal aerosol number and a distribution which is one-third of the normal amount at all heights. The calculated values show general qualitative agreement with the available experimental measurements. The position of the neutral points of the polarization in the principal plane is a sensitive indicator of the characteristics of the aerosol particles in the atmosphere, since it depends on the sign and value of the single scattered polarization for scattering angles around 20 deg and 160 deg for transmitted and reflected photons respectively.

Hitzfelder, S. J.↗
NASA NTRS1989

Microgravity processing of particulate reinforced metal matrix composites

The elimination of such gravity-related effects as buoyancy-driven sedimentation can yield more homogeneous microstructures in composite materials whose individual constituents have widely differing densities. A comparison of composite samples consisting of particulate ceramics in a nickel aluminide matrix solidified under gravity levels ranging from 0.01 to 1.8 G indicates that the G force normal to the growth direction plays a fundamental role in determining the distribution of the reinforcement in the matrix. Composites with extremely uniform microstructures can be produced by these methods.

Morel, Donald E.↗
NASA NTRS2005

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗
DOE OSTI2015

Phonon-phonon interactions: First principles theory

We present the details of a method to perform molecular-dynamics (MD) simulations without thermostat and with very small temperature fluctuations ±ΔT starting with MD step 1. It involves preparing the supercell at the time t = 0 in physically correct microstates using the eigenvectors of the dynamical matrix. Each initial microstate corresponds to a different distribution of kinetic and potential energies for each vibrational mode (the total energy of each microstate is the same). Averaging the MD runs over many initial microstates further reduces ΔT. The electronic states are obtained using first-principles theory (density-functional theory in periodic supercells). Three applications are discussed: the lifetime and decay of vibrational excitations, the isotope dependence of thermal conductivities, and the flow of heat at an interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Spectroscopy of states in 136 Ba using the 138 Ba(p,t) reaction

Background: The 136 Ba isotope is the daughter nucleus in 136 Xe $ββ$ decay. It also lies in a shape transitional region of the nuclear chart, making it a suitable candidate to test a variety of nuclear models. Purpose: To obtain spectroscopic information on states in 136 Ba, which will allow a better understanding of its low-lying structure. These data may prove useful to constrain future 136 Xe → 136 Ba neutrinoless $ββ$ decay matrix element calculations. Methods: Here, a 138 Ba(p, t) reaction was used to populate states in 136 Ba up to approximately 4.6 MeV in excitation energy. The tritons were detected using a high-resolution Q3D magnetic spectrograph. A distorted wave Born approximation analysis was performed for the measured triton angular distributions. Results: 102 excited states in 136 Ba were observed, out of which 52 are reported for the first time. Definite spin-parity assignments are made for 26 newly observed states, while previously ambiguous assignments for ten other states are resolved. Together with other available data, the results are used to determine level densities in 136 Ba. These were compared with theory predictions, obtained using shell model calculations with Hamiltonians previously used for 136 Xe neutrinoless $ββ$ decay matrix element evaluations. Conclusions: The shell model predicted level densities agree reasonably well for the two Hamiltonians. However the results for theory and experiment are found to agree only at lower energies, diverging from one another for the higher lying states, with the discrepancy increasing with energy. This is presumably because of lower production cross sections for a majority of the higher-lying predicted states and the experimental limitations in resolving a large number of nearly degenerate states predicted by the theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗
DOE OSTI2013

Optical and structural investigation of GSO gels made via sol gel method

Silica and doped silica materials that are produced by a solution gelation (sol gel) process create versatile materials that have a large range of applications. Coatings, thin films and fibers, nanomaterials, energy conservation materials and optical materials can be synthesized with this technique. Specifically, silica gels have many advantages for optical materials. The gels are transparent and possess a high homogeneity and density. In addition, silica materials can be doped with different compounds to enhance scintillating properties. Gadolinium Oxyorthosilicate (GSO) is a promising candidate as a host for rare earth ions for scintillating applications because of its high stability, density and fast decay time. When a sol gel is dried conventionally, it often shrinks or cracks, forming a xerogel. This work will first focus on forming a transparent GSO host xerogel with sol gel chemistry. Once a transparent GSO xerogel is formulated and characterized, Ce 3+ will be doped into the GSO network. Non-radiative energy transfer from the Gd matrix ions to the Ce 3+ ions occurs because the GSO crystals activated are close in energy to the Ce 3+ f → d transition. GSO: Ce 3+ crystals are currently produced by growing single crystals via Bridgeman and Czochralski methods, which are expensive. Optical ceramics offers a cheap alternative method for making scintillating materials. Structural and thermal characterization of the xerogel will be investigated using X-ray diffraction (XRD), Fourier-transformed infrared spectroscopy (FTIR), differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA). The materials’ optical properties will be examined by UV-Vis absorption spectroscopy. XRD will be utilized to determine the crystal structure. This will tell determine if the gel is amorphous or crystalline. FTIR was be used to identify functional groups in the gel. DSC and TGA will be used to thermally characterize the gel’s crystallization temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2019

The Mechanical and Corrosion Properties of Hot-Rolled TiB{sub 2}/Ni Composites

To quickly prepare TiB{sub 2}/Ni composites at low temperatures, a method combining spark plasma sintering and hot rolling was used. This paper presents the preparation, microstructural and mechanical properties, strengthening mechanism and corrosion properties of TiB{sub 2}/Ni composites in detail. X-ray diffraction analysis shows that the texture degree of the TiB{sub 2}/Ni composites after hot rolling greatly differed from that of the as-rolled pure Ni samples. The addition of TiB{sub 2} particles limited grain growth in the Ni matrix and led to grain refinement. Dislocation strengthening, grain refinement strengthening and load transfer strengthening contributed to a 30% increase in the tensile strength compared to that of the as-rolled pure Ni. In addition, the strengthening mechanism of the mechanical property reinforcement was analyzed. The electrochemical test results demonstrate that the as-rolled pure Ni exhibited a lower corrosion current density and corrosion resistance than the other as-rolled TiB{sub 2}/Ni composites fabricated in this paper. This work demonstrates a new approach for fabricating composites in less time at lower preparation temperatures than conventional methods.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constructing a Simulation Surrogate with Partially Observed Output

Gaussian process surrogates are a popular alternative to directly using computationally expensive simulation models. When the simulation output consists of many responses, dimension-reduction techniques are often employed to construct these surrogates. However, surrogate methods with dimension reduction generally rely on complete output training data. This article proposes a new Gaussian process surrogate method that permits the use of partially observed output while remaining computationally efficient. The new method involves the imputation of missing values and the adjustment of the covariance matrix used for Gaussian process inference. The resulting surrogate represents the available responses, disregards the missing responses, and provides meaningful uncertainty quantification. In conclusion, the proposed approach is shown to offer sharper inference than alternatives in a simulation study and a case study where an energy density functional model that frequently returns incomplete output is calibrated.

42 ENGINEERING↗
DOE OSTI2018

Multi-scale theory-assisted nano-engineering of plasmonic-organic hybrid electro-optic device performance

Multi-scale (correlated quantum and statistical mechanics) modeling methods have been advanced and employed to guide the improvement of organic electro-optic (OEO) materials, including by analyzing electric field poling induced electro-optic activity in nanoscopic plasmonic-organic hybrid (POH) waveguide devices. The analysis of in-device electro-optic activity emphasizes the importance of considering both the details of intermolecular interactions within organic electro-optic materials and interactions at interfaces between OEO materials and device architectures. Dramatic improvement in electro-optic device performance--including voltage-length performance, bandwidth, energy efficiency, and lower optical losses have been realized. These improvements are critical to applications in telecommunications, computing, sensor technology, and metrology. Multi-scale modeling methods illustrate the complexity of improving the electro-optic activity of organic materials, including the necessity of considering the trade-off between improving poling-induced acentric order through chromophore modification and the reduction of chromophore number density associated with such modification. Computational simulations also emphasize the importance of developing chromophore modifications that serve multiple purposes including matrix hardening for enhanced thermal and photochemical stability, control of matrix dimensionality, influence on material viscoelasticity, improvement of chromophore molecular hyperpolarizability, control of material dielectric permittivity and index of refraction properties, and control of material conductance. Consideration of new device architectures is critical to the implementation of chipscale integration of electronics and photonics and achieving the high bandwidths for applications such as next generation (e.g., 5G) telecommunications.

Engineering↗
DOE OSTI2020

Impact of raw material surface oxide removal on dual band infrared optical properties of As 2 Se 3 chalcogenide glass

The manufacturing of low loss chalcogenide glasses (ChGs) for optoelectronic applications is ultimately defined by the concentration of impurities present in starting materials or imparted via processing. We describe a rapid method for purifying metallic starting materials in As 2 Se 3 glass where oxide reduction is correlated to optical and physical properties. Specifically, As-O reduction enhances the glass’ dual-band optical transparency proportional to the extent (13-fold reduction) of oxide reduction, and is accompanied by a change in density and hardness associated with changes in matrix bonding. A significant modification of the glass’ index and LWIR Abbe number is reported highlighting the significant impact purification has on material dispersion control required in optical designs.

36 MATERIALS SCIENCE↗