SEARCH · Engineering Papers

Results for “sensitivity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant at or Near Delayed Critical [Abstract]

Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

The sensitivity of the forest carbon budget shifts across processes along with stand development and climate change

Abstract The future trajectory of atmospheric CO 2 concentration depends on the development of the terrestrial carbon sink, which in turn is influenced by forest dynamics under changing environmental conditions. An in‐depth understanding of model sensitivities and uncertainties in non‐steady‐state conditions is necessary for reliable and robust projections of forest development and under scenarios of global warming and CO 2 enrichment. Here, we systematically assessed if a biogeochemical process‐based model (3D‐CMCC‐CNR), which embeds similarities with many other vegetation models, applied in simulating net primary productivity (NPP) and standing woody biomass (SWB), maintained a consistent sensitivity to its 55 input parameters through time, during forest ageing and structuring as well as under climate change scenarios. Overall, the model applied at three contrasting European forests showed low sensitivity to the majority of its parameters. Interestingly, model sensitivity to parameters varied through the course of >100 yr of simulations. In particular, the model showed a large responsiveness to the allometric parameters used for initialize forest carbon and nitrogen pools early in forest simulation (i.e., for NPP up to ~37%, 256 g C·m −2 ·yr −1 and for SWB up to ~90%, 65 Mg C/ha, when compared to standard simulation), with this sensitivity decreasing sharply during forest development. At medium to longer time scales, and under climate change scenarios, the model became increasingly more sensitive to additional and/or different parameters controlling biomass accumulation and autotrophic respiration (i.e., for NPP up to ~30%, 167 g C·m −2 ·yr −1 and for SWB up to ~24%, 64 Mg C/ha, when compared to standard simulation). Interestingly, model outputs were shown to be more sensitive to parameters and processes controlling stand development rather than to climate change (i.e., warming and changes in atmospheric CO 2 concentration) itself although model sensitivities were generally higher under climate change scenarios. Our results suggest the need for sensitivity and uncertainty analyses that cover multiple temporal scales along forest developmental stages to better assess the potential of future forests to act as a global terrestrial carbon sink.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Regional climate change projections from NA-CORDEX and their relation to climate sensitivity

Abstract The climate sensitivity of global climate models (GCMs) strongly influences projected climate change due to increased atmospheric carbon dioxide. Reasonably, the climate sensitivity of a GCM may be expected to affect dynamically downscaled projections. However, there has been little examination of the effect of the climate sensitivity of GCMs on regional climate model (RCM) ensembles. Therefore, we present projections of temperature and precipitation from the ensemble of projections produced as a part of the North American branch of the international Coordinated Regional Downscaling Experiment (NA-CORDEX) in the context of their relationship to the climate sensitivity of their parent GCMs. NA-CORDEX simulations were produced at 50-km and 25-km resolutions with multiple RCMs which downscaled multiple GCMs that spanned nearly the full range of climate sensitivity available in the CMIP5 archive. We show that climate sensitivity is a very important source of spread in the NA-CORDEX ensemble, particularly for temperature. Temperature projections correlate with driving GCM climate sensitivity annually and seasonally across North America not only at a continental scale but also at a local-to-regional scale. Importantly, the spread in temperature projections would be reduced if only low, mid, or high climate sensitivity simulations were considered, or if only the ensemble mean were considered. Precipitation projections correlate with climate sensitivity, but only at a continental scale during the cold season, due to the increasing influence of other processes at finer scales. Additionally, it is shown that the RCMs do alter the projection space sampled by their driving GCMs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

High-density thermal sensitivity maps of the human body

