Search NASA⌕ Search

SEARCH · Search NASA

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Numerical solution of random singular integral equation appearing in crack problems

The solution of several elasticity problems, and particularly crack problems, can be reduced to the solution of one-dimensional singular integral equations with a Cauchy-type kernel or to a system of uncoupled singular integral equations. Here a method for the numerical solution of random singular integral equations of Cauchy type is presented. The solution technique involves a Chebyshev series approximation, the coefficients of which are the solutions of a system of random linear equations. This method is applied to the problem of periodic array of straight cracks inside an infinite isotropic elastic medium and subjected to a nonuniform pressure distribution along the crack edges. The statistical properties of the random solution are evaluated numerically, and the random solution is used to determine the values of the stress-intensity factors at the crack tips. The error, expressed as the difference between the mean of the random solution and the deterministic solution, is established. Values of stress-intensity factors at the crack tip for different random input functions are presented.

Sambandham, M.↗
DOE OSTI2013

Flow and axial dispersion in a sinusoidal-walled tube: Effects of inertial and unsteady flows

Dispersion in porous media flows has been the subject of much experimental, theoretical and numerical study. Here we consider a wavy-walled tube (a three-dimensional tube with sinusoidally-varying diameter) as a simplified conceptualization of flow in porous media, where constrictions represent pore throats and expansions pore bodies. A theoretical model for effective (macroscopic) longitudinal dispersion in this system has been developed by volume averaging the microscale velocity field. Direct numerical simulation using computational fluid dynamics (CFD) methods was used to compute velocity fields by solving the Navier-Stokes equations, and also to numerically solve the volume averaging closure problem, for a range of Reynolds numbers (Re) spanning the low-Re to inertial flow regimes, including one simulation at Re = 449 for which unsteady flow was observed. Dispersion values were computed using both the volume averaging solution and a random walk particle tracking method, and results of the two methods were shown to be consistent. Our results are compared to experimental measurements of dispersion in porous media and to previous theoretical results for the low-Re, Stokes flow regime. In the steady inertial regime we observe an power-law increase in effective longitudinal dispersion (DL) with Re, consistent with previous results. This rapid rate of increase is caused by trapping of solute in expansions due to flow separation (eddies). For the unsteady case (Re = 449), the rate of increase of DL with Re was smaller than that observed at lower Re. Velocity fluctuations in this regime lead to increased rates of solute mass transfer between the core flow and separated flow regions, thus diminishing the amount of tailing caused by solute trapping in eddies and thereby reducing longitudinal dispersion.

CFD↗
DOE OSTI2025

Randomized Preconditioned Solvers for Strong Constraint 4D-Var Data Assimilation

The Strong Constraint 4D Variational (SC-4DVAR) data assimilation method is widely used in climate and weather applications. SC-4DVAR involves solving a minimization problem to compute the maximum a posteriori estimate, which we tackle using the Gauss-Newton method. The computation of the descent direction is expensive since it involves the solution of a large-scale and potentially ill-conditioned linear system, solved using the preconditioned conjugate gradient (PCG) method. Here, to address this cost, we efficiently construct scalable preconditioners using three different randomization techniques, which all rely on a certain low-rank structure involving the Gauss-Newton Hessian. The proposed techniques come with theoretical guarantees on the condition number, and at the same time, are amenable to parallelization. We also develop an adaptive approach to estimate the sketch size and choose between the reuse or recomputation of the preconditioner. We demonstrate the performance and effectiveness of our methodology on two representative model problems—the Burgers and barotropic vorticity equation—showing a drastic reduction in both the number of PCG iterations and the number of Gauss-Newton Hessian products after including the preconditioner construction cost.

Gauss-Newton↗
DOE OSTI2016

A biorthogonal decomposition for the identification and simulation of non-stationary and non-Gaussian random fields

In this paper, a new method for the identification and simulation of non-Gaussian and non-stationary stochastic fields given a database is proposed. It is based on two successive biorthogonal decompositions aiming at representing spatio–temporal stochastic fields. The proposed double expansion allows to build the model even in the case of large-size problems by separating the time, space and random parts of the field. A Gaussian kernel estimator is used to simulate the high dimensional set of random variables appearing in the decomposition. The capability of the method to reproduce the non-stationary and non-Gaussian features of random phenomena is illustrated by applications to earthquakes (seismic ground motion) and sea states (wave heights).

