Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Models”

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 991 records · Page 55

Giga-z: A 100,000 OBJECT SUPERCONDUCTING SPECTROPHOTOMETER FOR LSST FOLLOW-UP

We simulate the performance of a new type of instrument, a Superconducting Multi-Object Spectrograph (SuperMOS), that uses microwave kinetic inductance detectors (MKIDs). MKIDs, a new detector technology, feature good quantum efficiency in the UVOIR, can count individual photons with microsecond timing accuracy, and, like X-ray calorimeters, determine their energy to several percent. The performance of Giga-z, a SuperMOS designed for wide field imaging follow-up observations, is evaluated using simulated observations of the COSMOS mock catalog with an array of 100,000 R{sub 423{sub nm}} = E/{Delta}E = 30 MKID pixels. We compare our results against a simultaneous simulation of LSST observations. In 3 yr on a dedicated 4 m class telescope, Giga-z could observe Almost-Equal-To 2 billion galaxies, yielding a low-resolution spectral energy distribution spanning 350-1350 nm for each; 1000 times the number measured with any currently proposed LSST spectroscopic follow-up, at a fraction of the cost and time. Giga-z would provide redshifts for galaxies up to z Almost-Equal-To 6 with magnitudes m{sub i} {approx}< 25, with accuracy {sigma}{sub {Delta}z/(1+z)} Almost-Equal-To 0.03 for the whole sample, and {sigma}{sub {Delta}z/(1+z)} Almost-Equal-To 0.007 for a select subset. We also find catastrophic failure rates and biases that are consistently lower than for LSST. The added constraint on dark energy parameters for WL + CMB by Giga-z using the FoMSWG default model is equivalent to multiplying the LSST Fisher matrix by a factor of {alpha} = 1.27 (w{sub p} ), 1.53 (w{sub a} ), or 1.98 ({Delta}{gamma}). This is equivalent to multiplying both the LSST coverage area and the training sets by {alpha} and reducing all systematics by a factor of 1/{radical}({alpha}), advantages that are robust to even more extreme models of intrinsic alignment.

79 ASTRONOMY AND ASTROPHYSICS↗

Bound states in the B -matrix formalism for the three-body scattering

In this work, we consider a model of relativistic three-body scattering with a bound state in the two-body subchannel. We show that the naiïve K-matrix type parametrization, here referred to as the B-matrix, has nonphysical singularities near the physical region. We show how to eliminate such singularities by using dispersion relations and also show how to reproduce unitarity relations by taking into account all relevant open channels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Reactive Transport Modeling of Shale–Fluid Interactions after Imbibition of Fracturing Fluids

Injection of hydraulic fracturing fluid (HFF) into shale formations for unconventional oil/gas production results in chemical reactions in the shale matrix. Our recent experimental study determined the depths to which different types of reactions between the shale matrix and the HFF extended. In the present study, we built continuum-scale reactive transport models to understand the coupling of chemical reactions and the transport of aqueous species in these shale–HFF systems. Calibration of the model with our previous experimental results reveals that it takes hours to months for the shale matrix to completely neutralize the imbibed acids, depending primarily upon the carbonate content of the shale. Both the HFF pH and pore pH affect the location of barite precipitation, resulting in unique barite precipitation profiles extending millimeters into calcite-rich Eagle Ford shale but only tens of micrometers into low-carbonate Marcellus shale. In addition, dissolved oxygen and extracted bitumen are key to reproducing the experimental observation of Fe(III) (oxyhydr)oxide formation in the shale matrix as a result of pyrite dissolution in the shales. A comparison between the modeling results of porosity in the present study to experimentally measured permeabilities in our previous study suggests that chemical reactions occurring at a greater depth than the observable reaction zone might have impacted permeability. Our model serves as a benchmark for efficiently modeling water–rock interactions in similar systems where bulk rock samples react with a solution in batch reactors. Important reactive transport processes were ascertained via modeling, which allows for quantitative prediction of shale–HFF interactions in shale matrices given the shale and HFF compositions.

58 GEOSCIENCES↗

Enhanced HLW Glass Property-Composition Models – Phase 1 (Final Report)