‘Personal comfort systems’ and thermally active clothing are able to warm and cool individual building occupants by transferring heat directly to and from their body surfaces. Such systems would ideally target local body surfaces with high temperature sensitivities. Such sensitivities have not been quantified in detail before. Here we report local thermal sensations and sensitivities for 318 local skin spots distributed over one side of the body, measured on a large number of subjects. Skin temperature changes were induced with a thermal probe 14 mm in diameter, and subjective thermal sensations were surveyed after 10 s. Our neutral base temperature was 31 °C and the spot stimulus was ±5 °C. Cool and warm sensitivities are seen to vary widely by body part. The foot, lower leg and upper chest are much less sensitive than average; in comparison, the cheek, neck back, and seat area are 2–3 times as sensitive to both cooling and warming stimuli. Every body part exhibits stronger sensitivity to cooling (1.3–1.6 times stronger) than to warming. Inter-personal differences and regional variance within body parts were observed to be 2–3 times greater than potential sex differences. These high-density thermal sensitivity maps with appended dataset provide the most comprehensive distributions of cold and warm sensitivity across the human body.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Regular sensitivity computation avoiding chaotic effects in particle-in-cell plasma methods

Particle-in-cell (PIC) simulation methods are attractive for representing species distribution functions in plasmas. However, as a model, they introduce uncertain parameters, and for quantifying their prediction uncertainty it is useful to be able to assess the sensitivity of a quantity-of-interest (QoI) to these parameters. Such sensitivity information is likewise useful for optimization. However, computing sensitivity for PIC methods is challenging due to the chaotic particle dynamics, and sensitivity techniques remain underdeveloped compared to those for Eulerian discretizations. This challenge is examined from a dual particle–continuum perspective that motivates a new sensitivity discretization. Two routes to sensitivity computation are presented and compared: a direct fully-Lagrangian particle-exact approach provides sensitivities of each particle trajectory, and a new particle-pdf discretization, which is formulated from a continuum perspective but discretized by particles to take the advantages of the same type of Lagrangian particle description leveraged by PIC methods. Since the sensitivity particles in this approach are only indirectly linked to the plasma-PIC particles, they can be positioned and weighted independently for efficiency and accuracy. The corresponding numerical algorithms are presented here in mathematical detail. The advantage of the particle-pdf approach in avoiding the spurious chaotic sensitivity of the particle-exact approach is demonstrated for Debye shielding and sheath configurations. In essence, the continuum perspective makes implicit the distinctness of the particles, which circumvents the Lyapunov instability of the N -body PIC system. The cost of the particle-pdf approach is comparable to the baseline PIC simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

SU-E-T-249: Determining the Sensitivity of Beam Profile Parameters for Detecting Energy Changes in Flattening Filter-Free Beams

Purpose: To identify the beam profile parameters that can be used to detect energy changes in a flattening filter-free photon beams. Methods: Flattening filter-free beam profiles (inline, crossline, and diagonals) were measured for multiple field sizes (25×25cm and 10×10cm) at 6MV on a clinical system (Truebeam, Varian Medical Systems Palo Alto CA). Profiles were acquired for baseline energy and detuned beams by changing the bending magnet current (BMC), above and below baseline. The following profile parameters were measured: flatness (off-axis ratio at 80% of field size), symmetry, uniformity, slope, and the off-axis ratio (OAR) at several off-axis distances. Tolerance values were determined from repeated measurements. Each parameter was evaluated for sensitivity to the induced beam changes, and the minimum detectable BMC change was calculated for each parameter by calculating the change in BMC that would Result in a change in the parameter above the measurement tolerance. Results: Tolerance values for the parameters were-Flatness≤0.1%; Symmetry≤0.4%; Uniformity≤0.01%; Slope≤ 0.001%/mm. The measurements made with a field size of 25cm and a depth of d=1.5cm showed the greatest sensitivity to bending magnet current variations. Uniformity had the highest sensitivity, able to detect a change in BMC of BMC=0.02A. The OARs and slope were sensitive to the magnitude and direction of BMC change. The sensitivity in the flatness parameter was BMC=0.04A; slope was sensitive to BMC=0.05A. The sensitivity decreased for OARs measured closer to central axis-BMC(8cm)=0.23A; BMC(5cm)=0.47A; BMC(2cm)=1.35A. Symmetry was not sensitive to changes in BMC. Conclusion: These tests allow for better QA of FFF beams by setting tolerance levels to beam parameter baseline values which reflect variations in machine calibration. Uniformity is most sensitive to BMC changes, while OARs provide information about magnitude and direction of miscalibration. Research funding provided by Varian Medical Systems. Dr. Sasa Mutic receives compensation for providing patient safety training services from Varian Medical Systems, the sponsor of this study.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

