Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 451 records · Page 25

Do We Need to Parameterize Mesoscale Convective Organization to Mitigate the MJO-Mean State Trade-Off?

Modifications in cumulus parameterizations that improve Madden-Julian oscillation (MJO) simulation tend to degrade the mean state—known as the MJO-mean state trade-off. The impacts of parameterizing mesoscale convective organization on the relationship between the MJO and the mean state simulation fidelity are examined. A series of experiments are made with a general circulation model that parameterizes the mesoscale convective organization and simulates well both the MJO and the mean state. In the control simulation, a prognostic nondimensional variable (Ω) represents the degree of convective organization. In order to examine the effect of the parameterized convective organization, fixed Ω values are imposed in a series of constrained experiments. The fixed-Ω simulations show a negative relationship between the MJO and the mean state simulation fidelity. The control simulation is deviated from the negative relationship, suggesting the parameterized mesoscale convective organization helps general circulation models to mitigate the MJO-mean state trade-off.

58 GEOSCIENCES↗

Sensitivity of Carbonaceous Aerosol Properties to the Implementation of a Dynamic Aging Parameterization in the Regional Climate Model RegCM

Freshly emitted soot is hydrophobic, but condensation of secondary aerosols and coagulation with other particles modify its hygroscopic optical properties. This conversion is referred to as “aerosol aging.” Many climate models represent this aging process with a fixed aging time scale, whereas in reality, it is a dynamic process that depends on environmental conditions. Here, we implement a dynamic aging parameterization scheme in the regional climate model RegCM4 in place of the fixed aging timescale of 1.15 days (~27.6 h) and examine its impact on the aerosol life cycle over the Indian subcontinent. The conversion from hydrophobic to hydrophilic aerosol is usually lower than 27.6 h over the entire landmass and lower than 10 h over the polluted Indo-Gangetic Basin (IGB), with seasonal variability. Due to the implementation of the dynamic aging scheme, the column burden and surface mass concentration of carbonaceous aerosols increase during the drier season (December–February) when washout is negligible. The burden is reduced during the wet season (June–September) due to a more efficient washout except over the IGB, where a reduction in precipitation as a result of radiative feedbacks increases the aerosol concentrations. Over the polluted IGB, surface dimming increases due to the dynamic aging scheme, with the top of the atmosphere forcing remaining mostly unchanged. As a result, atmospheric heating increases by at least 1.2 W/m 2 . Our results suggest that climate models should incorporate dynamic aging for a more realistic representation of aerosol simulations, especially in highly polluted regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of the Fitch Wind-Farm Wake Parameterization with Large-Eddy Simulations of Wakes Using the Weather Research and Forecasting Model

Abstract Wind-farm parameterizations in weather models can be used to predict both the power output and farm effects on the flow; however, their correctness has not been thoroughly assessed. We evaluate the wind-farm parameterization of the Weather Research and Forecasting Model with large-eddy simulations (LES) of the wake performed with the same model. We study the impact on the velocity and turbulence kinetic energy (TKE) of inflow velocity, roughness, resolution, number of turbines (one or two), and inversion height and strength. We compare the mesoscale with the LES by spatially averaging the LES within areas correspondent to the mesoscale horizontal spacing: one covering the turbine area and two downwind. We find an excellent agreement of the velocity within the turbine area between the two types of simulations. However, within the same area, we find the largest TKE discrepancies because in mesoscale simulations, the turbine-added TKE has to be highest at the turbine position to be advected downwind. Within the downwind areas, differences between velocities increase as the wake recovers faster in the LES, whereas for the TKE both types of simulations show similar levels. From the various configurations, the impact of inversion height and strength is small for these heights and inversion levels. The highest impact for the one-turbine simulations appears under the low-speed case due to the higher thrust, whereas the impact of resolution is low for the large-eddy simulations but high for the mesoscale simulations. Our findings demonstrate that higher-fidelity simulations are needed to validate wind-farm parameterizations.

17 WIND ENERGY↗

Optimization-based control using input convex neural networks

