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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Estimating electric field enhancement factors on an aircraft utilizing a small scale model: A method evaluation

A method for obtaining field enhancement factors at specific points on an aircraft utilizing a small scale model was evaluated by measuring several canonical shapes. Comparison of the form factors obtained by analytical means with measurements indicate that the experimental method has serious flaws. Errors of 200 to 300 percent were found between analytical values and measured values. As a result of the study, the analytical method is not recommended for calibration of field meters located on aircraft, and should not be relied upon in any application where the local spatial derivatives of the electric field on the model are large over the dimensions of the sensing probe.

Easterbrook, Calvin C.↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
NASA NTRS1978

Calculating parts factors for redundant systems

Method that is easily programmed simplifies calculation of parts factor. Individual module unreliabilities are computed as function of number of service intervals and service interval length. At each service interval, unreliability is sum of unreliabilities of replaced and original modules. It must be calculated for each module to obtain parts factor.

Derocher, W. L., Jr.↗
DOE OSTI2020

Embedding via the Exact Factorization Approach

We present a quantum electronic embedding method derived from the exact factorization approach to calculate static properties of a many-electron system. The method is exact in principle but the practical power lies in utilizing input from a low-level calculation on the entire system in a high-level method computed on a small fragment, as in other embedding methods. Here, the exact factorization approach defines an embedding Hamiltonian on the fragment. Further, various Hubbard models demonstrate that remarkably accurate ground-state energies are obtained over the full range of weak to strongly correlated systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Incorporating Physical Priors into Weakly Supervised Anomaly Detection

We propose a new machine-learning-based anomaly detection strategy for comparing data with a background-only reference (a form of weak supervision). The sensitivity of previous strategies degrades significantly when the signal is too rare or there are many unhelpful features. Our prior-assisted weak supervision (PAWS) method incorporates information from a class of signal models to significantly enhance the search sensitivity of weakly supervised approaches. As long as the true signal is in the prespecified class, PAWS matches the sensitivity of a dedicated, fully supervised method without specifying the exact parameters ahead of time. On the benchmark LHC Olympics anomaly detection dataset, our mix of semisupervised and weakly supervised learning is able to extend the sensitivity over previous methods by a factor of 10 in cross section. Furthermore, if we add irrelevant (noise) dimensions to the inputs, classical methods degrade by another factor of 10 in cross section while PAWS remains insensitive to noise. This new approach could be applied in a number of scenarios and pushes the frontier of sensitivity between completely model-agnostic approaches and fully model-specific searches.

artificial neural networks↗
NASA NTRS1987

Stress-intensity factor calculations using the boundary force method

The Boundary Force Method (BFM) was formulated for the three fundamental problems of elasticity: the stress boundary value problem, the displacement boundary value problem, and the mixed boundary value problem. Because the BFM is a form of an indirect boundary element method, only the boundaries of the region of interest are modeled. The elasticity solution for the stress distribution due to concentrated forces and a moment applied at an arbitrary point in a cracked infinite plate is used as the fundamental solution. Thus, unlike other boundary element methods, here the crack face need not be modeled as part of the boundary. The formulation of the BFM is described and the accuracy of the method is established by analyzing a center-cracked specimen subjected to mixed boundary conditions and a three-hole cracked configuration subjected to traction boundary conditions. The results obtained are in good agreement with accepted numerical solutions. The method is then used to generate stress-intensity solutions for two common cracked configurations: an edge crack emanating from a semi-elliptical notch, and an edge crack emanating from a V-notch. The BFM is a versatile technique that can be used to obtain very accurate stress intensity factors for complex crack configurations subjected to stress, displacement, or mixed boundary conditions. The method requires a minimal amount of modeling effort.

Tan, P. W.↗
NASA NTRS1974

Time-fixed rendezvous by impulse factoring with an intermediate timing constraint

