Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Response surface model”

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 289 records · Page 16

Ankle and hip postural strategies defined by joint torques

Previous studies have identified two discrete strategies for the control of posture in the sagittal plane based on EMG activations, body kinematics, and ground reaction forces. The ankle strategy was characterized by body sway resembling a single-segment-inverted pendulum and was elicited on flat support surfaces. In contrast, the hip strategy was characterized by body sway resembling a double-segment inverted pendulum divided at the hip and was elicited on short or compliant support surfaces. However, biomechanical optimization models have suggested that hip strategy should be observed in response to fast translations on a flat surface also, provided the feet are constrained to remain in contact with the floor and the knee is constrained to remain straight. The purpose of this study was to examine the experimental evidence for hip strategy in postural responses to backward translations of a flat support surface and to determine whether analyses of joint torques would provide evidence for two separate postural strategies. Normal subjects standing on a flat support surface were translated backward with a range of velocities from fast (55 cm/s) to slow (5 cm/s). EMG activations and joint kinematics showed pattern changes consistent with previous experimental descriptions of mixed hip and ankle strategy with increasing platform velocity. Joint torque analyses revealed the addition of a hip flexor torque to the ankle plantarflexor torque during fast translations. This finding indicates the addition of hip strategy to ankle strategy to produce a continuum of postural responses. Hip torque without accompanying ankle torque (pure hip strategy) was not observed. Although postural control strategies have previously been defined by how the body moves, we conclude that joint torques, which indicate how body movements are produced, are useful in defining postural control strategies. These results also illustrate how the biomechanics of the body can transform discrete control patterns into a continuum of postural corrections.

Non-NASA Center↗

Modeled Impact of Cirrus Cloud Increases Along Aircraft Flight Paths

The potential impact of contrails and alterations in the lifetime of background cirrus due to subsonic airplane water and aerosol emissions has been investigated in a set of experiments using the GISS GCM connected to a q-flux ocean. Cirrus clouds at a height of 12-15km, with an optical thickness of 0.33, were input to the model "x" percentage of clear-sky occasions along subsonic aircraft flight paths, where x is varied from .05% to 6%. Two types of experiments were performed: one with the percentage cirrus cloud increase independent of flight density, as long as a certain minimum density was exceeded; the other with the percentage related to the density of fuel expenditure. The overall climate impact was similar with the two approaches, due to the feedbacks of the climate system. Fifty years were run for eight such experiments, with the following conclusions based on the stable results from years 30-50 for each. The experiments show that adding cirrus to the upper troposphere results in a stabilization of the atmosphere, which leads to some decrease in cloud cover at levels below the insertion altitude. Considering then the total effect on upper level cloud cover (above 5 km altitude), the equilibrium global mean temperature response shows that altering high level clouds by 1% changes the global mean temperature by 0.43C. The response is highly linear (linear correlation coefficient of 0.996) for high cloud cover changes between 0. 1% and 5%. The effect is amplified in the Northern Hemisphere, more so with greater cloud cover change. The temperature effect maximizes around 10 km (at greater than 40C warming with a 4.8% increase in upper level clouds), again more so with greater warming. The high cloud cover change shows the flight path influence most clearly with the smallest warming magnitudes; with greater warming, the model feedbacks introduce a strong tropical response. Similarly, the surface temperature response is dominated by the feedbacks, and shows little geographical relationship to the high cloud input. Considering whether these effects would be observable, changing upper level cloud cover by as little as 0.4% produces warming greater than 2 standard deviations in the Microwave Sounding Unit (MSU) channels 4, 2 and 2r, in flight path regions and in the subtropics. Despite the simplified nature of these experiments, the results emphasize the sensitivity of the modeled climate to high level cloud cover changes, and thus the potential ability of aircraft to influence climate by altering clouds in the upper troposphere.

Rind, David↗

A singularity free approach to post glacial rebound calculations

