Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation 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 91 records · Page 5

A General Framework for Bounding Approximate Dynamic Programming Schemes

For years, there has been interest in approximation methods for solving dynamic programming problems, because of the inherent complexity in computing optimal solutions characterized by Bellman’s principle of optimality. A wide range of approximate dynamic programming (ADP) methods now exists. It is of great interest to guarantee that the performance of an ADP scheme be at least some known fraction, say ß , of optimal. This letter introduces a general approach to bounding the performance of ADP methods, in this sense, in the stochastic setting. The approach is based on new results for bounding greedy solutions in string optimization problems, where one has to choose a string (ordered set) of actions to maximize an objective function. This bounding technique is inspired by submodularity theory, but submodularity is not required for establishing bounds. Instead, the bounding is based on quantifying certain notions of curvature of string functions; the smaller the curvatures the better the bound. The key insight is that any ADP scheme is a greedy scheme for some surrogate string objective function that coincides in its optimal solution and value with those of the original optimal control problem. The ADP scheme then yields to the bounding technique mentioned above, and the curvatures of the surrogate objective determine the value ß of the bound. The surrogate objective and its curvatures depend on the specific ADP.

discrete event systems↗

Calculating Time-Integral Quantities in Depletion Calculations

A method referred to as tally nuclides is presented for accurately and efficiently calculating the time-step averages and integrals of any quantities that are weighted sums of atomic densities with constant weights during the step. The method allows all such quantities to be calculated simultaneously as a part of a single depletion solution with existing depletion algorithms. Some examples of the results that can be extracted include step-average atomic densities and macroscopic reaction rates, the total number of fissions during the step, and the amount of energy released during the step. Furthermore, the method should be applicable with several depletion algorithms, and the integrals or averages should be calculated with an accuracy comparable to that reached by the selected algorithm for end-of-step atomic densities. The accuracy of the method is demonstrated in depletion calculations using the Chebyshev rational approximation method. Here, we demonstrate how the ability to calculate energy release in depletion calculations can be used to determine the accuracy of the normalization in a constant-power burnup calculation during the calculation without a need for a reference solution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A new computational approach for the solutions of generalized pantograph-delay differential equations

In this study, a new computational approach is presented to solve the generalized pantograph-delay differential equations (PDDEs). The solutions obtained by our scheme represent by a linear combination of a special kind of basis functions, and can be deduced in a straightforward manner. Firstly, using the least squares approximation method and the Lagrange-multiplier method, the given PDDE is converted to a linear system of algebraic equations, and those unknown coefficients of the solution of the problem are determined by solving this linear system. Secondly, a PDDE related to the error function of the approximate solution is constructed based on the residual error function technique, and error estimation is presented for the suggested method. The convergence of the approximate solution is proved. Several numerical examples are given to demonstrate the accuracy and efficiency. Comparisons are made between the proposed method and other existing methods.

97 MATHEMATICS AND COMPUTING↗

Towards the Robust Small-Signal Stability Region of Power Systems Under Perturbations Such as Uncertain and Volatile Wind Generation

This paper addresses how to extend the concept of small-signal stability region (SSSR) to that of robust small-signal stability region (RSSSR), within which the system can remain stable even under perturbations caused by uncertain and volatile nodal injections, such as renewable generation. We first employ the structured perturbation theory to formulate the perturbations of nodal injections in state space. Both the intensity and the locations of perturbations can be taken into account. Then, we leverage the stability radius theory and structured singular value theory to define the RSSSR in parameter subspace, enabling a systematic analysis of small-signal stability of power systems under perturbations in a region-wise manner. The hyperplane-approximation method can be employed to construct a linear closed-form approximation of RSSSR boundaries. Case studies on the modified two-area system and New England 39-node system illustrate the new concept of RSSSR and its potential applications.

power system stability analysis↗

A Data-Driven Global Sensitivity Analysis Framework for Three-Phase Distribution System with PVs

Global sensitivity analysis (GSA) of distribution systems with respect to stochastic PV and load variations plays an important role in designing optimal voltage control schemes. This paper proposes a data-driven framework for GSA of distribution systems. In particular, two representative surrogate modeling-based approaches are developed, including the traditional Gaussian process-based and the analysis of variance (ANOVA) kernel ones. The key idea is to develop a surrogate model that captures the hidden global relationship between voltage and real and reactive power injections from the historical data. With the surrogate model, the Sobol indices can be conveniently calculated through either the sampling-based method or the analytical method to assess the global sensitivity of voltage to variations of PV and load power injections. The sampling-based method approximates the Sobol indices using Monte Carlo simulations while the analytical method calculates them by resorting to the ANOVA expansion framework. Comparison results with other model-based GSA methods on the unbalanced three-phase IEEE 37-bus and 123-bus distribution systems show that the proposed framework can achieve much higher computational efficiency with negligible loss of accuracy. The results on a real 240-node distribution system using actual smart meter data further validate the feasibility and scalability of the proposed framework.