In the present work, the previously compiled glass property-composition database for WTP HLW glasses has been updated with newly available data. The earlier evaluation of the baseline WTP models against the larger data sets, in terms of range of validity and prediction performance, has been used together with a review of the HLW glass data to identify gaps in the composition regions that need to be filled to support the development of expanded models. Subsequently, statistically-designed compositional matrices were developed to begin to fill these composition gaps. Crucible melts have been prepared for the matrix glasses and glass properties have been collected. Finally, preliminary models have been developed that are applicable to the expanded dataset. Additionally, the performance of the preliminary models have been reviewed to assess the relative advantages of developing one set of models for the entire HLW glass composition space or to developing separate models for discrete composition spaces. The longer-term objective of this multi-year task is to develop new glass property-composition models for PCT release, spinel crystallization (T 1% ), nepheline formation on CCC, TCLP Cd response, viscosity, and electrical conductivity to control glass compositions for HLW vitrification at Hanford.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Single discharge of the matrix source of negative hydrogen ions: Influence of the neutral particle dynamics

The study presents two-dimensional (2D) fluid-plasma-model description of a planar-coil inductively-driven discharge, considered as a single element of a matrix source of volume-produced negative hydrogen ions. Whereas the models developed up to now have been directed towards description of the charged particle behavior in the discharge, including that of the negative ions, this model stresses on the role of the neutral particle dynamics and of the surface processes in the formation of the discharge structure. The latter is discussed based on comparison of results obtained for discharges in a flowing gas and at a constant gas pressure as well as for different values of the coefficient of atom recombination on the walls. The conclusions are that the main plasma parameters – electron density and temperature and plasma potential – determining the gas discharge regime stay stable, regardless of changes in the redistribution of the densities of the neutral particles and of the positive ions. With regards to the volume production of the ions, which requires high density of (vibrationally excited) molecules, the impact on the degree of dissociation of the coefficient of atom recombination on the wall is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Influence of Hydraulic Fracturing on Carbon Storage Performance: Hydraulic Fracturing in Carbon Storage

We report that conventional principles of the design and operation of geologic carbon storage (GCS) require injecting CO 2 below the caprock fracturing pressure to ensure the integrity of the storage complex. In nonideal storage reservoirs with relatively low permeability, pressure buildup can lead to hydraulic fracturing of the reservoir and caprock. While the GCS community has generally viewed hydraulic fractures as a key risk to storage integrity, a carefully designed stimulation treatment under appropriate geologic conditions could provide improved injectivity while maintaining overall seal integrity. A vertically contained hydraulic fracture, either in the reservoir rock or extending a limited height into the caprock, provides an effective means to access reservoir volume far from the injection well. Employing a fully coupled numerical model of hydraulic fracturing, solid deformation, and matrix fluid flow, we study the enabling conditions, processes, and mechanisms of hydraulic fracturing during CO 2 injection. A hydraulic fracture's pressure-limiting behavior dictates that the near-well fluid pressure is only slightly higher than the fracturing pressure of the rock and is insensitive to injection rate and mechanical properties of the formation. Although a fracture contained solely within the reservoir rock with no caprock penetration, would be an ideal scenario, poroelastic principles dictate that sustaining such a fracture could lead to continuously increasing pressure until the caprock fractures. We also investigate the propagation pattern and injection pressure responses of a hydraulic fracture propagating in a caprock subjected to heterogeneous in situ stress. In conclusion, the results have important implications for the use of hydraulic fracturing as a tool for managing storage performance.

58 GEOSCIENCES↗

Segmentation of RDX and TNT in X‐Ray Computed Tomography Reconstructions of Melt‐Cast Explosives