58 GEOSCIENCES↗
DOE OSTI2013

Solution of stochastic media transport problems using a numerical quadrature-based method

We present a new conceptual framework for analyzing transport problems in random media. We decompose such problems into stratified subproblems according to the number of material pseudo-interfaces within realizations. For a given subproblem we assign pseudo-interface locations in each realization according to product quadrature rules, which allows us to deterministically generate a fixed number of realizations. Quadrature integration of the solutions of these realizations thus approximately solves each subproblem; the weighted superposition of solutions of the subproblems approximately solves the general stochastic media transport problem. We revisit some benchmark problems to determine the accuracy and efficiency of this approach in comparison to randomly generated realizations. We find that this method is very accurate and fast when the number of pseudo-interfaces in a problem is generally low, but that these advantages quickly degrade as the number of pseudo-interfaces increases. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Probing the astrophysical origin of high-energy cosmic-ray electrons with Monte Carlo simulation

High-energy cosmic-ray electrons reveal some remarkable spectral features, the most noteworthy of which is the rise in the positron fraction above 10 GeV. Due to strong energy loss during propagation, these particles can reach Earth only from nearby sources. Yet, the exact nature of these sources, which most likely manifest themselves in the observed anomalies, remains elusive. The many explanations put forward to resolve this case range from standard astrophysics to exotic physics. In this paper, we discuss the possible astrophysical origin of high-energy cosmic-ray electrons through a fully three-dimensional time-dependent Monte Carlo simulation. This approach, which takes advantage of the intrinsic random nature of cosmic-ray diffusive propagation, provides valuable information on the electron-by-electron fluctuations, making it particularly suitable for analyzing in depth the single-source scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Comparing the white dwarf cooling sequences in 47 Tuc and NGC 6397

Using deep Hubble Space Telescope imaging, color-magnitude diagrams are constructed for the globular clusters 47 Tuc and NGC 6397. As expected, because of its lower metal abundance, the main sequence of NGC 6397 lies well to the blue of that of 47 Tuc. A comparison of the white dwarf cooling sequences of the two clusters, however, demonstrates that these sequences are indistinguishable over most of their loci—a consequence of the settling out of heavy elements in the dense white dwarf atmosphere and the near equality of their masses. Lower quality data on M4 continues this trend to a third cluster whose metallicity is intermediate between these two. While the path of the white dwarfs in the color-magnitude diagram is nearly identical in 47 Tuc and NGC 6397, the numbers of white dwarfs along the path are not. This results from the relatively rapid relaxation in NGC 6397 compared to 47 Tuc and provides a cautionary note that simply counting objects in star clusters in random locations as a method of testing stellar evolutionary theory is likely dangerous unless dynamical considerations are included.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Quantification of Organ Motion During Chemoradiotherapy of Rectal Cancer Using Cone-Beam Computed Tomography

Purpose: There has been no previously published data related to the quantification of rectal motion using cone-beam computed tomography (CBCT) during standard conformal long-course chemoradiotherapy. The purpose of the present study was to quantify the interfractional changes in rectal movement and dimensions and rectal and bladder volume using CBCT and to quantify the bony anatomy displacements to calculate the margins required to account for systematic ({Sigma}) and random ({sigma}) setup errors. Methods and Materials: CBCT images were acquired from 16 patients on the first 3 days of treatment and weekly thereafter. The rectum and bladder were outlined on all CBCT images. The interfraction movement was measured using fixed bony landmarks as references to define the rectal location (upper, mid, and low), The maximal rectal diameter at the three rectal locations was also measured. The bony anatomy displacements were quantified, allowing the calculation of systematic ({Sigma}) and random ({sigma}) setup errors. Results: A total of 123 CBCT data sets were analyzed. Analysis of variance for standard deviation from planning scans showed that rectal anterior and lateral wall movement differed significantly by rectal location. Anterior and lateral rectal wall movements were larger in the mid and upper rectum compared with the low rectum. The posterior rectal wall movement did not change significantly with the rectal location. The rectal diameter changed more in the mid and upper than in the low rectum. No consistent relationship was found between the rectal and bladder volume and time, nor was a significant relationship found between the rectal volume and bladder volume. Conclusions: In the present study, the anterior and lateral rectal movement and rectal diameter were found to change most in the upper rectum, followed by the mid rectum, with the smallest changes seen in the low rectum. Asymmetric margins are warranted to ensure phase 2 coverage.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Peridynamic Models for Random Media Found by Coarse Graining

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Designing an Intrusion Detection for an Adjustable Speed Drive System Controlling a Critical Process

