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 793 records · Page 44

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

Development of a new casting method to fabricate U–Zr alloy containing minor actinides

Metal fuel slugs of U–Zr alloys for a sodium-cooled fast reactor (SFR) have conventionally been fabricated using an injection casting method. However, casting alloys containing volatile radioactive constituents, such as Am, are problematic in a conventional injection casting method. As an alternative fabrication method, low pressure gravity casting has been developed. Casting soundness, microstructural characteristics, alloying composition, density, and fuel losses were evaluated for the following as-cast fuel slugs: U–10 wt% Zr, U–10 wt% Zr–5 wt% RE, and U–10 wt% Zr–5 wt% RE–5 wt% Mn. The U and Zr contents were uniform throughout the matrix, and impurities such as oxyen, carbon, and nitrogen satisfied the specification of total impurities less than 2,000 ppm. The appearance of the fuel slugs was generally sound, and the internal integrity was shown to be satisfactory based on gamma-ray radiography. In a volatile surrogate casting test, the U–Zr–RE–Mn fuel slug showed that nearly all of the manganese was retained when casting was done under an inert atmosphere.

99 GENERAL AND MISCELLANEOUS↗

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.

Elder, Delwin L.↗

Witnessing Nonequilibrium Entanglement Dynamics in a Strongly Correlated Fermionic Chain

Many-body entanglement in condensed matter systems can be diagnosed from equilibrium response functions through the use of entanglement witnesses and operator-specific quantum bounds. Here, we investigate the applicability of this approach for detecting entangled states in quantum systems driven out of equilibrium. We use a multipartite entanglement witness, the quantum Fisher information, to study the dynamics of a paradigmatic fermion chain undergoing a time-dependent change of the Coulomb interaction. Our results show that the quantum Fisher information is able to witness distinct signatures of multipartite entanglement both near and far from equilibrium that are robust against decoherence. Here, we discuss implications of these findings for probing entanglement in light-driven quantum materials with time-resolved optical and x-ray scattering methods.

1-dimensional spin chains↗

Scalable, Infiltration-Free Ceramic Matrix Composite (SIF-CMC) manufacturing for Molten Salt Receiver: Investigation of particle-loaded preceramic resin

Due to their unique combination of high temperature, corrosion resistance, and toughness, Ceramic Matrix Composites (CMC) are the optimal materials to be used in the next generation of Concentrated Solar Power (CSP), which will operate at working temperatures up to 800 °C. However, today, CMC is too costly to be used in CSP applications due to its complex manufacturing process based on repeated infiltration-pyrolysis cycles. To solve this issue, PARC is developing a scalable, infiltration-free method (SIF) to manufacture CMC, which can reduce the production cost, increase batch consistency, and leverage all of the intrinsic advantages of current CMCs. The key to the SIF method is using a flowable preceramic resin with high loading of functional components. The formulation impregnates the fiber during curing and leads to high-density, low-shrinkage parts upon pyrolysis. In this paper, we demonstrate a proof-of-concept of the SIF process and report the relationships between key process parameters in the manufacturing process: functional material loading in the resin vs. resin’s flowability and morphology changes upon pyrolysis.

ceramic matrix composite, concentrated solar power↗

Atomic data and spectral line intensities for Fe XV

Electron impact collision strengths, energy levels, oscillator strengths, and spontaneous radiative decay rates are calculated for the n=4 and 5 complexes in Fe XV. We include in the calculations all the configurations in the n=3, 4, and 5 complexes, corresponding to 283 fine-structure levels. We also provide accurate A values for the n=3 complex. Collision strengths are calculated at five incident energies for all transitions: 12.1, 28.2, 50.4, 80.9 and 123.0 Ry above the threshold of each transition. An additional energy has been added, whose value is between 0.00075 and 0.40 Ry, depending on the transition. Calculations have been carried out using the Flexible Atomic Code and the distorted-wave approximation. Excitation rate coefficients are calculated as a function of electron temperature by assuming a Maxwellian electron velocity distribution. We combined the excitation rate coefficients and the radiative transition rates calculated in the present work with R-Matrix results available in the literature for the n=3 complex, to solve the statistical equilibrium equations for level populations at electron densities in the 10{sup 8}-10{sup 14} cm{sup -3} range and at an electron temperature of logT{sub e}(K)=6.3, corresponding to the maximum abundance of Fe XV. Spectral line intensities are calculated and compared with available observations. - Highlights: {yields} We calculate atomic data and transition rates to predict Fe XV line intensities. {yields} We provide data for the first time for many Fe XV lines observed in X-rays. {yields} These lines allow us to measure the physical parameters of solar and stellar coronae. {yields} Intensity ratios of lines emitted from n=3 and n=4 configurations do not agree with observations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

