Search NASA⌕ Search

SEARCH · Search NASA

Results for “methods - statistical”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A splitting Chebyshev collocation method for Schrödinger–Poisson system

We develop a splitting Chebyshev collocation (SCC) method for the time-dependent Schrödinger–Poisson (SP) system arising from theoretical analysis of quantum plasmas. By means of splitting technique in time, the time-dependant SP system is first reduced to uncoupled Schrödinger and Poisson equations at every time step. The space variables in Schrödinger and Poisson equations are next represented by high-order Chebyshev polynomials, and the resulting system are discretized by the spectral collocation method. Finally, matrix diagonalization technique is applied to solve the fully discretized system in one dimension, two dimensions and three dimensions, respectively. The newly proposed method not only achieves spectral accuracy in space but also reduces the computer-memory requirements and the computational time in comparison with conventional solver. Numerical results confirm the spectral accuracy and efficiency of this method, and indicate that the SCC method could be an efficient alternative method for simulating the dynamics of quantum plasmas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

On the Feynman-alpha method for reflected fissile assemblies

The Feynman-alpha method is a neutron noise technique that is used to estimate the prompt neutron period of fissile assemblies. The method and quantity are of widespread interest including in applications such as nuclear criticality safety, safeguards and nonproliferation, and stockpile stewardship; the prompt neutron period may also be used to infer the k eff multiplication factor. The Feynman-alpha method is predicated on time-correlated neutron detections that deviate from a Poisson random variable due to multiplication. Traditionally, such measurements are diagnosed with one-region point kinetics, but two-region models are required when the fissile assembly is reflected. This paper presents a derivation of the two-region point kinetics Feynman-alpha equations based on a double integration of the Rossi-alpha equations, develops novel propagation of measurement uncertainty, and validates the theory. Validation is achieved with organic scintillator measurements of weapons-grade plutonium reflected by various amounts of copper to achieve k eff values of 0.83–0.94 and prompt periods of 5–75 ns. The results demonstrate that Feynman-alpha measurements should use the two-region model instead of the one-region model. The simplified one-region model deviates from the validated two-region models by as much as 10% in the estimate of the prompt neutron period, and the two-region model reduces to the one-region model for small amounts of reflector. The Feynman-alpha estimates of the prompt neutron period are compared to those of the Rossi-alpha approach. The comparative results demonstrate that the Feynman-alpha method is more precise than the Rossi-alpha method and more accurate for k eff < 0.92, whereas the Rossi-alpha method is generally more accurate for higher multiplications. Here, the uncertainty propagation developed in this work should be used for all Feynman-alpha measurements and will therein improve fitting accuracy and appropriate precision estimates.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2021

An uncertainty quantification method relevant to material test reactors

Within material test reactor calculations, energy dependent flux and reaction rate uncertainties are typically not quantified when performing as-run analyses to determine the neutron field experienced by the experiment. When high fidelity Monte-Carlo codes are used in such analyses, straight forward methods to calculate output uncertainties are not available, instead expert opinion is used to postulate computational uncertainties. New methods to propagate uncertainties through these high fidelity simulations are available when sufficient computational power is available. A tool is developed here for sampling any part of an MCNP input from random distributions to determine output uncertainties based on those inputs. Another tool is developed to sample nuclear data cross-section in ACE format using multi-group nuclear data covariances. The Total Monte-Carlo Method and Gesellschaft für Anlagen-und Reaktorsicherheit method (GRS) are implemented and compared to one another as well as MCNP sensitivity and uncertainty calculations. The methods were applied to the Godiva critical sphere k-eigenvalue, the UAM pin-cell benchmark energy dependent flux and reaction rates, and the Advanced Test Reactor energy dependent flux within an experimental location. Furthermore, the two methods agree well, with GRS allowing for an order of magnitude speedup for reaction rate uncertainty calculations and several orders of magnitude for eigenvalue uncertainty calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

A propagation-based fault detection and discrimination method and the optimization of sensor deployment