IMRT QA: Selecting gamma criteria based on error detection sensitivity

Purpose: The gamma comparison is widely used to evaluate the agreement between measurements and treatment planning system calculations in patient-specific intensity modulated radiation therapy (IMRT) quality assurance (QA). However, recent publications have raised concerns about the lack of sensitivity when employing commonly used gamma criteria. Understanding the actual sensitivity of a wide range of different gamma criteria may allow the definition of more meaningful gamma criteria and tolerance limits in IMRT QA. We present a method that allows the quantitative determination of gamma criteria sensitivity to induced errors which can be applied to any unique combination of device, delivery technique, and software utilized in a specific clinic. Methods: A total of 21 DMLC IMRT QA measurements (ArcCHECK®, Sun Nuclear) were compared to QA plan calculations with induced errors. Three scenarios were studied: MU errors, multi-leaf collimator (MLC) errors, and the sensitivity of the gamma comparison to changes in penumbra width. Gamma comparisons were performed between measurements and error-induced calculations using a wide range of gamma criteria, resulting in a total of over 20 000 gamma comparisons. Gamma passing rates for each error class and case were graphed against error magnitude to create error curves in order to represent the range of missed errors in routine IMRT QA using 36 different gamma criteria. Results: This study demonstrates that systematic errors and case-specific errors can be detected by the error curve analysis. Depending on the location of the error curve peak (e.g., not centered about zero), 3%/3 mm threshold = 10% at 90% pixels passing may miss errors as large as 15% MU errors and ±1 cm random MLC errors for some cases. As the dose threshold parameter was increased for a given %Diff/distance-to-agreement (DTA) setting, error sensitivity was increased by up to a factor of two for select cases. This increased sensitivity with increasing dose threshold was consistent across all studied combinations of %Diff/DTA. Criteria such as 2%/3 mm and 3%/2 mm with a 50% threshold at 90% pixels passing are shown to be more appropriately sensitive without being overly strict. However, a broadening of the penumbra by as much as 5 mm in the beam configuration was difficult to detect with commonly used criteria, as well as with the previously mentioned criteria utilizing a threshold of 50%. Conclusions: We have introduced the error curve method, an analysis technique which allows the quantitative determination of gamma criteria sensitivity to induced errors. The application of the error curve method using DMLC IMRT plans measured on the ArcCHECK® device demonstrated that large errors can potentially be missed in IMRT QA with commonly used gamma criteria (e.g., 3%/3 mm, threshold = 10%, 90% pixels passing). Additionally, increasing the dose threshold value can offer dramatic increases in error sensitivity. This approach may allow the selection of more meaningful gamma criteria for IMRT QA and is straightforward to apply to other combinations of devices and treatment techniques.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SCALE Continuous-Energy Eigenvalue Sensitivity Coefficient Calculations

Sensitivity coefficients describe the fractional change in a system response that is induced by changes to system parameters and nuclear data. The Tools for Sensitivity and UNcertainty Analysis Methodology Implementation (TSUNAMI) code within the SCALE code system makes use of eigenvalue sensitivity coefficients for an extensive number of criticality safety applications, including quantifying the data-induced uncertainty in the eigenvalue of critical systems, assessing the neutronic similarity between different critical systems, and guiding nuclear data adjustment studies. The need to model geometrically complex systems with improved fidelity and the desire to extend TSUNAMI analysis to advanced applications has motivated the development of a methodology for calculating sensitivity coefficients in continuous-energy (CE) Monte Carlo applications. The Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) and Iterated Fission Probability (IFP) eigenvalue sensitivity methods were recently implemented in the CE-KENO framework of the SCALE code system to enable TSUNAMI-3D to perform eigenvalue sensitivity calculations using continuous-energy Monte Carlo methods. This work provides a detailed description of the theory behind the CLUTCH method and describes in detail its implementation. This work explores the improvements in eigenvalue sensitivity coefficient accuracy that can be gained through the use of continuous-energy sensitivity methods and also compares several sensitivity methods in terms of computational efficiency and memory requirements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2014