Input convex neural networks (ICNNs) are a family of deep learning models where the outputs are constructed to be convex functions of the inputs. By parameterizing system models using ICNNs, optimization-based control problems can be solved as convex optimization problems, leading to improved performance and robustness. This work proposes a novel framework where the control objective function and constraints are modelled using ICNNs. A case study of optimization-based control with output constraints is conducted on a process with input multiplicity and nonminimum phase behavior. The simulation results demonstrate improved economic yield compared with normal neural networks. Additionally, the input convexity formulation is compared with simple regularization techniques, and unique benefits such as improved data efficiency and robustness of the proposed formulation are shown. Finally, by explicitly incorporating prior knowledge about convexity, this framework provides a good balance between the universal approximation power of deep learning and computational feasibility required by control.

42 ENGINEERING↗

A cash flow model of an integrated industrial CCS-EOR project in a petrochemical corridor: A case study in Louisiana

Petroleum refineries and petrochemical plants are major CO 2 sources, however, they are also significant capital and employment assets that are unlikely to be replaced in the near term. As a result, nations and states that are interested in reducing the carbon intensity of their economies will need to find ways to reduce the emissions of their existing industrial capacity. Industrial carbon capture provides one potential mechanism for reducing the carbon intensity of existing industrial facilities, however, an economically feasible capture system requires that the captured CO 2 be integrated into a system of transport and storage with income generated either through tax credits, enhanced oil recovery (EOR), or both. Here, we present a cash-flow model of an integrated system with industrial capture, pipeline transport, and EOR, and we parameterize the model with data from Louisiana. Given a 50/bbl oil price, an integrated capture, transport and EOR system that uses ethylene oxide production, ammonia production, or natural gas processing as sources is predicted to have a net present value of about $500 million; hydrogen-based capture has a cash flow of -214 given the same assumptions. Further, we find that the recent 45Q Tax Credit expansion has a positive impact on the cash flow of the system but does not change the overall profitability of the systems under the specified assumptions such that without the tax credits natural gas processing, ammonia production and ethylene oxide production-based capture systems remain cost-effective, while hydrogen-based capture remains unprofitable with or without the tax credit.

42 ENGINEERING↗

A Decision-Relevant Factor-Fixing Framework: Application to Uncertainty Analysis of a High-Dimensional Water Quality Model

Factor Fixing (FF) is a common method for reducing the number of model parameters to lower computational cost. FF typically starts with distinguishing the insensitive parameters from the sensitive and pursues uncertainty quantification (UQ) on the resulting reduced-order model, fixing each insensitive parameter at a fixed value. There is a need, however, to expand such a common approach to consider the effects of decision choices in the FF-UQ procedure on metrics of interest. Therefore, to guide the use of FF and increase confidence in the resulting dimension-reduced model, we propose a new adaptive framework consisting of four principles: (a) re-parameterize the model first to reduce obvious non-identifiable parameter combinations, (b) focus on decision relevance especially with respect to errors in quantities of interest (QoI), (c) conduct adaptive evaluation and robustness assessment of errors in the QoI across FF choices as sample size increases, and (d) reconsider whether fixing is warranted. The framework is demonstrated on a spatially-distributed water quality model. The error in estimates of QoI caused by FF can be estimated using a Polynomial Chaos Expansion (PCE) surrogate model. Built with 70 model runs, the surrogate is computationally inexpensive to evaluate and can provide global sensitivity indices for free. For the selected catchment, just two factors may provide an acceptably accurate estimate of model uncertainty in the average annual load of Total Suspended Solids (TSS), suggesting that reducing the uncertainty in these two parameters is a priority for future work before undertaking further formal uncertainty quantification.

54 ENVIRONMENTAL SCIENCES↗

Analysis and Modeling of the End-to-End I/O Performance on OLCF's Titan Supercomputer

With the increase of scale and complexity seen in a variety of leadership-class scientific computation and simulation applications, it has become more important to understand their I/O performance characteristics. The user-observed performance is a combination of properties of how the application is using the HPC facility, as well as how others' use of the facility causes variability in the static machine capabilities. Our work leverages statistical analysis of I/O performance data gathered with fine time resolution over a full week from Titan supercomputer. Based on observed properties of the distribution of I/O latencies, we build a three-state hidden Markov model (HMM) to characterize the end-to-end I/O performance on Titan. We parameterize our model using part of the field-gathered I/O performance data and validate it against the rest. The validation results demonstrate that our model can capture the dynamics of end-to-end I/O performance on Titan accurately.