Calculating the post glacial response of a viscoelastic Earth model using the exponential decay normal mode technique leads to intrinsic singularities if viscosity varies continuously as a function of radius. We develop a numerical implementation of the Complex Real Fourier transform (CRFT) method as an accurate and stable procedure to avoid these singularities. Using CRFT, we investigate the response of a set of Maxwell Earth models to surface loading. We find that the effect of expanding a layered viscosity structure into a continuously varying structure is to destroy the modes associated with the boundary between layers. Horizontal motion is more sensitive than vertical motion to the viscosity structure just below the lithosphere. Horizontal motion is less sensitive to the viscosity of the lower mantle than the vertical motion is. When the viscosity increases at 670 km depth by a factor of about 60, the response of the lower mantle is close to its elastic limit. Any further increase of the viscosity contrast at 670 km depth or further increase of viscosity as a continuous function of depth starting from 670 km depth is unlikely to be resolved.

Fang, Ming↗

Data-model files associated with the manuscript "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS)

This data package contains the model inputs and outputs used in "The Effects of Spatial and Temporal Resolution of Gridded Meteorological Forcing on Watershed Hydrological Responses" (Shuai et al., 2022 HESS). The data.zip file contains the data used to drive the model simulations. The model.zip file contains the XML input file for ATS. The notebook.zip file contains the Jupyter notebooks for pre- and post- processing model results. The figures.zip file contains the raw figures associated with the manuscript.Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (10s of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow, and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) -- Daily Surface Weather and Climatological Summaries (Daymet), Parameter-elevation Regressions on Independent Slopes Model (PRISM), and North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30% of the watershed area). However, the 4 km resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications on model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

Deficient precipitation sensitivity to Sahel land surface forcings among CMIP5 models

Abstract The overall performance of the simulated seasonal precipitation response to local terrestrial forcings, namely vegetation abundance and soil moisture, in the Sahel among the Coupled Model Intercomparison Project Phase Five (CMIP5) Earth System Models (ESMs) is systematically investigated and compared with its observational counterpart using a multivariate statistical method. The observed seasonal precipitation response is evaluated against a large ensemble of observational, reanalysis, and satellite data sets to provide quantification of uncertainties. The behaviour of models with and without a Dynamic Global Vegetation Model (DGVM) component is also explored, along with the mechanisms responsible for terrestrial feedback on rainfall. In general, the CMIP5 models can reasonably capture the seasonal evolution of Sahel precipitation and soil moisture, albeit with wet biases during the pre‐monsoon period and dry biases during the peak monsoon period. The non‐DGVM ESMs simulate comparable leaf area indices (LAIs) with observations, while DGVM‐enabled ESMs simulate too much year‐round LAI. The variance of precipitation that is attributed to oceanic forcings in CMIP5 is comparable with observations; however, the variance of precipitation that is attributed to terrestrial forcings is smaller in CMIP5 models than observed, especially for non‐DGVM ESMs. CMIP5 models, especially those without DGVMs, undervalue precipitation's observed response strength to soil moisture anomalies. In both observations and CMIP5 models, none of the atmospheric variables show significant responses to direct vegetation forcing, except for the response in transpiration. Although vegetation has minimal direct effect on the atmospheric state, it can affect the atmosphere by modifying soil moisture and transpiration rate indirectly, which helps explain the more realistic simulation of rainfall in DGVM‐enabled ESMs than non‐DGVM ESMs. Coupling of an ESM to a DGVM is critical in generating reasonable land–atmosphere feedback and examining future ecological and climatic changes over the Sahel.

54 ENVIRONMENTAL SCIENCES↗

Root Foraging Alters Global Patterns of Ecosystem Legacy From Climate Perturbations