Climate Sensitivity in the Anthropocene

Climate sensitivity in its most basic form is defined as the equilibrium change in global surface temperature that occurs in response to a climate forcing, or externally imposed perturbation of the planetary energy balance. Within this general definition, several specific forms of climate sensitivity exist that differ in terms of the types of climate feedbacks they include. Based on evidence from Earth's history, we suggest here that the relevant form of climate sensitivity in the Anthropocene (e.g. from which to base future greenhouse gas (GHG) stabilization targets) is the Earth system sensitivity including fast feedbacks from changes in water vapour, natural aerosols, clouds and sea ice, slower surface albedo feedbacks from changes in continental ice sheets and vegetation, and climate-GHG feedbacks from changes in natural (land and ocean) carbon sinks. Traditionally, only fast feedbacks have been considered (with the other feedbacks either ignored or treated as forcing), which has led to estimates of the climate sensitivity for doubled CO2 concentrations of about 3 C. The 2×CO2 Earth system sensitivity is higher than this, being approx. 4-6 C if the ice sheet/vegetation albedo feedback is included in addition to the fast feedbacks, and higher still if climate-GHG feedbacks are also included. The inclusion of climate-GHG feedbacks due to changes in the natural carbon sinks has the advantage of more directly linking anthropogenic GHG emissions with the ensuing global temperature increase, thus providing a truer indication of the climate sensitivity to human perturbations. The Earth system climate sensitivity is difficult to quantify due to the lack of palaeo-analogues for the present-day anthropogenic forcing, and the fact that ice sheet and climate-GHG feedbacks have yet to become globally significant in the Anthropocene. Furthermore, current models are unable to adequately simulate the physics of ice sheet decay and certain aspects of the natural carbon and nitrogen cycles. Obtaining quantitative estimates of the Earth system sensitivity is therefore a high priority for future work.

Anthropocene↗

DOE OSTI2016

Adjoint-Based Sensitivity and Uncertainty Analysis for Density and Composition: A User's Guide

The ability to perform sensitivity analyses using adjoint-based sensitivity theory has existed for decades . This paper provides guidance on how adjoint sensitivity methods can be used to predict the effect of material density and composition uncertainties in critical experiments, including when these uncertain parameters are correlated or constrained. Two widely used Monte Carlo codes, MCNP6 and SCALE, are both capable of computing isotopic sensitivities in continuous energy and angle. Additionally, Perko et al. have shown how individual isotope sensitivities, easily computed using adjoint methods, can be combined to compute constrained sensitivities that may be used in the uncertainty analysis. This paper provides details on how the codes are used to compute sensitivities and how the sensitivities are used in an uncertainty analysis. Constrained sensitivities are computed in a simple example problem. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Silencing of secretory clusterin sensitizes NSCLC cells to V-ATPase inhibitors by downregulating survivin

Highlights: • V-ATPase inhibitors stimulate sCLU protein expression. • Knockdown of sCLU enhances NSCLC cell sensitivity to V-ATPase inhibitors. • The enhanced cell sensitivity is associated with decreased survivin protein. • Inhibition of PI3K/AKT/mTOR sensitizes NSCLC cells to V-ATPase inhibitors. Secretory clusterin (sCLU) is a stress-associated protein that confers resistance to therapy when overexpressed. In this study, we observed that the V-ATPase inhibitors bafilomycin A1 and concanamycin A significantly stimulated sCLU protein expression. Knockdown of sCLU with siRNA sensitized non-small cell lung cancer (NSCLC) cells to bafilomycin A1, suggesting that sCLU expression renders cells resistant to V-ATPase inhibitors. The dual PI3K/AKT and mTOR inhibitor BEZ235 suppressed sCLU expression and enhanced cell sensitivity induced by bafilomycin A1. Notably, sCLU knockdown further decreased the expression of the survivin protein by bafilomycin A1, and the ectopic expression of survivin alleviated the cell sensitivity by bafilomycin A1 and sCLU depletion, suggesting that increased sensitivity to sCLU depletion in the cells with V-ATPase inhibitors is due, at least in part, to the down-regulation of survivin. Taken together, we demonstrated that the depletion of sCLU expression enhances the sensitivity of NSCLC cells to V-ATPase inhibitors by decreasing survivin expression. Inhibition of the PI3K/AKT/mTOR pathway enhances the sensitivity of NSCLC cells to V-ATPase inhibitors, leading to decreased sCLU and survivin expression. Thus, we suggest that a combination of PI3K/AKT/mTOR inhibitors with V-ATPase inhibitors might be an effective approach for NSCLC treatment.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

