Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Modular evaluation method for subsurface activities (MEMSA). A novel approach for integrating social acceptance in a permit decision-making process for subsurface activities

We investigate how the decision support system ‘Modular Evaluation Method Subsurface Activities’ (MEMSA) can help facilitate an informed decision-making process for permit applications of subsurface activities. To this end, we analyze the extent the MEMSA approach allows for a dialogue between stakeholders in a transparent manner. We use the exploration permit for the underground gas storage facility at the Pieterburen salt dome (Netherlands) as a case study. The results suggest that the MEMSA approach is flexible enough to adjust to changing conditions. Furthermore, MEMSA provides a novel way for identifying structural problems and possible solutions in permit decision-making processes for subsurface activities, on the basis of the sensitivity analysis of intermediate rankings. We suggest that the planned size of an activity should already be specified in the exploration phase, because this would allow for a more efficient use of the subsurface as a whole. We conclude that the host community should be involved to a greater extent and in an early phase of the permit decision-making process, for example, already during the initial analysis of the project area of a subsurface activity. We suggest that strategic national policy goals are to be re-evaluated on a regular basis, in the form of a strategic vision for the subsurface, to account for timing discrepancies between the realization of activities and policy deadlines, because this discrepancy can have a large impact on the necessity and therefore acceptance of a subsurface activity.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Process-informed adsorbent design guidelines for direct air capture

Direct air capture using solid adsorbents is a proven technology critical to reducing our net greenhouse gas emissions to zero and beyond. Currently, academic research into the technology mainly focuses on the development of new adsorbents. However, there is a discord between the adsorbent design and process performance. Many materials scientists focus on maximising metrics such as the CO 2 capacity of their adsorbent. Here, we combine detailed process modelling, machine learning, and extensive global sensitivity analysis, which entails varying all of the model parameters together, on a direct air capture process to show that the dry CO 2 adsorption capacity does not influence process performance for an amine-functionalised adsorbent operating in a temperature vacuum swing adsorption (TVSA) process, while it is important in a steam-assisted TVSA (S-TVSA) process. In fact, adsorption kinetics, density, and thermal conductivity are all critical attributes to obtaining a low energy penalty and reduced costs. The analysis also highlights the importance of heat transfer, directing process engineers to (alternative) adsorber designs that maximise this. By an in-depth evaluation of how process performance indicators are affected by materials properties and process operating parameters, this work provides guidance to both material scientists and process engineers towards the design of a “unicorn adsorbent” and intensified DAC processes. This will improve the performance of solid adsorbent direct air capture and help drive down the costs of this vital technology to avert the worst impacts of climate change.

42 ENGINEERING↗
DOE OSTI2011

Advanced Test Reactor Core Modeling Update Project Annual Report for Fiscal Year 2011

Legacy computational reactor physics software tools and protocols currently used for support of Advanced Test Reactor (ATR) core fuel management and safety assurance and, to some extent, experiment management are obsolete, inconsistent with the state of modern nuclear engineering practice, and are becoming increasingly difficult to properly verify and validate (V&V). Furthermore, the legacy staff knowledge required for application of these tools and protocols from the 1960s and 1970s is rapidly being lost due to staff turnover and retirements. In 2009 the Idaho National Laboratory (INL) initiated a focused effort to address this situation through the introduction of modern high-fidelity computational software and protocols, with appropriate V&V, within the next 3-4 years via the ATR Core Modeling and Simulation and V&V Update (or 'Core Modeling Update') Project. This aggressive computational and experimental campaign will have a broad strategic impact on the operation of the ATR, both in terms of improved computational efficiency and accuracy for support of ongoing DOE programs as well as in terms of national and international recognition of the ATR National Scientific User Facility (NSUF). The ATR Core Modeling Update Project, targeted for full implementation in phase with the anticipated ATR Core Internals Changeout (CIC) in the 2014 time frame, began during the last quarter of Fiscal Year 2009, and has just completed its first full year. Key accomplishments so far have encompassed both computational as well as experimental work. A new suite of stochastic and deterministic transport theory based reactor physics codes and their supporting nuclear data libraries (SCALE, KENO-6, HELIOS, NEWT, and ATTILA) have been installed at the INL under various permanent sitewide license agreements and corresponding baseline models of the ATR and ATRC are now operational, demonstrating the basic feasibility of these code packages for their intended purpose. Furthermore, a capability for rigorous sensitivity analysis and uncertainty quantification based on the TSUNAMI system is being implemented and initial computational results have been obtained. This capability will have many applications in 2011 and beyond as a tool for understanding the margins of uncertainty in the new models as well as for validation experiment design and interpretation. Finally we note that although full implementation of the new computational models and protocols will extend over a period 3-4 years as noted above, interim applications in the much nearer term have already been demonstrated. In particular, these demonstrations included an analysis that was useful for understanding the cause of some issues in December 2009 that were triggered by a larger than acceptable discrepancy between the measured excess core reactivity and a calculated value that was based on the legacy computational methods. As the Modeling Update project proceeds we anticipate further such interim, informal, applications in parallel with formal qualification of the system under the applicable INL Quality Assurance procedures and standards.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2012