ABSTRACT Three‐dimensional mesoscale characterization of heterogeneous melt‐cast high explosives is challenging because of the difficulty differentiating binder from explosive crystals: two functionally different materials which are typically similar in density by design. Here, we report an algorithm which can differentiate hexahydro‐1,3,5‐trinitro‐1,3,5‐triazine (RDX) from 2,4,6‐trinitrotoluene (TNT) in x‐ray computed tomography (CT) volumes with tens of microns resolution. This method allows us to quantify RDX/TNT content, porosity, and RDX domain size. We calibrated the segmentation algorithm using simulated x‐ray CT volumes containing object models of RDX crystals within a TNT matrix. We then segmented and analyzed CT data for Composition B (Comp B), a 60/40 RDX/TNT mixture, and Cyclotol, a 75/25 RDX/TNT mixture. We examined melt‐cast samples fabricated with 100% theoretical maximum density (TMD) and 85% TMD. For the 100% TMD Comp B and Cyclotol samples, the RDX content values calculated by segmentation were 3% and 9% lower, respectively, than the values measured by high‐performance liquid chromatography on material from the same synthesis lots. This result is consistent with the expected underreporting of RDX content resulting from x‐ray CT resolution limits on RDX particles with diameters smaller than 25 µm. The 85% TMD samples were less accurately segmented with our algorithm due to the confounding presence of voids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrodynamic Limit of Multiple SLE

Recently del Monaco and Schleißinger addressed an interesting problem whether one can take the limit of multiple Schramm–Loewner evolution (SLE) as the number of slits N goes to infinity. When the N slits grow from points on the real line ℝ in a simultaneous way and go to infinity within the upper half plane ℍ, an ordinary differential equation describing time evolution of the conformal map g{sub t}(z) was derived in the N→∞ limit, which is coupled with a complex Burgers equation in the inviscid limit. It is well known that the complex Burgers equation governs the hydrodynamic limit of the Dyson model defined on ℝ studied in random matrix theory, and when all particles start from the origin, the solution of this Burgers equation is given by the Stieltjes transformation of the measure which follows a time-dependent version of Wigner’s semicircle law. In the present paper, first we study the hydrodynamic limit of the multiple SLE in the case that all slits start from the origin. We show that the time-dependent version of Wigner’s semicircle law determines the time evolution of the SLE hull, K{sub t}⊂ℍ∪ℝ, in this hydrodynamic limit. Next we consider the situation such that a half number of the slits start from a>0 and another half of slits start from −a<0, and determine the multiple SLE in the hydrodynamic limit. After reporting these exact solutions, we will discuss the universal long-term behavior of the multiple SLE and its hull K{sub t} in the hydrodynamic limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Crashworthiness of a Composite Wing Section: Numerical Investigation of the Bird Strike Phenomenon by Using a Coupled Eulerian–Lagrangian Approach

Crashworthiness is defined as the capability of a structure to guarantee its occupants safety during a crash event. In the present work, the crashworthiness of a composite wing section, subjected to a bird strike event, has been investigated. Indeed, the mechanical behavior of the impacted wing section has been inspected by means of the FE code Abaqus/Explicit. The impact phenomenon has been numerically simulated with a coupled Eulerian–Lagrangian (CEL) approach. The impacted structure has been modeled considering a Lagrangian formulation, while a Eulerian one has been used to model the bird. Moreover, Eulerian–Lagrangian contacts are used for transfer the loads arising from the impact event to the Lagrangian structure. The proposed numerical model, capable to take into account the matrix and fiber traction and compression damages, uses a 3D finite element approach to better predict the deformations and the damages onset and propagation during the impact event. The numerical results, in terms of structural deformation and damage onset and propagation, have been obtained by considering different impact locations and different impact angles.

36 MATERIALS SCIENCE↗

Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds

The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Separation of the 1 + /1 - parity doublet in 20 Ne

The (J, T) = (1, 1) parity doublet in 20 Ne at 11.26 MeV is a good candidate to study parity violation in nuclei. However, its energy splitting is known with insufficient accuracy for quantitative estimates of parity violating effects. To improve on this unsatisfactory situation, nuclear resonance fluorescence experiments using linearly and circularly polarized γ-ray beams were used to determine the energy difference of the parity doublet ΔE = E(1 - ) -E(1 + ) = 3.2(±0.7) stat ($^{+0.6}_{-1.2}$) sys keV and the ratio of their integrated cross sections I$^{(+)}_{s,0}$/I$^{(-)}_{s,0}$ = 29(±3) stat ($^{+14}_{-7}$) sys . Shell-model calculations predict a parity-violating matrix element having a value in the range 0.46–0.83 eV for the parity doublet. The small energy difference of the parity doublet makes 20 Ne an excellent candidate to study parity violation in nuclear excitations.