SU-E-T-291: Sensitivity of a Simple 2D EPID in Vivo Dosimetry

Purpose: As radiotherapy (RT) increases in complexity, so does motivation for in vivo dosimetry (IVD), which may detect errors such as: setup, beam shaping and dose delivered. We have recently developed an easy-toimplement method for two-dimensional IVD based on images taken with the electronic portal imaging device (EPID) in cine mode during treatment. The purpose of this work is to characterize its sensitivity to possible RT delivery errors. Methods: We introduced a series of modifications to a simple RT field (10×10, 100MU, 300RR, 20cm homogeneous phantom) to simulate errors. These modifications included multi-leaf collimator (MLC) position, number of MUs, and collimator angle. We quantified the sensitivity to inhomogeneities by inserting variable amounts of solid lung and bone. Finally we delivered realistic fields to an anthropomorphic phantom to estimate sensitivity to gantry angle and setup errors. Results: Our EPIDIVD is sensitive to MLC positioning errors of 1mm and 3mm in the closed and open directions respectively, and to 3% MU variations. Sensitivity to collimator angle depends on field shape irregularity; in the case of a 10x10 field, we are sensitive to errors of 0.8°. The sensitivity to inhomogeneities is limited by the nature of MV imaging: approximately 1% signal change is noted when switching 5cm of water to equal amounts of bone or lung. This suggests that the EPID-IVD is likely not sensitive to small setup or gantry angle errors, as confirmed by anthropomorphic tests. Conclusion: We have characterized a simple method of 2D dose reconstruction at isocenter depth inside the patient, which is sensitive to possible RT delivery errors. This method may be useful as a secondary safety check, to prevent large errors from being carried on to following fractions, and to record delivered dose. By using readily available hardware, it is easily implemented and may prove especially useful in centers with limited resources.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Validation of Jezebel Reactivity Coefficients and Sensitivity Analysis

Nuclear data validation is often performed today using criticality measurements. The gold standard for criticality measurements is the International Criticality Safety Benchmark Experiment Project (ICSBEP). The validation specifically focuses on the effective multiplication factor ($k_{eff}$). $K_{eff}$ is a relatively easy parameter to infer and has reduced uncertainty due to being at or above critical. However, while $k_{eff}$ is the most documented parameter and its uncertainties and sensitivities have been evaluated in great detail, it cannot be used as a standalone metric to determine inaccuracies in nuclear data (e.g., cross section data, PFNS, nu), which is based on theory, physics, and differential measurements. The Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) project aims to identify compensating errors in specific isotope nuclear data by optimally designing experiments that are, or are not sensitive to a suite of measurement parameters beyond $k_{eff}$. By identifying parameters that are sensitive to each other, oppositely sensitive, or have substantial magnitude differences in sensitivity, experiments can be designed to constrain questionable nuclear data. One sensitivity that is of particular interest to this project includes the sensitivity of reactivity coefficients. Reactivity coefficients compare reactivity, which is related to $k_{eff}$ at two different states therefore being sensitive to small changes in the system. The most common type of reactivity coefficient measurements is comparison to void for a small sample within the assembly. It is key that the sample sizes are small enough to not affect the flux of the full system. Reactivity coefficients were evaluated for many early experiments to better understand transport corrected cross sections. In fact, ICSBEP includes reactivity coefficient results as “Supplemental Measurements” in appendices for a handful of older benchmarks. One of those benchmarks is Jezebel, the bare Pu critical assembly. This paper compares new simulations of reactivity coefficients for Jezebel, and explores the sensitivity of reactivity coefficients to small changes in nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

