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 1,783 records · Page 99

Damage Prediction Models for Advanced Materials and Composites

In the present study, the assessment and evaluation of various acoustic tile designs were conducted using three-dimensional finite element analysis, which included static analysis, thermal analysis and modal analysis of integral and non-integral tile design options. Various benchmark specimens for acoustic tile designs, including CMC integral T-joint and notched CMC plate, were tested in both room and elevated temperature environment. Various candidate ceramic matrix composite materials were used in the numerical modeling and experimental study. The research effort in this program evolved from numerical modeling and concept design to a combined numerical analysis and experimental study. Many subjects associated with the design and performance of the acoustic tile in jet engine exhaust nozzle have been investigated.

Xie, Ming↗

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↗

Elasto-plastic analysis of interface layers for fiber reinforced metal matrix composites

The mismatch in coefficients of thermal expansion (CTE) of fiber and matrix in metal matrix composites reinforced with ceramic fibers induces high thermal stresses in the matrix. Elasto-plastic analyses - with different degrees of simplification and modelization - show that an interface layer with a sufficiently high CTE can reduce the tensile hoop stress in the matrix substantially.

Doghri, I.↗

How to Train Your Jacobian: Least-Squares System Identification for Space-Based Coronagraphy

Stellar coronagraphs use closed-loop focal-plane wavefront sensing and control algorithms to create high-contrast dark zones suitable for imaging exoplanets and exozodiacal dust clouds around nearby stars. Model-based algorithms are susceptible to model mismatch, wherein a departure of the coronagraph's true optical characteristics from the assumed model causes reduced control loop performance. Here, we describe a simple technique for empirically tuning the wavefront control Jacobian matrix using applied deformable mirror commands and observed images. This mitigates model mismatch and recovers near-optimal control loop performance.

coronagraphy↗

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↗

Load redistribution considerations in the fracture of ceramic matrix composites

Using a macroscopic viewpoint, composite laminae are homogeneous orthotropic solids whose directional strengths are random variables. Incorporation of these random variable strengths into failure models, either interactive or noninteractive, allows for the evaluation of the lamina reliability under a given stress state. Using a noninteractive criterion for demonstration purposes, laminate reliabilities are calculated assuming previously established load sharing rules for the redistribution of load as the failure of laminae occur. The matrix cracking predicted by ACK theory is modeled to allow a loss of stiffness in the fiber direction. The subsequent failure in the fiber direction is controlled by a modified bundle theory. Results are compared with previous models which did not permit separate consideration of matrix cracking, as well as to results obtained from experimental data. The effects of variations from the ideal physical geometry which is normally used to depict the matrix cracking are also studied.

Thomas, David J.↗

Understanding Variability in Mechanical Behavior of Adhesively Bonded Sandwich Joints due to Stochastic Material Damage

A mesoscale progressive damage analysis was applied to the analysis of recently tested adhesively bonded sandwich panels, in which ten nominally identical tension loaded panels failed in three different failure modes, with peakloads ranging from 13.0 –16.1 kips. A parametric finite element model for progressive damage analysis of an adhesively bonded joint (ABJ) considering matrix and fiber damage, and delamination cracks in unidirectional and fabric composite sections of the joint, as well as core crushing, and adhesive fracture was developed. The model was developed using Abaqus/Explicit and utilized continuum damage mechanics material models for intralaminar damage in the facesheet and doubler plies, cohesive elements for delamination and adhesive debonding, and discretized cells considering plasticity to model crushing of the aluminum honeycomb core. Analysis predictions obtained with a pristine model agreed well with the average experimental peak load, strain in the joint, and post-test damage states. Then, various methods to incorporate stochastic material variability into ABJ model were evaluated. These methods included simulation of pre-existing matrix damage, pre-existing adhesive porosity, and uniform modification of adhesive properties. These material and manufacturing defects changed the damage mode predictions in some cases. The analysis of stochastic strength due to the considered defects revealed unintuitive damage mode interaction in the ABJ, explaining some structure-property relationships observed experimentally.

Richard Larson↗

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↗