Industrial processes can be affected by faults having a serious impact on operation when not promptly detected and diagnosed. Here in this paper, a propagation-based fault detection and discrimination(PFDD) method is proposed to develop a strategy for fault diagnosis while in the design phase of a system. The PFDD method constructs the system model using the Integrated System Fault Analysis(ISFA) technique. Based on the system model, the propagation of hardware and software faults are simulated qualitatively. Given the results of the simulation, the process by which a fault propagates can be characterized using the qualitative features of system variables including the deviation of the system variables from their expected values, the variation of the system variables over time, and the order in which each variable is influenced during the propagation of the fault. The strategy by which a fault can be detected and discriminated is defined using those features. The PFDD method supports the detection and discrimination of faults in both steady states and transient states. Based on the PFDD method, the optimization of sensor deployment in a system is discussed. A brute force algorithm is developed to examine the system’s capability at diagnosing faults and the cost of sensor deployment for all possible configurations of sensors. The optimal sensor deployment strategy can be derived accordingly. However, the brute force method is only applicable to small-scale systems due to its high computational cost. A genetic algorithm is used to optimize sensor deployment in large-scale systems. The PFDD and sensor deployment optimization methods are applied to the Experimental Breeder Reactor II (EBR-II) for verification.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Particulate Fuel Modeling of MC 2 -3 using Iterative Local Spatial Self-shielding Method

We report a new spatial self-shielding method for particulate fuels has been developed based on disadvantage factors and implemented in the MC 2 -3 code. This method named the iterative local spatial self-shielding (ILSS) method considers the shadowing effect of randomly distributed particles on spatial self-shielding in particles through a homogenized composition region added outside the particle of interest at the center. The self-shielded cross sections of the central particle are determined iteratively since they are used in determining the cross sections of the homogenized composition region. The ILSS method was verified for infinite stochastic medium problems of single and multiple types of particles, VHTR unit cell problems, and HTTR assembly problems. The verification test results show that the ILSS method accurately predicts the stochastic particle shadowing effect and reaction rates in particles, whereas the regular array model and the stochastic collision probability method underpredict the particle shadowing effect and overestimate reaction rates in particles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A Hybrid Domain Overlapping Method for Coupling System Thermal Hydraulics and CFD Codes

A hybrid multiscale coupling methodology based on a domain overlapping approach has been developed for coupling System Thermal Hydraulics (STH) and Computational Fluid Dynamics (CFD) codes. The method has been implemented between the modern STH code SAM and the CFD code NekRS, using the coupling tool Cardinal. The coupling aims to extend the STH code's applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional mixing. Two coupling strategies are implemented and compared: a hybrid domain overlapping method and the conventional domain decomposition method. Here, the strategies are applied to two closed-loop applications, and the present method shows superior stability behavior when compared to the domain decomposition method. Then, the present coupling method is validated against experimental data from a double T-junction experiment. The present STH/CFD coupling shows improved agreement with experimental data when compared to STH standalone simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

On the convergence of the fixed point method for solving neutron transport alpha eigenvalue problems

It was shown that the Fixed Point Method (also known as the Rayleigh Quotient Method) is several times faster than the Critical Search Method for solving neutron transport alpha eigenvalue problems. It was also shown that the Fixed Point Method is able to determine the alpha eigenvalues of sub-critical systems that are beyond the reach of the Critical Search Method. Despite these significant advances, the Fixed Point Method remains an unproven algorithm. Here, this report provides a proof.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Sizing ramping reserve using probabilistic solar forecasts: A data-driven method

Ramping products have been introduced or proposed in several U.S. power markets to mitigate the impact of load and renewable uncertainties on market efficiency and reliability. Current methods often rely on historical data to estimate the requirements of ramping products and fail to take into account the effects of the latest weather conditions and their uncertainties, which could lead to overly conservative or insufficient requirements. This study proposes a k-nearest-neighbor-based method to give weather-informed estimates of ramping needs based on short-term probabilistic solar irradiance forecasts. Forecasts from multiple sites are employed in conjunction with principal component analysis to derive numerical classifiers to characterize system-level weather conditions. In addition, we develop a data-driven method to optimize the model parameters in a rolling-forward manner. By using real-world data from the California Independent System Operator, we design two metrics to evaluate method performance: 1) frequency of shortage and 2) oversupply of ramping product. Our proposed method presents advantages in comparison with the baseline and a set of benchmark methods: without compromising system reliability, it reduces system ramping requirements by up to 25%, therefore improving both system reliability and economics.

