Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 793 records · Page 44

The zonally averaged transport characteristics of the atmosphere as determined by a general circulation model

Two dimensional modeling has become an established technique for the simulation of the global structure of trace constituents. Such models are simpler to formulate and cheaper to operate than three dimensional general circulation models, while avoiding some of the gross simplifications of one dimensional models. Nevertheless, the parameterization of eddy fluxes required in a 2-D model is not a trivial problem. This fact has apparently led some to interpret the shortcomings of existing 2-D models as indicating that the parameterization procedure is wrong in principle. There are grounds to believe that these shortcomings result primarily from incorrect implementations of the predictions of eddy transport theory and that a properly based parameterization may provide a good basis for atmospheric modeling. The existence of these GCM-derived coefficients affords an unprecedented opportunity to test the validity of the flux-gradient parameterization. To this end, a zonally averaged (2-D) model was developed, using these coefficients in the transport parameterization. Results from this model for a number of contrived tracer experiments were compared with the parent GCM. The generally good agreement substantially validates the flus-gradient parameterization, and thus the basic principle of 2-D modeling.

Plumb, R. A.↗

Multiple Quasi-Equilibria of the ITCZ and the Origin of Monsoon Onset: Rotational ITCZ Attractors - Part 2

Chao's numerical and theoretical work on multiple quasi-equilibria of the intertropical convergence zone (ITCZ) and the origin of monsoon onset is extended to solve two additional puzzles. One is the highly nonlinear dependence on latitude of the "force" acting on the ITCZ due to earth's rotation, which makes the multiple quasi-equilibria of the ITCZ and monsoon onset possible. The other is the dramatic difference in such dependence when different cumulus parameterization schemes are used in a model. Such a difference can lead to a switch between a single ITCZ at the equator and a double ITCZ, when a different cumulus parameterization scheme is used. Sometimes one of the double ITCZ can diminish and only the other remain, but still this can mean different latitudinal locations for the single ITCZ. A single idea based on two off-equator attractors for the ITCZ, due to earth's rotation and symmetric with respect to the equator, and the dependence of the strength and size of these attractors on the cumulus parameterization scheme solves both puzzles. The origin of these rotational attractors, explained in Part I, is further discussed. The "force" acting on the ITCZ due to earth's rotation is the sum of the "forces" of the two attractors. Each attractor exerts on the ITCZ a "force" of simple shape in latitude; but the sum gives a shape highly varying in latitude. Also the strength and the domain of influence of each attractor vary, when change is made in the cumulus parameterization. This gives rise to the high sensitivity of the "force" shape to cumulus parameterization. Numerical results, of experiments using Goddard's GEOS general circulation model, supporting this idea are presented. It is also found that the model results are sensitive to changes outside of the cumulus parameterization. The significance of this study to El Nino forecast and to tropical forecast in general is discussed.

Chao, Winston C.↗

Numerical Study of the Role of Shallow Convection in Moisture Transport and Climate

The objective of this investigation was to study the role of shallow convection on the regional water cycle of the Mississippi and Little Washita Basins of the Southern Great Plains (SGP) using a 3-D mesoscale model, the PSU/NCAR MM5. The underlying premise of the project was that current modeling of regional-scale climate and moisture cycles over the continents is deficient without adequate treatment of shallow convection. At the beginning of the study, it was hypothesized that an improved treatment of the regional water cycle can be achieved by using a 3-D mesoscale numerical model having high-quality parameterizations for the key physical processes controlling the water cycle. These included a detailed land-surface parameterization (the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) sub-model of Wetzel and Boone), an advanced boundary-layer parameterization (the 1.5-order turbulent kinetic energy (TKE) predictive scheme of Shafran et al.), and a more complete shallow convection parameterization (the hybrid-closure scheme of Deng et al.) than are available in most current models. PLACE is a product of researchers working at NASA's Goddard Space Flight Center in Greenbelt, MD. The TKE and shallow-convection schemes are the result of model development at Penn State. The long-range goal is to develop an integrated suite of physical sub-models that can be used for regional and perhaps global climate studies of the water budget. Therefore, the work plan focused on integrating, improving, and testing these parameterizations in the MM5 and applying them to study water-cycle processes over the SGP. These schemes have been tested extensively through the course of this study and the latter two have been improved significantly as a consequence.