Behavior of an heterogeneous annular FBR core during an unprotected loss of flow accident: Analysis of the primary phase with SAS-SFR

In the framework of a substantial improvement on FBR core safety connected to the development of a new Gen IV reactor type, heterogeneous core with innovative features are being carefully analyzed in France since 2009. At EDF R and D, the main goal is to understand whether a strong reduction of the Na-void worth - possibly attempting a negative value - allows a significant improvement of the core behavior during an unprotected loss of flow accident. Also, the physical behavior of such a core is of interest, before and beyond the (possible) onset of Na boiling. Hence, a cutting-edge heterogeneous design, featuring an annular shape, a Na-plena with a B{sub 4}C plate and a stepwise modulation of fissile core heights, was developed at EDF by means of the SDDS methodology, with a total Na-void worth of -1 $. The behavior of such a core during the primary phase of a severe accident, initiated by an unprotected loss of flow, is analyzed by means of the SAS-SFR code. This study is carried-out at KIT and EDF, in the framework of a scientific collaboration on innovative FBR severe accident analyses. The results show that the reduction of the Na-void worth is very effective, but is not sufficient alone to avoid Na-boiling and, hence, to prevent the core from entering into the primary phase of a severe accident. Nevertheless, the grace time up to boiling onset is greatly enhanced in comparison to a more traditional homogeneous core design, and only an extremely low fraction of the fuel (<0.1%) enters into melting at the end of this phase. A sensitivity analysis shows that, due to the inherent neutronic characteristics of such a core, the gagging scheme plays a major role on the core behavior: indeed, an improved 4-zones gagging scheme, associated with an enhanced control rod drive line expansion feed-back effect, finally prevents the core from entering into sodium boiling. This major conclusion highlights both the progress already accomplished and the need for more detailed future analyses, particularly concerning: the neutronic burn-up scheme, the modeling of the diagrid effect and the control rod drive line expansion feed-backs, as well as the primary/secondary systems thermal-hydraulics behavior. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Development of control quality factor for HVAC control loop performance assessment I—Methodology (ASHRAE RP-1587)

This article is the second paper from the ASHRAE research project RP-1587, focusing on the methodology of obtaining a control quality factor (CQF). This article presents a development of two CQFs for assessing the heating, ventilation, and air-conditioning (HVAC) control loop performance. Both CQFs are able to detect whether the control loop is able to maintain the setpoint and identify the loop's ability to handle disturbances. In addition, the reversal behaviors are assessed as well. The first CQF (CQF-Harris) is proposed based on the normalized Harris index using the recursive least squares method. This recursive least squares method is selected because of its computational efficiency compared with the maximum likelihood estimation method. The second CQF (CQF-EWMA) is based on the exponentially weighted moving average of the error ratio. The assessment scale of excellent, good, fair, bad, and failed, which indicates the quality of the HVAC control loops, is established as well. The sensitivity analysis for both CQFs is also conducted, and it provides insights on choosing the appropriate parameters to compute such CQFs. Such parameters include the sampling frequency, the length of the moving window, and the variance of the unmeasured disturbance. The field evaluations and tests of the proposed CQFs for simulated control loops and real control loops can be found in the companion paper with the title 'Development of Control Quality Factor for HVAC Control Loop Performance Assessment - III: Field Testing and Results (ASHRAE RP-1587).'

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2012

