Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “forward sensitivity analysis”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Extended forward sensitivity analysis of one-dimensional isothermal flow

Sensitivity analysis and uncertainty quantification is an important part of nuclear safety analysis. In this work, forward sensitivity analysis is used to compute solution sensitivities on 1-D fluid flow equations typical of those found in system level codes. Time step sensitivity analysis is included as a method for determining the accumulated error from time discretization. The ability to quantify numerical error arising from the time discretization is a unique and important feature of this method. By knowing the relative sensitivity of time step with other physical parameters, the simulation is allowed to run at optimized time steps without affecting the confidence of the physical parameter sensitivity results. The time step forward sensitivity analysis method can also replace the traditional time step convergence studies that are a key part of code verification with much less computational cost. One well-defined benchmark problem with manufactured solutions is utilized to verify the method; another test isothermal flow problem is used to demonstrate the extended forward sensitivity analysis process. Through these sample problems, the paper shows the feasibility and potential of using the forward sensitivity analysis method to quantify uncertainty in input parameters and time step size for a 1-D system-level thermal-hydraulic safety code. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Extended Forward Sensitivity Analysis for Uncertainty Quantification

This paper presents the extended forward sensitivity analysis as a method to help uncertainty qualification. By including time step and potentially spatial step as special sensitivity parameters, the forward sensitivity method is extended as one method to quantify numerical errors. Note that by integrating local truncation errors over the whole system through the forward sensitivity analysis process, the generated time step and spatial step sensitivity information reflect global numerical errors. The discretization errors can be systematically compared against uncertainties due to other physical parameters. This extension makes the forward sensitivity method a much more powerful tool to help uncertainty qualification. By knowing the relative sensitivity of time and space steps with other interested physical parameters, the simulation is allowed to run at optimized time and space steps without affecting the confidence of the physical parameter sensitivity results. The time and space steps forward sensitivity analysis method can also replace the traditional time step and grid convergence study with much less computational cost. Two well-defined benchmark problems with manufactured solutions are utilized to demonstrate the method.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2011

Extended Forward Sensitivity Analysis for Uncertainty Quantification

Verification and validation (V&V) are playing more important roles to quantify uncertainties and realize high fidelity simulations in engineering system analyses, such as transients happened in a complex nuclear reactor system. Traditional V&V in the reactor system analysis focused more on the validation part or did not differentiate verification and validation. The traditional approach to uncertainty quantification is based on a 'black box' approach. The simulation tool is treated as an unknown signal generator, a distribution of inputs according to assumed probability density functions is sent in and the distribution of the outputs is measured and correlated back to the original input distribution. The 'black box' method mixes numerical errors with all other uncertainties. It is also not efficient to perform sensitivity analysis. Contrary to the 'black box' method, a more efficient sensitivity approach can take advantage of intimate knowledge of the simulation code. In these types of approaches equations for the propagation of uncertainty are constructed and the sensitivities are directly solved for as variables in the simulation. This paper presents the forward sensitivity analysis as a method to help uncertainty qualification. By including time step and potentially spatial step as special sensitivity parameters, the forward sensitivity method is extended as one method to quantify numerical errors. Note that by integrating local truncation errors over the whole system through the forward sensitivity analysis process, the generated time step and spatial step sensitivity information reflect global numerical errors. The discretization errors can be systematically compared against uncertainties due to other physical parameters. This extension makes the forward sensitivity method a much more powerful tool to help uncertainty qualification. By knowing the relative sensitivity of time and space steps with other interested physical parameters, the simulation is allowed to run at optimized time and space steps without affecting the confidence of the physical parameter sensitivity results. The time and space steps forward sensitivity analysis method can also replace the traditional time step and grid convergence study with much less computational cost. Several well defined benchmark problems with manufactured solutions are utilized to demonstrate the extended forward sensitivity analysis method. All the physical solutions, parameter sensitivity solutions, even the time step sensitivity in one case, have analytical forms, which allows the verification to be performed in the strictest sense.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

On the dual advantage of placing observations through forward sensitivity analysis