(U) SENSMG: First-Order Sensitivities of Neutron Reaction Rates, Reaction-Rate Ratios, Leakage, k eff , α , and Subcritical Multiplication Using PARTISN

SENSMG is a tool for computing first-order sensitivities of neutron reaction rates, reaction-rate ratios, leakage, k eff , α, and subcritical multiplication using the PARTISN multigroup discrete-ordinates code. SENSMG computes sensitivities to all of the transport cross sections and data (total, fission, a nu, chi, and all scattering moments), two edit cross sections (absorption and capture), and the density for every nuclide and energy group. It also computes sensitivities to the mass density for every material and derivatives with respect to all interface locations and outer boundaries. It computes sensitivities to user specified reactions whose cross sections are available in a user-supplied NJOY output file. The tool can be used for one-dimensional spherical and slab (r) and two-dimensional cylindrical (r-z) geometries. The tool can be used for fixed-source and eigenvalue problems. For most responses, the tool implements Generalized Perturbation Theory (GPT) as discussed by Williams and Stacey. The tool is thus limited to computing sensitivities only for GPT-allowable responses. For subcritical multiplication, the tool implements sensitivities derived by O’Brien and Clark. SENSMG has a similar role as the old SWANLAKE (Ref. 8), FORSS (Ref. 9), and SENSIT (Ref. 10) codes. It has capabilities similar to those of SUSD3D (Refs. 11 and 12), which also uses PARTISN. Section II of this report describes the theory behind adjoint-based sensitivities, gives the equations that SENSMG solves, and defines the sensitivities that are output. Section III describes the user interface, including the input file and command line options. Section IV describes the output. Section V gives some notes about the coding that may be of interest. Section VI presents some sample problems and discusses verification, which is ongoing. Section VII lists needs and ideas for future work. Appendix A lists most of the input files whose results are presented in Sec. VI. Appendix B provides some useful details on one of the cross-section libraries that SENSMG supports.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

Numerical modeling and sensitivity analysis of seawater intrusion in a dual-permeability coastal karst aquifer with conduit networks

Long-distance seawater intrusion has been widely observed through the subsurface conduit system in coastal karst aquifers as a source of groundwater contaminant. In this study, seawater intrusion in a dual-permeability karst aquifer with conduit networks is studied by the two-dimensional density-dependent flow and transport SEAWAT model. Local and global sensitivity analyses are used to evaluate the impacts of boundary conditions and hydrological characteristics on modeling seawater intrusion in a karst aquifer, including hydraulic conductivity, effective porosity, specific storage, and dispersivity of the conduit network and of the porous medium. The local sensitivity analysis evaluates the parameters' sensitivities for modeling seawater intrusion, specifically in the Woodville Karst Plain (WKP). A more comprehensive interpretation of parameter sensitivities, including the nonlinear relationship between simulations and parameters, and/or parameter interactions, is addressed in the global sensitivity analysis. The conduit parameters and boundary conditions are important to the simulations in the porous medium because of the dynamical exchanges between the two systems. The sensitivity study indicates that salinity and head simulations in the karst features, such as the conduit system and submarine springs, are critical for understanding seawater intrusion in a coastal karst aquifer. The evaluation of hydraulic conductivity sensitivity in the continuum SEAWAT model may be biased since the conduit flow velocity is not accurately calculated by Darcy's equation as a function of head difference and hydraulic conductivity. In addition, dispersivity is no longer an important parameter in an advection-dominated karst aquifer with a conduit system, compared to the sensitivity results in a porous medium aquifer. In the end, the extents of seawater intrusion are quantitatively evaluated and measured under different scenarios with the variabilities of important parameters identified from sensitivity results, including salinity at the submarine spring with rainfall recharge, sea level rise, and a longer simulation time under an extended low rainfall period.

58 GEOSCIENCES↗

NASA NTRS1988