Wan, Lipeng↗

Review of Wind–Wave Coupling Models for Large-Eddy Simulation of the Marine Atmospheric Boundary Layer

Abstract We present a review of existing wind–wave coupling models and parameterizations used for large-eddy simulation of the marine atmospheric boundary layer. The models are classified into two main categories: (i) the wave-phase-averaged, sea surface–roughness models and (ii) the wave-phase-resolved models. Both categories are discussed from their implementation, validity, and computational efficiency viewpoints, with emphasis given on their applicability in offshore wind energy problems. In addition to the various models discussed, a review of laboratory-scale and field-measurement databases is presented thereafter. The majority of the presented data have been gathered over many decades of studying air–sea interaction phenomena, with the most recent ones compiled to reflect an offshore wind energy perspective. Both provide valuable data for model validation. We also discuss the modeling knowledge gaps and computational challenges ahead.

54 ENVIRONMENTAL SCIENCES↗

Reduced-Order Parameterized Short-Circuit Current Model of Inverter-Interfaced Distributed Generators

The short-circuit response of inverter-interfaced distributed generators (IIDGs) is not adequately represented in many conventional protection studies. This paper presents an in-depth analysis of IIDG behavior during grid faults and proposes a more accurate reduced-order parameterized short-circuit current (RPSC) model of inverters. Typical inverter components are thoroughly investigated to identify those that play dominant roles during faults. The paper shows that the current limiter is the dominant factor for the steady-state fault-current of inverter; while the inverter filter along with the severity of the voltage disturbance largely determine the initial transient spike of inverter fault current. The proposed model is low-order and can be used in large scale simulations. The parameters of the proposed RPSC model can be extracted from laboratory experiments without requiring proprietary manufacturer information. The proposed fault-current model is analogous to the well-known synchronous machine model that segregates the inverter fault current into subtransient, transient, and steady-state fault currents. Finally, experimental and simulation tests are presented to validate the model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗

Colloid-Facilitated Transport of 137 Cs in Fracture-Fill Material. Experiments and Modeling

In this study, we demonstrate how a combination of batch sorption/desorption experiments and column transport experiments were used to effectively parameterize a model describing the colloid-facilitated transport of Cs in the Grimsel granodiorite/FFM system. Cs partition coefficient estimates onto both the colloids and the stationary media obtained from the batch experiments were used as initial estimates of partition coefficients in the column experiments, and then the column experiment results were used to obtain refined estimates of the number of different sorption sites and the adsorption and desorption rate constants of the sites. The desorption portion of the column breakthrough curves highlighted the importance of accounting for adsorption-desorption hysteresis (or a very nonlinear adsorption isotherm) of the Cs on the FFM in the model, and this portion of the breakthrough curves also dictated that there be at least two different types of sorption sites on the FFM. In the end, the two-site model parameters estimated from the column experiments provided excellent matches to the batch adsorption/desorption data, which provided a measure of assurance in the validity of the model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Diagnosing Near-Surface Model Errors with Candidate Physics Parameterization Schemes for the Multiphysics Rapid Refresh Forecast System (RRFS) Ensemble during Winter over the Northeastern United States and Southern Great Plains

Abstract During the winter of 2020/21 an ensemble of FV3-LAM forecasts was produced over the contiguous United States for the Winter Weather Experiment using five physics suites. These forecasts are evaluated with the goal of optimizing physics parameterizations within the future operational Rapid Refresh Forecast System (RRFS) in the Unified Forecast System (UFS) realm and for selecting suitable physics suites for a multiphysics RRFS ensemble. The five physics suites have different combinations of land surface models (LSMs), planetary boundary layer (PBL) parameterizations, and surface layer schemes, chosen from those used in current and possible future operational systems and likely to be supported in the operational UFS. Full-season evaluation reveals a persistent near-surface cold bias in the U.S. Northeast from one suite and a nighttime warm bias in the southern Great Plains in another suite, while other suites have smaller biases. A representative case is chosen to diagnose the cause for each of these biases using sensitivity simulations with different physics combinations or modified parameters and verified with additional mesonet observations. The cold bias in the Northeast is attributed to aspects of the Noah-MP LSM over snow cover, where Noah-MP simulates lower soil water content, and thus lower thermal conductivity than other LSMs, leading to less upward ground heat flux during nighttime and consequently lower surface temperature. The nighttime warm bias found in the southern Great Plains is attributed to overestimation of vertical mixing in the K -profile-based eddy-diffusivity mass-flux (K-EDMF) PBL scheme and insufficient land–atmospheric coupling from the GFS surface layer scheme over short vegetation. A few key parameters driving these systematic biases are identified.