14 SOLAR ENERGY↗
DOE OSTI2022

Assessing mass balance-based inverse modeling methods via a pseudo-observation test to constrain NO x emissions over South Korea

This study constrains NO x emissions by utilizing mass balance–based inverse modeling methods and examines the optimization of inverse modeling conditions through a pseudo-observation test for South Korea, which has complex topography. Here, we applied mass balance–based inverse modeling methods such as Basic Mass Balance (BMB), Finite Difference Mass Balance (FDMB), and Iterative Finite Difference Mass Balance (IFDMB). We performed a numerical simulation of air quality using the Community Multi-scale Air Quality (CMAQ) model to calculate the NO 2 column density required for the inverse modeling. The pseudo-observation test was performed according to season, modeling resolution, and regridding methodology of satellite observation data to identify various conditions while applying the inverse modeling for South Korea. Comparing the inverse modeling results from the BMB, FDMB, and IFDMB methods, IFDMB was the most effective method in constraining NO x emissions in the South Korean region since it minimized smearing effects (i.e., transport-induced errors) through iterative calculations. The accuracy of the constrained NO x emissions using mass balance–based inversions in South Korea was the highest in the summer due to the minimized smearing effects, and the 9 km resolution modeling was the most efficient for inverse modeling of the region. In addition, the results of inverse modeling varied depending on the regridding methods, implying the importance of using a suitable regridding method and modeling resolution. This study used pseudo-observations, but the inversions are expected to be applied based on actual satellite data in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

A new comprehensive method for detection of livestock-related pathogenic viruses using a target enrichment system

Highlights: • Development of a comprehensive method to detect viruses using SureSelect system. • A method by which almost the full length of the viral genome can be determined without virus isolation. • Widely available methods for comprehensive analysis of virus genomic information. We tested usefulness of a target enrichment system SureSelect, a comprehensive viral nucleic acid detection method, for rapid identification of viral pathogens in feces samples of cattle, pigs and goats. This system enriches nucleic acids of target viruses in clinical/field samples by using a library of biotinylated RNAs with sequences complementary to the target viruses. The enriched nucleic acids are amplified by PCR and subjected to next generation sequencing to identify the target viruses. In many samples, SureSelect target enrichment method increased efficiencies for detection of the viruses listed in the biotinylated RNA library. Furthermore, this method enabled us to determine nearly full-length genome sequence of porcine parainfluenza virus 1 and greatly increased Breadth, a value indicating the ratio of the mapping consensus length in the reference genome, in pig samples. Our data showed usefulness of SureSelect target enrichment system for comprehensive analysis of genomic information of various viruses in field samples.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

An updated Lagrangian discontinuous Galerkin hydrodynamic method for gas dynamics

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for gas dynamics. The new method evolves conserved unknowns in the current configuration, which obviates the Jacobi matrix that maps the element in a reference coordinate system or the initial coordinate system to the current configuration. The density, momentum, and total energy (ρ, ρu, E ) are approximated with conservative higher-order Taylor expansions over the element and are limited toward a piecewise constant field near discontinuities using a limiter. Two new limiting methods are presented for enforcing the bounds on the primitive variables of density, velocity, and specific internal energy (ρ, u, e). The nodal velocity, and the corresponding forces, are calculated by solving an approximate Riemann problem at the element nodes. An explicit second-order method is used to temporally advance the solution. This new Lagrangian DG hydrodynamic method conserves mass, momentum, and total energy. 1D Cartesian coordinates test problem results are presented to demonstrate the accuracy and convergence order of the new DG method with the new limiters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A higher-order Lagrangian discontinuous Galerkin hydrodynamic method for solid dynamics