The four-dimensional variational data assimilation methodology for assimilating noisy observations into a deterministic model has been the workhorse in forecasting centres for over three decades. While this method provides a computationally efficient framework for dynamic data assimilation, it is largely silent on the important question concerning the minimum number and placement of observations. To answer this question, we demonstrate the dual advantage of placing the observations where the square of the sensitivity of the model solution with respect to the unknown control variables, called forward sensitivities, attains its maximum. By following this approach, we can force the observability Gramian to be of full rank, which in turn guarantees efficient recovery of the optimal values of the control variables, which is the first of the two advantages of this strategy. We further show that the proposed strategy of placing observations has another inherent optimality: the square of the sensitivity of the optimal estimates of the control with respect to the observations (used to obtain these estimates) attains its minimum value, a second advantage that is a direct consequence of the above strategy for placing observations. Our analytical framework and numerical experiments on linear and nonlinear systems confirm the effectiveness of our proposed strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Forward sensitivity analysis and mode dependent control for closure modeling of Galerkin systems

Model reduction by projection-based approaches is often associated with losing some of the important features that contribute towards the dynamics of the retained scales. As a result, a mismatch occurs between the predicted trajectories of the original system and the truncated one. We put forth a framework to apply a continuous time control signal in the latent space of the reduced order model (ROM) to account for the effect of truncation. We set the control input using parameterized models by following energy transfer principles. Our methodology relies on observing the system behavior in the physical space and using the projection operator to restrict the feedback signal into the latent space. Then, we leverage the forward sensitivity method (FSM) to derive relationships between the feedback and the desired mode-dependent control. We test the performance of the proposed approach using two test cases, corresponding to viscous Burgers and vortex merger problems at high Reynolds number. Results show that the ROM trajectory with the applied FSM control closely matches its target values in both the data-dense and data-sparse regimes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Differential methods for assessing sensitivity in biological models

Differential sensitivity analysis is indispensable in fitting parameters, understanding uncertainty, and forecasting the results of both thought and lab experiments. Although there are many methods currently available for performing differential sensitivity analysis of biological models, it can be difficult to determine which method is best suited for a particular model. In this paper, we explain a variety of differential sensitivity methods and assess their value in some typical biological models. First, we explain the mathematical basis for three numerical methods: adjoint sensitivity analysis, complex perturbation sensitivity analysis, and forward mode sensitivity analysis. We then carry out four instructive case studies. (a) The CARRGO model for tumor-immune interaction highlights the additional information that differential sensitivity analysis provides beyond traditional naive sensitivity methods, (b) the deterministic SIR model demonstrates the value of using second-order sensitivity in refining model predictions, (c) the stochastic SIR model shows how differential sensitivity can be attacked in stochastic modeling, and (d) a discrete birth-death-migration model illustrates how the complex perturbation method of differential sensitivity can be generalized to a broader range of biological models. Finally, we compare the speed, accuracy, and ease of use of these methods. We find that forward mode automatic differentiation has the quickest computational time, while the complex perturbation method is the simplest to implement and the most generalizable.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Efficient Adjoint Computation of Hybrid Systems of Differential Algebraic Equations with Applications in Power Systems

Sensitivity analysis is an important tool to describe power system dynamic behavior in response to parameter variations. It is a central component in preventive and corrective control applications. The existing approaches for sensitivity calculations, namely, finite-difference and forward sensitivity analysis, require a computational effort that increases linearly with the number of sensitivity parameters. In this work, we investigate, implement, and test a discrete adjoint sensitivity approach whose computational effort is effectively independent of the number of sensitivity parameters. The proposed approach is highly efficient for calculating trajectory sensitivities of larger systems and is consistent, within machine precision, with the function whose sensitivity we are seeking. This is an essential feature for use in optimization applications. Moreover, our approach includes a consistent treatment of systems with switching, such as DC exciters, by deriving and implementing the adjoint jump conditions that arise from state and time-dependent discontinuities. The accuracy and the computational efficiency of the proposed approach are demonstrated in comparison with the forward sensitivity analysis approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Discrete Adjoint Sensitivity Analysis of Hybrid Dynamical Systems With Switching [Discrete Adjoint Sensitivity Analysis of Hybrid Dynamical Systems]