Meteorology & Atmospheric Sciences↗

Composition of the core from gallium metal–silicate partitioning experiments

We present gallium concentration (normalized to CI chondrites) in the mantle is at the same level as that of lithophile elements with similar volatility, implying that there must be little to no gallium in Earth's core. Metal-silicate partitioning experiments, however, have shown that gallium is a moderately siderophile element and should be therefore depleted in the mantle by core formation. Moreover, gallium concentrations in the mantle (4 ppm) are too high to be only brought by the late veneer; and neither pressure, nor temperature, nor silicate composition has a large enough effect on gallium partitioning to make it lithophile. We therefore systematically investigated the effect of core composition (light element content) on the partitioning of gallium by carrying out metal–silicate partitioning experiments in a piston–cylinder press at 2 GPa between 1673 K and 2073 K. Four light elements (Si, O, S, C) were considered, and their effect was found to be sufficiently strong to make gallium lithophile. The partitioning of gallium was then modeled and parameterized as a function of pressure, temperature, redox and core composition. A continuous core formation model was used to track the evolution of gallium partitioning during core formation, for various magma ocean depths, geotherms, core light element contents, and magma ocean composition (redox) during accretion. The only model for which the final gallium concentration in the silicate Earth matched the observed value is the one involving a light-element rich core equilibrating in a FeO-rich deep magma ocean (>1300 km) with a final pressure of at least 50 GPa. More specifically, the incorporation of S and C in the core provided successful models only for concentrations that lie far beyond their allowable cosmochemical or geophysical limits, whereas realistic O and Si amounts (less than 5 wt.%) in the core provided successful models for magma oceans deeper that 1300 km. In conclusion, these results offer a strong argument for an O- and Si-rich core, formed in a deep terrestrial magma ocean, along with oxidizing conditions.

58 GEOSCIENCES↗

Tropical Cirrus in Global Storm-Resolving Models: 1. Role of Deep Convection

Pervasive cirrus clouds in the upper troposphere and tropical tropopause layer (TTL) influence the climate by altering the top-of-atmosphere radiation balance and stratospheric water vapor budget. These cirrus are often associated with deep convection, which global climate models must parameterize and struggle to accurately simulate. By comparing high-resolution global storm-resolving models from the Dynamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) intercomparison that explicitly simulate deep convection to satellite observations, we assess how well these models simulate deep convection, convectively generated cirrus, and deep convective injection of water into the TTL over representative tropical land and ocean regions. The DYAMOND models simulate deep convective precipitation, organization, and cloud structure fairly well over land and ocean regions, but with clear intermodel differences. All models produce frequent overshooting convection whose strongest updrafts humidify the TTL and are its main source of frozen water. Intermodel differences in cloud properties and convective injection exceed differences between land and ocean regions in each model. We argue that, with further improvements, global storm-resolving models can better represent tropical cirrus and deep convection in present and future climates than coarser-resolution climate models. To realize this potential, they must use available observations to perfect their ice microphysics and dynamical flow solvers.

54 ENVIRONMENTAL SCIENCES↗

Intensified Soil Moisture Extremes Decrease Soil Organic Carbon Decomposition: A Mechanistic Modeling Analysis