Adaptive Algebraic Smoothers

methods based on Local Sensitivity Analysis (LSA). The method can be used in the context of geometric and algebraic multigrid methods for constructing smoothers, and in the context of Krylov methods for constructing block preconditioners. It is suitable for both constant and variable coecient problems. Furthermore, the method can be applied to systems arising from both scalar and coupled system partial differential equations (PDEs), as well as linear systems that do not arise from PDEs. The simplicity of the method will allow it to be easily incorporated into existing multigrid and Krylov solvers while providing a powerful tool for adaptively constructing methods tuned to a problem.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

A Mathematical Model of Breast Tumor Progression Based on Immune Infiltration

Breast cancer is the most prominent type of cancer among women. Understanding the microenvironment of breast cancer and the interactions between cells and cytokines will lead to better treatment approaches for patients. In this study, we developed a data-driven mathematical model to investigate the dynamics of key cells and cytokines involved in breast cancer development. We used gene expression profiles of tumors to estimate the relative abundance of each immune cell and group patients based on their immune patterns. Dynamical results show the complex interplay between cells and molecules, and sensitivity analysis emphasizes the direct effects of macrophages and adipocytes on cancer cell growth. In addition, we observed the dual effect of IFN-γ on cancer proliferation, either through direct inhibition of cancer cells or by increasing the cytotoxicity of CD8+ T-cells.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Acoustic Tomography of the Atmosphere: A Large-Eddy Simulation Sensitivity Study

Accurate measurement of atmospheric turbulent fluctuations is critical for understanding environmental dynamics and improving models in applications such as wind energy. Advanced remote sensing technologies are essential for capturing instantaneous velocity and temperature fluctuations. Acoustic tomography (AT) offers a promising approach that utilizes sound travel times between an array of transducers to reconstruct turbulence fields. This study presents a systematic evaluation of the time-dependent stochastic inversion (TDSI) algorithm for AT using synthetic travel-time measurements derived from large-eddy simulation (LES) fields under both neutral and convective atmospheric boundary-layer conditions. Unlike prior work that relied on field observations or idealized fields, the LES framework provides a ground-truth atmospheric state, enabling quantitative assessment of TDSI retrieval reliability, sensitivity to travel-time measurement noise, and dependence on covariance model parameters and temporal data integration. A detailed sensitivity analysis was conducted to determine the best-fit model parameters, identify the tolerance thresholds for parameter mismatch, and establish a maximum spatial resolution. The TDSI algorithm successfully reconstructed large-scale velocity and temperature fluctuations with root mean square errors ( RMSE s) below 0.35 m/s and 0.12 K, respectively. Spectral analysis established a maximum spatial resolution of approximately 1.4 m, and reconstructions remained robust for travel-time measurement uncertainties up to 0.002 s. These findings provide critical insights into the operational limits of TDSI and inform future applications of AT for atmospheric turbulence characterization and system design.

17 WIND ENERGY↗
DOE OSTI2012

Appendix I3-1 to Wind HUI Initiative 1: AWST-WindNET-Phase 1 Final Report