Sensitivity analysis is an important tool for describing power system dynamic behavior in response to parameter variations. It is a central component in preventive and corrective control applications. The existing approaches for sensitivity calculations, namely, finite-difference and forward sensitivity analysis, require a computational effort that increases linearly with the number of sensitivity parameters. In this paper, we investigate, implement, and test a discrete adjoint sensitivity approach whose computational effort is effectively independent of the number of sensitivity parameters. The proposed approach is highly efficient for calculating sensitivities of larger systems and is consistent, within machine precision, with the function whose sensitivity we are seeking. This is an essential feature for use in optimization applications. Moreover, our approach includes a consistent treatment of systems with switching, such as dc exciters, by deriving and implementing the adjoint jump conditions that arise from state-dependent and time-dependent switchings. The accuracy and the computational efficiency of the proposed approach are demonstrated in comparison with the forward sensitivity analysis approach. In conclusion, this paper focuses primarily on the power system dynamics, but the approach is general and can be applied to hybrid dynamical systems in a broader range of fields.

42 ENGINEERING↗
DOE OSTI2012

Generalized Perturbation Theory-Free Sensitivity Analysis for Eigenvalue Problems

Generalized perturbation theory (GPT) has been recognized as the most computationally efficient approach for performing sensitivity analysis for models with many input parameters, which renders forward sensitivity analysis computationally overwhelming. In critical systems, GPT involves the solution of the adjoint form of the eigenvalue problem with a response-dependent fixed source. Although conceptually simple to implement, most neutronics codes that can solve the adjoint eigenvalue problem do not have a GPT capability unless envisioned during code development. We introduce in this manuscript a reduced-order modeling approach based on subspace methods that requires the solution of the fundamental adjoint equations but allows the generation of response sensitivities without the need to set up GPT equations, and that provides an estimate of the error resulting from the reduction. Moreover, the new approach solves the eigenvalue problem independently of the number or type of responses. This allows for an efficient computation of sensitivities when many responses are required. This paper introduces the theory and implementation details of the GPT-free approach and describes how the errors could be estimated as part of the analysis. The applicability is demonstrated by estimating the variations in the flux distribution everywhere in the phase space of a fast critical sphere and a high-temperature gas-cooled reactor prismatic lattice. The variations generated by the GPT-free approach are benchmarked to the exact variations generated by direct forward perturbations.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2012

Use of Forward Sensititvity Analysis Method to Improve Code Scaling, Applicability, and Uncertainty (CSAU) Methodology

Since the Code Scaling, Applicability, and Uncertainty (CSAU) methodology was proposed about two decades ago, it has been widely used for new reactor designs and existing LWRs power uprates. In spite of these huge successes, CSAU has been criticized for the need of further improvement, focusing on two main issues - lack of objectiveness and high cost. With the effort to develop next generation safety analysis codes, new opportunities appear to take advantage of new numerical methods, better physical models, and modern uncertainty qualification methods. Forward sensitivity analysis (FSA) directly solves the partial differential equations for parameter sensitivities. Moreover, our work shows that time and space steps can be treated as special sensitivity parameters so that numerical errors can be directly compared with physical uncertainties. When the FSA is implemented in a new advanced system analysis code, CSAU could be significantly improved by quantifying numerical errors and allowing a quantitative PIRT (Q-PIRT) to reduce subjective judgement and improving efficiency. This paper will review the issues related to the current CSAU implementations, introduce FSA, show a simple example to perform FSA, and discuss potential improvements on CSAU with FSA. Finally, the general research direction and requirements to use FSA in an advanced system analysis code will be discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Analytical Sensitivity Analysis of a Spent Nuclear Fuel Cask