Seaman, Nelson L.↗

Pion, Kaon, Proton and Antiproton Production in Proton-Proton Collisions

Inclusive pion, kaon, proton, and antiproton production from proton-proton collisions is studied at a variety of proton energies. Various available parameterizations of Lorentz-invariant differential cross sections as a function of transverse momentum and rapidity are compared with experimental data. The Badhwar and Alper parameterizations are moderately satisfactory for charged pion production. The Badhwar parameterization provides the best fit for charged kaon production. For proton production, the Alper parameterization is best, and for antiproton production the Carey parameterization works best. However, no parameterization is able to fully account for all the data.

Norbury, John W.↗

Relativistic Three-Dimensional Lippmann-Schwinger Cross Sections for Space Radiation Applications

Radiation transport codes require accurate nuclear cross sections to compute particle fluences inside shielding materials. The Tripathi semi-empirical reaction cross section, which includes over 60 parameters tuned to nucleon-nucleus (NA) and nucleus-nucleus (AA) data, has been used in many of the world’s best-known transport codes. Although this parameterization fits well to reaction cross section data, the predictive capability of any parameterization is questionable when it is used beyond the range of the data to which it was tuned. Using uncertainty analysis, it is shown that a relativistic Three- Dimensional Lippmann-Schwinger (LS3D) equation model based on Multiple Scatter- ing Theory (MST) that uses 5 parameterizations—3 fundamental parameterizations to nucleon-nucleon (NN) data and 2 nuclear charge density parameterizations—predicts NA and AA reaction cross sections as well as the Tripathi cross section parameterization for reactions in which the kinetic energy of the projectile in the laboratory frame (TLab) is greater than 220 MeV/n. The relativistic LS3D model has the additional advantage of being able to predict highly accurate total and elastic cross sections. Consequently, it is recommended that the relativistic LS3D model be used for space radiation applications in which TLab > 220 MeV/n.

C M Werneth↗

Observational Constraints on Warm Cloud Microphysical Processes Using Machine Learning and Optimization Techniques

We introduce new parameterizations for autoconversion and accretion rates that greatly improve representation of the growth processes of warm rain. The new parameterizations capitalize on machine‐learning and optimization techniques and are constrained by in situ cloud probe measurements from the recent Atmospheric Radiation Measurement Program field campaign at Azores. The uncertainty in the new estimates of autoconversion and accretion rates is about 15% and 5%, respectively, outperforming existing parameterizations. Our results confirm that cloud and drizzle water content are the most important factors for determining accretion rates. However, for autoconversion, in addition to cloud water content and droplet number concentration, we discovered a key role of drizzle number concentration that is missing in current parameterizations. The robust relation between autoconversion rate and drizzle number concentration is surprising but real, and furthermore supported by theory. Thus, drizzle number concentration should be considered in parameterizations for improved representation of the autoconversion process.

54 ENVIRONMENTAL SCIENCES↗

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗

Local structure graph models with higher-order dependence

Local structure graph models (LSGMs) describe random graphs and networks as a Markov random field (MRF)—each graph edge has a specified conditional distribution dependent on explicit neighbourhoods of other graph edges. Centered parameterizations of LSGMs allow for direct control and interpretation of parameters for large- and small-scale structures (e.g., marginal means vs. dependence). Here, we extend this parameterization to account for triples of dependent edges and illustrate the importance of centered parameterizations for incorporating covariates and interpreting parameters. Using a MRF framework, common exponential random graph models are also shown to induce conditional distributions without centered parameterizations and thereby have undesirable features. This work attempts to advance graph models through conditional model specifications with modern parameterizations, covariates and higher-order dependencies.

97 MATHEMATICS AND COMPUTING↗

Optimization and Evaluation of Stochastic Unified Convection Using Single-Column Model Simulations at Multiple Observation Sites

