A Fast Quasi-Static Time Series Simulation Method for PV Smart Inverters with Var Control using Linear Sensitivity Model.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
Understanding the impact of distributed photovoltaic (PV) resources on various elements of the distribution feeder is imperative for their cost effective integration. A year-long quasi-static time series (QSTS) simulation at 1-second granularity is often necessary to fully study these impacts. However, the significant computational burden associated with running QSTS simulations is a major challenge to their adoption. In this paper, we propose a fast scalable QSTS simulation algorithm that is based on a linear sensitivity model for estimating voltage-related PV impact metrics of a three-phase unbalanced, nonradial distribution system with various discrete step control elements including tap changing transformers and capacitor banks. Here, the algorithm relies on computing voltage sensitivities while taking into account all the effects of discrete controllable elements in the circuit. Consequently, the proposed sensitivity model can accurately estimate the state of controllers at each time step and the number of control actions throughout the year. For the test case of a real distribution feeder with 2969 buses (5469 nodes), 6 load/PV time series power profiles, and 9 voltage regulating elements including controller delays, the proposed algorithm demonstrates a dramatic time reduction, more than 180 times faster than traditional QSTS techniques.
Abstract. Multi-meter sea level rise (SLR) is thought to be possible within the next few centuries, with most of the uncertainty originating from the Antarctic land ice contribution. One source of uncertainty relates to the ice sheet model initialization. Since ice sheets have a long response time (compared to other Earth system components such as the atmosphere), ice sheet model initialization methods can have significant impacts on how the ice sheet responds to future forcings. To assess this, we generated 25 different ice sheet spin-ups, using the Community Ice Sheet Model (CISM) at a 4 km resolution. During each spin-up, we varied two key parameters known to impact the sensitivity of the ice sheet to future forcing: one related to the sensitivity of the ice shelf melt rate to ocean thermal forcing (TF) and the other related to the basal friction. The spin-ups all nudge toward observed thickness and enforce a no-advance calving criterion, such that all final spin-up states resemble observations but differ in their melt and friction parameter settings. Each spin-up was then forced with future ocean thermal forcings from 13 different CMIP6 models under the Shared Socioeconomic Pathway (SSP)5-8.5 emissions scenario and modern climatological surface mass balance data. Our results show that the effects of the ice sheet and ocean parameter settings used during the spin-up are capable of impacting multi-century future SLR predictions by as much as 2 m. By the end of this century, the effects of these choices are more modest, but still significant, with differences of up to 0.2 m of SLR. We have identified a combined ocean and ice parameter space that leads to widespread mass loss within 500 years (low friction and high melt rate sensitivity). To explore temperature thresholds, we also ran a synthetically forced CISM ensemble that is focused on the Amundsen region only. Given certain ocean and ice parameter choices, Amundsen mass loss can be triggered with thermal forcing anomalies between 1.5 and 2 ∘C relative to the spin-up. Our results emphasize the critical importance of considering ice sheet and ocean parameter choices during spin-up for SLR predictions and suggest the importance of including glacial isostatic adjustment in ice sheet simulations.
We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.
Abstract. Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO2 injection, seasonality of SO2 and black carbon (BC) particulate emissions, and the assumed fraction of SO2 emissions that is injected into the atmosphere as particulate phase sulfate (SO4) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO2, SO4, and BC, with a particularly large relative variation for SO2, which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of −0.35 W m−2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8 d), surface SO2 concentration (up to 59 % decrease), and surface sulfate concentration (up to 23 % increase). Emitting SO2 at height consistently increased SO2 and SO4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO4 emission fraction are accurately and consistently represented in global models.
Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO 2 injection, seasonality of SO 2 and black carbon (BC) particulate emissions, and the assumed fraction of SO 2 emissions that is injected into the atmosphere as particulate phase sulfate (SO 4 ) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO 2 , SO 4 , and BC, with a particularly large relative variation for SO 2 , which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO 2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of -0.35?W?m-2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8?d), surface SO 2 concentration (up to 59?% decrease), and surface sulfate concentration (up to 23?% increase). Emitting SO 2 at height consistently increased SO 2 and SO 4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO 4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO 4 emission fraction are accurately and consistently represented in global models.
Fatigue life knockdown due to surface roughness induces local stress concentrations. Physical fatigue experiments are crucial to qualify materials for use but are limited in their capability to examine the mechanisms underlying the fatigue crack formation driving force due to these stress concentrations. Microstructure-sensitive computational models can aid in this regard and are employed in this work to examine the effects of surface roughness in Al 7075-T6. Digital microstructure models are subjected to crystal plasticity finite element method (CPFEM) simulations to examine the combined effects of microstructure and grain scale asperities on the driving force for fatigue crack formation. Following elastic–plastic shakedown, mesoscale volume-averaged fatigue indicator parameters (FIPs) are computed within fatigue damage process zones of grains. An increase in the intensity of realistic surface roughness profiles corresponds accordingly to larger FIPs but it is difficult to precisely interpret the mechanisms and effects of this intensification. We thus parametrically investigate surface asperities that couple with microstructure and lack of constraint on slip at the free surface to describe microstructure-sensitive surface roughness knockdown effects on fatigue resistance. The depth of a single valley (i.e., notch) due to surface roughness is found to be more detrimental to fatigue resistance than the radius (i.e., notch width, not to be confused with notch acuity). Finally, we also find that different crystallographic textures, grain morphologies, and number of sampled grains have a limited effect on the depth of influence of these notches.
In this work we studied the sensitivity of different computational fluid dynamics (CFD) numerical solutions on turbulence model, domain discretization mesh type, and inlet boundary conditions. The case under examination is a turbulent one-phase flow (water) through a scaled nuclear pressurized water reactor (PWR) fuel assembly. The main feature of this geometrical configuration is the presence of a spacer grid with mixing vanes. The design of a spacer grid with mixing vanes plays a key role for the core heat transfer efficiency, increasing the coolant heat transfer coefficient. The main reason for this is the introduction of the top mixing vanes, which consist of an ordered repetitive succession of wings that generate secondary swirling flows. The increase of heat transfer efficiency also increases the critical heat flux (CHF), making it possible to increase the design and operational limits, far from boiling crises (departure of nucleate boiling in case of a PWR). The main focus of this work was to analyze the simulation capabilities of Reynolds-averaged Navier-Stokes (RANS) models (a cubic non-linear k-epsilon model formulation) for this complex flow geometry. The CFD code used was CD-Adapco STAR-CCM+.
The eastern coastal area of the Korean Peninsula, Yeongdong region, often experiences severe snowfall during the winter season. We studied here the development mechanism of a snowfall case occurred in February 2014 using the Weather Research and Forecasting (WRF) model at the convection-permitting resolution of 1 km. The model reasonably captured the observed snowfall, and three factors were discovered as important contributors to this snowfall case: a temperature contrast between land and ocean, a high-pressure system over the northeastern part of the Korean Peninsula, and the mountainous terrain. Cold and dry air blown from the high-pressure system gains heat and moisture when it passes over the East Sea enhancing the low-level instability. When the air mass reaches the coastline, updraft is generated by the density difference with the continental atmosphere and by the mountain slope. Further sensitivity experiments were conducted by adjusting strategies of nesting down (offline and online). Sensitivity experiments consist of three offline nesting experiments with different boundary condition update frequencies in 6, 3 and 1 h and one online nesting experiment that updates the boundary condition every integration time step. The major snowfall along the mountain range located in the middle of domain does not change significantly regardless of the nesting frequency. The offline nesting causes unphysical discontinuities of precipitation along the north and east sides of lateral boundary where wind blows in, which are caused by the absence of condensed hydrometeor variables in the lateral boundary condition. The domain nesting strategy tested in this study has dramatically impacted on the upper troposphere structure in the model: artificial fluctuations at the upper troposphere in the offline nesting simulations were significantly alleviated in the online nesting simulation.
Offshore wind is an emerging industry in the with rapidly changing technology advancements and an accelerating global deployment. Policy demands and interest in the United States has grown, but the long-term future for offshore wind remains uncertain. This study is designed to evaluate offshore wind's potential role in the future U.S. energy system. The analysis uses long-term power system models to assess a wide-range of power system possibilities. Specific objectives of the analysis are to: (1) identify the determinants of offshore wind deployment under different scenarios and (2) understand the impacts of substantial levels of offshore wind deployment on the power system. We also highlight future research needs for more robust understanding the opportunities and challenges to U.S. offshore wind development.
Carbon tetrachloride (CCl 4 ) is an ozone-depleting substance, which is controlled by the Montreal Protocol and for which the atmospheric abundance is decreasing. But, the current observed rate of this decrease is known to be slower than expected based on reported CCl 4 emissions and its estimated overall atmospheric lifetime. Here we use a three-dimensional (3-D) chemical transport model to investigate the impact on its predicted decay of uncertainties in the rates at which CCl 4 is removed from the atmosphere by photolysis, by ocean uptake and by degradation in soils. The largest sink is atmospheric photolysis (74 % of total), but a reported 10 % uncertainty in its combined photolysis cross section and quantum yield has only a modest impact on the modelled rate of CCl 4 decay. This is partly due to the limiting effect of the rate of transport of CCl 4 from the main tropospheric reservoir to the stratosphere, where photolytic loss occurs. The model suggests large interannual variability in the magnitude of this stratospheric photolysis sink caused by variations in transport. The impact of uncertainty in the minor soil sink (9 % of total) is also relatively small. In contrast, the model shows that uncertainty in ocean loss (17 % of total) has the largest impact on modelled CCl 4 decay due to its sizeable contribution to CCl 4 loss and large lifetime uncertainty range (147 to 241 years). Furthermore, with an assumed CCl 4 emission rate of 39 Gg year -1 , the reference simulation with the best estimate of loss processes still underestimates the observed CCl 4 (overestimates the decay) over the past 2 decades but to a smaller extent than previous studies. Changes to the rate of CCl 4 loss processes, in line with known uncertainties, could bring the model into agreement with in situ surface and remote-sensing measurements, as could an increase in emissions to around 47 Gg year -1 . Further progress in constraining the CCl 4 budget is partly limited by systematic biases between observational datasets. For example, surface observations from the National Oceanic and Atmospheric Administration (NOAA) network are larger than from the Advanced Global Atmospheric Gases Experiment (AGAGE) network but have shown a steeper decreasing trend over the past 2 decades. The observed differences imply a difference in emissions which is significant relative to uncertainties in the magnitudes of the CCl 4 sinks.
Biomass-burning smoke drives large uncertainty in climate projections of the Earth's radiative balance. This is due to the chemical and physical evolution of smoke and its impact on clouds and radiation. Here we focus on the southeastern Atlantic region and its inflow of African biomass-burning smoke during August 2017. We evaluate smoke properties and processes in two coupled earth-system models, the Energy Exascale Earth System Model (E3SM) and Community Earth System Model (CESM). These are compared against in situ aircraft observations from two field campaigns, ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) and CLoud–Aerosol–Radiation Interaction and Forcing: Year 2017 (CLARIFY-2017). Observations reveal an increase and subsequent decrease in smoke mean diameter, and a steady decrease in the mass ratio of organic aerosol (OA) to black carbon aerosol (BC) (OA:BC) over 4–12 days of aging, neither captured by the base models. Implementation of a photolytic loss scheme for secondary organic aerosol (SOA)—as a proxy for other heterogeneous volatilization chemistry—and a ∼1-day conversion for primary OA to SOA significantly improves the representation of this loss. In the boundary layer, both models show dimethyl sulfide driving a large increase in the sulfate aerosol mass fraction from the free troposphere, which is consistent with observations. Finally, models tend to underpredict cloud droplet number concentration partially due to weak modeled turbulent updraft strength, and model performance improves when the parameterized turbulent updraft strength is increased substantially. These results are expected to provide insights into future model development to reduce climate model uncertainties.
Abstract The subpolar North Atlantic (SPNA) shows contrasting responses in two sensitivity experiments with increased stratospheric aerosols, offering insight into the physical processes that may impact the Atlantic meridional overturning circulation (AMOC) in a warmer climate. In one, the upper ocean becomes warm and salty, but in the other it becomes cold and fresh. The changes are accompanied by diverging AMOC responses. The first experiment strengthens the AMOC, opposing the weakening trend in the reference simulation. The second experiment shows a much smaller impact. Both simulations use the Community Earth System Model with the Whole Atmosphere Community Climate Model component (CESM-WACCM) but differ in model versions and stratospheric aerosol specifications. Despite both experiments using similar approaches to increase stratospheric aerosols to counteract the rising global temperature, the contrasting SPNA and AMOC responses indicate a considerable dependency on model physics, climate states, and model responses to forcings. This study focuses on examining the physical processes involved with the impact of stratospheric aerosols on the SPNA salinity changes and their potential connections with the AMOC and the Arctic. We find that in both cases, increased stratospheric aerosols act to enhance the SPNA upper-ocean salinity by reducing freshwater export from the Arctic, which is closely tied to the Arctic sea ice changes. The impact on AMOC is primarily through the thermal component of the surface buoyancy fluxes, with negligible contributions from the freshwater component. These experiments shed light on the physical processes that dictate the important connections between the SPNA, the Arctic, the AMOC, and their subsequent feedbacks on the climate system.
Predicting the effects of microstructure on high-temperature creep responses of steel components is critical to minimize the risks of failure and maximize economic viability in the energy sector. In this research, a recently developed advanced mechanistic constitutive model is employed to study the effect of microstructure on the creep responses and to rationalize experimentally reported variability in the performances of grade-91 alloy. Hundreds of experimental creep tests from the literature are used to assess the predictability of the model. The model proposed is shown to accurately predict the steady-state creep rates and also the less commonly considered, yet important, primary-to-secondary transients for a wide range of creep conditions. Further, the model is capable of extrapolating the creep response of grade 91 alloy even outside of the calibration regime. Using this model, the roles of initial microstructure described in terms of grain size, dislocation density, and precipitate content, on the creep responses are investigated. The effect of initial microstructure on the steady-state creep rates and creep-strain is found to vary with stress and temperature. Grain size plays a significant role in the low-stress creep, where the diffusional creep is dominant. On the other hand, the effects of dislocation densities and precipitate content are significant in the dislocation-plasticity-dominated high-stress and temperature regimes. Furthermore, by comparing model predictions against a large experimental creep database, it is found that variability in the creep responses can be rationalized on the basis of differences in the initial microstructure. Overall, this work provides a robust pathway for microstructure-aware engineering design, verification, and validation of metallic components for energy applications.
Significant particle-in-cell code development has enabled simulations of power flow in multi-MA accelerators to include the desorption of surface contaminants, their ionization into surface plasmas, and the impact of these plasmas on efficiency. The simulations base desorption on an Arrhenius equation, whose most significant unknown is the surface contaminant inventory. The sensitivity of power-flow simulations to this inventory is studied here using Sandia National Laboratories' Z accelerator with a 7-nH MagLIF load [Phys. Plasmas 17, 056303 (2010)]. Simulations are conducted in 3D cylindrical coordinates for the current-adder, or “convolute,” region of Z and in 2D for the final feed only. Simulated contaminant inventories are varied from 1 to 32 monolayers (MLs) in 2D, and 2 to 4 ML in 3D. The results reveal sensitivities to the local ratio of E/B. The high B-field, low E-field region near the short-circuit load is insensitive to the contaminant inventory, where assumed values of 4–32 ML change the load current by ≤ 2%, and agree with experiment to within 2% at peak current. A 1-ML value is the outlier, increasing the load current by 5%, but still within measurement uncertainty. In contrast, the relatively higher E-field, lower B-field convolute region has slower contaminant desorption and higher-magnitude E-field penetration of the surface plasmas. The current loss in the convolute region does increase with contaminant inventory. The loss assuming 4 ML is 12% larger than for 2 ML, with 4 ML being the better match to experiment.
Abstract not provided.
Abstract not provided.