A method is presented for factoring a two-impulse orbital transfer into a three- or four-impulse transfer which solves the rendezvous problem and satisfies an intermediate timing constraint. Both the time of rendezvous and the intermediate time of a alinement are formulated as any element of a finite sequence of times. These times are integer multiples of a constant plus an additive constant. The rendezvous condition is an equality constraint, whereas the intermediate alinement is an inequality constraint. The two timing constraints are satisfied by factoring the impulses into collinear parts that vectorially sum to the original impulse and by varying the resultant period differences and the number of revolutions in each orbit. Five different types of solutions arise by considering factoring either or both of the two impulses into two or three parts with a limit for four total impulses. The impulse-factoring technique may be applied to any two-impulse transfer which has distinct orbital periods.

Green, R. N.↗
NASA NTRS1974

Modification of an impulse-factoring orbital transfer technique to account for orbit determination and maneuver execution errors

A method has previously been developed to satisfy terminal rendezvous and intermediate timing constraints for planetary missions involving orbital operations. The method uses impulse factoring in which a two-impulse transfer is divided into three or four impulses which add one or two intermediate orbits. The periods of the intermediate orbits and the number of revolutions in each orbit are varied to satisfy timing constraints. Techniques are developed to retarget the orbital transfer in the presence of orbit-determination and maneuver-execution errors. Sample results indicate that the nominal transfer can be retargeted with little change in either the magnitude (Delta V) or location of the individual impulses. Additonally, the total Delta V required for the retargeted transfer is little different from that required for the nominal transfer. A digital computer program developed to implement the techniques is described.

Kibler, J. F.↗
DOE OSTI2023

Viewfactor and Raytracing for AgriPV Modeling

View factor models are used in due diligence software to calculate rear irradiance for bifacial modules. An intermediate step in this calculation is the irradiance at the ground level, which can be leveraged for evaluating the Photo Active Radiation available for crops in Agrivoltaic setups. This paper presents the metrics and modifications to the model for ground irradiance study with the view factor approach, compares it to the raytracing method, and validates it with field measurements of ground irradiance. It is found that for the clearances, row-to-row setups, and tilts studied, the view factor method matches with raytracing results within 2% MBD. The comparison is performed for the nine most common agriPV configurations using high-performance computing and the NSRDB database for the whole US, with results and data made available open-source on the InSPIRE AgriPV website.

AgriPV↗
DOE OSTI2025

Evaluation of data driven low-rank matrix factorization for accelerated solutions of the Vlasov equation

Low-rank methods have shown success in accelerating simulations of a collisionless plasma described by the Vlasov equation, but still rely on computationally costly linear algebra every time step. We propose a data-driven factorization method using artificial neural networks, specifically with convolutional layer architecture, that trains on existing simulation data. At inference time, the model outputs a low-rank decomposition of the distribution field of the charged particles, and we demonstrate that this step is faster than the standard linear algebra technique. Numerical experiments show that the method achieves comparable reconstruction accuracy for interpolation tasks, generalizing to unseen test data in a manner beyond just memorizing training data; patterns in factorization also inherently followed the same numerical trend as those within algebraic methods (e.g., truncated singular-value decomposition). However, when training on the first 70% of a time-series data and testing on the remaining 30%, the method fails to meaningfully extrapolate. Despite this limiting result, the technique may have benefits for simulations in a statistical steady-state or otherwise showing temporal stability. These results suggest that while the model offers a computationally efficient alternative for datasets with temporal stability, its current formulation is best suited for interpolation rather than for predicting future states in time-evolving systems. This study thus lays the groundwork for further refinement of neural network-based approaches to low-rank matrix factorization in high-dimensional plasma simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sullivan Process near Threshold and the Pion Gravitational Form Factors

We propose a novel method to experimentally access the gravitational form factors of the charged pion 𝜋 + through the Sullivan process in electron-proton scattering. We demonstrate that the cross sections of 𝐽/𝜓 photoproduction and 𝜙 electroproduction near the respective thresholds are dominated by the gluon gravitational form factor of the pion to next-to-leading order in perturbative QCD. We predict cross sections for the Electron-Ion Collider and the Jefferson Lab experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1995

Comparison of Implicit Schemes for the Incompressible Navier-Stokes Equations