New method for analysis of X-ray computed tomography scans of TRISO fuel forms

X-ray computed tomography (XCT) is a powerful characterization tool that has been used for examination of both tristructural-isotropic (TRISO) coated particles and fuel forms. Rapid, nondestructive imaging of complete specimens by XCT provides significant characterization advantages in both speed and comprehensiveness over traditional materialographic methods. This paper discusses a method that has been developed to analyze XCT images of pebble fuel forms containing ~18,000 TRISO particles each. The XCT images are three-dimensional with contrast based on relative x-ray transparency. In the case of TRISO fuel forms, XCT image intensity distinguishes between the high-density kernel, the medium-density SiC, and the low-density matrix and carbon layers. Image analysis software has been developed to process XCT data sets to extract spatial distribution information, including particle nearest-neighbor distances, local packing fraction, and fuel-free zone thickness in TRISO fuel forms. Finally, although this type of high-resolution XCT analysis is likely not feasible as a quality control measurement for commercial TRISO fuel, it is highly useful in generating feedback for the development of fuel form fabrication processes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Communication: Towards ab initio self-energy embedding theory in quantum chemistry

The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green’s function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THEORY OF A QUODON GAS WITH APPLICATION TO PRECIPITATION KINETICS IN SOLIDS UNDER IRRADIATION

Rate theory of the radiation-induced precipitation in solids is modified with account of non-equilibrium fluctuations driven by the “gas” of lattice solitons (a.k.a. “quodons”) produced by irradiation. According to quantitative estimations, a steady-state density of the quodon gas under sufficiently intense irradiation can be comparable to the density of classical phonon gas. The modified rate theory is applied to modelling of copper precipitation in FeCu binary alloys under electron irradiation. In contrast to the classical rate theory, which disagrees strongly with experimental data on all precipitation parameters, the modified rate theory describes quite well both the evolution of precipitates and the matrix concentration of copper measured by different methods.

Dubinko, Volodymyr↗

The trends and challenges of fiber reinforced additive manufacturing

In the last few years, utilizing fiber reinforced additive manufacturing (FRAM)-based components in several industries has become quite popular. Compared to conventional AM technologies, FRAM offered complementary solutions to their needs. In general, fibers have been traditionally used in many manufacturing processes for various reasons. However, using conventional methods, there are obstacles in obtaining the desired complex geometries and low setup costs. AM offers possible avoidance of these limitations. Shape complexity, infill density, and manufacturing lead times are no longer barriers. Bridging AM with fiber reinforced materials offers a vast opportunity for lightweight and strong parts. Depending on the affinity, fibers with different structures can be mixed with different matrix materials and, thus, create stronger parts with improved mechanical properties. Process parameters like raster angle, infill speed, layer thickness, and nozzle temperature also strongly impact physical properties of FRAM products and are considered carefully. FRAM-based components are used in many industries such as aerospace, motorsports, and biomedicine, where the weight, strength, and complexity of parts are critical. Hence, numerous industrial companies and research facilities are investigating the implementation and adaptation of FRAM to their requirements. Studies are generally conducted on new materials, new FRAM technologies, the effect of fiber orientations and fraction on the performance of parts, improving the printing parameters, and other subjects. Furthermore, this study reports the current trends and challenges that FRAM is bringing to AM ecosystem.

36 MATERIALS SCIENCE↗

Nonclassical effects in a highly nonlinear generalized homogeneous Dicke model

An intensity dependent nonlinear coupling model of N two-level atoms (generalized Dicke model) interacting dispersively with a bimodal cavity field via two-photon transitions is investigated in a scenario where the rotating wave approximation is assumed. The model becomes homogeneous in the sense that the spin transition frequency is the same for all atoms and the coupling constants emerging from the collective interactions of the atomic system with the cavity field depend only on the particular radiation field mode. This allows us to represent the Dicke Hamiltonian entirely in terms of the total angular momentum J. It is assumed that, initially, the atomic system and the field are in a disentangled state where the field modes are in Glauber coherent states and the atomic system is a superposition of states |JM> (Dicke states). The model is numerically tested against simulations of normal squeezing variance of the field, squeezing factors based on the Heisenberg uncertainty principle, along with the statistical properties of the light leading to the possible production of nonclassical effects, such as degree of second-order coherence in the modes, degree of intermode correlation, as well as violation of the Cauchy-Schwartz inequality. Analytical expression of the total density operator matrix elements at t>0 shows the present nonlinear model to be strongly entangled, which is reflected in the time evolution of the linear entropy, where the superposition states are reduced to statistical mixtures. Thus, the present generalized Dicke model does not preserve the modulus of the Bloch vector. The computations, performed in the weak coupling and strong field limits, were conducted via second-order Dyson perturbative expansion of the time evolution operator matrix elements for the totality of the angular momentum states of the atomic system. - Highlights: > The present highly nonlinear generalized Dicke model involves two-photon transitions. > Quantum coherences significantly affect the dynamics of the atomic system. > At t>0 the time-dependent density operator matrix becomes strongly entangled. > Linear entropy shows that the superposition states turn into statistical mixtures. > The present Dicke model does not preserve the modulus of the Bloch vector.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical homogenization of concrete microstructures without explicit meshes