Results of an integrated structure-control law design sensitivity analysis

Next generation air and space vehicle designs are driven by increased performance requirements, demanding a high level of design integration between traditionally separate design disciplines. Interdisciplinary analysis capabilities have been developed, for aeroservoelastic aircraft and large flexible spacecraft control for instance, but the requisite integrated design methods are only beginning to be developed. One integrated design method which has received attention is based on hierarchal problem decompositions, optimization, and design sensitivity analyses. This paper highlights a design sensitivity analysis method for Linear Quadratic Cost, Gaussian (LQG) optimal control laws, which predicts change in the optimal control law due to changes in fixed problem parameters using analytical sensitivity equations. Numerical results of a design sensitivity analysis for a realistic aeroservoelastic aircraft example are presented. In this example, the sensitivity of the optimally controlled aircraft's response to various problem formulation and physical aircraft parameters is determined. These results are used to predict the aircraft's new optimally controlled response if the parameter was to have some other nominal value during the control law design process. The sensitivity results are validated by recomputing the optimal control law for discrete variations in parameters, computing the new actual aircraft response, and comparing with the predicted response. These results show an improvement in sensitivity accuracy for integrated design purposes over methods which do not include changess in the optimal control law. Use of the analytical LQG sensitivity expressions is also shown to be more efficient that finite difference methods for the computation of the equivalent sensitivity information.

Gilbert, Michael G.↗

NASA NTRS1994

The sensitivity to parametric variation in direct minimization techniques

Solutions of some objective analysis techniques are known to depend upon the subjective values of internal parameters. The change in the solution per change in the parameter is the sensitivity. Parameters with low sensitivity can be varied with large increments during preliminary searches for near-optimal parameter values. Only terms with high sensitivity must be thoroughly investigated once the parameters are determined to be close to optimal. Both absolute and relative sensitivities are discussed and a sensitivity-based definition of the solution uncertainty is proposed. The sensitivity of direct minimization analysis to parametric variation is evaluated using a set of 'response functions' that characterize different aspects of the solution. It is shown that solutions of direct minimization techniques have low absolute sensitivity. Two examples are used to illustrate the usefulness of the technique. Both involved measurements of air-sea quantities (e.g., wind stress and latent heat flux) from a variety of data sources using a direct minimization technique. The examples demonstrate that sensitivity analysis is capable of quantifying regional sensitivities as well as indicating the magnitude and relationship between the various parameters.

Meyers, S. D.↗

DOE OSTI2015

Improved Sensitivity Relations in State Constrained Optimal Control

Sensitivity relations in optimal control provide an interpretation of the costate trajectory and the Hamiltonian, evaluated along an optimal trajectory, in terms of gradients of the value function. While sensitivity relations are a straightforward consequence of standard transversality conditions for state constraint free optimal control problems formulated in terms of control-dependent differential equations with smooth data, their verification for problems with either pathwise state constraints, nonsmooth data, or for problems where the dynamic constraint takes the form of a differential inclusion, requires careful analysis. In this paper we establish validity of both ‘full’ and ‘partial’ sensitivity relations for an adjoint state of the maximum principle, for optimal control problems with pathwise state constraints, where the underlying control system is described by a differential inclusion. The partial sensitivity relation interprets the costate in terms of partial Clarke subgradients of the value function with respect to the state variable, while the full sensitivity relation interprets the couple, comprising the costate and Hamiltonian, as the Clarke subgradient of the value function with respect to both time and state variables. These relations are distinct because, for nonsmooth data, the partial Clarke subdifferential does not coincide with the projection of the (full) Clarke subdifferential on the relevant coordinate space. We show for the first time (even for problems without state constraints) that a costate trajectory can be chosen to satisfy the partial and full sensitivity relations simultaneously. The partial sensitivity relation in this paper is new for state constraint problems, while the full sensitivity relation improves on earlier results in the literature (for optimal control problems formulated in terms of Lipschitz continuous multifunctions), because a less restrictive inward pointing hypothesis is invoked in the proof, and because it is validated for a stronger set of necessary conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