In this article, we address the cyber-security problem of industrial control systems (ICSs) when their sensor measurements may be compromised due to an attacker who has intercepted those measurements via a network. We introduce a general-purpose method “Dynamic Watermarking (DW)” to detect potential cyber-intrusions on speed sensor measurements within industrial control systems, which deploy an adjustable speed drive (ASD) to control a critical process. The DW method is injecting a random private low-amplitude signal with a zero mean Gaussian distribution, “watermark”, into one of the input phase voltages powering the ASD system. The watermark signal propagates through the system including pulse width modulation (PWM) power conversion stage and motor, then ultimately appears in the speed sensor measurements. By deploying two statistical DW tests with two proper thresholds, the system can detect potential cyber-intrusions or unobservable cyber-attacks such as replay attacks and false data injection attacks (FDIA). The DW method tested on a laboratory-scale ASD system experimentally to protect the system against cyber-intrusions. This system, powered by a commercial PWM drive operating at 208 V, 3-phase, and 3.7 kW, served as our experimental platform.

42 ENGINEERING↗
DOE OSTI2013

Cubature on Wiener Space: Pathwise Convergence

Cubature on Wiener space (Lyons and Victoir in Proc. R. Soc. Lond. A 460(2041):169-198, 2004) provides a powerful alternative to Monte Carlo simulation for the integration of certain functionals on Wiener space. More specifically, and in the language of mathematical finance, cubature allows for fast computation of European option prices in generic diffusion models.We give a random walk interpretation of cubature and similar (e.g. the Ninomiya-Victoir) weak approximation schemes. By using rough path analysis, we are able to establish weak convergence for general path-dependent option prices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ –1 . Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

An atomistic methodology of energy release rate for graphene at nanoscale

Graphene is a single layer of carbon atoms packed into a honeycomb architecture, serving as a fundamental building block for electric devices. Understanding the fracture mechanism of graphene under various conditions is crucial for tailoring the electrical and mechanical properties of graphene-based devices at atomic scale. Although most of the fracture mechanics concepts, such as stress intensity factors, are not applicable in molecular dynamics simulation, energy release rate still remains to be a feasible and crucial physical quantity to characterize the fracture mechanical property of materials at nanoscale. This work introduces an atomistic simulation methodology, based on the energy release rate, as a tool to unveil the fracture mechanism of graphene at nanoscale. This methodology can be easily extended to any atomistic material system. We have investigated both opening mode and mixed mode at different temperatures. Simulation results show that the critical energy release rate of graphene is independent of initial crack length at low temperature. Graphene with inclined pre-crack possesses higher fracture strength and fracture deformation but smaller critical energy release rate compared with the graphene with vertical pre-crack. Owing to its anisotropy, graphene with armchair chirality always has greater critical energy release rate than graphene with zigzag chirality. The increase of temperature leads to the reduction of fracture strength, fracture deformation, and the critical energy release rate of graphene. Also, higher temperature brings higher randomness of energy release rate of graphene under a variety of predefined crack lengths. The energy release rate is independent of the strain rate as long as the strain rate is small enough.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

Nonlinear analyses of composite aerospace structures in sonic fatigue

The primary research effort of this project is the development of analytical methods for the prediction of nonlinear random response of composite aerospace structures subjected to combined acoustic and thermal loads. The progress, accomplishments, and future plans of three random response research topics are discussed, namely acoustics-structure interactions using boundary/finite element methods, nonlinear vibrations of beams and composite plates under harmonic and random excitations, and numerical simulation of the nonlinear response of composite plates under combined thermal and acoustic loading.

Mei, Chuh↗
DOE OSTI2023

A novel data gaps filling method for solar PV output forecasting