We present a new multidimensional high-order Lagrangian discontinuous Galerkin (DG) hydrodynamic method that supports hypoelastic and hyperelastic strength models for simulating solid dynamics with higher-order elements. We also present new one-dimensional test problems that have an analytic solution corresponding to a hyperelastic–plastic wave. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated high-order Taylor series polynomials. The stress fields are represented using nodal quantities. The constitutive models used to calculate the deviatoric stress are either a hypoelastic–plastic, infinitesimal strain hyperelastic–plastic, or finite strain hyperelastic–plastic model. These constitutive models require new methods for calculating high-order polynomials for the velocity gradient and deformation gradient in an element. The plasticity associated with the strength model is determined using a radial return method to with a J 2 yield criterion and perfect plasticity. The temporal evolution of the governing equations is achieved with the total variation diminishing Runge–Kutta (TVD RK) time integration method. A diverse suite of 1D and 2D test problems are calculated. Here, the new 1D piston test problems, which have analytic solutions for each elastic–plastic model, are presented and calculated to demonstrate the stability and formal accuracy of the various models with the new Lagrangian DG method. 2D test problems are calculated to demonstrate the stability and robustness of the new Lagrangian DG method on multidimensional problems with high-order elements, which have faces that can bend.

42 ENGINEERING↗
DOE OSTI2021

Parallel exponential time differencing methods for geophysical flow simulations

Two ocean models are considered for geophysical flow simulations: the multilayer shallow water equations and the multilayer primitive equations. For the former, we investigate the parallel performance of exponential time differencing (ETD) methods, including exponential Rosenbrock–Euler, ETD2wave, and B-ETD2wave. For the latter, we take advantage of the splitting of barotropic and baroclinic modes and propose a new two-level method in which an ETD method is applied to solve the fast barotropic mode. Furthermore, these methods could improve the computational efficiency of numerical simulations because ETD methods allow for much larger time step sizes than traditional explicit time-stepping techniques that are commonly used in existing computational ocean models. Several standard benchmark tests for ocean modeling are performed and comparison of the numerical results demonstrates a great potential of applying the parallel ETD methods for simulating real-world geophysical flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

A regularized deconvolution method for turbulent closure modeling in implicitly filtered large-eddy simulation

Turbulence closure models that are based on a regularized deconvolution method (RDM) are proposed for application to implicitly filtered Large-Eddy Simulation (LES). This method reconstructs sub-grid scale contributions by approximately recovering filtered quantities through an optimization approach. Physical principles on the boundedness and conservation of reconstructed quantities are preserved by constraining the optimization problem. Further, to account for effects of scales that are not resolved in implicitly filtered LES, a reconstruction method is proposed that combines sub-grid projection and energy-similarity to enrich unresolved scales. With this method, closure models are developed for turbulent scalar fluxes and turbulence–chemistry interaction that are intrinsically consistent with the LES formulation. These closure models are evaluated in application to a piloted turbulent jet flame with inhomogeneous inlet composition, in which combustion is represented by a two-scalar manifold model. The RDM-formulation is benchmarked against simulations obtained using a presumed PDF method and eddy diffusivity closure model. In partially-premixed flame regions, it is observed that RDM improves the prediction of turbulent mixing, and faster grid convergence is observed for results with RDM compared to simulations employing a presumed PDF method.

42 ENGINEERING↗
DOE OSTI2017

System dynamic simulation: A new method in social impact assessment (SIA)

Many complex social questions are difficult to address adequately with conventional methods and techniques, due to the complicated dynamics, and hard to quantify social processes. Despite these difficulties researchers and practitioners have attempted to use conventional methods not only in evaluative modes but also in predictive modes to inform decision making. The effectiveness of SIAs would be increased if they were used to support the project design processes. This requires deliberate use of lessons from retrospective assessments to inform predictive assessments. Social simulations may be a useful tool for developing a predictive SIA method. There have been limited attempts to develop computer simulations that allow social impacts to be explored and understood before implementing development projects. In light of this argument, this paper aims to introduce system dynamic (SD) simulation as a new predictive SIA method in large development projects. We propose the potential value of the SD approach to simulate social impacts of development projects. We use data from the SIA of Gareh-Bygone floodwater spreading project to illustrate the potential of SD simulation in SIA. It was concluded that in comparison to traditional SIA methods SD simulation can integrate quantitative and qualitative inputs from different sources and methods and provides a more effective and dynamic assessment of social impacts for development projects. We recommend future research to investigate the full potential of SD in SIA in comparing different situations and scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A method for automated modeling of plasma shape control response in tokamaks

