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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Framework for Robust Long-Term Voltage Stability of Distribution Systems

Power injection uncertainties in distribution power grids, which are mostly induced by aggressive introduction of intermittent renewable sources, may drive the system away from normal operating regimes and potentially lead to the loss of long-term voltage stability (LTVS). Naturally, there is an ever increasing need for a tool for assessing the LTVS of a distribution system. This paper presents a fast and reliable tool for constructing inner approximations of LTVS regions in multidimensional injection space such that every point in our constructed region is guaranteed to be solvable. Numerical simulations demonstrate that our approach outperforms all existing inner approximation methods in most cases. Furthermore, the constructed regions are shown to cover substantial fractions of the true voltage stability region. In conclusion, the paper will later discuss a number of important applications of the proposed technique, including fast screening for viable injection changes, constructing an effective solvability index and rigorously certified loadability limits.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Preliminary Study on TRISO Fuel Cross Section Generation

Cross section self-shielding methodologies for TRISO fuel were assessed to provide accurate multigroup cross sections for a high-fidelity reactor physics code so that the code is able to accurately model and simulate advanced reactors with TRISO fuel. Initially, the two existing methodologies (the SCALE method and the Sanchez-Pomraning method) were studied and implemented to MC2-3 for detailed performance tests. Additionally, a new spatial self-shielding method, named the iterative local spatial self-shielding (ILSS) method, for particulate fuels was developed based on the disadvantage factor and implemented to MC2-3 as well. The new method approximately accounts for the effect of randomly distributed particles on the particle shadowing effect using a homogenized compact region surrounding a particle of interest at the center. The self-shielded cross sections of the particle at the center are determined iteratively since the cross sections of the homogenized compact region are calculated using them. For the energy range above 100 keV where the fuel-to-moderator ratio is more important than the random distribution of particles, a single particle unit-cell model is used by preserving the average amount of moderator per fuel particle in the system. The three self-shielding methods implemented in MC2-3 were tested using numerical benchmark problems made based on fuel compact problems of a prismatic-type very high temperature reactor. Test results indicated that the ILSS method produced slightly better results than the SCALE and Sanchez-Pomraning methods, compared to the Serpent-2 Monte Carlo results obtained with 25 independent random particle configurations. The SCALE and Sanchez-Pomraning methods tend to underestimate the heterogeneity effect by 150 and 100 pcm, respectively, while the new ILSS method overestimates the heterogeneity effect by 70 pcm. In future, the new self-shielding method will be extended to perform pebble calculations and compare results with those from the SCALE and Sanchez-Pomraning methods. Furthermore, the new method will be optimized for practical applications to on-the-fly resonance treatment for lattice or whole-core calculations for advanced reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

WE-AB-209-07: Explicit and Convex Optimization of Plan Quality Metrics in Intensity-Modulated Radiation Therapy Treatment Planning

Purpose: To formulate objective functions of a multicriteria fluence map optimization model that correlate well with plan quality metrics, and to solve this multicriteria model by convex approximation. Methods: In this study, objectives of a multicriteria model are formulated to explicitly either minimize or maximize a dose-at-volume measure. Given the widespread agreement that dose-at-volume levels play important roles in plan quality assessment, these objectives correlate well with plan quality metrics. This is in contrast to the conventional objectives, which are to maximize clinical goal achievement by relating to deviations from given dose-at-volume thresholds: while balancing the new objectives means explicitly balancing dose-at-volume levels, balancing the conventional objectives effectively means balancing deviations. Constituted by the inherently non-convex dose-at-volume measure, the new objectives are approximated by the convex mean-tail-dose measure (CVaR measure), yielding a convex approximation of the multicriteria model. Results: Advantages of using the convex approximation are investigated through juxtaposition with the conventional objectives in a computational study of two patient cases. Clinical goals of each case respectively point out three ROI dose-at-volume measures to be considered for plan quality assessment. This is translated in the convex approximation into minimizing three mean-tail-dose measures. Evaluations of the three ROI dose-at-volume measures on Pareto optimal plans are used to represent plan quality of the Pareto sets. Besides providing increased accuracy in terms of feasibility of solutions, the convex approximation generates Pareto sets with overall improved plan quality. In one case, the Pareto set generated by the convex approximation entirely dominates that generated with the conventional objectives. Conclusion: The initial computational study indicates that the convex approximation outperforms the conventional objectives in aspects of accuracy and plan quality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Spatial Correlations of the Poisson Model for Radiation Transport