We extend the previously developed stochastic unified convection scheme (UNICON) for shallow convection to deep convection by parameterizing the impact of mesoscale organized flow on updraft properties. The extended stochastic UNICON parameterizes thermodynamic properties of updrafts at the near-surface as a multivariate Gaussian distribution, where the variances of the distribution are the summation of variances from non-organized turbulence and mesoscale organized flow. The distribution of updraft radius is parameterized as a power-law distribution with a scale break which is parameterized as a linear function of the strength of mesoscale organized flow. The proposed parameterization is validated using a series of large-eddy simulations of deep convection. The free parameters introduced in the formulation of stochastic UNICON are optimized using 10 cases of single-column model simulations over the ocean. Stochastic UNICON with the optimized parameters significantly reduces the biases of thermodynamic profiles and surface precipitation rates simulated in the original UNICON for tropical convection cases. The simulation of the variation in anomalies of temperature and moisture associated with the Madden-Julian oscillation is also improved. The overall improvements in simulated thermodynamic profiles are found to be due to the increased heating and drying tendencies by convective processes in stochastic UNICON. An additional simulation of an idealized deep convection case shows that stochastic UNICON produces enhanced cloud variabilities with dependency on updraft radius, indicating its ability to represent the coexistence of shallow and deep convection.

54 ENVIRONMENTAL SCIENCES↗

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗

Formulation of autoconversion and drop spectra shape in shallow cumulus clouds

Two-moment autoconversion parameterizations as compared to accretion parameterizations exhibit significant errors suggesting that additional moments are needed to increase their accuracy. We develop a three-moment autoconversion parameterization using output from an LES model with size-resolved microphysics. Adding the third moment decreases the errors of parameterization and improves precipitation prediction. However, the errors are still significantly larger than errors of accretion rate. An analysis of the cloud drop size distributions (DSDs) in the simulated tropical convective cloud system reveals that most of DSDs have a significant fraction of cloud liquid water content (qc) in the mid-size droplet range (radii from 20 to 40 microns). Our data indicates that more than 30% of DSDs have over half of qc contained in the mid-size range and about 60% of spectra have, at least, one third of qc in this range. Even when the rain/drizzle mode is small (radar reflectivity Z < -10 dBZ), there is a significant number of spectra in which fraction of qc in the mid-size range is as large as 60%. These DSDs are more complex than the frequently used single-mode Gamma or Log-normal distributions, which usually do not have a considerable tail extending to the mid-size range and can be defined by three microphysical moments. The need to define DSDs by more than three moments explains the large errors in the three-moment autoconversion parameterization. The limitation of three-parameter Gamma or Log-normal distributions should be kept in mind when applying them in precipitating shallow Cu clouds.

Kogan, Yefim↗

Resolved Gravity Waves in the Extratropical Stratosphere: Effect of Horizontal Resolution Increase from O (10) to O (1) km

Global ECMWF IFS simulations with horizontal grid spacings of 1, 4, and 9 km are used to assess gravity wave forcing (GWF) in the extratropical stratosphere. Results with important implications for GWF parameterizations at high and intermediate resolutions are presented. A doubling in the zonal-mean resolved GWF is observed when the horizontal resolution is increased from 9 to 1 km. Small-scale gravity waves with horizontal wavelengths < 100 km dominate this increase. Over most regions, excluding the polar night jet in the Antarctic spring, the total (resolved + parameterized) GWF at 9 km (4 km) is underestimated by up to 30% (15%). This implies that the parameterization of GWF is still required at 9 and 4 km horizontal resolutions. Despite the small land area in the Southern Hemisphere (SH), the resolved orographic and nonorographic GWF contribute equally to the total GWF in the SH at 1 km resolution. This is not reflected in the partitioning of the parameterized GWF, which has a significantly larger nonorographic contribution at 9 km. As a result, a zonal-mean momentum budget analysis reveals that the total GWF contributes one-third of SH springtime polar vortex deceleration at 1 km, whereas the contribution is as much as 50% at 9 km. In conclusion, this suggests that a rebalancing of the parameterized nonorographic and orographic GWF is required.

54 ENVIRONMENTAL SCIENCES↗

Improving GCM Representation of Convective Cloud Microphysics by Using ARM Raman Lidar and Cloud Radar Observations