The response of terrestrial ecosystems to climate perturbations typically persist longer than the timescale of the forcing, a phenomenon broadly referred to as legacy. Understanding the strength of legacy is critical for predicting ecosystem sensitivity to climate extremes and the extent that persistent changes in surface-atmosphere exchange might feedback onto climate. The cause of ecosystem legacy has been associated with myriad factors such as changes in aboveground biomass, however, few studies have tested how changes in the depth distribution of fine roots in response to perturbation might alter an ecosystem's recovery time. We explore this question using an Earth System Model that includes a dynamic root module where vegetation can forage for water and nutrients by altering their root profiles. The global simulations presented here show that in response to water stress events most ecosystems deepen their root profiles. In semi-arid ecosystems, this response expedites recovery (i.e., less legacy) relative to simulations without dynamics roots because access to deeper water pools after the initial event remains favorable. In wetter ecosystems, the development of deeper root profiles slows down the recovery timescale (i.e., more legacy) because the deeper root profile reduces access to nutrients and is thus unfavorable. Here, the simulations show that while the inclusion of dynamic roots might only minimally affect global patterns of Gross Primary Productivity and Evapotranspiration, the shift in root profile alters the timescale of recovery. Studies interested in the sustained response of land surfaces fluxes to climate disturbances should consider models that include dynamic root capability.

54 ENVIRONMENTAL SCIENCES↗

The effects of spatial and temporal resolution of gridded meteorological forcing on watershed hydrological responses

Abstract. Meteorological forcing plays a critical role in accurately simulating the watershed hydrological cycle. With the advancement of high-performance computing and the development of integrated watershed models, simulating the watershed hydrological cycle at high temporal (hourly to daily) and spatial resolution (tens of meters) has become efficient and computationally affordable. These hyperresolution watershed models require high resolution of meteorological forcing as model input to ensure the fidelity and accuracy of simulated responses. In this study, we utilized the Advanced Terrestrial Simulator (ATS), an integrated watershed model, to simulate surface and subsurface flow and land surface processes using unstructured meshes at the Coal Creek Watershed near Crested Butte (Colorado). We compared simulated watershed hydrologic responses including streamflow and distributed variables such as evapotranspiration, snow water equivalent (SWE), and groundwater table driven by three publicly available, gridded meteorological forcings (GMFs) – Daily Surface Weather and Climatological Summaries (Daymet), the Parameter-elevation Regressions on Independent Slopes Model (PRISM), and the North American Land Data Assimilation System (NLDAS). By comparing various spatial resolutions (ranging from 400 m to 4 km) of PRISM, the simulated streamflow only becomes marginally worse when spatial resolution of meteorological forcing is coarsened to 4 km (or 30 % of the watershed area). However, the 4 km-resolution has much worse performance than finer resolution in spatially distributed variables such as SWE. Using the temporally disaggregated PRISM, we compared models forced by different temporal resolutions (hourly to daily), and sub-daily resolution preserves the dynamic watershed responses (e.g., diurnal fluctuation of streamflow) that are absent in results forced by daily resolution. Conversely, the simulated streamflow shows better performance using daily resolution compared to that using sub-daily resolution. Our findings suggest that the choice of GMF and its spatiotemporal resolution depends on the quantity of interest and its spatial and temporal scale, which may have important implications for model calibration and watershed management decisions.

54 ENVIRONMENTAL SCIENCES↗

Analysis of PICA-NuSil at the Hypersonic Materials Environmental Test System (HyMETS)

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the stardust mission and is the baseline material for missions to Mars. The phenolic phase of PICA is friable, and therefore the surface is treated with a siloxane resin (NuSil) to mitigate the shedding of particulate matter. Historically, material response models have not accounted for the presence of this coating. Therefore, a test campaign was executed at the Hypersonic Materials Environmental Test System (HyMETS) to investigate the performance of NuSil under multiple heating rates and gas composition. Sphere-cone models were instrumented with thermocouples to measure the in-depth temperature response, and surface temperatures were monitored with a pyrometer. Furthermore, emission spectrometers were used to collect spectra from species emitting at the surface of the test article and the post-shock stagnation region. Post-test analysis suggests that the coating decomposes according to multiple stages. Inspection of silicon emission data reveals a rapid rise and decline in intensity within the first couple of seconds of testing, which is attributed to pyrolysis of NuSil. Pyrolysis is followed by a stagnation in silicon emission and surface temperature which are attributed to the formation of a thin layer of oxidation-resistant silicon oxycarbide. Moreover, the presence of the silicon oxycarbide layer significantly reduces the in-depth thermal response when compared to the baseline material. Eventually, the silicon oxycarbide layer decomposes via high-temperature reactions, which expose the underlying char layer to a highly reactive flow.