This report is an appendix to the Hawaii WindHUI efforts to develop and operationalize short-term wind forecasting and wind ramp event forecasting capabilities. The report summarizes the WindNET Phase 1 efforts on the Big Island of Hawaii and includes descriptions of modeling methodologies, use of field validation data, results and recommendations. The objective of the WindNET project was to investigate the improvement that could be obtained in short-term wind power forecasting for wind generation facilities operating on the island grids operated by Hawaiian Electric Companies through the use of atmospheric sensors deployed at targeted locations. WindNET is envisioned as a multiphase project that will address the short-term wind forecasting issues of all of the wind generation facilities on the all of the Hawaiian Electric Companies' island grid systems. The first phase of the WindNET effort (referred to as WindNET-1) was focused on the wind generation facilities on the Big Island of Hawaii. With complex terrain and marine environment, emphasis was on improving the 0 to 6 hour forecasts of wind power ramps and periods of wind variability, with a particular interest in the intra-hour (0-1 hour) look-ahead period. The WindNET project was built upon a foundation that was constructed with the results from a previously completed observation targeting study for the Big Island that was conducted as part of a project supported by the National Renewable Energy Laboratory (NREL) and interactions with the western utilities. The observational targeting study provided guidance on which variables to measure and at what locations to get the most improvement in forecast performance at a target forecast site. The recommendations of the observation targeting study were based on the application two techniques: (1) an objective method called ensemble sensitivity analysis (ESA) (Ancell and Hakim, 2007; Torn and Hakim, 2008; Zack et al, 2010); and (2) a subjective method based on a diagnostic analysis of large ramp events. The analysis was completed for both the wind farm on the southern tip of the Big Island and on the northern tip of the island. The WindNET project was designed to also deploy sensors to validate the Big Island observational targeting study and enhance operator's understanding of predominate causes of wind variability conditions at the wind facilities. Compromises had to be made with the results from the observation targeting study to accommodate project resource limitations, availability of suitable sites, and other factors. To focus efforts, field sensor deployment activities focused on the wind facility on the southern point of Big Island.

14 SOLAR ENERGY↗
DOE OSTI2017

Life Cycle Assessment to support the quantification of the environmental impacts of an event

In recent years, several tools have been used to define and quantify the environmental impacts associated with an event; however, a lack of uniform approaches for conducting environmental evaluations has been revealed. The aim of this paper is to evaluate whether the Life Cycle Assessment methodology, which is rarely applied to an event, can be an appropriate tool for calculating the environmental impacts associated with the assembly, disassembly, and use phase of an event analysing in particular the components and the displays used to establish the exhibits. The aim is also to include the issues reported by ISO 20121:2012 involving the interested parties that can be monitored but also affected by the event owner, namely the event organiser, the workforce and the supply chain. A small event held in Northern Italy was selected as the subject of the research. The results obtained show that the main contributors are energy consumption for lighting and heating and the use of aluminium materials, such as bars for supporting the spotlights, carpet and the electronic equipment. A sensitivity analysis for estimating the effects of the impact assessment method chosen has also been conducted and an uncertainty analysis has been performed using the Monte Carlo technique. This study highlighted the importance of the energy consumed by heating and lighting on the environmental implications, and indicated that the preparation and assembly should always be considered when quantifying the environmental profile of an event. - Highlights: • LCA methodology, developed for products and services, is applied to an event. • A small event held in Northern Italy is analysed. • The main contributors are energy consumption and the use of aluminium and carpet. • Exhibition site preparation can have important environmental implications. • This study demonstrates the importance of the assembly, disassembly and use phase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Time-frequency energy density precipitation method for time-of-flight extraction of narrowband Lamb wave detection signals

