Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization 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 163 records · Page 9

Coupling Fractional Thermal Runaway Calorimetry (FTRC) Results with Statistical Analysis Methods

Important factors to consider when optimizing Lithium-ion (Li-ion) battery thermal management systems, from a battery safety perspective, include the total energy released during a single cell thermal runaway (TR) event and the fraction of the total energy that is released through the cell casing versus through the ejecta material. A fractional thermal runaway calorimeter (FTRC) was used to examine the TR behavior of the following 18650-format Li-ion cells: (1) Molicel 18650-J, (2) LG 18650-HG2, (3) Samsung 18650-30Q, (4) 3.35 Ah LG 18650 test cell, (5) Sony 18650-VC7, and (6) LG 18650-MJ1. These cells were intentionally selected to represent a variety of manufacturers, chemistries, capacities, and safety features. Because of the large number of random and non-random variables associated with the FTRC results, a linear regression model, developed using engineering statistical methods, was utilized to provide final interpretation of the TR behavior for each cell type and design variable combination.

Walker, William↗

Comparison of Resource Requirements for a Wind Tunnel Test Designed with Conventional vs. Modern Design of Experiments Methods

The factors that determine data volume requirements in a typical wind tunnel test are identified. It is suggested that productivity in wind tunnel testing can be enhanced by managing the inference error risk associated with evaluating residuals in a response surface modeling experiment. The relationship between minimum data volume requirements and the factors upon which they depend is described and certain simplifications to this relationship are realized when specific model adequacy criteria are adopted. The question of response model residual evaluation is treated and certain practical aspects of response surface modeling are considered, including inference subspace truncation. A wind tunnel test plan developed by using the Modern Design of Experiments illustrates the advantages of an early estimate of data volume requirements. Comparisons are made with a representative One Factor At a Time (OFAT) wind tunnel test matrix developed to evaluate a surface to air missile.

DeLoach, Richard↗

A Method of Moments Wide Band Adaptive Rational Interpolation Method for High-Quality Factor Resonant Cavities

A new adaptive rational interpolation method is proposed to obtain the wideband frequency response of a resonant cavity simulated with the method of moments (MoM). This interpolation method uses both the Loewner matrix to construct a rational expression for the solution vector of MoM’s matrix system and an error estimator generated by the solution vectors and their derivatives. This error estimator is implemented in the adaptive procedure to gain a minimum set of frequencies and solution vectors required in the interpolation. The resulting set of frequencies and solution vectors is applied to interpolate other system variables, such as shielding effectiveness and input impedance. Here numerical results of a slotted cylindrical cavity supporting high-quality factor resonances are presented, showing that the new rational interpolation method is accurate and efficient in interpolating the complicated resonant response of the solution vector functions.

42 ENGINEERING↗

Precision Measurement of the Neutron Magnetic Form Factor via the Ratio Method at Jefferson Lab Hall A

Protons and neutrons, collectively known as nucleons, are composed of quarks and gluons. The Sachs electromagnetic form factors encode information about the spatial distributions of charge and magnetization in the nucleon, particularly at low momentum transfer. In particular, the neutron magnetic form factor (GMn) provides crucial information about the distribution of magnetization inside the neutron and helps constrain theoretical models of nucleon structure. Quasi-elastic electron scattering from deuterium was measured up to Q^2=13.5 GeV^2 using the Super BigBite Spectrometer in Hall A at Jefferson Lab. In this work, the neutron magnetic form factor GMn was extracted at Q^2 = 3.0 GeV^2 and Q^2=4.5 GeV^2 using the Ratio Method. These results represent a subset of the full dataset collected in this experiment, which extended to significantly higher Q^2. The extracted GMn values agree with the existing global fit within approximately two standard deviations at Q^2=3.0 and show excellent agreement at Q^2=4.5. The measurements achieved systematic uncertainties of about 2% and statistical uncertainties below 0.5%, among the most precise determinations of GMn at these kinematics. These results demonstrate the robustness of the experimental technique and provide an important validation point for future extractions at higher Q^2, where data remain scarce. In addition, the GRINCH heavy gas Cherenkov detector—a key component of the experimental apparatus—was commissioned and achieved an electron detection efficiency of approximately 97%, supporting reliable particle identification. Together, the analysis presented here advances both our understanding of nucleon structure and the validation of the experimental methods and instrumentation used to access it.