Characterizing the relationship between bulk physical properties and mixing in randomly heterogeneous media is a central challenge across many areas of science and engineering. A benchmark model for such studies is the Poisson model, a random tessellation of space by a Poisson process of hyperplanes. In radiation transport studies, the lack of exact expressions for the Poisson model’s spatial multipoint functions has led to approximate methods being used, introducing unquantified sources of error. Here, we recently introduced an exact solution for the Poisson model’s multipoint functions and closely related conditional probability functions (CPFs), providing a new opportunity to understand and reduce these sources of error. In this paper, we enable a more rigorous investigation of radiation transport in stochastic media by applying the recently introduced exact solution for the Poisson model’s CPFs. This paper consists of three main contributions. First, we introduce a unified framework for CPFs of the Poisson model, encompassing the recently introduced exact CPFs as well as the previously introduced atomic mix, nearest-neighbor, and combination CPFs. This framework also includes existing pruning techniques for the approximate CPFs, such as angular exclusion, as well as a novel form of angular exclusion suitable for the exact CPFs. Second, we use the exact CPFs to characterize the spatial regions where each approximate three-point CPF is most accurate, thereby explaining the observed hierarchy of accuracy among the approximate models. Finally, we evaluate material transmittance, reflectance, and flux in a three-dimensional test problem using conditional point sampling, demonstrating the relationship between CPF accuracy and transport simulation accuracy.

Poisson model↗
DOE OSTI2020

27 AlNMR chemical shift of Al(OH)$^{-}_{4}$ calculated from first principles: Assessment of error cancellation in chemically distinct reference and target systems

Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H 2 O)$^{3+}_{6}$, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)$^{-}_{4}$ requires an understanding of how these errors cancel, in order to define the limits of accurately predicting 27 Al chemical shielding in Al(OH)$^{-}_{4}$. In this work, we estimate the absolute shielding of the 27 Al nucleus in Al(OH)$^{-}_{4}$ at the coupled cluster level (515.1± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab-initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H 2 O)$^{3+}_{6}$ and Al(OH)$^{-}_{4}$ ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the 27 Al chemical shift of Al(OH)$^{-}_{4}$ from first principles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Optimal control of wave energy converters

A wave energy converter and method for extracting energy from water waves maximizes the energy extraction per cycle by estimating an excitation force of heave wave motion on the buoy, computing a control force from the estimated excitation force using a dynamic model, and applying the computed control force to the buoy to extract energy from the heave wave motion. Analysis and numerical simulations demonstrate that the optimal control of a heave wave energy converter is, in general, in the form of a bang-singular-bang control; in which the optimal control at a given time can be either in the singular arc mode or in the bang-bang mode. The excitation force and its derivatives at the current time can be obtained through an estimator, for example, using measurements of pressures on the surface of the buoy in addition to measurements of the buoy position. A main advantage of this approximation method is the ease of obtaining accurate measurements for pressure on the buoy surface and for buoy position, compared to wave elevation measurements.

16 TIDAL AND WAVE POWER↗
DOE OSTI2016

Physical optics simulations with PHASE for SwissFEL beamlines

PHASE is a software tool for physical optics simulation based on the stationary phase approximation method. The code is under continuous development since about 20 years and has been used for instance for fundamental studies and ray tracing of various beamlines at the Swiss Light Source. Along with the planning for SwissFEL a new hard X-ray free electron laser under construction, new features have been added to permit practical performance predictions including diffraction effects which emerge with the fully coherent source. We present the application of the package on the example of the ARAMIS 1 beamline at SwissFEL. The X-ray pulse calculated with GENESIS and given as an electrical field distribution has been propagated through the beamline to the sample position. We demonstrate the new features of PHASE like the treatment of measured figure errors, apertures and coatings of the mirrors and the application of Fourier optics propagators for free space propagation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Shape evolution of neutron-rich 106,108,110 Mo isotopes in the triaxial degree of freedom