The time-of-flight of the Lamb wave provides an important basis for defect evaluation in metal plates and is the input signal for Lamb wave tomographic imaging. However, the time-of-flight can be difficult to acquire because of the Lamb wave dispersion characteristics. This work proposes a time-frequency energy density precipitation method to accurately extract the time-of-flight of narrowband Lamb wave detection signals in metal plates. In the proposed method, a discrete short-time Fourier transform is performed on the narrowband Lamb wave detection signals to obtain the corresponding discrete time-frequency energy density distribution. The energy density values at the center frequency for all discrete time points are then calculated by linear interpolation. Next, the time-domain energy density curve focused on that center frequency is precipitated by least squares fitting of the calculated energy density values. Finally, the peak times of the energy density curve obtained relative to the initial pulse signal are extracted as the time-of-flight for the narrowband Lamb wave detection signals. An experimental platform is established for time-of-flight extraction of narrowband Lamb wave detection signals, and sensitivity analysis of the proposed time-frequency energy density precipitation method is performed in terms of propagation distance, dispersion characteristics, center frequency, and plate thickness. For comparison, the widely used Hilbert–Huang transform method is also implemented for time-of-flight extraction. The results show that the time-frequency energy density precipitation method can accurately extract the time-of-flight with relative error of <1% and thus can act as a universal time-of-flight extraction method for narrowband Lamb wave detection signals.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Rainfall erosivity in Slovenia: Sensitivity estimation and trend detection

Highlights: • Rainfall erosivity for 10 rainfall stations located in Slovenia was evaluated. • Rainfall erosivity trends are mostly positive while mostly statistically insignificant. • The maximum rainfall erosivity is obtained under Alpine climate. • Monthly rainfall erosivity was computed. • Sensitivity analyses were carried out. Slovenia is one of the EU countries with the largest values and largest amounts of variability in rainfall erosivity, with maximum annual values exceeding 10,000 MJ mm ha{sup -1} h{sup -1} yr{sup -1}. Five-minute rainfall data was analysed from 10 Slovenian rainfall stations with data-length availability longer than 25 years with a maximum data length of 69 years and a total data-station length equal to 443 years. Trends in the rainfall erosivity R-factor were detected for four different sub-samples using monthly, half-year, and annual rainfall erosivity values. The results indicate that rainfall erosivity trends for the selected Slovenian stations are mostly statistically insignificant, with the selected significance level of 0.05. However, a larger share of identified trends are positive than negative. The maximum annual rainfall erosivity values were obtained for one specific mountain station. Furthermore, a sensitivity analysis regarding the rainfall erosivity factor R calculation showed that the rainfall threshold parameter (12.7 mm) that is used to remove the small-magnitude rainfall events in order to reduce the computational burden can attribute up to 10% of the average annual R-values in cases where this threshold is not used. Other parameters have, on average, a smaller impact on the calculated rainfall erosivity. Furthermore, the application of local kinetic energy equations resulted in, on average, about 20% higher annual rainfall erosivity values compared to the equation that is proposed by the Revised Universal Soil Loss Equation (RUSLE) manual and was not developed specifically for this region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Environmental performance of household waste management in Europe - An example of 7 countries

Highlights: • Attributional LCA of household waste management in 7 EU countries. • Country-specific technologies and waste composition are crucial for impacts. • Proper waste management is a benefit to society. • Waste hierarchy focus needs to change from quantity to quality. • Need for overall improvement in waste data. - Abstract: An attributional life cycle assessment (LCA) of the management of 1 ton of household waste was conducted in accordance with ISO 14044:2006 and the ILCD Handbook for seven European countries, namely Germany, Denmark, France, UK, Italy, Poland and Greece, representing different household waste compositions, waste management practices, technologies, and energy systems. National data were collected from a range of sources regarding household waste composition, household sorting efficiency, collection, waste treatments, recycling, electricity and heat composition, and technological efficiencies. The objective was to quantify the environmental performance in the different countries, in order to analyze the sources of the main environmental impacts and national differences which affect the results. In most of the seven countries, household waste management provides environmental benefits when considering the benefits of recycling of materials and recovering and utilization of energy. Environmental benefits come from paper recycling and, to a lesser extent, the recycling of metals and glass. Waste-to-energy plants can lead to an environmental load (as in France) or a saving (Germany and Denmark), depending mainly on the composition of the energy being substituted. Sensitivity analysis and a data quality assessment identified a range of critical parameters, suggesting from where better data should be obtained. The study concluded that household waste management is environmentally the best in European countries with a minimum reliance on landfilling, also induced by the implementation of the Waste Hierarchy, though environmental performance does not correlate clearly with the rate of material recycling. From an environmental point of view, this calls for a change in the waste management paradigm, with less focus on where the waste is routed and more of a focus on the quality and utilization of recovered materials and energy.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2013