This study proposes a modified gaps filling method, expanding the column mean imputation method and evaluated using randomly generated missing values comprising 5%, 10%, 15%, and 20% of the original data on power output. The XGBoost algorithm was implemented as a forecasting model using the original and processed datasets and two sources of solar radiation data, namely, Shortwave Radiation (SWR) from Advanced Himawari Imager 8 (AHI-8) and Surface Solar Radiation Downward (SSRD) from ERA5 global reanalysis data. Further, the accuracy of the two sets of forecasted power output was evaluated using Root Mean Square Error (RMSE) and Mean Absolute Error (MAE). Results show that by applying the proposed gap filling method and using SWR in forecasting solar photovoltaic (PV) output, the improvement in the RMSE and MAE values range from 12.52% to 24.30% and from 21.10% to 31.31%, respectively. Meanwhile, using SSRD, the improvement in the RMSE values range from 14.01% to 28.54% and MAE values from 22.39% to 35.53%. To further evaluate the accuracy of the proposed gap-filling method, the proposed method could be validated using different datasets and other forecasting methods. Future studies could also consider applying the said method to datasets with data gaps higher than 20%.

Energy & Fuels↗
DOE OSTI2014

SU-E-T-549: A Combinatorial Optimization Approach to Treatment Planning with Non-Uniform Fractions in Intensity Modulated Proton Therapy

Purpose: Non-uniform fractionation, i.e. delivering distinct dose distributions in two subsequent fractions, can potentially improve outcomes by increasing biological dose to the target without increasing dose to healthy tissues. This is possible if both fractions deliver a similar dose to normal tissues (exploit the fractionation effect) but high single fraction doses to subvolumes of the target (hypofractionation). Optimization of such treatment plans can be formulated using biological equivalent dose (BED), but leads to intractable nonconvex optimization problems. We introduce a novel optimization approach to address this challenge. Methods: We first optimize a reference IMPT plan using standard techniques that delivers a homogeneous target dose in both fractions. The method then divides the pencil beams into two sets, which are assigned to either fraction one or fraction two. The total intensity of each pencil beam, and therefore the physical dose, remains unchanged compared to the reference plan. The objectives are to maximize the mean BED in the target and to minimize the mean BED in normal tissues, which is a quadratic function of the pencil beam weights. The optimal reassignment of pencil beams to one of the two fractions is formulated as a binary quadratic optimization problem. A near-optimal solution to this problem can be obtained by convex relaxation and randomized rounding. Results: The method is demonstrated for a large arteriovenous malformation (AVM) case treated in two fractions. The algorithm yields a treatment plan, which delivers a high dose to parts of the AVM in one of the fractions, but similar doses in both fractions to the normal brain tissue adjacent to the AVM. Using the approach, the mean BED in the target was increased by approximately 10% compared to what would have been possible with a uniform reference plan for the same normal tissue mean BED.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2005

Modeling for Ultrasonic Health Monitoring of Foams with Embedded Sensors

In this report analytical and numerical methods are proposed to estimate the effective elastic properties of regular and random open-cell foams. The methods are based on the principle of minimum energy and on structural beam models. The analytical solutions are obtained using symbolic processing software. The microstructure of the random foam is simulated using Voronoi tessellation together with a rate-dependent random close-packing algorithm. The statistics of the geometrical properties of random foams corresponding to different packing fractions have been studied. The effects of the packing fraction on elastic properties of the foams have been investigated by decomposing the compliance into bending and axial compliance components. It is shown that the bending compliance increases and the axial compliance decreases when the packing fraction increases. Keywords: Foam; Elastic properties; Finite element; Randomness

Wang, L.↗
NASA NTRS1973

Investigation of Vibrational Energy Transfer in Connected Structures

The results are reported of an analytical and experimental investigation on the vibrational energy transfer between connected substructures under random excitation. In the analytical area, the basic foundation and assumptions of the statistical energy analysis (SEA) method, a major tool in random response analysis of structures, were examined and reviewed. A new SEA formulation based on the strong coupling condition of the substructures was carried out and presented. Also presented were the results of vibration energy transfer study based on the wave equations applied to connected structures. In the experimental phase, three simple structural models were fabricated and tested. Additional tests were performed on selected substructures which formed parts of the test models. The test results were presented and evaluated against the analytical data.

Hwang, C.↗