Life management of electric hydro or nuclear power plants requires to estimate long-term concrete properties on facilities, for obvious safety and serviceability reasons. Decades-old structures are foreseen to be operational for several more decades. As a large number of different concrete formulations are found in EDF facilities, empirical models based on many experiments cannot be an option for a large fleet of power plant buildings. To build predictive models, homogenization techniques offer an appealing alternative. To properly upscale creep, especially at long term, a rather precise description of the microstructure is required. However, the complexity of the morphology of concrete poses several challenges. In particular, concrete is formulated to maximize the packing density of the granular skeleton, leading to aggregates spanning several length scales with small inter particle spacings. Thus, explicit meshing of realistic concrete microstructures is either out of reach of current meshing algorithms or would produce a number of degrees of freedom far higher than the current generic FEM codes capabilities. This paper proposes a method to deal with complex matrix-inclusions microstructures such as the ones encountered at the mortar or concrete scales, without explicitly meshing them. The microstructure is superimposed to an independent mesh, which is a regular Cartesian grid. This inevitably yields so called 'gray elements', spanning across multiple phases. As the reliability of the estimate of the effective properties highly depends on the behavior affected to these gray elements, special attention is paid to them. As far as the question of the solvers is concerned, generic FEM codes are found to lack efficiency: they cannot reach high enough levels of discretization with classical free meshes, and they do not take advantage of the regular structure of the mesh. Thus, a specific finite differences/finite volumes solver has been developed. At first, generic off-the-shelf linear system solvers were used. To further improve the efficiency in terms of memory requirements, specific variants of the preconditioned conjugate gradient were implemented. This allowed to homogenize the conductivity of a concrete-like microstructure using more than 10{sup 9} degrees of freedom on a rather common hardware for 2010 (a PC embedding 48 GB of RAM). Taking benefit of the properties of the regular Cartesian grid we have also investigated a multi-level method to improve the CPU efficiency of the code.

36 MATERIALS SCIENCE↗

Copper-polydopamine composite derived from bioinspired polymer coating

Metal matrix composites with nanocarbon phases, such carbon nanotube (CNT) and graphene, have shown potentials to achieve improved mechanical, thermal, and electrical properties. However, incorporation of these nanocarbons into the metal matrix usually involves complicated processes. Here, this study explored a new processing method to fabricate copper (Cu) matrix composite by coating Cu powder particles with nanometer-thick polydopamine (PDA) thin films and sintering of the powder compacts. For sintering temperatures between 300°C and 750°C, the Cu-PDA composite samples showed higher electrical conductivity and thermal conductivity than the uncoated Cu samples, which is likely related to the higher mass densities of the composite samples. After being sintered at 950°C, the thermal conductivity of the Cu-PDA sample was approximately 12% higher than the Cu sample, while the electrical conductivity did not show significant difference. On the other hand, Knoop micro-hardness values were comparable between the Cu-PDA and Cu samples sintered at the same temperatures.

36 MATERIALS SCIENCE↗

Renormalized Singles Green’s Function for Quasi-Particle Calculations beyond the G 0 W 0 Approximation

Quasi-particle energies and band gaps in particular are critical for investigating novel materials. Commonly used density functional approximations (DFAs) systematically underestimate band gaps, and GW approximation is the established method of choice for good accuracy and reliability. However, G 0 W 0 has some undesired dependence on the DFA, while self-consistent GW (scGW) is expensive and not consistent in accuracy improvement. Here a simple and efficient G RS W 0 approach has been developed: a subspace diagonalization of the Hartree-Fock (HF) Hamiltonian with the DFA density matrix provides the new reference Green's function GRS that incorporates the effect of all single excitation contributions to the self-energy, thereby essentially eliminating the starting-point dependence. Calculations for molecules and large band gap solids demonstrate the significant improvement over G 0 W 0 and greatly reduced dependence on the initial DFA. GRS W0 approach also improve results for other bulks over G 0 W 0 but to a lesser extent, which could be due to the limitations in current implementation for bulks. The results demonstrate that to achieve good accuracy, it is not necessary to use hybrid DFA, which is expensive for bulks. Lastly, this work should be greatly significant in making GW a more robust approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design of a multi-functional gel polymer electrolyte with a 3D compact stacked polymer micro-sphere matrix for high-performance lithium metal batteries