Satnik, Maria [College of William and Mary, Willia↗

The Intercalibration of Geostationary Visible Imagers Using Operational Hyperspectral SCIAMACHY Radiances

Spectral band differences between sensors can complicate the process of intercalibration of a visible sensor against a reference sensor. This can be best addressed by using a hyperspectral reference sensor whenever possible because they can be used to accurately mitigate the band differences. This paper demonstrates the feasibility of using operational Scanning Imaging Absorption Spectrometer for Atmospheric Cartography (SCIAMACHY) large-footprint hyperspectral radiances to calibrate geostationary Earth-observing (GEO) sensors. Near simultaneous nadir overpass measurements were used to compare the temporal calibration of SCIAMACHY with Aqua Moderate Resolution Imaging Spectroradiometer band radiances, which were found to be consistent to within 0.44% over seven years. An operational SCIAMACHY/GEO ray-matching technique was presented, along with enhancements to improve radiance pair sampling. These enhancements did not bias the underlying intercalibration and provided enough sampling to allow up to monthly monitoring of the GEO sensor degradation. The results of the SCIAMACHY/GEO intercalibration were compared with other operational four-year Meteosat-9 0.65-μm calibration coefficients and were found to be within 1% of the gain, and more importantly, it had one of the lowest temporal standard errors of all the methods. This is more than likely that the GEO spectral response function could be directly applied to the SCIAMACHY radiances, whereas the other operational methods inferred a spectral correction factor. This method allows the validation of the spectral corrections required by other methods.

Doelling, David R.↗

Local self-interaction correction method with a simple scaling factor

A recently proposed local self-interaction correction (LSIC) method [Zope et al., J. Chem. Phys., 2019, 151, 214108] when applied to the simplest local density approximation provides a significant improvement over standard Perdew–Zunger SIC (PZSIC) for both equilibrium properties such as total or atomization energies as well as properties involving stretched bond such as barrier heights. The method uses an iso-orbital indicator to identify the single-electron regions. To demonstrate the LSIC method, Zope et al. used the ratio z σ of von Weizsäcker τ W σ and total kinetic energy densities τ σ , (z σ = τ W σ /τ σ ) as a scaling factor to scale the self-interaction correction. The present work further explores the LSIC method using a ratio of orbital and spin densities as a simpler scaling factor in place of the ratio of kinetic energy densities. We compute a wide array of both, equilibrium and non-equilibrium properties using LSIC and orbital scaling methods using this simple scaling factor and compare them with previously reported results. Our study shows that LSIC with the simple scaling factor performs better than PZSIC, with results comparable to those obtained by LSIC( z σ ) for most properties, but has slightly larger errors than LSIC( z σ ). Furthermore, we study the binding energies of small water clusters using both scaling factors. Our results show that LSIC with zσ has limitations in predicting the cluster binding energies of weakly bonded systems due to the inability of z σ to distinguish weakly bonded regions from slowly varying density regions. LSIC when used with the density ratio as a scaling factor, on the other hand, provides a good description of water cluster binding energies, thus highlighting the appropriate choice of the iso-orbital indicator.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the equivalent Dancoff factor cell method performing slowing down calculations on the fly for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for pin-cell, unit-cell, and fuel assembly problems with various resonance self-shielding conditions based on very high temperature reactor, high temperature test reactor, and Empire benchmark cores, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also showed that the methodology worked well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The use of parabolic variations and the direct determination of stress intensity factors using the BIE method

Two advances in the numerical techniques of utilizing the BIE method are presented. The boundary unknowns are represented by parabolas over each interval which are integrated in closed form. These integrals are listed for easy use. For problems involving crack tip singularities, these singularities are included in the boundary integrals so that the stress intensity factor becomes just one more unknown in the set of boundary unknowns thus avoiding the uncertainties of plotting and extrapolating techniques. The method is applied to the problems of a notched beam in tension and bending, with excellent results.

Mendelson, A.↗

The use of parabolic variations and the direct determination of stress intensity factors using the BIE method

Two advances in the numerical techniques of utilizing the BIE method are presented. The boundary unknowns are represented by parabolas over each interval which are integrated in closed form. These integrals are listed for easy use. For problems involving crack tip singularities, these singularities are included in the boundary integrals so that the stress intensity factor becomes just one more unknown in the set of boundary unknowns thus avoiding the uncertainties of plotting and extrapolating techniques. The method is applied to the problems of a notched beam in tension and bending, with excellent results.

Mendelson, A.↗

Shuttle user analysis (study 2.2). Volume 4: Standardized subsystem modules analysis

The capability to analyze payloads constructed of standardized modules was provided for the planning of future mission models. An inventory of standardized module designs previously obtained was used as a starting point. Some of the conclusions and recommendations are: (1) the two growth factor synthesis methods provide logical configurations for satellite type selection; (2) the recommended method is the one that determines the growth factor as a function of the baseline subsystem weight, since it provides a larger growth factor for small subsystem weights and results in a greater overkill due to standardization; (3) the method that is not recommended is the one that depends upon a subsystem similarity selection, since care must be used in the subsystem similarity selection; (4) it is recommended that the application of standardized subsystem factors be limited to satellites with baseline dry weights between about 700 and 6,500 lbs; and (5) the standardized satellite design approach applies to satellites maintainable in orbit or retrieved for ground maintenance.

Source record↗

A Practical Approach to Uncertainty Quantification Using Probability Boxes

To date, while the use of CFD for aerospace vehicle design and development is prevalent, the documentation of uncertainties associated with the simulations are rare. Instead, the current state-of-the-art relies heavily on the experience of the CFD practitioner to estimate the uncertainty associated with their simulations through simple sensitivity studies or subject matter expertise. This practice will have to be replaced with a formal uncertainty quantification (UQ) process if CFD is to play an expanded role in the research and engineering design community, test and evaluation community, and ultimately certification for flight. Accounting for uncertainties in a formal manner is a tedious process. Moreover, the typical CFD practitioner is not likely to be familiar with formal UQ methods. These factors have prevented the adoption of UQ methods in the engineering design and development cycle. This presentation will outline a credible approach to UQ using Probability Boxes that is straightforward to apply, and can readily be automated using existing UQ tool sets such as the DAKOTA packaged developed at Sandia. The added expense incurred when moving away from a deterministic CFD process to a stochastic one that captures uncertainties to enable risk-informed decision making will be discussed, as well as effective ways to reduce the computational costs.

Uncertainty Quantification↗

Sparse Approximate Multifrontal Factorization with Composite Compression Methods

This article presents a fast and approximate multifrontal solver for large sparse linear systems. In a recent work by Liu et al., we showed the efficiency of a multifrontal solver leveraging the butterfly algorithm and its hierarchical matrix extension, HODBF (hierarchical off-diagonal butterfly) compression to compress large frontal matrices. The resulting multifrontal solver can attain quasi-linear computation and memory complexity when applied to sparse linear systems arising from spatial discretization of high-frequency wave equations. To further reduce the overall number of operations and especially the factorization memory usage to scale to larger problem sizes, in this article we develop a composite multifrontal solver that employs the HODBF format for large-sized fronts, a reduced-memory version of the nonhierarchical block low-rank format for medium-sized fronts, and a lossy compression format for small-sized fronts. This allows us to solve sparse linear systems of dimension up to 2.7 × larger than before and leads to a memory consumption that is reduced by 70% while ensuring the same execution time. The code is made publicly available in GitHub.

97 MATHEMATICS AND COMPUTING↗

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

Here, the exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron–ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron–nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The factorization of large composite numbers on the MPP

The continued fraction method for factoring large integers (CFRAC) was an ideal algorithm to be implemented on a massively parallel computer such as the Massively Parallel Processor (MPP). After much effort, the first 60 digit number was factored on the MPP using about 6 1/2 hours of array time. Although this result added about 10 digits to the size number that could be factored using CFRAC on a serial machine, it was already badly beaten by the implementation of Davis and Holdridge on the CRAY-1 using the quadratic sieve, an algorithm which is clearly superior to CFRAC for large numbers. An algorithm is illustrated which is ideally suited to the single instruction multiple data (SIMD) massively parallel architecture and some of the modifications which were needed in order to make the parallel implementation effective and efficient are described.

Mckurdy, Kathy J.↗

Improvement of the AC orbit bump method of local kick factor measurement using reference bunches

Here we have improved the AC orbit bump technique of local kick factor measurement. This technique, developed at NSLS-II a few years ago, is based on local sine-wave excitation of the electron beam by four fast correctors adjacent to the impedance location. By implementing reference bunches and the time gate function of beam position monitors, we are now able to measure the AC orbits of bunches with different intensities simultaneously. This improvement allows us to obtain the impedance-driven orbit distortion in a single measurement, significantly improving the measurement accuracy and reducing the measurement time.

36 MATERIALS SCIENCE↗