Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Matrix”

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 451 records · Page 25

Backscattering by nonspherical particles - A review of methods and suggested new approaches

Scattering of electromagnetic radiation near the backward direction is more sensitive to particle shape than scattering near the forward direction. Mie theory is therefore of dubious applicability to predicting backscattering by atmospheric particles known to be irregular or to inverting measurements on such particles. An irregular particle is one with an uncertain shape. In the face of uncertainty one must adopt a statistical approach in which scattering properties of ensembles are determined. To obtain ensemble averages, a basis is needed for averaging over a set of electromagnetic microstates. Ensemble averages based on the Rayleigh theory for small ellipsoids and on the T matrix method for spheroids agree better with measurements than Mie theory does. The coupled-dipole method also provides a basis for ensemble averaging. This method also leads to a simple physical interpretation of why backscattering is so sensitive to particle shape and can be used to calculate scattering by one- and two-dimensional analogs to three-dimensional irregular particles.

Bohren, Craig F.↗

Resolved resonance region evaluations of n+ 206,207,208 Pb for fast spectrum applications

Resolved resonance region evaluations of the major isotopes of natural lead, 206 Pb, 207 Pb, 208 Pb, have been per formed to support the development of Generation IV reactors. In this study, validation of nuclear data for lead fast reactors was performed with simulations of shielding benchmarks, integral critical benchmarks, and quasi-differential scattering measurements. Sensitivity analyses of these systems showed that elastic scattering reactions above 100 keV were the dominant reactions driving system performance. The resolved resonance regions (RRRs) of the lead isotopes extend past 100 keV, making the RRR an ideal starting point to evaluate lead cross sections. Since the R-matrix requires knowledge of bound, distant, and observed resonances, it was necessary to evaluate from 10 -5 eV up to the respective limit of the RRR. The 208 Pb RRR evaluation was extended to 1.5 MeV in order to obtain resonance parameters used to calculate new elastic scattering angular distributions up to 1.5 MeV. Resonance parameter uncertainties and covariance were generated using the R-matrix code SAMMY. The new RRR parameters show a direct improvement to the scattering kernel below 1.5 MeV which in turn greatly improves fast critical experiments over ENDF/B-VIII.0.

206Pb↗

Driving rapid atomic order in MnAl via low-magnitude magnetic field annealing

Application of a mild (60 mT), uniform magnetic field during short-term thermal treatment of kinetically retained, atomically disordered (paramagnetic) ε-MnAl was found to deliver a significant ~50 % increase in the formation of L1 0 atomically ordered (ferromagnetic) τ-MnAl product phase, compared to that produced by conventional (i.e., zero-field) annealing under identical thermal conditions. The magnetic field, applied in a passive closed-circuit configuration during annealing, induced significant changes in the structural, magnetic, and phase evolution of the material. Computational results based on electronic structure calculations demonstrate that the effective magnetic susceptibility of τ-MnAl is sensitive to the orientation, rather than the magnitude, of an applied magnetic field in the vicinity of the Curie temperature. The uniaxial magnetocrystalline anisotropy of the L1 0 structure is proposed to act as a filter for selective propagation of the population of τ-MnAl variants that are favorably aligned with the applied field. In this manner, crystallographic “gridlock” is alleviated that would otherwise arise from the coexistence of multiple, energetically equivalent τ-phase variants within the parent ε-phase matrix. These results confirm that static, low-magnitude magnetic field annealing is able to accelerate L1 0 atomic ordering in the MnAl system and likely can exert similar influences in relevant magnetic systems, facilitating efficient tailoring of structure-sensitive magnetic properties for the manufacture of magnetic materials.

36 MATERIALS SCIENCE↗

The influence of random packed column parameters on the liquid holdup and interfacial area

Abstract Carbon dioxide capture via solvent absorption in packed columns has emerged as a potential technology to mitigate coal‐fired power plant CO 2 emissions. Parameters, including packing types, solvent properties, and operating conditions, could potentially affect the packed column CO 2 capture efficiency. To understand the importance of those parameters and help packed column optimization, a design of experiments (DoEs) method was proposed to generate input parameter matrix. Combined with multiphase computational fluid dynamics (CFD), the random packed column parameter influence on the liquid holdup and interfacial area can be efficiently investigated. Surrogate‐based sensitivity analysis shows that the solvent flow rate and contact angle are key factors dictating liquid holdup and interfacial area. On the other hand, solvent viscosity has a marginal impact on the interfacial area. The sensitivity scores were calculated for each input parameter to guide the selection of dimensionless numbers for the liquid holdup and interfacial area correlation development.