14 SOLAR ENERGY↗

Conditional Point Sampling: A Monte Carlo Method for Radiation Transport in Stochastic Media.

Current methods for stochastic media transport are either computationally expensive or, by nature, approximate. Moreover, none of the well-developed, benchmarked approximate methods can compute the variance caused by the stochastic mixing, a quantity especially important to safety calculations. Therefore, we derive and apply a new conditional probability function (CPF) for use in the recently developed stochastic media transport algorithm Conditional Point Sampling (CoPS), which 1) leverages the full intra-particle memory of CoPS to yield errorless computation of stochastic media outputs in 1D, binary, Markovian-mixed media, and 2) leverages the full inter-particle memory of CoPS and the recently developed Embedded Variance Deconvolution method to yield computation of the variance in transport outputs caused by stochastic material mixing. Numerical results demonstrate errorless stochastic media transport as compared to reference benchmark solutions with the new CPF for this class of stochastic mixing as well as the ability to compute the variance caused by the stochastic mixing via CoPS. Using previously derived, non-errorless CPFs, CoPS is further found to be more accurate than the atomic mix approximation, Chord Length Sampling (CLS), and most of memory-enhanced versions of CLS surveyed. In addition, we study the compounding behavior of CPF error as a function of cohort size (where a cohort is a group of histories that share intra-particle memory) and recommend that small cohorts be used when computing the variance in transport outputs caused by stochastic mixing.

61 RADIATION PROTECTION AND DOSIMETRY↗

Numeric Modified Adomian Decomposition Method for Power System Simulations

This paper investigates the applicability of numeric Wazwaz El Sayed modified Adomian Decomposition Method (WES-ADM) for time domain simulation of power systems. WESADM is a numerical method based on a modified Adomian decomposition (ADM) technique. WES-ADM is a numerical approximation method for the solution of nonlinear ordinary differential equations. The non-linear terms in the differential equations are approximated using Adomian polynomials. In this paper WES-ADM is applied to time domain simulations of multimachine power systems. WECC 3-generator, 9-bus system and IEEE 10-generator, 39-bus system have been used to test the applicability of the approach. Several fault scenarios have been tested. It has been found that the proposed approach is faster than the trapezoidal method with comparable accuracy.

Adomian Decomposition↗

A Low-Rank Solver for the Navier--Stokes Equations with Uncertain Viscosity

In this work, we study an iterative low-rank approximation method for the solution of the steady-state stochastic Navier--Stokes equations with uncertain viscosity. The method is based on linearization schemes using Picard and Newton iterations and stochastic finite element discretizations of the linearized problems. For computing the low-rank approximate solution, we adapt the nonlinear iterations to an inexact and low-rank variant, where the solution of the linear system at each nonlinear step is approximated by a quantity of low rank. This is achieved by using a tensor variant of the GMRES method as a solver for the linear systems. We explore the inexact low-rank nonlinear iteration with a set of benchmark problems, using a model of flow over an obstacle, under various configurations characterizing the statistical features of the uncertain viscosity, and we demonstrate its effectiveness by extensive numerical experiments.

97 MATHEMATICS AND COMPUTING↗

A simple method of measuring the effective SRS coefficient in single-mode optical fibres and its applicability limits

This paper describes techniques for measuring the SRS coefficient in a wide spectral range, including the region of small Stokes shifts. A simple, approximate method is proposed for evaluating the SRS coefficient near a gain peak. Spectral dependences of the SRS coefficient are presented for various telecom fibres. (nonlinear optical phenomena)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SolarPILOT: A power tower solar field layout and characterization tool