We report lithium metal batteries (LMBs) are considered as promising high energy density batteries. However, they still suffer from poor cyclability due to the instability of solid electrolyte interphases (SEIs) and Li dendrite issues. Herein, a novel 3D densely packed multifunctional crosslinked gel polymer electrolyte (PHGPE) has been synthesized by in situ copolymerization of pentaerythritol tetraacrylate and 2-hydroxyethyl acrylate in a liquid electrolyte via a thermal initiation method. Due to intermolecular hydrogen bonding in the frameworks, the PHGPE exhibits a 3D compact-stacked and crosslinked micro-sphere structure, which benefits the construction of homogeneous Li + flux, enabling uniform Li plating/stripping and formation of robust SEIs onto the lithium electrode. Benefiting from the above advantages, symmetric Li cells in the PHGPE deliver excellent cyclabilities of 6000 h at a current density of 1 mA cm -2 (1 mA h cm -2 ) and 7700 h at 0.5 mA cm -2 (2 mA h cm -2 ). Moreover, the PHGPE has such a wide electrochemical window that it is compatible with both high-voltage LiNi 0.8 Co 0.1 Mn 0.1 O 2 and sulfur cathodes to realize long cycle life. Additionally, a LiFePO 4 based pouch size LMB also presents promising performance. Therefore, this study of the PHGPE has opened a new avenue for high-performance LMBs with excellent versatility and practicality.

25 ENERGY STORAGE↗

Predicting the properties of the lead alloys from DFT calculations

We provide qualitative results for the physical properties of the lead alloys at atomic scale by using DFT calculations. Our approach is based on the two assumptions: (i) the geometric structure of lead atoms provides a matrix where the alloying elements can take their positions in the structure as substitutions and (ii) there is a small probability of a direct interaction between the alloying elements, thus the interactions of each alloying element may be approximated by the interactions to the lead matrix. DFT calculations are used to investigate the interaction between several types of impurities and the lead matrix for low concentrations of the alloying element. We report results such as the enthalpy of formation, charge transfer and mechanical stress induced by the impurities in the lead matrix; these results can be used as qualitative guide in tuning the physico-chemical properties of the lead alloys.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lattice QCD exploration of parton pseudo-distribution functions

Here, we demonstrate a new method of extracting parton distributions from lattice calculations. The starting idea is to treat the generic equal-time matrix element ${\cal M} (Pz_3, z_3^2)$ as a function of the Ioffe time $\nu = Pz_3$ and the distance $z_3$. The next step is to divide ${\cal M} (Pz_3, z_3^2)$ by the rest-frame density ${\cal M} (0, z_3^2)$. Our lattice calculation shows a linear exponential $z_3$-dependence in the rest-frame function, expected from the $Z(z_3^2)$ factor generated by the gauge link. Still, we observe that the ratio ${\cal M} (Pz_3 , z_3^2)/{\cal M} (0, z_3^2)$ has a Gaussian-type behavior with respect to $z_3$ for 6 values of $P$ used in the calculation. This means that $Z(z_3^2)$ factor was canceled in the ratio. When plotted as a function of $\nu$ and $z_3$, the data are very close to $z_3$-independent functions. This phenomenon corresponds to factorization of the $x$- and $k_\perp$-dependence for the TMD ${\cal F} (x, k_\perp^2)$. For small $z_3 \leq 4a$, the residual $z_3$-dependence is explained by perturbative evolution, with $\alpha_s/\pi =0.1$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Microporous structure in silicone polymers with polylactic acid and polyethylene glycol microspheres: Synthesis and characterization

Porosity in polymer parts could improve mechanical properties while reducing weight and density. Individual pores must be small enough that bulk mechanical properties are retained but significant enough to reduce the density of the material. Filler materials can be incorporated into polymer matrices and removed to create voids. Here, this paper focuses on creating spherical pores in silicone matrices by inserting and removing core materials or fillers during post-cure. A method was developed to incorporate pore formers into an uncured polymer resin and remove the pore former materials from the matrix. Polyethylene glycol (PEG) and polylactic acid (PLA) microspheres were created and incorporated into silicone matrices as pore formers. Microspheres with controllable size distributions were formed via an emulsion formation. Microspheres are removed from the matrix by solvent extraction and calcination. Thermogravimetric Analysis (TGA), Dynamic Mechanical Analysis (DMA), and optical imaging characterize and determine the successful creation of micropores in the polymer matrices.

36 MATERIALS SCIENCE↗