42 ENGINEERING↗

Spontaneous Raman–LIF–CO–OH measurements of species concentration in turbulent spray flames

This paper presents new measurements of species concentrations, temperature and mixture fraction in selected regions of a turbulent ethanol spray flame. The line-Raman–LIF–COsingle bondOH setup developed at the Sandia's Combustion Research Facility is utilised to probe regions of a spray flame where laser breakdown of liquid droplets is avoided and the remaining interferences can be corrected. The spray flame is stabilised on the piloted Sydney needle spray burner, where axial translation of the liquid injecting needle in the air-blast stream can transition the spray from dilute to dense. The solution to obtaining successful measurements is found to be multifaceted and includes: the appropriate selection of flame conditions; high sensitivity of the Raman detection system permitting reduced laser energies; development of a pre-processing algorithm to reject strong droplet interferences; and application of the hybrid matrix inversion method combined with wavelet denoising to account for interference corrections and noise at the very low signal levels obtained. Unique and necessary for the successful measurements reported in this paper, a pre-processing algorithm is outlined that removes data points corrupted with strong interferences from droplets. These interferences arise from a range of sources, but the most intense are due to the laser interaction with surrounding mist or liquid fragments, such that measurements near the jet centreline are corrupted and hence discarded. Reliable measurements of mixture fraction, temperature obtained from the sum of the species number densities, and species mole fractions are reported for regions in the flames sufficiently far from the centreline. The paper demonstrates the feasibility of the judicious use of Raman scattering in turbulent spray flames, the results of which will be extremely useful for validating numerical simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗

The Lyman- α forest from LBGs: First 3D correlation measurement with DESI and prospects for cosmology

The Lyman-α (Lyα) forest is a key tracer of large-scale structure at redshifts z > 2, traditionally studied using the spectra of luminous but relatively rare quasars. In this work, we explore the viability of using the fainter yet significantly more abundant Lyman Break Galaxies (LBGs) as alternative background sources for Lyα forest studies. We analyze 4,151 Lyα forest skewers extracted from LBG spectra obtained in the DESI pilot surveys conducted in the COSMOS and XMM-LSS fields. From this dataset, we present the first measurement of the Lyα forest auto-correlation function derived exclusively from LBG spectra, probing comoving separations up to 48 h -1 Mpc at an effective redshift of z eff = 2.70. The measured LBG Lyα forest auto-correlation is consistent with that derived from DESI DR2 quasar Lyα forest spectra at a comparable redshift, validating the use of LBGs as reliable background sources for Lyα forest analyses. In addition, we measure the cross-correlation between the LBG Lyα forest and the positions of 13,362 galaxies, demonstrating that this observable serves as a sensitive diagnostic for assessing the precision and accuracy of galaxy redshift estimates, and for identifying and correcting systematic offsets. Finally, using both synthetic LBG spectra and Fisher matrix forecasts, we show that a future wide-area survey covering ∼5,000 deg 2 , targeting 1,000 LBGs per square degree at signal-to-noise levels comparable to our sample, could enable LBG-based Lyα forest baryon acoustic oscillation (BAO) measurements with expected uncertainties of σ α ISO = 0.4% (isotropic) and σ α AP = 1.3% (Alcock-Paczynski). This performance is further enhanced when combining the BAO analysis with a Lyα forest Full Shape (FS) approach, yielding a predicted uncertainty of σ α ISO FS = 0.6%. These results open a new avenue for precision cosmology at high redshift using the Lyα forest in dense LBG samples.

Lyman alpha forest↗

Searching for parity violation in SDSS DR16 Lyman-α forest data