Earth system models have predicted that there will be more frequent and severe precipitation and drought events in terrestrial ecosystems. Microbially mediated decomposition of soil organic carbon (SOC) tends to increase as soils wet and decrease as soils dry. However, the long-term SOC change under intensified moisture extremes remains poorly known as it depends on the frequency and intensity of soil drying and wetting. In this study, we explored long-term SOC dynamics under scenarios of alternating drying-wetting cycles using the Microbial-ENzyme Decomposition model, a mechanistic microbial model. The model was parameterized with 11 years of observations from a temperate deciduous broadleaf forest site, showing satisfactory model performance in both model calibration (R 2 = 0.67) and validation (R 2 = 0.69) against heterotrophic respiration. We then used the model to simulate the long-term SOC dynamics under five scenarios of alternating drying-wetting cycles with different frequencies and severities over a period of 100 years. Results showed that the changes in active microbial biomass C and the corresponding turnover rates of SOC pools were more sensitive to soil drying than soil wetting. As a result, the cumulative soil carbon emission from microbial respiration decreased by 433.7 g C m -2 after the 100-year simulation in the highest frequency and intensity moisture scenario, but was not significantly affected by the lowest frequency and intensity scenario. This study emphasizes the nonlinear response of SOC decomposition to soil moisture changes, which causes decreased decomposition by microbes under drying that is, not compensated by increased decomposition under wetting conditions.

58 GEOSCIENCES↗

Real-Time Ice Nuclei Measurements at Oliktok Point Field Campaign Report

The ability of atmospheric particles (aerosols) to serve as heterogeneous ice nuclei (IN) for formation of ice and mixed-phase clouds is one of the least understood microphysical processes resulting in large uncertainties in climate modeling. As stated in the Intergovernmental Panel on Climate Change 2007 report, “The effects of aerosol particles on heterogeneous ice formation are currently insufficiently understood and present another level of challenge for both observations and modeling.” IN can strongly affect cloud radiative properties and precipitation formation. Only a small fraction (often less than 1%) of all aerosol particles act as atmospheric IN. Therefore, it becomes important to identify the active fraction of aerosol particles under a variety of environmental conditions. Proposed IN measurement techniques provide IN concentrations as a function of temperature and humidity conditions, and this unique data set can be analyzed to develop or constrain IN parameterizations that models rely on to calculate primary ice formation. Currently, the processes of ice formation within models are poorly represented, such that climate models often inaccurately partition cloud water between liquid and solid ice. Model studies have shown that uncertainties in cloud phase partitioning induce large uncertainties in the impact of mixedphase clouds on climate change. Arctic clouds composed of ice particles and supercooled liquid droplets at temperatures warmer than -37°C are important for the Arctic’s climate because of their extensive horizontal coverage over long periods, and their impact on the surface radiative balance and thus ice sheet coverage. The complexity of microphysical processes within these clouds and the sensitivity of precipitation to the details of microphysics are not fully understood. In addition, representation of these processes in a variety of cloud models continues to be challenging. Arctic clouds provide unique conditions to understand heterogeneous ice nucleation mechanisms, particularly immersion freezing, and will be the initial focus of the ice nucleation group. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM)-North Slope of Alaska (NSA) site is well equipped with other aerosol measurement (number concentration, size) instruments, and a relationship among long-term IN measurements and aerosol physical and chemical properties under different environmental conditions can be identified and established. The long-term IN measurements will be critical for understanding the ice nucleation process in response to the seasonal variations of aerosol and developing robust ice nucleation parameterizations for cloud and climate models. In the future, IN measurements could also be made at other ARM-NSA sites, specifically, Oliktok Point.

Arctic processes (meteorology), Oliktok Point, aer↗

An analysis of the operational GFS simplified Arakawa Schubert parameterization within a WRF framework: A Hurricane Sandy (2012) long‐term track forecast perspective

Abstract Hurricane Sandy (2012) is known as an incredibly destructive storm and one defined meteorologically by its large size, and its significant forecast track spreads among various operational models roughly 1 week before landfall. While the operational European Centre for Medium‐Range Weather Forecasts model accurately depicted a northeastern United States landfall, the Global Forecasting System (GFS) model consistently forecast a track toward the North Atlantic Ocean. Using a Weather Research and Forecasting (WRF) model framework, Bassill ( ) suggested that these differences were primarily a function of differences between the two models' cumulus parameterization (CP). This study also uses a WRF model framework to examine the simplified Arakawa Schubert CP used in the GFS model. It is found that increasing the deep convective entrainment coefficient produces more realistic forecast tracks for forecasts initialized roughly 1 week before landfall. This occurs through a reorientation of the precipitation (and associated latent heating) around Sandy during a critical time period in which it was interacting with a series of upper troughs to its west and northwest. Reorienting the latent heating reshapes the upper tropospheric steering pattern toward the one that is more negatively tilted and consistent with observations.

Bassill, Nick P.↗