Background: Neutron-rich nuclei with mass number between 100 and 110 attract much attention, since several kinds of shapes, such as spherical, prolate, oblate, and triaxial shapes, are predicted. In particular, for neutron-rich Mo isotopes, different models predict different magnitudes and rigidity of triaxial deformation. Previous interpretations of experimental results based solely on low-lying $2^+_2$ states are insufficient to distinguish between the rigid triaxial shape, $\gamma$ vibration, or $\gamma$-soft rotor. Purpose: The shape evolution of 106 Mo, 108 Mo, and 110 Mo is investigated through their $2^+_1$-state lifetimes, decay-branching ratios of the $2^+_2$ state, and energies of the low-lying collective excited states with $K^π = 0^+, 2^+$, and $4^+$. Method: $\beta$-delayed $\gamma$-ray spectroscopy was employed for neutron-rich Nb and Zr isotopes produced at the RIKEN RI Beam Factory to populate excited states in 106 Mo, 108 Mo, and 110 Mo. The EUroball-RIKEN Cluster Array was used for high-resolution $\gamma$-ray detection and lifetimes of the $2^+_1$ states were determined using the UK fast-timing array of LaBr 3 (Ce) detectors. Results: New $\gamma$-ray transitions and levels are reported, including newly assigned $0^+_2$ states in 108,110 Mo. Quadrupole deformations were obtained for 106,108,110 Mo from their $2^+_1$ energies and lifetimes. The $\beta$-delayed neutron-emission probabilities of 108 Nb and 110 Nb were determined by examining the $\gamma$ rays of their respective daughter decays. Conclusions: In this work, the even-odd energy staggering in the $2^+_2$ band was compared with typical patterns of the $\gamma$-vibrational band, rigid triaxial rotor, and $\gamma$-soft rotor. The very small even-odd staggering of 106 Mo, 108 Mo, and 110 Mo favors a $\gamma$-vibrational band assignment. The kinematic moment of inertia for the $2^+_2$ band showed a trend similar to the ground-state band, which is as expected for the $\gamma$-vibrational band. Beyond-mean-field calculations employing the constrained Hartree-Fock-Bogoliubov and local quasiparticle-random-phase approximation method using the SLy5 + T interaction reproduced the ground and $2^+_2$ bands in 106 Mo and 108 Mo. The collective wave functions are consistent with the interpretation of the $2^+_2$ band as the $\gamma$-vibrational band of the prolate shape. However, the staggering pattern observed in 110 Mo differs from the one suggested in the calculations which predict a $\gamma$-soft rotor. There was no experimental indication of the oblate shape or the $\gamma$-soft rotor predicted in these Mo isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Electric radiation mapping of silver/zinc oxide nanoantennas by using electron holography

In this work, we report the fabrication of self-assembled zinc oxide nanorods grown on pentagonal faces of silver nanowires by using microwaves irradiation. The nanostructures resemble a hierarchal nanoantenna and were used to study the far and near field electrical metal-semiconductor behavior from the electrical radiation pattern resulting from the phase map reconstruction obtained using off-axis electron holography. As a comparison, we use electric numerical approximations methods for a finite number of ZnO nanorods on the Ag nanowires and show that the electric radiation intensities maps match closely the experimental results obtained with electron holography. The time evolution of the radiation pattern as generated from the nanostructure was recorded under in-situ radio frequency signal stimulation, in which the generated electrical source amplitude and frequency were varied from 0 to 5 V and from 1 to 10 MHz, respectively. The phase maps obtained from electron holography show the change in the distribution of the electric radiation pattern for individual nanoantennas. The mapping of this electrical behavior is of the utmost importance to gain a complete understanding for the metal-semiconductor (Ag/ZnO) heterojunction that will help to show the mechanism through which these receiving/transmitting structures behave at nanoscale level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Uncertainty in the visibility mask of a survey and its effects on the clustering of biased tracers

The forecasted accuracy of upcoming surveys of large-scale structure cannot be achieved without a proper quantification of the error induced by foreground removal (or other systematics like 0-point photometry offset). Because these errors are highly correlated on the sky, their influence is expected to be especially important at very large scales, at and beyond the first Baryonic Acoustic Oscillation (BAO). In this work we quantify how the uncertainty in the visibility mask of a survey, that gives the survey depth in a specific sky area, influences the measured power spectrum of a sample of tracers of the density field and its covariance matrix. We start from a very large set of 10,000 catalogs of dark matter (DM) halos in periodic cosmological boxes, produced with the PINOCCHIO approximate method. To make an analytic approach feasible, we assume luminosity-independent halo bias and an idealized geometry for the visibility mask, that is constant in square tiles of physical length l ; this should be interpreted as the projection, at the observation redshift, of the angular correlation scale of the foreground residuals. We find that the power spectrum of these biased tracers can be expressed as the sum of a cosmological term, a mask term and a term involving their convolution. The mask and convolution terms scale like P ∝ l {sup 2}σ {sub A} {sup 2}, where σ {sub A} {sup 2} is the variance of the uncertainty on the visibility mask. With l = 30−100 Mpc/ h and σ {sub A} = 5−20%, the mask term can be significant at k ∼ 0.01−0.1 h /Mpc, and the convolution term can amount to ∼ 1−10% of the total. The influence of mask uncertainty on power spectrum covariance is more complicated: the coupling of the convolution term with the other two gives rise to several mixed terms, that we quantify by difference using the mock catalogs. These are found to be of the same order of the mask covariance, and to introduce non-diagonal terms at large scales. As a consequence, the power spectrum covariance matrix cannot be expressed as the sum of a cosmological and of a mask term. More realistic settings (realistic foregrounds, luminosity-dependent bias) make the analytical approach not feasible, and the problem requires on the one hand usage of extended sets of mock catalogs, on the other hand detailed knowledge of the correlations among errors in the visibility masks. Our results lie down the theoretical bases to quantify the impact that uncertainties in the mask calibration have on the derivation of cosmological constraints from large spectroscopic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Tailoring the Electronic Structures and Spectral Properties of ZnO with Irradiation Defects Generated Under Intense Electronic Excitation: A Combined Experimental and DFT Approach