The four-point correlation function is the lowest order correlation function for scalar fields that can be used to probe statistical parity invariance in an isotropic universe. There are intriguing claims of detection of parity violation in the 4-point function of BOSS galaxy clustering data. We apply the same estimator to the public SDSS Data Release 16 Lyman-α forest data. Lyman-α forest data probes a different redshift range and is sensitive to a different density regime using a completely different technique. A detection would therefore be a strong indication of new physics. We identify an accurate covariance matrix as a crucial impediment to performing this measurement accurately, consistent with existing literature on galaxy 4-point function. Here, we discuss several approaches to estimating the covariance matrix, several of which produce spurious detection. Using a robust, but very suboptimal, covariance matrix derived from subsample bootstrapping, we find no evidence for parity violation.

79 ASTRONOMY AND ASTROPHYSICS↗

Self-assembled complex micelle phase stability in ABA-type triblock copolymers and related core-shell bottlebrushes

We report the synthesis and temperature-dependent morphologies of a series of polylactide--poly (-decalactone)--polylactide (LDL) triblock copolymers with M n = 16.0 – 18.1 kg / mol and volume fractions f L = 0.27 – 0.31 and associated core-shell bottlebrush (csBB) polymers, which derive from enchaining LDL triblocks through a polymerizable midchain functionality. While the LDL triblocks form micellar Frank-Kasper A15 and phases due to the conformational asymmetry of this monomer pair, the csBB morphologies sensitively depend on the backbone degree of polymerization ( N bb ). At low N bb values, micellar Frank-Kasper phases with the brush backbone situated in the matrix domain are stable, albeit with a modest reduction in the mean interfacial curvature evidenced by a to A15 order-to-order transition. However, larger N bb values drive csBBs to form hexagonally packed cylinders phases. This N bb -dependent phase behavior is rationalized in terms of a star-to-bottlebrush transition. At low N bb values, the csBBs are akin to star polymers with pointlike junctions that can support complex micelle packings. As N bb increases, the csBBs adopt cylindrical molecular geometries with extended backbones situated in the matrix domain that prefer hexagonally packed cylinders morphologies. Published by the American Physical Society 2024

Materials Science↗

Code Coverage Status of the ARC Code PERSENT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. PERSENT fulfills the role of generating reactivity coefficients for a given time point of a REBUS calculation usable in a point kinetics based safety analysis capability. PERSENT also provides a sensitivity coefficient capability on eigenvalue, reactivity worth, and several other key coefficients that are used in the follow-on safety analysis. Given a co-variance matrix, PERSENT can carry out the uncertainty quantification to indicate the amount of error in the reactivity coefficients derived from the errors in the cross section measurements. With continued improvement of computational resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. PERSENT is part of the ARC code system and is built around DIF3D-VARIANT and the flux solution it provides. The purpose of the present work is to identify a set of test problems for PERSENT and assess the code coverage of PERSENT for those test problems. PERSENT treats the DIF3D executable as an external executable and thus the code coverage considerations only need to focus on the PERSENT source code and only a fraction of the connected modules in the existing ARC software library. The goal is to document what parts of the existing PERSENT code are touched by the set of test problems and which are not. Because the verification work done on PERSENT was focused on the most common uses of PERSENT for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that users of PERSENT rely upon. The code coverage analysis of PERSENT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of PERSENT. The detailed coverage tables are given for each submodule of PERSENT. Most of the uncovered parts/files could be easily ignored because they are either for error message and debugging output or not needed by PERSENT today. Only a few uncovered parts of PERSENT deserve extending the verification test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Influence of convection on microstructure

The primary motivation for this research has been to determine the cause for space processing altering the microstructure of some eutectics, especially the MnBi-Bi eutectic. Prior experimental research at Grumman and here showed that the microstructure of MnBi-Bi eutectic is twice as fine when solidified in space or in a magnetic field, is uninfluenced by interfacial temperature gradient, adjusts very quickly to changes in freezing rate, and becomes coarser when spin-up/spin-down (accelerated crucible rotation technique) is used during solidification. Theoretical work at Clarkson predicted that buoyancy driven convection on earth could not account for the two fold change in fiber spacing caused by solidification in space. However, a lamellar structure with a planar interface was assumed, and the Soret effect was not included in the analysis. Experimental work at Clarkson showed that the interface is not planar, and that MnBi fibers project out in front of the Bi matrix on the order of one fiber diameter. Originally four primary hypotheses were to be tested under this current grant: (1) a fibrous microstructure is much more sensitive to convection than a lamellar microstructure, which was assumed in our prior theoretical treatment; (2) an interface with one phase projecting out into the melt is much more sensitive to convection than a planar interface, which was assumed in our prior theoretical treatment; (3) the Soret effect is much more important in the absence of convection and has a sufficiently large influence on microstructure that its action can explain the flight results; and (4) the microstructure is much more sensitive to convection when the composition of the bulk melt is off eutectic. As reported previously, we have learned that while a fibrous microstructure and a non-planar interface are more sensitive to convection than a lamellar microstructure with a planar interface, the influence of convection remains too small to explain the flight and magnetic field results. Similarly addition of the Soret effect does not explain the flight and magnetic field results.

