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 487 records · Page 27

Design and Testing of an Automated System using Thermochromatic Liquid Crystals to Determine Local Heat Transfer Coefficients for an Impinging Jet

Using thermochromatic liquid crystal to measure surface temperature, an automated transient method with time-varying free-stream temperature is developed to determine local heat transfer coefficients. By allowing the free-stream temperature to vary with time, the need for complicated mechanical components to achieve a step temperature change is eliminated, and by using the thermochromatic liquid crystals as temperature indicators, the labor intensive task of installing many thermocouples is omitted. Bias associated with human perception of the transition of the thermochromatic liquid crystal is eliminated by using a high speed digital camera and a computer. The method is validated by comparisons with results obtained by the steady-state method for a circular Jet impinging on a flat plate. Several factors affecting the accuracy of the method are evaluated.

Tan, Benjamin↗

On the Effectiveness of Neural Operators at Zero-Shot Weather Downscaling [SWR-25-20]

Code repository for the experiments performed in the paper: On the Effectiveness of Neural Operators at Zero-Shot Weather Downscaling (https://doi.org/10.1017/eds.2025.11) Overall, our work investigates the zero-shot downscaling potential of neural operators. To summarize, our contributions are: 1. We provide a comparative analysis based on two challenging weather downscaling problems, between various neural operator and non-neural-operator methods with large upsampling factors (e.g., 8x and 15x) and fine grid resolutions (e.g., 2 km × 2 km wind speed). 2. We examine whether neural operator layers provide unique advantages when testing downscaling models on upsampling factors higher than those seen during training, i.e., zero-shot downscaling. Our results instead show the surprising success of an approach that combines a powerful transformer-based model with a parameter-free interpolation step at zero-shot weather downscaling. 3. We find that this Swin-Transformer-based approach mostly outperforms all neural operator models in terms of average error metrics, whereas an enhanced super-resolution generative adversarial network (ESRGAN)-based approach is better than most models in capturing the physics of the system, and suggests their use in future work as strong baselines. However, these approaches still do not capture variations at smaller spatial scales well, including the physical characteristics of turbulence in the HR data. This suggests a potential for improvement in transformer or GAN-based methods and neural-operator-based methods for zero-shot weather downscaling.

Sinha, Saumya [National Renewable Energy Laborator↗

A Measurement of the Neutron Electromagnetic Form Factor Ratio from a Rosenbluth Technique with Simultaneous Detection of Neutrons and Protons

The internal structure of protons and neutrons provides insight into both the dynamical behavior of the constitute quarks and gluons, and emergent properties of the nucleons (such as mass, spin, and electromagnetic distributions). Elastic electron-nucleon scattering can probe the elastic electromagnetic form factors of the nucleon. The electric and magnetic form factors, respectively, encode information about the internal charge and magnetization distributions within the nucleon. Precision data for these form factors, over a broad range of the four-momentum transfer squared, Q^2, can benchmark theoretical models describing the strong interaction of nuclear physics. The Super BigBite Spectrometer (SBS) program in Hall A at Jefferson Lab, is a series of high-precision experiments which seek to significantly extend the Q^2 reach of previous data for the nucleon electromagnetic form factors. The first two experiments of this program are known as G_M^n and the neutron Two Photon Exchange (nTPE) and the data were collected from October 2021 to February 2022. Both experiments were conducted with the simultaneous measurement of D(e,e'n) and D(e,e'p) reactions for quasi-elastic electron-deuteron scattering. The scattered electrons were detected in the BigBite Spectrometer, which features multiple large-acceptance Gas Electron Multiplier (GEM) detectors. The Super BigBite Spectrometer provided simultaneous detection of scattered nucleons, and utilized a large acceptance dipole magnet and Hadron Calorimeter (HCal). The G_M^n experiment provides precision measurements of the neutron magnetic form factor, via the ratio method, over a Q^2 range of 3.0 to 13.5 (GeV/c)2. From this data analysis, preliminary values for G_M^n/µ_n G_D are extracted. For Q^2=4.48 (GeV/c)2 we find G_M^n/µ_n G_D=0.9546±0.0132 and for Q^2=4.476 (GeV/c)2 we find G_M^n/µ_n G_D=0.9563±0.0110. These preliminary G_M^n/µ_n G_D values are more precise than existing world data in this Q^2 regime and are consistent with the most recent parameterization of the G_M^n/µ_n G_D world data. The nTPE experiment provides a first measurement of the neutron Rosenbluth Slope and seeks to quantify the two-photon exchange(TPE) contribution to elastic electron-neutron scattering at a fixed Q^2=4.5 (GeV/c)2 with two different beam energies and scattering angle values. For data of the proton form factor ratio, µ_p G_E^p/G_M^p, significant discrepancies exist between values obtained from Rosenbluth Separation and polarization transfer measurement, particularly at large Q^2, and TPE contributions are thought to resolve this discrepancy. The impacts of TPE contributions have not yet been experimentally established for the neutron. From the data analysis presented in this dissertation, a preliminary result for the neutron Rosenbluth Slope is found as S^n=(G_E^n )^2/t_n (G_M^n )^2=0.0916±0.0476 for Q^2=4.48 (GeV/c)2. This value of the neutron Rosenbluth Slope is consistent with the world data extrapolation and the absence of large TPE corrections.

Wertz, Ezekiel [Thomas Jefferson National Accelera↗

Efficient computer algebra algorithms for polynomial matrices in control design

The theory of polynomial matrices plays a key role in the design and analysis of multi-input multi-output control and communications systems using frequency domain methods. Examples include coprime factorizations of transfer functions, cannonical realizations from matrix fraction descriptions, and the transfer function design of feedback compensators. Typically, such problems abstract in a natural way to the need to solve systems of Diophantine equations or systems of linear equations over polynomials. These and other problems involving polynomial matrices can in turn be reduced to polynomial matrix triangularization procedures, a result which is not surprising given the importance of matrix triangularization techniques in numerical linear algebra. Matrices with entries from a field and Gaussian elimination play a fundamental role in understanding the triangularization process. In the case of polynomial matrices, matrices with entries from a ring for which Gaussian elimination is not defined and triangularization is accomplished by what is quite properly called Euclidean elimination. Unfortunately, the numerical stability and sensitivity issues which accompany floating point approaches to Euclidean elimination are not very well understood. New algorithms are presented which circumvent entirely such numerical issues through the use of exact, symbolic methods in computer algebra. The use of such error-free algorithms guarantees that the results are accurate to within the precision of the model data--the best that can be hoped for. Care must be taken in the design of such algorithms due to the phenomenon of intermediate expressions swell.

Baras, J. S.↗

Excavation Volume Growth Factors at Environmental Remediation Sites - 20318

The U.S. Army Corps of Engineers (USACE) manages numerous environmental remediation projects in accordance with the Comprehensive Environmental Response, Compensation, and Liability Act, otherwise known as CERCLA. These include projects for the Formerly Utilized Sites Remedial Action Program (FUSRAP). Utilizing the CERCLA framework to address these projects requires detailed cost estimates and a significant portion of costs are associated with the volume of material to be excavated and that material's final disposition. Accordingly, understanding the expected volumes and uncertainty associated with those volumes is vital to project cost estimating. These estimates are completed at various stages of the CERCLA process, thus the project team has varying degrees of information available from which to develop these estimates depending on how far along the project is. Understanding the factors contributing to volume growth (increase in volume of material from that originally estimated) is a crucial element to determining volume and cost uncertainty. This paper presents a case study of a FUSRAP site which involves excavation and offsite disposal of contaminated soils from complex commercial, industrial, and residential properties. Original volume estimates were compared to post excavation volumes on a property by property basis. Reasons for differences are determined along with magnitude of impacts to cost estimating. For both residential and commercial property types, volume estimating factors, uncertainty considerations, and methods to account for each are discussed with emphasis on volume growth considerations. Cost estimating model input considerations are discussed with specific emphasis on accounting for volume uncertainty at the Feasibility Study (FS) phase of CERCLA. Underestimating or overestimating costs at the FS phase has repercussions under CERCLA that may impact project close out or completion schedule. Discussions presented will be helpful to both reviewers and preparers of cost estimates involving excavations. Additionally, project planners may be able to use the information to plan for volume growth contingencies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Storage and computationally efficient permutations of factorized covariance and square-root information matrices

A unified method to permute vector-stored upper-triangular diagonal factorized covariance (UD) and vector stored upper-triangular square-root information filter (SRIF) arrays is presented. The method involves cyclical permutation of the rows and columns of the arrays and retriangularization with appropriate square-root-free fast Givens rotations or elementary slow Givens reflections. A minimal amount of computation is performed and only one scratch vector of size N is required, where N is the column dimension of the arrays. To make the method efficient for large SRIF arrays on a virtual memory machine, three additional scratch vectors each of size N are used to avoid expensive paging faults. The method discussed is compared with the methods and routines of Bierman's Estimation Subroutine Library (ESL).

Muellerschoen, R. J.↗

Stress Analysis of Columns and Beam Columns by the Photoelastic Method

Principles of similarity and other factors in the design of models for photoelastic testing are discussed. Some approximate theoretical equations, useful in the analysis of results obtained from photoelastic tests are derived. Examples of the use of photoelastic techniques and the analysis of results as applied to uniform and tapered beam columns, circular rings, and statically indeterminate frames, are given. It is concluded that this method is an effective tool for the analysis of structures in which column action is present, particularly in tapered beam columns, and in statically indeterminate structures in which the distribution of loads in the structures is influenced by bending moments due to axial loads in one or more members.

BEAM-COLUMNS - STRESSES↗

An iterative method for obtaining the optimum lightning location on a spherical surface

A brief introduction to the basic principles of an eigen method used to obtain the optimum source location of lightning is presented. The location of the optimum source is obtained by using multiple direction finders (DF's) on a spherical surface. An improvement of this method, which takes the distance of source-DF's as a constant, is presented. It is pointed out that using a weight factor of signal strength is not the most ideal method because of the inexact inverse signal strength-distance relation and the inaccurate signal amplitude. An iterative calculation method is presented using the distance from the source to the DF as a weight factor. This improved method has higher accuracy and needs only a little more calculation time. Some computer simulations for a 4DF system are presented to show the improvement of location through use of the iterative method.

Chao, Gao↗

Methods for analysis of cracks in three-dimensional solids

Analytical and numerical methods evaluating the stress-intensity factors for three-dimensional cracks in solids are presented, with reference to fatigue failure in aerospace structures. The exact solutions for embedded elliptical and circular cracks in infinite solids, and the approximate methods, including the finite-element, the boundary-integral equation, the line-spring models, and the mixed methods are discussed. Among the mixed methods, the superposition of analytical and finite element methods, the stress-difference, the discretization-error, the alternating, and the finite element-alternating methods are reviewed. Comparison of the stress-intensity factor solutions for some three-dimensional crack configurations showed good agreement. Thus, the choice of a particular method in evaluating the stress-intensity factor is limited only to the availability of resources and computer programs.

Raju, I. S.↗

Addressing Human and Organizational Factors in Nuclear Industry Modernization: An Operationally Focused Approach to Process and Methodology

Utility owners and operators of commercial nuclear power plants in the United States (U.S.) are and will be modernizing their nuclear power plants by performing a digital transformation involving design of an integrated set of systems that together enable a technology centric operating plant. The Plant Modernization Pathway of the U.S. Department of Energy Light Water Reactor Sustainability Program has a strategic action plan that lays the groundwork for a digital transformation of the nuclear industry. The model for this transformation is an advanced concept of operations, with an end point vision, “To achieve the maximum aggregate benefit enabled by this digital transformation.” To achieve this, the digital infrastructure for a nuclear plant must be designed as an integrated set of systems that together enable a technology centric operating model. The digital transformation process obviously needs to involve technology considerations and systems engineering, but it also needs to include human and organizational expertise. Thus, human and organizational factors, including sociotechnical systems methods and techniques (e.g., Cognitive Systems Engineering, Systems Theoretic Accident Modeling and Processes, human systems integration, and Macroergonomics) need to be considered for digital transformation projects in order to effectively integrate human and organizational expertise efforts into the new work system that results from nuclear power plant digital modernization. That is, the work system is the basic unit of sociotechnical systems analysis and contains three components: personnel, technical, and organization and management. These components should be jointly optimized with respect to the interdependence of systems performance criteria of effectiveness, efficiency and safety. Joint optimization can be achieved through the application of three human and organization functions: knowledge representation, knowledge elicitation, and cross-functional integration. This report provides a strategic framework for effective integration of human and organizational expertise within nuclear power plant digital modernization efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A machine learning decision criterion for reducing scan time for hyperspectral neutron computed tomography systems

We present the first machine learning-based autonomous hyperspectral neutron computed tomography experiment performed at the Spallation Neutron Source. Hyperspectral neutron computed tomography allows the characterization of samples by enabling the reconstruction of crystallographic information and elemental/isotopic composition of objects relevant to materials science. High quality reconstructions using traditional algorithms such as the filtered back projection require a high signal-to-noise ratio across a wide wavelength range combined with a large number of projections. This results in scan times of several days to acquire hundreds of hyperspectral projections, during which end users have minimal feedback. To address these challenges, a golden ratio scanning protocol combined with model-based image reconstruction algorithms have been proposed. This novel approach enables high quality real-time reconstructions from streaming experimental data, thus providing feedback to users, while requiring fewer yet a fixed number of projections compared to the filtered back projection method. In this paper, we propose a novel machine learning criterion that can terminate a streaming neutron tomography scan once sufficient information is obtained based on the current set of measurements. Our decision criterion uses a quality score which combines a reference-free image quality metric computed using a pre-trained deep neural network with a metric that measures differences between consecutive reconstructions. The results show that our method can reduce the measurement time by approximately a factor of five compared to a baseline method based on filtered back projection for the samples we studied while automatically terminating the scans.

97 MATHEMATICS AND COMPUTING↗

Cryogenic single-port calibration for superconducting microwave resonator measurements

Superconducting circuit testing and materials loss characterization requires robust and reliable methods for the extraction of internal and coupling quality factors of microwave resonators. Additionally, a common method, imposed by limitations on the device design or experimental configuration, is the single-port reflection geometry, i.e. reflection-mode. However, impedance mismatches in cryogenic systems must be accounted for through calibration of the measurement chain while it is at low temperatures. In this paper, we demonstrate a data-based, single-port calibration using commercial microwave standards and a vector network analyzer with samples at millikelvin temperature in a dilution refrigerator, making this method useful for measurements of quantum phenomena. Finally, we cross reference our data-based, single-port calibration and reflection measurement with over-coupled 2D- and 3D-resonators against well established two-port techniques corroborating the validity of our method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient phase-factor evaluation in quantum signal processing

Quantum signal processing (QSP) is a powerful quantum algorithm to exactly implement matrix polynomials on quantum computers. Asymptotic analysis of quantum algorithms based on QSP has shown that asymptotically optimal results can in principle be obtained for a range of tasks, such as Hamiltonian simulation and the quantum linear system problem. A further benefit of QSP is that it uses a minimal number of ancilla qubits, which facilitates its implementation on near-to-intermediate term quantum architectures. However, there is so far no classically stable algorithm allowing computation of the phase factors that are needed to build QSP circuits. Existing methods require the use of variable precision arithmetic and can only be applied to polynomials of a relatively low degree. We present here an optimization-based method that can accurately compute the phase factors using standard double precision arithmetic operations. We demonstrate the performance of this approach with applications to Hamiltonian simulation, eigenvalue filtering, and quantum linear system problems. Furthermore, our numerical results show that the optimization algorithm can find phase factors to accurately approximate polynomials of a degree larger than 10000 with errors below 10 -12 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hanford Site Composite Analysis Data Package: Exposure Scenarios and Radionuclide Specific Dose Conversion Factors

This data package summarizes the exposure assumptions, equations, and methods used to calculate radionuclide-specific unit dose factors and the radiological doses for both groundwater and atmospheric pathways as a part of the updated Hanford Site Composite Analysis (CA). An All Pathways Representative Person exposure scenario is considered to evaluate exposure via both groundwater and atmospheric transport pathways. The radiological dose assessments for both groundwater and atmospheric pathways are included in the performance assessments (PAs) for various Waste Management Areas at the Hanford Site. This data package calculates exposure route-specific and total unit dose factors for CA-specific radionuclides of concern based on the exposure assumptions used in the CA and PAs. This data package also presents the results and comparison of the radionuclide-specific unit dose factors based on the exposure assumptions used in the revised Hanford Site Composite Analysis and Hanford Site performance assessments for low-level waste disposal facilities.

61 RADIATION PROTECTION AND DOSIMETRY↗

Hanford Site Composite Analysis Data Package: Exposure Scenarios and Radionuclide Specific Dose Conversion Factors.

This data package summarizes the exposure assumptions, equations, and methods used to calculate radionuclide-specific unit dose factors and the radiological doses for both groundwater and atmospheric pathways as a part of the revised Hanford Site Composite Analysis. An All-Pathways Representative Person exposure scenario is considered to evaluate exposure via both groundwater and atmospheric transport pathways. The radiological dose assessments for both groundwater and atmospheric pathways are included in the performance assessments for various Waste Management Areas at the Hanford Site. This data package calculates exposure route-specific and total unit dose factors for composite-analysis-specific radionuclides of concern based on the exposure assumptions used in the composite analysis and performance assessments. This data package presents the results and comparison of the radionuclide-specific unit dose factors based on the exposure assumptions used in the revised composite analysis and various performance assessments.

61 RADIATION PROTECTION AND DOSIMETRY↗

Vibration Analysis Reduces Computer Time

Improved calculation method promises to reduce computer time by factor of 10 for vibration analysis of complex structures. Method originally developed to predict mechanical response of spacecraft to maneuvers in various cases in which same booster is used with different payloads. Method is more exact than conventional normal-mode methods; it involves no initial approximations of assumptions. Calculation cost is reduced without loss of accuracy.

Engels, R. C.↗

Three-dimensional finite-element analysis of chevron-notched fracture specimens

Stress-intensity factors and load-line displacements were calculated for chevron-notched bar and rod fracture specimens using a three-dimensional finite-element analysis. Both specimens were subjected to simulated wedge loading (either uniform applied displacement or uniform applied load). The chevron-notch sides and crack front were assumed to be straight. Crack-length-to-specimen width ratios (a/w) ranged from 0.4 to 0.7. The width-to-thickness ratio (w/B) was 1.45 or 2. The bar specimens had a height-to-width ratio of 0.435 or 0.5. Finite-element models were composed of singularity elements around the crack front and 8-noded isoparametric elements elsewhere. The models had about 11,000 degrees of freedom. Stress-intensity factors were calculated by using a nodal-force method for distribution along the crack front and by using a compliance method for average values. The stress intensity factors and load-line displacements are presented and compared with experimental solutions from the literature. The stress intensity factors and load-line displacements were about 2.5 and 5 percent lower than the reported experimental values, respectively.

Raju, I. S.↗

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