Here, the relationship between the composition of the internal defect states, spectral properties, and correlated electronic structures of wurtzite zinc oxide (ZnO) crystals under 645 MeV Xe 35+ irradiation is systematically investigated, employing experimental characterizations combined with first-principle calculations. Based on the ion irradiation-induced thermal expansion and relaxation processes, the high concentration of vacancy/interstitial defects produced from the transient disordered phase in molten track states trigger photoelectric changes, as follows: i) the generation of internal defect states effectively reduces the intrinsic bandgap (3.25 eV → 2.66 eV); ii) a large number of defective active sites inhibits the recombination between electron–hole pairs, causing dark conductance and photoconductance to increase with increasing damage levels until optimal fluence is achieved. Based on the density functional theory (DFT) with the GGA + U (GGA = generalized gradient approximation) method, the defective models associated with the different electronic structures, density of states, formation energy, and the nature of the chemical bonding are established. The narrowing of the bandgap observed experimentally and the enhancement of carrier concentration originating from the internal electron defect states are qualitatively verified, therefore laying the foundation for designing future nanoscale photoelectronic devices and microelectronics applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Real corrections to Higgs boson pair production at NNLO in the large top quark mass limit

In this paper we consider the next-to-next-to-leading order total cross section of Higgs boson pair production in the large top quark mass limit and compute four expansion terms in 1 $ {m}_t^2 $ . To this end, we analytically compute the real-virtual and double-real contributions to the total cross section and combine them with the existing virtual contribution. Good convergence is observed below the top quark threshold, which makes our results a valuable input for approximation methods which aim for next-to-next-to-leading order corrections over the whole kinematic range. We present details on various steps of our calculation; in particular, we provide results for three- and four-particle phase-space master integrals and describe in detail the evaluation of the collinear counterterms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Effects of La Doping and Zn or O Vacancy on the Magnetic Property of ZnO

Research on the magnetic origin and mechanism of La-doped ZnO system has been controversial. To solve this problem, we studied the effects of La doping and Zn or O vacancy on the magnetic property of ZnO using geometry optimization and energy calculation on the basis of first-principle generalized gradient approximation method of density functional theory. The electronic structures and magnetic properties of undoped ZnO, La-doped ZnO system, and systems with coexisting La doping and Zn or O vacancies were calculated. Results showed that the systems of La-doped ZnO and the coexistence of La doping and O vacancy are nonmagnetic. In addition, the system of La replacing Zn and one Zn vacancy exhibits long-range orderly ferromagnetism, and the Curie temperature of the doping system can achieve room temperature. The magnetism source of the systems with La doping and Zn vacancy coexisting in ZnO demonstrates strong hybrid coupling electron exchange effects existing among Zn-4s, Zn-3p, O-2p, and La-5s orbits that are nearest to the Zn vacancy. The next nearest distance between doping and Zn vacancy leads to the lowest formation energy and highest stability under the same doping condition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Constitutive Modeling of the Flow Stress of GCr15 Continuous Casting Bloom in the Heavy Reduction Process

