Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Randomized 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 1,333 records · Page 74

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Disorder Enhanced Thermalization in Interacting Many-Particle System

We introduce an extension of the non-equilibrium dynamical mean field theory to incorporate the effects of static random disorder in the dynamics of a many-particle system by integrating out different disorder configurations resulting in an effective time-dependent density-density interaction. We use this method to study the non-equilibrium transient dynamics of a system described by the Fermi Anderson-Hubbard model following an interaction and disorder quench. The method recovers the solution of the disorder-free case for which the system exhibits qualitatively distinct dynamical behaviors in the weak-coupling (prethermalization) and strong-coupling regimes (collapse-and-revival oscillations). However, we find that weak random disorder promotes thermalization. In the weak coupling regime, the jump in the quasiparticle weight in the prethermal regime is suppressed by random disorder while in the strong-coupling regime, random disorder reduces the amplitude of the quasiparticle weight oscillations. These results highlight the importance of disorder in the dynamics of realistic many-particle systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A comparison of acoustic and random vibration testing of spacecraft

Before selecting the most suitable testing technique in a particular case, it has to be considered which of the two methods duplicates better the distribution of vibration responses to be expected during the actual launch phase. An investigation is conducted of the case of a large, low density spacecraft launched by a vehicle which produces significant acoustic noise and transmits relatively little vibration through the mechanical path provided by the vehicle structure in the area of the spacecraft interface. Criteria for vibration-acoustic equivalence are discussed along with tests conducted with the Synchronous Meteorological Satellite, the Radio Astronomy Explorer, and the Orbiting Geophysical Spacecraft. The investigation indicates that the vibration test is generally a poor substitute for an acoustic study.-

Bangs, W. F.↗

Response of nonlinear panels to random loads

Lightweight aircraft structures exposed to a high intensity noise environment can fatigue prematurely if adequate consideration is not given to the problem. Design methods and design criteria for sonic fatigue prevention were developed based on analytical and experimental techniques. Most of the analytical work was based upon small deflection or linear structural theory which did not agree with the experimental results. A large deflection geometrical nonlinearity was incorporated into the analysis methods for determining the structural response to high intensity noise. The Karman-Herrmann large deflection equations with a single mode Galerkin approximation, and the method of equivalent linearization were used to predict mean square amplitude, mean square stresses, and nonlinear frequency at various acoustic loadings for rectangular panels. Both simply supported and clamped support conditions with immovable or movable inplane edges are considered. Comparisons with experimental results are presented.

Mei, C.↗

Optimum seismic structural design based on random vibration and fuzzy graded damages

This paper presents the fuzzy dynamical reliability and failure probability as well as the basic principles and the analytical method of loss assessment for nonlinear seismic steel structures. Also presented is the optimization formulation and a numerical example for double objectives, initial construction cost and expected failure loss, and dynamical reliability constraints. The earthquake ground motion is based on a stationary filtered non-white noise and the fuzzy damage grade is described by damage index.

Cheng, Franklin Y.↗

Calculating Launch Vehicle Flight Performance Reserve

This paper addresses different methods for determining the amount of extra propellant (flight performance reserve or FPR) that is necessary to reach orbit with a high probability of success. One approach involves assuming that the various influential parameters are independent and that the result behaves as a Gaussian. Alternatively, probabilistic models may be used to determine the vehicle and environmental models that will be available (estimated) for a launch day go/no go decision. High-fidelity closed-loop Monte Carlo simulation determines the amount of propellant used with each random combination of parameters that are still unknown at the time of launch. Using the results of the Monte Carlo simulation, several methods were used to calculate the FPR. The final chosen solution involves determining distributions for the pertinent outputs and running a separate Monte Carlo simulation to obtain a best estimate of the required FPR. This result differs from the result obtained using the other methods sufficiently that the higher fidelity is warranted.

Hanson, John M.↗

Predicting Large‐Scale Systematic Missing Pipe Attributes in Water Distribution Networks