Parity doublet↗

Comparative environmental and human health evaluations of thermolysis and solvolysis recycling technologies of carbon fiber reinforced polymer waste

Highlights: • Performed life cycle impact assessment (LCIA) of pyrolysis for CFRP waste recycling. • Performed LCIA of solvolysis using supercritical water (SCW) for CFRP recycling. • Performed comparative LCIA of pyrolysis vs SCW-based solvolysis for CFRP recycling. • Closed current knowledge gap about sustainability of SCW-based solvolysis. • Compared LCIA of pyrolyzer’s heating using natural gas in lieu of electricity. - Abstract: This quantitative research aims to compare environmental and human health impacts associated with two recycling technologies of CFRP waste. The ‘baseline’ recycling technology is the conventional thermolysis process via pyrolysis and the ‘alternative’ recycling technology is an emerging chemical treatment via solvolysis using supercritical water (SCW) to digest the thermoset matrix. Two Gate-to-Gate recycling models are developed using GaBi LCA platform. The selected functional unit (FU) is 1 kg CFRP waste and the geographical boundary of this comparative LCIA is defined to be within the U.S. The results of this comparative assessment brought to light new insights about the environmental and human health impacts of CFRP waste recycling via solvolysis using SCW and, therefore, helped close a gap in the current state of knowledge about sustainability of SCW-based solvolysis as compared to pyrolysis. Two research questions are posed to identify whether solvolysis recycling offers more environmental and human health gains relative to the conventional pyrolysis recycling. These research questions lay the basis for formulating two null hypotheses (H{sub 0,1} and H{sub 0,2}) and their associated research hypotheses (H{sub 1,1} and H{sub 1,2}). LCIA results interpretation included ‘base case’ scenarios, ‘sensitivity studies,’ and ‘scenarios analysis.’ The results revealed that: (a) recycling via solvolysis using SCW exhibits no gains in environmental and human health impacts relative to those impacts associated with recycling via pyrolysis and (b) use of natural gas in lieu of electricity for pyrolyzer’s heating reduces the environmental and human health impacts by 37% (lowest) and up to 95.7% (highest). It is recommended that on-going experimental efforts that focus only on identifying the best solvent for solvolysis-based recycling should also consider quantification of the energy intensity as well as environmental and human health impacts of the proposed solvents.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Matrix Approach to Accelerate Spin-Up of CLM5

Numerical models have been developed to investigate and understand responses of biogeochemical cycle to global changes. Steady state, when a system is in dynamic equilibrium, is generally required to initialize these model simulations. However, the spin-up process that is used to achieve steady state pose a great burden to computational resources, limiting the efficiency of global modeling analysis on biogeochemical cycles. This study introduces a new Semi-Analytical Spin-Up (SASU) to tackle this grand challenge. We applied SASU to Community Land Model version 5 and examined its computational efficiency and accuracy. At the Brazil site, SASU is computationally 7 times more efficient than (or saved up to 86% computational cost in comparison with) the traditional native dynamics (ND) spin-up to reach the same steady state. Globally, SASU is computationally 8 times more efficient than the accelerated decomposition spin-up and 50 times more efficient than ND. In summary, SASU achieves the highest computational efficiency for spin-up on site and globally in comparison with other spin-up methods. It is generalizable to wide biogeochemical models and thus makes computationally costly studies (e.g., parameter perturbation ensemble analysis and data assimilation) possible for a better understanding of biogeochemical cycle under climate change.

54 ENVIRONMENTAL SCIENCES↗

Classical molecular dynamics simulation of electronically non-adiabatic processes