Wilcox, William R.↗

Preconditioned domain decomposition scheme for three-dimensional aerodynamic sensitivity analysis

A preconditioned domain decomposition scheme is introduced for the solution of the 3D aerodynamic sensitivity equation. This scheme uses the iterative GMRES procedure to solve the effective sensitivity equation of the boundary-interface cells in the sensitivity analysis domain-decomposition scheme. Excluding the dense matrices and the effect of cross terms between boundary-interfaces is found to produce an efficient preconditioning matrix.

Eleshaky, Mohammed E.↗

X-ray Mapping of Terrestrial and Extraterrestrial Materials Using the Electron Microprobe

Lunar samples returned from the Apollo program motivated development of the Bence-Albee algorithm for the rapid and accurate analysis of lunar materials, and established interlaboratory comparability through its common use. In the analysis of mineral and rock fragments it became necessary to combine micro- and macroscopic analysis by coupling electron-probe microanalysis (EPMA) with automated stage point counting. A coarse grid that included several thousand points was used, and initially wavelength-dispersive (WDS) and later energydispersive (EDS) data were acquired at discrete stage points using approx. 5 sec count times. A approx 50 micrometer beam diameter was used for WDS and up to 500 micrometer beam diameter for EDS analysis. Average analyses of discretely sampled phases were coupled with the point count data to calculate the bulk composition using matrix algebra. Use of a defocused beam resulted in a contribution from multiple phases to each analytical point, and the analytical data were deconvolved relative to end-member phase chemistry on the fly. Impressive agreement was obtained between WDS and EDS measurements as well as comparison with bulk chemistry obtained by other methods. In the 30 years since these methods were developed, significant improvements in EPMA automation and computer processing have taken place. Digital beam control allows routine collection of x-ray maps by EDS, and stage mapping for WDS is conducted continuously at slew speed and incrementally by sampling at discrete points. Digital pulse processing in EDS systems has significantly increased the throughput for EDS mapping, and the ongoing development of Si-drift detector systems promises mapping capabilities rivaling WDS systems. Spectrum imaging allows a data cube of EDS spectra to be acquired and sophisticated processing of the original data is possible using matrix algebra techniques. The study of lunar and meteoritic materials includes the need to conveniently: (1) Characterize the sample at microscopic and macroscopic scales with relatively high sensitivity, (2) Determine the modal abundance of minerals, and (3) Identify and relocate discrete features of interest in terms of size and chemistry. The coupled substitution of cations in minerals can result in significant variation in mineral chemistry, but at similar average Z, leading to poor backscattered-electron (BSE) contrast discrimination of mineralogy. It is necessary to discriminate phase chemistry at both the trace element level and the major element level. To date, the WDS of microprobe systems is preferred for mapping due to high throughput and the ability to obtain the necessary intensity to discriminate phases at both trace and major element concentrations. It is desirable to produce fully quantitative compositional maps of geological materials, which requires the acquisition of k-ratio maps that are background and dead-time corrected, and which have been corrected by phi(delta z> or an equivalent algorithm at each pixel. To date, turnkey systems do not allow the acquisition of k-ratio maps and the rigorous correction in this manner. X-ray maps of a chondrule from the Ourique meteorite, and a comb-layered xenolith from the San Francisco volcanic field, have been analyzed and processed to extract phase information. The Ourique meteorite presents a challenge due to relatively low BSE contrast, and has been studied using spectrum imaging. X-ray maps for Si, Mg, and FeK(alpha) were used to produce RGB images. The xenolith sample contains sector-zoned augite, olivine, plagioclase, and basaltic glass. X-ray maps were processed using Lispix and ImageJ software to produce mineral phase maps. The x-ray maps for Mg, Ca, and Ti were used with traceback to generate binary images that were converted to RGB images. These approaches are successful in discriminating phases, but it is desirable to achieve the methods that were used on lunar samples 30 years ago on current microprobe systems. Curnt research includes x-ray mapping analysis of the Dalgety Downs chondrite by micro x-ray fluorescence and spectrum imaging, in collaboration with Kenny Witherspoon of IXRF Systems and Dale Newbury of NIST.

Carpenter, P.↗

Simplified micromechanical equations for thermal residual stress analysis of coated fiber composites

The fabrication of metal matrix composites poses unique problems to the materials engineer. The large thermal expansion coefficient (CTE) mismatch between the fiber and matrix leads to high tensile residual stresses at the fiber/matrix (F/M) interface which could lead to premature matrix cracking during cooldown. Fiber coatings could be used to reduce thermal residual stresses. A simple closed form analysis, based on a three phase composite cylinder model, was developed to calculate thermal residual stresses in a fiber/interphase/matrix system. Parametric studies showed that the tensile thermal residual stresses at the F/M interface were very sensitive to the CTE and thickness of the interphase layer. The modulus of the layer had only a moderate effect on tensile residual stresses. For a silicon carbide titanium aluminide composite, the tangential stresses were 20 to 30 pct. larger than the axial stresses, over a wide range of interphase layer properties, indicating a tendency to form radial matrix cracks during cooldown. Guidelines for the selection of appropriate material properties of the fiber coating were also derived in order to minimize thermal residual stresses in the matrix during fabrication.

Naik, R. A.↗

Search for Majorana neutrinos in same-sign WW scattering events from pp collisions at $\sqrt{s}=13$ TeV

A search for Majorana neutrinos in same-sign WW scattering events is presented. The analysis uses $\sqrt{s}=13$ TeV proton–proton collision data with an integrated luminosity of 140 fb -1 recorded during 2015–2018 by the ATLAS detector at the Large Hadron Collider. The analysis targets final states including exactly two same-sign muons and at least two hadronic jets well separated in rapidity. The modelling of the main backgrounds, from Standard Model same-sign WW scattering and WZ production, is constrained with data in dedicated signal-depleted control regions. The distribution of the transverse momentum of the second-hardest muon is used to search for signals originating from a heavy Majorana neutrino with a mass between 50 GeV and 20 TeV. No significant excess is observed over the background expectation. The results are interpreted in a benchmark scenario of the Phenomenological Type-I Seesaw model. In addition, the sensitivity to the Weinberg operator is investigated. Upper limits at the 95% confidence level are placed on the squared muon-neutrino–heavy-neutrino mass-mixing matrix element |V μN | 2 as a function of the heavy Majorana neutrino’s mass m N , and on the effective μμ Majorana neutrino mass |m μμ |.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Engineering applications of heuristic multilevel optimization methods

Some engineering applications of heuristic multilevel optimization methods are presented and the discussion focuses on the dependency matrix that indicates the relationship between problem functions and variables. Coordination of the subproblem optimizations is shown to be typically achieved through the use of exact or approximate sensitivity analysis. Areas for further development are identified.

Barthelemy, Jean-Francois M.↗

Engineering applications of heuristic multilevel optimization methods

Some engineering applications of heuristic multilevel optimization methods are presented and the discussion focuses on the dependency matrix that indicates the relationship between problem functions and variables. Coordination of the subproblem optimizations is shown to be typically achieved through the use of exact or approximate sensitivity analysis. Areas for further development are identified.

Barthelemy, Jean-Francois M.↗

Engineering applications of heuristic multilevel optimization methods

Some engineering applications of heuristic multilevel optimization methods are presented and the discussion focuses on the dependency matrix that indicates the relationship between problem functions and variables. Coordination of the subproblem optimizations is shown to be typically achieved through the use of exact or approximate sensitivity analysis. Areas for further development are identified.

Barthelemy, Jean-Francois M.↗