This paper develops and demonstrates a new SolarPower tower Integrated Layout and Optimization Tool (SolarPILOT). The tool uses the analytical flux image Hermite series approximation originally implemented in the DELSOL3 software developed by Sandia National Laboratory in the early 1980s. By applying the analytical model to individual heliostat images rather than large groups or zones of heliostats, SolarPILOT can characterize a wide variety of heliostat field layouts. The individual heliostat modeling approach increases computational expense in comparison with DELSOL3, so SolarPILOT implements a number of improvements to the analytical approximation method to improve model accuracy and computational efficiency. This paper discusses several of these methods, including dynamic heliostat grouping to reduce the expense of intercept factor evaluation, approximation of annual productivity with a subset of time steps throughout the year, polygon clipping to accurately calculate inter-heliostat shadowing and blocking, receiver and tower geometry optimization, and a trigonometric image transformation technique that ensures model accuracy for small heliostats. SolarPILOT also integrates a Monte-Carlo ray tracing engine (SolTrace), providing improved receiver optical modeling capability, a user-friendly front end for geometry definition, and side-by-side validation of the analytical algorithms.

14 SOLAR ENERGY↗

Exact analytical solution of shear-induced flexural vibration of functionally graded piezoelectric beam

The priority of this paper is to obtain the exact analytical solution for free flexural vibration of FGPM beam actuated using the d{sub 15} effect. In piezoelectric actuators, the potential use of d{sub 15} effect has been of particular interest for engineering applications since shear piezoelectric coefficient d15 is much higher than the other piezoelectric coupling constants d{sub 31} and d{sub 33}. The applications of shear actuators are to induce and control the flexural vibrations of beams and plates. In this study, a modified Timoshenko beam theory is used where electric potential is assumed to vary sinusoidaly along the thickness direction. The material properties are assumed to be graded across the thickness in accordance with power law distribution. Hamilton's principle is employed to obtain the equations of motion along with the associated boundary conditions for FGPM beams. Exact analytical solution is derived thus obtained equations of motion. Results for clamped-clamped and clamped-free boundary conditions are presented. The presented result and method shell serve as benchmark for comparing the results obtained from the other approximate methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficient prediction of time- and angle-resolved photoemission spectroscopy measurements on a nonequilibrium BCS superconductor

We study how time- and angle-resolved photoemission (tr-ARPES) reveals the dynamics of BCS-type, s -wave superconducting systems with time-varying order parameters. Approximate methods are discussed, based on previous approaches to either optical conductivity or quantum dot transport, to enable computationally efficient prediction of photoemission spectra. One use of such predictions is to enable extraction of the underlying order parameter dynamics from experimental data, which is topical given the rapidly growing use of tr-ARPES in studying unconventional superconductivity. The methods considered model the two-time lesser Green's functions with an approximated lesser self-energy that describes relaxation by coupling of the system to two types of baths. The approach primarily used here also takes into consideration the relaxation of the excited states into equilibrium by explicitly including the level-broadening of the retarded and advanced Green's functions. Here, we present equilibrium and nonequilibrium calculations of tr-ARPES spectrum from our model and discuss the signatures of different types of superconducting dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Marginal unbiased score expansion and application to CMB lensing

Here, we present the marginal unbiased score expansion (MUSE) method, an algorithm for generic high-dimensional hierarchical Bayesian inference. MUSE performs approximate marginalization over arbitrary non-Gaussian latent parameter spaces, yielding Gaussianized asymptotically unbiased and near-optimal constraints on global parameters of interest. It is computationally much cheaper than exact alternatives like Hamiltonian Monte Carlo (HMC), excelling on funnel problems which challenge HMC, and does not require any problem-specific user supervision like other approximate methods such as variational inference or many simulation-based inference methods. MUSE makes possible the first joint Bayesian estimation of the delensed Cosmic Microwave Background (CMB) power spectrum and gravitational lensing potential power spectrum, demonstrated here on a simulated data set as large as the upcoming South Pole Telescope 3G 1500 deg 2 survey, corresponding to a latent dimensionality of ~6 million and of order 100 global bandpower parameters. On a subset of the problem where an exact but more expensive HMC solution is feasible, we verify that MUSE yields nearly optimal results. We also demonstrate that existing spectrum-based forecasting tools which ignore pixel-masking underestimate predicted error bars by only ~10%. This method is a promising path forward for fast lensing and delensing analyses which will be necessary for future CMB experiments such as SPT-3G, Simons Observatory, or CMB-S4, and can complement or supersede existing HMC approaches. The success of MUSE on this challenging problem strengthens its case as a generic procedure for a broad class of high-dimensional inference problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Late-time small body disruptions for planetary defense