Here, we report nuclear science and engineering is a field increasingly dominated by computational studies resulting from increasingly powerful computational tools. As a result, analytical studies, which previously pioneered nuclear engineering, are increasingly viewed as secondary or unnecessary. However, analytical solutions to reduced-fidelity models can provide important information concerning the underlying physics of a problem and aid in guiding computational studies. Similarly, there is increased interest in sensitivity analysis studies. These studies commonly use computational tools. However, providing a complementary sensitivity study of relevant analytical models can lead to a deeper analysis of a problem. This work provides the analytical sensitivity analysis of the one-dimensional (1D) cylindrical mono-energetic neutron diffusion equation using the forward sensitivity analysis procedure (FSAP) developed by Cacuci. Further, these results are applied to a reduced-fidelity model of a spent nuclear fuel cask, demonstrating how computational analysis might be improved with a complementary analytic sensitivity analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRSDate not supplied

Assimilation of Active MW and Radar Observations in the NWP Models

Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. This talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗
NASA NTRSDate not supplied

Enhancing Radiative Transfer Models for Optimized Assimilation of Microwave and Radar Observations

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗
NASA NTRSDate not supplied

Optimizing Assimilation of Microwave and Radar Observations in the NWP Models

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗
NASA NTRS2023

Developing a Radar Signal Simulator for the Community Radiative Transfer Model

Active radar instruments provide vertically resolved clouds and precipitation measurements that cannot be provided by the passive instruments. These active measurements are not conventionally assimilated into the data assimilation systems because of the lack of fast forward radiative transfer models and also difficulties in the error modelling of the measurements. This paper describes the development, evaluation, and sensitivity analysis for a forward radar model implemented in the Community Radiative Transfer Model (CRTM). The scattering properties required by the forward model are provided by the hydrometeor lookup tables that were generated using the discrete dipole approximation. The model is able to calculate both the reflectivity and the attenuated reflectivity for any given radar instrument at any given zenith angles as long as CRTM instrument specific coefficients are available. The evaluation using CloudSat measurements shows a very good agreement between the simulations and measurements as long as the input profiles of hydrometeors are consistent with the measured reflectivity profiles. Major sources contributing to the differences between the measured and simulated reflectivities are input hydrometeor profiles, scattering lookup tables, lack of melting layer in the forward model, CRTM scattering solvers, and attenuation calculations. In addition to the forward model, both Tangent Linear and Adjoint of the model are also implemented and tested within CRTM. These components may be required by some data assimilation systems for the assimilation of radar measurements.

radar↗
NASA NTRSDate not supplied

Assimilation of Spaceborne Microwave and Radar Observations

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. The Community Radiative Transfer Model (CRTM) is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk optical properties of hydrometeors in the form of lookup tables in order to perform all-sky RT calculations. However, the current cloud scattering lookup tables in CRTM assume spherical shapes for all frozen hydrometeors, whereas actual clouds contain frozen particles with diverse shapes. The first part of this talk presents the implementation and validation of a comprehensive Discrete Dipole Approximation (DDA) cloud scattering database into CRTM, specifically targeting microwave frequencies. The DDA technique proves effective in simulating the optical properties of non-spherical hydrometeors in the microwave region. The original DDA database assumes total random orientation in calculating single scattering properties. To generate the required mass scattering parameters for CRTM, the single scattering properties and water content dependent particle size distributions are used. The evaluation of results involved a collocated dataset comprising short-term forecasts from the Integrated Forecast System of the European Center for Medium-Range Weather Forecasts and satellite microwave data. The findings demonstrate that the DDA lookup tables significantly reduce discrepancies between simulated and observed values when compared to the Mie tables. Passive instruments lack the ability to provide vertically resolved measurements of clouds and precipitation, which can be obtained by active radar instruments. However, incorporating these active measurements into data assimilation systems presents challenges due to the absence of fast forward radiative transfer models and difficulties in error modeling. The second part of the talk focuses on the development, evaluation, and sensitivity analysis of a forward radar model integrated into CRTM. The forward radar model utilizes scattering properties obtained from hydrometeor lookup tables generated using the discrete dipole approximation. By utilizing CRTM instrument-specific coefficients, the model can calculate both reflectivity and attenuated reflectivity for any given radar instrument and zenith angles. Evaluation using CloudSat measurements demonstrates a strong agreement between simulations and observations when the input profiles of hydrometeors align with the measured reflectivity profiles.

Isaac Moradi↗