Here, a simple method commonly used in modeling of tokamak conductor and plasma response for plasma control has been extended to allow more general modeling of tokamak PF conductor circuits. This extension combines a previously existing method for integrating plasma response and tokamak circuits with a more capable circuit modeling technique derived from the theory that underlies present-day automatic circuit simulation codes. This combined method can incorporate any combination of inductors (e.g., coils), resistors, capacitors, voltage sources, and current sources, while retaining the same model structure and ease of interpretation as the simple method. For many tokamaks, use of this extended method would not be the preferred approach, but for certain devices or for specialized models of a device, it is useful and occasionally essential. For example, the method has enabled fast and automated modeling of the complicated and daily-changing circuit topology used in the DIII-D tokamak. This automated modeling is essential to the support of controller design for the frequently changing requirements for DIII-D plasmas. It is likely to also find application in the workflows for new device designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A high-order multiscale finite-element method for time-domain elastic wave modeling in strongly heterogeneous media

Efficient and accurate numerical methods for elastic wave modeling in complex media have many important applications. However, it is fairly challenging to model elastic wave propagation in strongly heterogeneous media with high computational efficiency and high-order accuracy simultaneously. We develop a novel high-order multiscale finite-element method to model elastic wave propagation in strongly heterogeneous media in the time domain. The most important feature of our method is a generalization of standard multiscale finite element method by using high-order multiscale finite-element basis functions to capture the fine-scale heterogeneities on the coarse mesh, in contrast to conventional finite-element basis functions that are merely determined by the order of polynomials. These multiscale basis functions leads to a system matrix with significantly reduced dimension, thus enable us to solve the elastic wave equation on the coarse mesh with high-order accuracy and very low computational time cost. In conclusion, we use 2D and 3D numerical examples to demonstrate the superior efficiency and accuracy of our new modeling method compared with the conventional spectral-element method.

58 GEOSCIENCES↗
DOE OSTI2017

A robust bi-orthogonal/dynamically-orthogonal method using the covariance pseudo-inverse with application to stochastic flow problems

We develop a new robust methodology for the stochastic Navier–Stokes equations based on the dynamically-orthogonal (DO) and bi-orthogonal (BO) methods . Both approaches are variants of a generalized Karhunen–Loève (KL) expansion in which both the stochastic coefficients and the spatial basis evolve according to system dynamics, hence, capturing the low-dimensional structure of the solution. The DO and BO formulations are mathematically equivalent , but they exhibit computationally complimentary properties. Specifically, the BO formulation may fail due to crossing of the eigenvalues of the covariance matrix, while both BO and DO become unstable when there is a high condition number of the covariance matrix or zero eigenvalues. To this end, we combine the two methods into a robust hybrid framework and in addition we employ a pseudo-inverse technique to invert the covariance matrix. The robustness of the proposed method stems from addressing the following issues in the DO/BO formulation: (i) eigenvalue crossing: we resolve the issue of eigenvalue crossing in the BO formulation by switching to the DO near eigenvalue crossing using the equivalence theorem and switching back to BO when the distance between eigenvalues is larger than a threshold value; (ii) ill-conditioned covariance matrix: we utilize a pseudo-inverse strategy to invert the covariance matrix; (iii) adaptivity: we utilize an adaptive strategy to add/remove modes to resolve the covariance matrix up to a threshold value. In particular, we introduce a soft-threshold criterion to allow the system to adapt to the newly added/removed mode and therefore avoid repetitive and unnecessary mode addition/removal. When the total variance approaches zero, we show that the DO/BO formulation becomes equivalent to the evolution equation of the Optimally Time-Dependent modes . We demonstrate the capability of the proposed methodology with several numerical examples, namely (i) stochastic Burgers equation: we analyze the performance of the method in the presence of eigenvalue crossing and zero eigenvalues; (ii) stochastic Kovasznay flow: we examine the method in the presence of a singular covariance matrix; and (iii) we examine the adaptivity of the method for an incompressible flow over a cylinder where for large stochastic forcing thirteen DO/BO modes are active.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