Here, diverting hazardous small bodies on impact trajectories with the Earth can in some circumstances be impossible without risking disrupting them. Disruption is a much more difficult planetary defense scenario to assess, being linked both to the response of the body to shock loading and the much more complicated gravitational dynamics of the fragments in the solar system relative to pure deflection scenarios. In this work we present a new simulation suite built on N-body gravitational methods that solves fragment orbits in the full gravitational system without recourse to more approximate methods. We assess the accuracy of our simulations and the simplifying assumptions we adopt to make the system tractable, and then discuss in more detail several specific, plausible planetary defense scenarios based on real close encounters. We find that disruption can be a very effective planetary defense strategy even for very late (sub-year) interventions, and should be considered an effective backup strategy should preferred methods, which require long warning times, fail.

42 ENGINEERING↗

Ion energy-angle distribution functions at the plasma-material interface in oblique magnetic fields

The ion energy-angle distribution (IEAD) at the wall of a magnetized plasma is of fundamental importance for the determination of the material processes occurring at the plasma-material interface, comprising secondary emissions and material sputtering. Here, we present a numerical characterization of the IEAD at the wall of a weakly collisional magnetized plasma with the magnetic field inclined at an arbitrary angle with respect to the wall. The analysis has been done using two different techniques: (1) a fluid-Monte Carlo method, and (2) particle-in-cell simulations, the former offering a fast but approximate method for the determination of the IEADs, the latter giving a computationally intensive but self-consistent treatment of the plasma behavior from the quasi-neutral region to the material boundary. The two models predict similar IEADs, whose similarities and differences are discussed. Data are presented for magnetic fields inclined at angles from normal to grazing incidence (0°–85°). We show the scaling factors of the average and peak ion energy and trends of the pitch angle at the wall as a function of the magnetic angle, for use in the correlation of fluid plasma models to material models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Application of Variance Reduction Techniques in RMC Burnup and Thermal-Hydraulic Coupled Calculation

In reactor analog simulation, the Monte Carlo method has relatively high fidelity since it can reduce calculation errors resulting from using approximation methods by describing complex geometry and neutron energy spectra. The Monte Carlo method is, in essence, a stochastic simulation of the neutron transport process, and statistical uncertainty of all the computed results must be provided. In large-scale reactor computation, this method can accurately estimate integral quantities such as the effective multiplication factor (K{sub eff}), and it can limit the uncertainty within an acceptable range via reasonably selecting simulation particles. Nevertheless, it faces great challenges in estimating local quantities like neutron flux and fission power, with uncertainty of estimators beyond acceptance. To decrease the uncertainty of local quantities by increasing simulation particles may cost much time and reduce efficiency. The main reasons are as follows. Firstly, the scale of reactors is so large that common pressurized water reactors have tens of thousands of lattice cells; secondly, uneven distribution of neutron flux and fission power causes uneven distribution of uncertainty. Smith K. puts forward the 95/95 principle to estimate whether uncertainty calculated with Monte Carlo Code satisfies the requirement. According to the principle, with a confidence level of 95%, statistical errors of power density in more than 95% of the regions should be less than 1%. To meet the requirement of the principle, the uniform fission site (UFS) method is raised currently. UFS aims at redistributing fission sources to uniformize the fission source sampling, which can make statistical variance distribution more uniform. This method, which is easy to use, will not obviously increase computation time. Many studies show that UFS has a good variance reduction effect in Monte Carlo criticality calculation, but less attention is paid to its performance in burnup calculation and thermal-hydraulic coupled calculation. Uncertainty arising from the Monte Carlo method has a large impact on burnup calculation and thermal-hydraulic coupled calculation. Variance propagates between burnup steps and between iterative steps of coupled calculation, which causes asymmetry of results of the physical symmetrical regions and instability of coupled iterative convergence errors. To further verify the effect of UFS, this study uses it in the RMC (the Reactor Monte Carlo code) burnup calculation module and the neutronics/thermohydraulics coupled calculation module, and conducts two-dimensional pressurized water reactor burnup calculation and BEAVRS benchmark thermalhydraulic coupled calculation. In this way, it verifies the improvement effect of UFS on variance propagation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Photoionization of Atomic Systems Using the Random-Phase Approximation Including Relativistic Interactions

Approximation methods are unavoidable in solving a many-electron problem. One of the most successful approximations is the random-phase approximation (RPA). Miron Amusia showed that it can be used successfully to describe atomic photoionization processes of many-electron atomic systems. In this article, the historical reasons behind the term “random-phase approximation” are revisited. A brief introduction to the relativistic RPA (RRPA) developed by Walter Johnson and colleagues is provided and some of its illustrative applications are presented.

74 ATOMIC AND MOLECULAR PHYSICS↗

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