The overall objective of this Department of Energy (DOE) Atmospheric System Research (ASR) funded project is to improve the representation of convective cloud microphysics in global climate models (GCMs) and check it by comparing the model simulations with observations. We have derived the cloud ice water content by synthesizing ARM Raman lidar (RL) and cloud radar observations at the Atmospheric Radiation Measurement (ARM) sites. Noting that the simulated anvil clouds in terms of their macro- and micro-physical properties are sensitive to the parameterization of convective microphysical processes, observed ice water content in anvil clouds provides a useful constraint on these parameterizations. We have improved the convective microphysics parameterization scheme by (1) considering sedimentation for cloud ice crystals that do not fall in the original scheme, (2) applying a new terminal velocity parameterization that depends on the environmental conditions for convective snow, (3) adding a new hydrometeor category, “rimed ice,” to the original four-class (cloud liquid, cloud ice, rain, and snow) scheme, and (4) allowing convective clouds to detrain snow particles into stratiform clouds. We have examined the impact of improved convective cloud microphysics parameterization on the simulated global climate from GCM simulations.

54 ENVIRONMENTAL SCIENCES↗

Long-term variability in immersion-mode marine ice-nucleating particles from climate model simulations and observations

Abstract. Ice-nucleating particles (INPs) in the Southern Ocean (SO) atmosphere have significant impacts on cloud radiative and microphysical properties. Yet, INP prediction skill in climate models remains poorly understood, in part because of the lack of long-term measurements. Here we show, for the first time, how model-simulated INP concentrations compare with year-round INP measurements during the Macquarie Island Cloud Radiation Experiment (MICRE) campaign from 2017–2018. We simulate immersion-mode INP concentrations using the Energy Exascale Earth System Model version 1 (E3SMv1) by combining simulated aerosols with recently developed deterministic INP parameterizations and the native classical nucleation theory (CNT) for mineral dust in E3SMv1. Because MICRE did not collect aerosol measurements of super-micron particles, which are more effective ice nucleators, we evaluate the model's aerosol fields at other high-latitude sites using long-term in situ observations of dust and sea spray aerosol. We find that the model underestimates dust and overestimates sea spray aerosol concentrations by 1 to 2 orders of magnitude for most of the high-latitude sites in the Southern Hemisphere. We next compare predicted INP concentrations with concentrations of INPs collected on filter samples (typically for 2 or 3 d) and processed offline using the Colorado State University ice spectrometer (IS) in immersion freezing mode. We find that when deterministic parameterizations for both dust and sea spray INPs are used, simulated INPs are within a factor of 10 of observed INPs more than 60 % of the time during summer. Our results also indicate that the E3SM's current treatment of mineral dust immersion freezing in the SO is impacted by compensating biases – an underprediction of dust amount was compensated by an overprediction of its effectiveness as INPs. We also perform idealized droplet freezing experiments to quantify the implications of the time-dependent behavior assumed by the E3SM's CNT-parameterization and compare with the ice spectrometer observations. We find that the E3SM CNT 10 s diagnostic used in this study is a reasonable approximation of the exact formulation of CNT, when applied to ice spectrometer measurements in low-INP conditions similar to Macquarie Island. However, the linearized 10 s diagnostic underestimates the exact formula by an order of magnitude or more in places with high-INP conditions like the Sahara. Overall, our findings suggest that it is important to correct the biases in E3SM's simulated dust life cycle and update E3SM's INP parameterizations. INP prediction errors of 2 to 3 orders of magnitude can have considerable impacts on the simulated cloud and radiative properties in global climate models. On comparing INP concentrations during MICRE against ship-based campaigns, Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) and Antarctic Circumnavigation Expedition (ACE), we find that INPs from the latter are significantly higher only in regions closer to Macquarie Island. This alludes to the fact that physical, chemical and biological processes affecting INP concentrations as stimulated by the island could be partly responsible for the high INP concentrations observed at Macquarie Island during the MICRE campaign. Therefore, improvements to both aerosol simulation and INP parameterizations are required to adequately simulate INPs and their cloud impacts in E3SM. It will be helpful to include a parallel measurement of the size-resolved aerosol composition and explore opportunities for long-term measurement platforms in future field campaigns studying INP sources in remote marine regions.

54 ENVIRONMENTAL SCIENCES↗

Emulating aerosol optics with randomly generated neural networks