Water distribution network (WDN) models are an essential tool used by water utilities for hydraulic analysis. Unfortunately, missing data and insufficient resources often make creating and maintaining these models unfeasible. Existing methods to address missing pipe properties, like sequential imputation for missing values and reconstruction using graph metrics, are designed to accommodate random patterns of missing information and require a significant percentage of the system's attributes to be known. However, these data completeness assumptions do not always align with real‐world scenarios where large sections of the WDN model have missing data. To address this challenge, this study proposes a data‐driven approach for estimating pipe diameter when considering different spatial patterns and degrees of data completeness (i.e., 0%–90%). Using data from 16 WDNs in Kentucky, this study compares the use of machine learning (ML) using topological and geospatial features against an existing deterministic approach. Results demonstrate that WDN models with pipe diameters predicted by the proposed ML method had comparable hydraulic performance to the ground truth models. Moreover, results showed that ML method performance varies between WDNs of differing topological classification. Insights from this study help advance the ability to leverage partial data to create and maintain WDN models amid uncertainty and inadequate resources.

Poff, Jason W. [Oregon State Univ., Corvallis, OR ↗

A new response surface approach for structural reliability analysis

This paper describes a new approach for computing structural reliability by post-processing previously computed probabilistic results for stress and strength. The objective is to provide an accurate method whereby independent probabilistic analyses for stress and strength functions can be performed independently and combined at a later time to compute probability of failure. The method provides a capability for testing different strength measures without the need for re-computing the probabilistic stress response. The proposed approach takes full account of the basic random variables effecting both stress and strength, and the failure region in the variable space identified during separate stress/strength probabilistic analyses. A simple closed-form example and a more complex analysis of a turbine blade subject to creep rupture is used to illustrate the method.

Thacker, B. H.↗

Reliability Sensitivity Analysis and Design Optimization of Composite Structures Based on Response Surface Methodology

This report discusses the development and application of two alternative strategies in the form of global and sequential local response surface (RS) techniques for the solution of reliability-based optimization (RBO) problems. The problem of a thin-walled composite circular cylinder under axial buckling instability is used as a demonstrative example. In this case, the global technique uses a single second-order RS model to estimate the axial buckling load over the entire feasible design space (FDS) whereas the local technique uses multiple first-order RS models with each applied to a small subregion of FDS. Alternative methods for the calculation of unknown coefficients in each RS model are explored prior to the solution of the optimization problem. The example RBO problem is formulated as a function of 23 uncorrelated random variables that include material properties, thickness and orientation angle of each ply, cylinder diameter and length, as well as the applied load. The mean values of the 8 ply thicknesses are treated as independent design variables. While the coefficients of variation of all random variables are held fixed, the standard deviations of ply thicknesses can vary during the optimization process as a result of changes in the design variables. The structural reliability analysis is based on the first-order reliability method with reliability index treated as the design constraint. In addition to the probabilistic sensitivity analysis of reliability index, the results of the RBO problem are presented for different combinations of cylinder length and diameter and laminate ply patterns. The two strategies are found to produce similar results in terms of accuracy with the sequential local RS technique having a considerably better computational efficiency.

Rais-Rohani, Masoud↗

Approach for Input Uncertainty Propagation and Robust Design in CFD Using Sensitivity Derivatives

An implementation of the approximate statistical moment method for uncertainty propagation and robust optimization for quasi 3-D Euler CFD code is presented. Given uncertainties in statistically independent, random, normally distributed input variables, first- and second-order statistical moment procedures are performed to approximate the uncertainty in the CFD output. Efficient calculation of both first- and second-order sensitivity derivatives is required. In order to assess the validity of the approximations, these moments are compared with statistical moments generated through Monte Carlo simulations. The uncertainties in the CFD input variables are also incorporated into a robust optimization procedure. For this optimization, statistical moments involving first-order sensitivity derivatives appear in the objective function and system constraints. Second-order sensitivity derivatives are used in a gradient-based search to successfully execute a robust optimization. The approximate methods used throughout the analyses are found to be valid when considering robustness about input parameter mean values.

Putko, Michele M.↗

Standard random number generation for MBASIC

A machine-independent algorithm is presented and analyzed for generating pseudorandom numbers suitable for the standard MBASIC system. The algorithm used is the polynomial congruential or linear recurrence modulo 2 method. Numbers, formed as nonoverlapping adjacent 28-bit words taken from the bit stream produced by the formula a sub m + 532 = a sub m + 37 + a sub m (modulo 2), do not repeat within the projected age of the solar system, show no ensemble correlation, exhibit uniform distribution of adjacent numbers up to 19 dimensions, and do not deviate from random runs-up and runs-down behavior.

Tausworthe, R. C.↗

Statistical analysis of secondary particle distributions in relativistic nucleus-nucleus collisions

The use is described of several statistical techniques to characterize structure in the angular distributions of secondary particles from nucleus-nucleus collisions in the energy range 24 to 61 GeV/nucleon. The objective of this work was to determine whether there are correlations between emitted particle intensity and angle that may be used to support the existence of the quark gluon plasma. The techniques include chi-square null hypothesis tests, the method of discrete Fourier transform analysis, and fluctuation analysis. We have also used the method of composite unit vectors to test for azimuthal asymmetry in a data set of 63 JACEE-3 events. Each method is presented in a manner that provides the reader with some practical detail regarding its application. Of those events with relatively high statistics, Fe approaches 0 at 55 GeV/nucleon was found to possess an azimuthal distribution with a highly non-random structure. No evidence of non-statistical fluctuations was found in the pseudo-rapidity distributions of the events studied. It is seen that the most effective application of these methods relies upon the availability of many events or single events that possess very high multiplicities.

Mcguire, Stephen C.↗

Data‐Driven Insights into Rare Earth Mineralization: Machine Learning Applications Using Functional Material Synthesis Data

Understanding rare‐earth element (REE) mineralization mechanisms is essential for developing efficient separation strategies. Although the geochemical pathways that generate REE deposits are qualitatively known, quantitative links between specific conditions and mineralization outcomes remain limited. Herein, the repurpose laboratory REE hydrothermal synthesis data—originally collected for functional‐materials fabrication—as a surrogate for studying mineralization with data‐driven methods. The compiled 1,200+ hydrothermal reaction records and trained three machine‐learning models—K‐nearest neighbors (KNN), random forest (RF), and extreme gradient boosting (XGB)—to predict product elements and phases from precursors, additives, reaction conditions, and engineered features. Validation shows XGB achieves the highest accuracy. Feature importance indicates thermodynamic properties of cations and anions dominate model decisions. Correlations reveal positive relationships among precursor concentration, reaction time, pH, and temperature, consistent with classical crystallization behavior. XGB‐based regressors are built to predict crystallization temperature and pH from precursor/product attributes. Performance is strongest when similar training examples exist, while accuracy declines for underrepresented reactions, notably REE carbonates and heavy‐REE systems. Overall, the study shows that functional‐materials datasets can illuminate REE mineralization and provide priors for exploration and processing. Expanding datasets with less‐studied chemistries and conditions will improve generality and support deposit discovery and more efficient REE recovery.

feature importance analysis↗

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗

Long-term thermal stability and calibration of Type-II fiber Bragg grating array inscribed in radiation-hardened fibers

This paper investigates the long-term thermal stability of Type-II fiber Bragg grating (FBG) arrays, inscribed by femtosecond laser in radiation-hardened fiber, for potential applications as multiplexed sensors in high-temperature energy systems. The thermal stability of FBG sensors was assessed through 16 thermal cycles from room temperature (RT) to 750 ℃ about two months, involving 100 FBG sensors. The results show that the absolute temperature drift of FBG sensors can be reduced to less than 0.4 pm/day after 54 h thermal annealing process at a constant temperature of 800 ℃. As temperature sensors, the FBGs demonstrated stable performance, achieving a standard deviation (STD) of 1.8 pm (corresponding to a temperature resolution of 0.118 ℃) post-annealing. Repeated thermal cycles revealed a random drift of 2.3 pm in the FBG wavelength at RT. Polynomial fitting was explored as a calibration method to convert FBG wavelength shifts into absolute temperature measurements. By optimizing calibration temperature points (RT, 200 ℃, 400 ℃, and 750 ℃), the study shows that cubic polynomial calibration using four points yields an average R2 of 0.9997 and an RMSE of 3.58 ℃ across the entire temperature range (RT to 750 ℃). This approach represents an 11.53-fold improvement over empirical slope calibration and a 1.34-fold improvement over four-point piecewise fitting. The findings indicate that Type-II FBGs inscribed in radiation-hardened fibers can function as accurate temperature sensors, with performance on par with or exceeding that of thermocouples. With their multiplexing capability, robust signal transmission over long lead cables, and immunity to electromagnetic interference, FBG sensors offer a promising alternative to traditional electronic sensors for energy system monitoring.

Dominguez-Ontiveros, Elvis [ORNL] (ORCID:000000018↗

Sensitivity and uncertainty of the IFR-1 BISON benchmark

The fuel performance code BISON is being used to evaluate metallic fuel for a new fast-spectrum test reactor called the Versatile Test Reactor, which is being considered for adoption by the US Department of Energy. To quantify the accuracy of BISON predictions, researchers at Oak Ridge National Laboratory have been developing a series of benchmarks based on legacy metallic fuel experiments. As part of this effort, the sensitivity of BISON predictions to variations in model inputs and the uncertainties associated with BISON predictions must be established. This paper summarizes efforts to perform a comprehensive sensitivity analysis (SA) and uncertainty quantification (UQ) on a benchmark based on the IFR-1 experiment.For the SA, at least one input was chosen from every BISON model and physics module used in the benchmark. The inputs were varied individually in a series of BISON simulations. Here, the resulting variations in benchmark predictions were normalized to calculate sensitivities. These sensitivities were then used to inform input selections for the UQ.The UQ was performed using the Monte Carlo UQ method. A literature review was conducted to estimate uncertainty distributions for the selected inputs, and values were sampled randomly from each distribution in a series of BISON simulations. Variations in the benchmark predictions were used to estimate uncertainty distributions and confidence intervals. It was found that nearly 100% of the benchmark predictions matched the corresponding legacy values within the confidence intervals. However, this is at least partially because the confidence intervals associated with benchmark predictions were wide. The uncertainty contributions of assumptions in the benchmark, experimental uncertainties, and BISON models were quantified. Some analysis was performed to identify inputs that contributed to the uncertainties. Finally, recommendations are made for future benchmark and future BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Operator Relaxation and the Optimal Depth of Classical Shadows

Classical shadows are a powerful method for learning many properties of quantum states in a sample-efficient manner, by making use of randomized measurements. Here we study the sample complexity of learning the expectation value of Pauli operators via “shallow shadows,” a recently proposed version of classical shadows in which the randomization step is effected by a local unitary circuit of variable depth t. Here we show that the shadow norm (the quantity controlling the sample complexity) is expressed in terms of properties of the Heisenberg time evolution of operators under the randomizing (“twirling”) circuit—namely the evolution of the weight distribution characterizing the number of sites on which an operator acts nontrivially. For spatially contiguous Pauli operators of weight k, this entails a competition between two processes: operator spreading (whereby the support of an operator grows over time, increasing its weight) and operator relaxation (whereby the bulk of the operator develops an equilibrium density of identity operators, decreasing its weight). From this simple picture we derive (i) an upper bound on the shadow norm which, for depth t~log⁡(k), guarantees an exponential gain in sample complexity over the t=0 protocol in any spatial dimension, and (ii) quantitative results in one dimension within a mean-field approximation, including a universal subleading correction to the optimal depth, found to be in excellent agreement with infinite matrix product state numerical simulations. Our Letter connects fundamental ideas in quantum many-body dynamics to applications in quantum information science, and paves the way to highly optimized protocols for learning different properties of quantum states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