PICA-NuSil↗

Carbon-phosphorus cycle models overestimate CO 2 enrichment response in a mature Eucalyptus forest

The importance of phosphorus (P) in regulating ecosystem responses to climate change has fostered P-cycle implementation in land surface models, but their CO 2 effects predictions have not been evaluated against measurements. Here, we perform a data-driven model evaluation where simulations of eight widely used P-enabled models were confronted with observations from a long-term free-air CO 2 enrichment experiment in a mature, P-limited Eucalyptus forest. We show that most models predicted the correct sign and magnitude of the CO 2 effect on ecosystem carbon (C) sequestration, but they generally overestimated the effects on plant C uptake and growth. We identify leaf-to-canopy scaling of photosynthesis, plant tissue stoichiometry, plant belowground C allocation, and the subsequent consequences for plant-microbial interaction as key areas in which models of ecosystem C-P interaction can be improved. Together, this data-model intercomparison reveals data-driven insights into the performance and functionality of P-enabled models and adds to the existing evidence that the global CO 2 -driven carbon sink is overestimated by models.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Changes in the Arctic Shortwave Cloud Radiative Effects

The shortwave cloud radiative effect (SWCRE) is important on the Arctic surface radiation budget and the major source of inter-model spread in predictions of Arctic climate. To better understand the individual contributions of various radiative processes to changes in SWCRE, the paper presents the use of the extended APRP (Atmospheric Radiative Perturbation Potential) method. This involves adding the absorptivity for the upward beam and considering differences in reflectivity between upward and downward beams, as well as analyzing the cloud masking effect resulting from changes in surface albedo in more detail. Using data from the CMIP5 and CMIP6 climate models, the study decomposes the SWCRE over the Arctic surface and analyzes inter-model differences in quadrupled CO2 simulations. The study takes into account the fact that the response of SWCRE to Arctic warming is influenced by changes in surface albedo, cloud amount, and cloud microphysics. Results show that in the sunlight season, the reduction in surface albedo associated with sea ice loss is directly linked to strong negative SWCRE, which explains the considerable model discrepancy. Arctic clouds can hinder the positive surface albedo feedback by changing the albedo in two ways: (1) decreasing incoming shortwave radiation due to cloud reflection and (2) by decreasing the shortwave reaching the surface after being reflected by clouds. In addition, increased (decreased) cloud amount and cloud liquid water are shown to be less (more) incoming shortwave fluxes at the surface, but not dominating factors to the Arctic surface radiation budget and its inter-model variation. Overall, the extended APRP method offers a useful tool for analyzing the complex interactions between clouds and radiative process, reasonably decomposes the individual SWCRE responses at the Arctic surface, and emphasizes that considering not only the cloud amount or its properties, but also surface albedo change is critical for the prediction of SWCRE on the Arctic surface.

Arctic cloud↗

Rocket Engine Turbine Blade Surface Pressure Distributions Experiment and Computations

Understanding the unsteady aspects of turbine rotor flow fields is critical to successful future turbine designs. A technology program was conducted at NASA's Marshall Space Flight Center to increase the understanding of unsteady environments for rocket engine turbines. The experimental program involved instrumenting turbine rotor blades with miniature surface mounted high frequency response pressure transducers. The turbine model was then tested to measure the unsteady pressures on the rotor blades. The data obtained from the experimental program is unique in two respects. First, much more unsteady data was obtained (several minutes per set point) than has been possible in the past. Also, an extensive steady performance database existed for the turbine model. This allowed an evaluation of the effect of the on-blade instrumentation on the turbine's performance. A three-dimensional unsteady Navier-Stokes analysis was also used to blindly predict the unsteady flow field in the turbine at the design operating conditions and at +15 degrees relative incidence to the first-stage rotor. The predicted time-averaged and unsteady pressure distributions show good agreement with the experimental data. This unique data set, the lessons learned for acquiring this type of data, and the improvements made to the data analysis and prediction tools are contributing significantly to current Space Launch Initiative turbine airflow test and blade surface pressure prediction efforts.