For a computational flow simulation tool to be useful in a design environment, it must be very robust and efficient. To develop such a tool for incompressible flow applications, a number of different implicit schemes are compared for several two-dimensional flow problems in the current study. The schemes include Point-Jacobi relaxation, Gauss-Seidel line relaxation, incomplete lower-upper decomposition, and the generalized minimum residual method preconditioned with each of the three other schemes. The efficiency of the schemes is measured in terms of the computing time required to obtain a steady-state solution for the laminar flow over a backward-facing step, the flow over a NACA 4412 airfoil, and the flow over a three-element airfoil using overset grids. The flow solver used in the study is the INS2D code that solves the incompressible Navier-Stokes equations using the method of artificial compressibility and upwind differencing of the convective terms. The results show that the generalized minimum residual method preconditioned with the incomplete lower-upper factorization outperforms all other methods by at least a factor of 2.

Rogers, Stuart E.↗
DOE OSTI2020

Analytic model of electron transport through and over non-linear barriers

Tunneling barriers are an essential component of electron sources, sensors, detectors, and vacuum nanoelectronics, and a pivotal factor in their performance, but the barriers themselves routinely depart from the analytic models used to model their behavior. Here, a new formalism is developed to analytically and accurately model emission through and over barriers associated with depletion layers, nanotip barriers, and MIM structures. The transmission probability for depletion layer and metal-insulator-metal (MIM) and metaloxide-semiconductor (MOS) barriers is accurately modeled as the electron energy exceeds the barrier height using approaches designed for rapid implementation demanded by simulation codes, and extensible to general barriers. The models supersede conventional thermal and field models in depletion and MIM/MOS barrier studies. Thermal-field methods are used to treat the transmission probability, and shape factor methods to treat the tunneling factor. Analytic formulae for current density are obtained. The methods ease device simulation and characterization of current-voltage relations for emerging technologically interesting barriers with better accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1973

Estimation of peak winds from hourly observations

Two closely related methods to obtain estimates of the hourly peak wind at Cape Kennedy were compared by statistical tests. The methods evaluated the Monin-Obukhov stability length and the standard deviation of the hourly observed wind speed, so as to augment the latter quantity by F standard deviations. F is an optimized factor. A third method utilizing an optimized gust factor was also applied to the hourly wind. The latter procedure estimated 2952 peak winds with an rms error of 2.81 knots, an accuracy which was not surpassed by the other methods. Peak ground wind speed data were developed for use in space shuttle design operation analyses.

Graves, M. E.↗
DOE OSTI2022

Resonance form factors from finite-volume correlation functions with the external field method

A novel method for the extraction of form factors of unstable particles on the lattice is proposed. The approach is based on the study of two-particle scattering in a static, spatially periodic external field by using a generalization of the Lüscher method in the presence of such a field. It is shown that the resonance form factor is given by the derivative of the resonance pole position in the complex plane with respect to the coupling constant to the external field. Unlike the standard approach, this proposal does not suffer from problems caused by the presence of the triangle diagram.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1979

Alternating direction implicit methods for parabolic equations with a mixed derivative

Alternating direction implicit (ADI) schemes for two-dimensional parabolic equations with a mixed derivative are constructed by using the class of all A sub 0-stable linear two-step methods in conjunction with the method of approximation factorization. The mixed derivative is treated with an explicit two-step method which is compatible with an implicit A sub 0-stable method. The parameter space for which the resulting ADI schemes are second order accurate and unconditionally stable is determined. Some numerical examples are given.

Beam, R. M.↗
NASA NTRS1980

Comparison of methods for prediction of transition by stability analysis

The paper compares several methods of transition prediction of linear stability analysis. The spectral stability analysis code SALLY is used to analyze flows over laminar flow control wings. It is shown that transition prediction by the envelope method and a new modified wave packet method are comparable in reliability, but that the envelope method is more efficient computationally. This is based on the results which show that the wave packet method provides N factors which are at best as consistent as those of the envelope method; since the wave packet method is at least 3 times as expensive as the envelope method, the latter is recommended for engineering design calculations.

Malik, M. R.↗