Evolution of US DOE Performance Assessments Over 20 Years - 13597

Performance assessments (PAs) have been used for many years for the analysis of post-closure hazards associated with a radioactive waste disposal facility and to provide a reasonable expectation of the ability of the site and facility design to meet objectives for the protection of members of the public and the environment. The use of PA to support decision-making for LLW disposal facilities has been mandated in United States Department of Energy (US DOE) directives governing radioactive waste management since 1988 (currently DOE Order 435.1, Radioactive Waste Management). Prior to that time, PAs were also used in a less formal role. Over the past 20+ years, the US DOE approach to conduct, review and apply PAs has evolved into an efficient, rigorous and mature process that includes specific requirements for continuous improvement and independent reviews. The PA process has evolved through refinement of a graded and iterative approach designed to help focus efforts on those aspects of the problem expected to have the greatest influence on the decision being made. Many of the evolutionary changes to the PA process are linked to the refinement of the PA maintenance concept that has proven to be an important element of US DOE PA requirements in the context of supporting decision-making for safe disposal of LLW. The PA maintenance concept is central to the evolution of the graded and iterative philosophy and has helped to drive the evolution of PAs from a deterministic compliance calculation into a systematic approach that helps to focus on critical aspects of the disposal system in a manner designed to provide a more informed basis for decision-making throughout the life of a disposal facility (e.g., monitoring, research and testing, waste acceptance criteria, design improvements, data collection, model refinements). A significant evolution in PA modeling has been associated with improved use of uncertainty and sensitivity analysis techniques to support efficient implementation of the graded and iterative approach. Rather than attempt to exactly predict the migration of radionuclides in a disposal unit, the best PAs have evolved into tools that provide a range of results to guide decision-makers in planning the most efficient, cost effective, and safe disposal of radionuclides. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2012

Improvement design study on steam generator of MHR-50/100 aiming higher safety level after water ingress accident

Mitsubishi Heavy Industries, Ltd. (MHI) has been studying on MHI original High Temperature Gas cooled Reactor (HTGR), namely MHR-50/100, for commercialization with supported by JAEA. In the heat transfer system, steam generator (SG) is one of the most important components because it should be imposed a function of heat transfer from reactor power to steam turbine system and maintaining a nuclear grade boundary. Then we especially focused an effort of a design study on the SG having robustness against water ingress accident based on our design experience of PWR, FBR and HTGR. In this study, we carried out a sensitivity analysis from the view point of economic and plant efficiency. As a result, the SG design parameter of helium inlet/outlet temperature of 750 deg. C/300 deg. C, a side-by-side layout and one unit of SG attached to a reactor were selected. In the next, a design improvement of SG was carried out from the view point of securing the level of inherent safety without reliance on active steam dump system during water ingress accident considering the situation of the Fukushima nuclear power plant disaster on March 11, 2011. Finally, according to above basic design requirement to SG, we performed a conceptual design on adapting themes of SG structure improvement. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Computational assessment of smooth and rough parameter dependence of statistics in chaotic dynamical systems

An assumption of smooth response to small parameter changes, of statistics or long-time averages of a chaotic system, is generally made in the field of sensitivity analysis, and the parametric derivatives of statistical quantities are critically used in science and engineering. In this paper, we propose a numerical procedure to assess the differentiability of statistics with respect to parameters in chaotic systems. Here, we numerically show that the existence of the derivative depends on the Lebesgue integrability of a certain density gradient function, which we define as the derivative of logarithmic SRB density along the unstable manifold. We develop a recursive formula for the density gradient that can be efficiently computed along trajectories, and demonstrate its use in determining the differentiability of statistics. Our numerical procedure is illustrated on low-dimensional chaotic systems whose statistics exhibit both smooth and rough regions in parameter space.

42 ENGINEERING↗