Hudson, Susan T.↗

Effects of Longwave Cloud Radiative Forcing Anomalies on the Atmospheric Response to Equatorial Pacific Sea Surface Temperature Anomalies

The latest version of the National Center for Atmospheric Research community climate model (CCM2) has been used to investigate cloud radiative forcing (CRF) anomalies associated with equatorial Pacific sea surface temperature (SST) anomalies and the effects of the longwave CRF (LWCRF) anomalies on the atmospheric response to the SST anomalies. The SST anomalies cause large CRF anomalies, both longwave and shortwave, as well as latent heat anomalies at low latitudes on a global scale. The relative magnitude of the simulated longwave and shortwave CRF anomalies is consistent with the result of the Earth Radiation Budget Experiment (ERBE), implying that cloud height and cloud radiative properties such as emissivity and reflectivity are well simulated by the model. The LWCRF anomaly strongly enhances the precipitation anomaly in the whole tropical belt. The positive (negative) LWCRF anomaly warms (cools) the troposphere and destabilizes (stabilizes) the upper troposphere. The LWCRF anomaly enhances the Southern Oscillation and the related Walker circulation anomaly. The effects of the LWCRF anomaly are essential to the northern hemispheric extratropical circulation anomaly, the Pacific/North American pattern.

Chen, M.↗

Mathematical Framework Underlying the In Situ Electrochemical Diagnosis of Adsorbed Intermediates Formed during Redox Reactions at Electrode Surfaces

Previously, we have presented an electrochemical technique wherein an electroactive tracer species is employed to probe the rate-limiting factors governing redox reactions at an electrode surface. In this technique, the electrode is first held potentiostatically to facilitate a redox process (step 1), and then the potential is released to open circuit conditions (step 2) so as to monitor the time-dependent re-equilibration of the electrode potential in the presence of the tracer. The time-dependent potential response in step 2 has been shown to contain information about diffusion—limited or desorption—limited steps, enabling in situ probing of the electrochemistry at the electrode surface during step 1. In the present contribution, a theoretical model governing the transient response in step 2 is developed for two limiting cases: diffusion—limited and desorption—limited recovery of the electrode potential. Mathematical modeling shows that, during re-equilibration, the step 2 potential transient corresponding to a case where step 1 involves surface adsorbed species which undergo desorption in step 2 exhibits a much longer time constant than that when re-equilibration occurs under diffusion limitations. The mathematical framework presented herein provides a sound fundamental basis for applying the aforementioned technique to studying adsorption-desorption processes during electrochemistry. Also, technique limitations are presented in light of the modeling findings.

25 ENERGY STORAGE↗

Responses of the Tropical Pacific to Wind Forcing as Observed by Spaceborne Sensors and Simulated by an Ocean General Circulation Model

In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific. The surface wind fields derived from observations by the spaceborne special sensor microwave imager (SSM/I) and from the operational products of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared. When SSM/I winds are used to force a primitive-equation ocean general circulation model (OGCM), they produce 3 C more surface cooling than ECMWF winds for the eastern equatorial Pacific during the cool phase of an El Nino-Southern Oscillation event. The stronger cooling by SSM/I winds is in good agreement with measurements at the moored buoys and observations by the advanced very high resolution radiometer, indicating that SSM/I winds are superior to ECMWF winds in forcing the tropical ocean. In comparison with measurements from buoys, tide gauges, and the Geosat altimeter, the OGCM simulates the temporal variations of temperature, steric, and sea level changes with reasonable realism when forced with the satellite winds. There are discrepancies between model simulations and observations that are common to both wind forcing fields, one of which is the simulation of zonal currents; they could be attributed to model deficiencies. By examining model simulations under two winds, vertical heat advection and uplifting of the thermocline are found to be the dominant factors in the anomalous cooling of the ocean mixed layer.