According to the calculation results of a 3D thermomechanical-coupled finite-element (FE) model of GCr15 bearing steel bloom during a heavy reduction (HR) process, the variation ranges in the strain rate and strain under HR were described. In addition, the hot deformation behavior of the GCr15 bearing steel was studied over the temperature range from 1023 K to 1573 K (750 °C to 1300 °C) with strain rates of 0.001, 0.01, and 0.1 s{sup −1} in single-pass thermosimulation compression experiments. To ensure the accuracy of the constitutive model, the temperature range was divided into two temperature intervals according to the fully austenitic temperature of GCr15 steel [1173 K (900 °C)]. Two sets of material parameters for the constitutive model were derived based on the true stress–strain curves of the two temperature intervals. A flow stress constitutive model was established using a revised Arrhenius-type constitutive equation, which considers the relationships among the material parameters and true strain. This equation describes dynamic softening during hot compression processes. Considering the effect of glide and climb on the deformation mechanism, the Arrhenius-type constitutive equation was modified by a physically based approach. This model is the most accurate over the temperatures ranging from 1173 K to 1573 K (900 °C to 1300 °C) under HR deformation conditions (ignoring the range from 1273 K to 1573 K (1000 °C to 1300 °C) with a strain rate of 0.1 s{sup −1}). To ensure the convergence of the FE calculation, an approximated method was used to estimate the flow stress at temperatures greater than 1573 K (1300 °C).

36 MATERIALS SCIENCE↗
DOE OSTI2018

Exact Bayesian computation using H-functions

In the Bayesian Statistics, the use of conjugate prior distributions allows to obtain exact posterior quantities, avoiding the use of approximate methods. However, there are only a few conjugate structures available, which limits the use in practical problems. In this work we propose a general procedure to obtain the posterior distribution in an exact form in non-conjugate location parameter models with known (possibly, heteroskedastic) variance. The theory is based on special functions, specifically the H-function, which embraces most of the probability distributions. We express explicitly in a computable form all the Bayesian tools for posterior inferences, that is the posterior and the posterior predictive distributions and their moments, as well as the cumulative posterior distribution. As an illustration we consider a problem in astronomy of estimating the distance to the Large Magellanic Cloud, the largest galaxy now orbiting the Milky Way Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

ONIX: An open-source depletion code

Open Source software enables innovative, community-based software development. ONIX brings this concept to the field of depletion calculations. It is an open-source depletion software to be used for nuclear reactor simulations, for fissile material production analysis as well as for nuclear arms control applications. ONIX provides a module to solve the depletion equation using a Chebyshev Rational Approximation Method. For the generation of one-group cross sections, it includes a coupling interface for the open-source neutron transport code, OpenMC, as well as a module to read pre-computed values in a stand-alone mode. ONIX has special features to optimize nuclear data libraries, to update isomeric branching ratio during burnup, and to support automation of simulations for nuclear archaeology. In conclusion, ONIX has been validated against results from numerical and experimental benchmarks, and its results agree with other methods within expected error ranges.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Cost-optimal evaluation of centralized and distributed microgrid topologies considering voltage constraints

Optimal design of hybrid renewable mini-grids requires both economic and power quality considerations. Existing modeling approaches address these considerations via separate or loosely coupled models. Here, we extend REopt—a techno-economic optimization model developed at the National Renewable Energy Laboratory—to consider both within a single model. REopt formulates the design problem as a mixed-integer linear program that solves for a site's optimal technology mix, sizing, and operation to minimize life cycle cost. REopt has traditionally assumed a single node system. In the work presented here, we expand the REopt platform to consider multiple connected nodes with associated voltage constraints. In order to do this, we model power flow using a fixed-point linear approximation method. Additionally, we then use the model to explore design considerations of mini-grids in Sub-Saharan Africa. Specifically, we evaluate under what combinations of transmission line distance and capacity it is technically viable and economically preferable to build multiple isolated mini-grids versus an interconnected, centralized system.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Non -degenerate marginal-likelihood calibration with application to quantum characterization

Here, we propose a marginal likelihood strategy within the Kennedy-O’Hagan (KOH) Bayesian framework, where a Gaussian process (GP) models the discrepancy between a physical system and its simulator. Our approach introduces a novel marginalized likelihood by integrating out the degenerate eigenspace of the covariance matrix, rather than approximating the original likelihood. Unlike approximation methods that compromise accuracy for computational efficiency, our method defines an exact likelihood—distinct from the original but preserving all relevant information. This formulation achieves computational efficiency and stability, even for large datasets where the covariance matrix nears degeneracy. Applied to the characterization of a superconducting quantum device at Lawrence Livermore National Laboratory, the approach enhances the predictive accuracy of the Lindblad master equations for modeling Ramsey measurement data by effectively quantifying uncertainties consistent with the quantum data.

general physics↗