Both classical and quantum mechanics (as well as hybrids thereof, i.e., semiclassical approaches) find widespread use in simulating dynamical processes in molecular systems. For large chemical systems, however, which involve potential energy surfaces (PES) of general/arbitrary form, it is usually the case that only classical molecular dynamics (MD) approaches are feasible, and their use is thus ubiquitous nowadays, at least for chemical processes involving dynamics on a single PES (i.e., within a single Born-Oppenheimer electronic state). Our paper reviews recent developments in an approach which extends standard classical MD methods to the treatment of electronically non-adiabatic processes, i.e., those that involve transitions between different electronic states. This approach treats nuclear and electronic degrees of freedom (DOF) equivalently (i.e., by classical mechanics, thereby retaining the simplicity of standard MD), and provides "quantization" of the electronic states through a symmetrical quasi-classical (SQC) windowing model. The approach is seen to be capable of treating extreme regimes of strong and weak coupling between the electronic states, as well as accurately describing coherence effects in the electronic DOF (including the de-coherence of such effects caused by coupling to the nuclear DOF). A survey of recent applications is presented to illustrate the performance of the approach. We also described is a newly developed variation on the original SQC model (found universally superior to the original) and a general extension of the SQC model to obtain the full electronic density matrix (at no additional cost/complexity).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Amplitude analysis of resonant production in three pions

We present some results on the analysis of three pion resonances. The analyses are motivated by the recent release of the largest data set on diffractively produced three pions by the COMPASS collaboration. We construct reaction amplitudes that satisfy fundamental $S$-matrix principles, which allows the use of models that have physical constraints to be used in fitting data. The models are motivated by the isobar model that satisfy unitarity constraints. The model consist of a Deck production amplitude with which final state interactions are constrained by unitarity. We employ the isobar model where two of the pions form a quasi-stable particle. The analysis is performed in the high-energy, single Regge limit. We specifically discuss the examples of the three pion $J^{PC}=2^{-+}$ resonance in the $\rho\pi$ and $f_2\pi$ channels.

Jackura, Andrew↗

Nonadiabatic anharmonic electron transfer

The effect of an inner sphere, local mode vibration on an electron transfer is modeled using the nonadiabatic transition probability (rate) expression together with both the anharmonic Morse and the harmonic oscillator potential. For an anharmonic inner sphere mode, a variational analysis uses harmonic oscillator basis functions to overcome the difficulties evaluating Morse-model Franck-Condon overlap factors. Individual matrix elements are computed with the use of new, fast, robust, and flexible recurrence relations. The analysis therefore readily addresses changes in frequency and/or displacement of oscillator minimums in the different electron transfer states. Direct summation of the individual Boltzmann weighted Franck-Condon contributions avoids the limitations inherent in the use of the familiar high-temperature, Gaussian form of the rate constant. The effect of harmonic versus anharmonic inner sphere modes on the electron transfer is readily seen, especially in the exoergic, inverted region. The behavior of the transition probability can also be displayed as a surface for all temperatures and values of the driving force/exoergicity {Delta}=-{Delta}G. The temperature insensitivity of the transfer rate is clearly seen when the exoergicity equals the collective reorganization energy ({Delta}={Lambda}{sub s}) along a maximum ln (w) vs. {Delta} ridge of the surface. The surface also reveals additional regions for {Delta} where ln (w) appears to be insensitive to temperature, or effectively activationless, for some kinds of inner sphere contributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhanced ionic diffusion in ionomer-filled nanopores

Coarse-grained simulations in the united-atom-model approximation are used to investigate confinement-induced morphological changes in Nafion-like ionomers. The system we study models a cylindrical pore in a hydrophobic matrix of supporting material with pore diameters that vary from 0.7 to 3.96 nm. Simulation results indicate a strong dependence of the equilibrium ionomer structures both on the pore diameter and on the sulfonate concentration in the pore. In the case of larger pores, the ionic clustering has the shape of a branched wire-like network oriented parallel to the pore axis. In the case of narrow pores, the ionic clusters occupy the pore center and exhibit strong density modulations both along the pore axis and across the pore diameter. The calculated diffusion coefficients for the ions indicate a sharp increase within the narrow pores. This finding is explained by ballistic-type ionic motion at shorter times and by the collective motion of ions in hydrophilic clusters. The influence of the hydrophobic walls on the distribution of ions and solvent molecules is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