Liu, W. Timothy↗

Risks of Synchronized Low Yields Are Underestimated in Climate and Crop Model Projections

Simultaneous harvest failures across major crop-producing regions are a threat to global food security. Concurrent weather extremes driven by a strongly meandering jet stream could trigger such events, but so far this has not been quantified. Specifically, the ability of state-of-the art crop and climate models to adequately reproduce such high impact events is a crucial component for estimating risks to global food security. Here we find an increased likelihood of concurrent low yields during summers featuring meandering jets in observations and models. While climate models accurately simulate atmospheric patterns, associated surface weather anomalies and negative effects on crop responses are mostly underestimated in bias-adjusted simulations. Given the identified model biases, future assessments of regional and concurrent crop losses from meandering jet states remain highly uncertain. Our results suggest that model-blind spots for such high-impact but deeply-uncertain hazards have to be anticipated and accounted for in meaningful climate risk assessments.

Atmospheric dynamics↗

A Combined CFD/Material Response Analysis of 3MDCP Arcjet Experiments

The Mars Sample Return Earth Entry System (MSR-EES) project has selected 3-D Woven Mid-Density Carbon Phenolic (3MDCP) as the baseline thermal protection system (TPS) material for the capsule that will return Martian soil samples to Earth sometime in the 2030’s. A series of experiments in the NASA Ames AHF and IHF arcjets will be used to characterize the performance of 3MDCP and to develop and refine material response models for it. The test objectives for these experiments include obtaining in-depth and surface temperature data, mass loss and surface recession measurements, and char depth measurements. To provided pre- and post-test support for the arcjet experiments, a combined CFD/material response was performed using the DPLR and Icarus code. Pre-and post-test analysis on the 3MDCP arcjet experiments was performed using the DPLR CFD and Icarus material response codes. The codes were not tightly coupled, but boundary condition data required by Icarus (e.g., heat transfer coefficients and surface pressures), were extracted from the CFD solutions. This poster shows comparisons between DPLR/Icarus and the experimental data taken during arcjet experiment AHF-348 where 4-inch diameter iso-q models made from 3MDCP were tested in the 12-inch nozzle of the AHF arcjet. This poster also includes a top-level overview of the CFD/material response arcjet simulation process. Best-practices for running CFD solutions of arcjet experiments including how to determine the inflow conditions for the CFD solution from the arc heater settings are discussed. The required inputs to the material response solver are presented, as well as a discussion of how they can be extracted from the CFD solution.

thermal protection systems↗

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Green's Function Model Intercomparison Project (GFMIP) Protocol

The atmospheric Green's function method is a technique for modeling the response of the atmosphere to changes in the spatial field of surface temperature. While early studies applied this method to changes in atmospheric circulation, it has also become an important tool to understand changes in radiative feedbacks due to evolving patterns of warming, a phenomenon called the “pattern effect.” To better study this method, this paper presents a protocol for creating atmospheric Green's functions to serve as the basis for a model intercomparison project, GFMIP. The protocol has been developed using a series of sensitivity tests performed with the HadAM3 atmosphere-only general circulation model, along with existing and new simulations from other models. Our preliminary results have uncovered nonlinearities in the response of the atmosphere to surface temperature changes, including an asymmetrical response to warming versus cooling patch perturbations, and a change in the dependence of the response on the magnitude and size of the patches. These nonlinearities suggest that the pattern effect may depend on the heterogeneity of warming as well as its location. These experiments have also revealed tradeoffs in experimental design between patch size, perturbation strength, and the length of control and patch simulations. The protocol chosen on the basis of these experiments balances scientific utility with the simulation time and setup required by the Green's function approach. Running these simulations will further our understanding of many aspects of atmospheric response, from the pattern effect and radiative feedbacks to changes in circulation, cloudiness, and precipitation.

54 ENVIRONMENTAL SCIENCES↗