Abstract. Atmospheric aerosols have a substantial impact on climate and remain one of the largest sources of uncertainty in climate prediction. Accurate representation of their direct radiative effects is a crucial component of modern climate models. However, direct computation of the radiative properties of aerosol populations is far too computationally expensive to perform in a climate model, so optical properties are typically approximated using a parameterization. This work develops artificial neural networks (ANNs) capable of replacing the current aerosol optics parameterization used in the Energy Exascale Earth System Model (E3SM). A large training dataset is generated by using Mie code to directly compute the optical properties of a range of atmospheric aerosol populations given a large variety of particle sizes, wavelengths, and refractive indices. Optimal neural architectures for shortwave and longwave bands are identified by evaluating ANNs with randomly generated wirings. Randomly generated deep ANNs are able to outperform conventional multilayer-perceptron-style architectures with comparable parameter counts. Finally, the ANN-based parameterization produces significantly more accurate bulk aerosol optical properties than the current parameterization when compared with direct Mie calculations using mean absolute error. The success of this approach makes possible the future inclusion of much more sophisticated representations of aerosol optics in climate models that cannot be captured by extension of the existing parameterization scheme and also demonstrates the potential of random-wiring-based neural architecture search in future applications in the Earth sciences.

54 ENVIRONMENTAL SCIENCES↗

A simplified, efficient approach to hybrid wind and solar plant site optimization

Abstract. Wind plant layout optimization is a difficult, complex problem with a large number of variables and many local minima. Layout optimization only becomes more difficult with the addition of solar generation. In this paper, we propose a parameterized approach to wind and solar hybrid power plant layout optimization that greatly reduces problem dimensionality while guaranteeing that the generated layouts have a desirable regular structure. Thus far, hybrid power plant optimization research has focused on system sizing. We go beyond sizing and present a practical approach to optimizing the physical layout of a wind–solar hybrid power plant. We argue that the evolution strategy class of derivative-free optimization methods is well-suited to the parameterized hybrid layout problem, and we demonstrate how hard layout constraints (e.g., placement restrictions) can be transformed into soft constraints that are amenable to optimization using evolution strategies. Next, we present experimental results on four test sites, demonstrating the viability, reliability, and effectiveness of the parameterized evolution strategy approach for generating optimized hybrid plant layouts. Completing the tool kit for parameterized layout generation, we include a brief tutorial describing how the parameterized evolutionary approach can be inspected, understood, and debugged when applied to hybrid plant layouts.

14 SOLAR ENERGY↗

The role of convective model choice in calculating the climate impact of doubling CO2

The role of the parameterization of vertical convection in calculating the climate impact of doubling CO2 is assessed using both one-dimensional radiative-convective vertical models and in the latitude-dependent Hadley-baroclinic model of Lindzen and Farrell (1980). Both the conventional 6.5 K/km and the moist-adiabat adjustments are compared with a physically-based, cumulus-type parameterization. The model with parameterized cumulus convection has much less sensitivity than the 6.5 K/km adjustment model at low latitudes, a result that can be to some extent imitiated by the moist-adiabat adjustment model. However, when averaged over the globe, the use of the cumulus-type parameterization in a climate model reduces sensitivity only approximately 34% relative to models using 6.5 K/km convective adjustment. Interestingly, the use of the cumulus-type parameterization appears to eliminate the possibility of a runaway greenhouse.

Lindzen, R. S.↗

Predictive Compensator Optimization for Head Tracking Lag in Virtual Environments

We examined the perceptual impact of plant noise parameterization for Kalman Filter predictive compensation of time delays intrinsic to head tracked virtual environments (VEs). Subjects were tested in their ability to discriminate between the VE system's minimum latency and conditions in which artificially added latency was then predictively compensated back to the system minimum. Two head tracking predictors were parameterized off-line according to cost functions that minimized prediction errors in (1) rotation, and (2) rotation projected into translational displacement with emphasis on higher frequency human operator noise. These predictors were compared with a parameterization obtained from the VE literature for cost function (1). Results from 12 subjects showed that both parameterization type and amount of compensated latency affected discrimination. Analysis of the head motion used in the parameterizations and the subsequent discriminability results suggest that higher frequency predictor artifacts are contributory cues for discriminating the presence of predictive compensation.

Adelstein, Barnard D.↗