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 361 records · Page 20

Scrutinizing entrainment and mass flux closures in shallow cumulus parameterizations using cloud-radar observations and large-eddy simulation (Final Report)

Representing warm, shallow clouds is a challenge for Global Climate Models, which leads to uncertainty in model projections of future climate. Our study was particularly concerned with how climate models represent the effects of shallow, non-precipitating cumulus clouds. Climate models must make certain assumptions when representing low clouds such as how much dry air from outside the cloud is mixed into the warm, moist cloud updraft. This process is called “entrainment” and tends to weaken the cloud updraft and dry out the cloud, compared to a theoretical cloud with no mixing. This project was a joint observational and modeling study to derive estimates of entrainment from field observations and use high-resolution modeling to evaluate the assumptions used in climate model representations of shallow cumulus. Our research took advantage of a wide variety of Department of Energy measurements from fixed climate-monitoring locations over the Southern Great Plains of the U.S. and the in the Eastern North Atlantic, along with short-term field projects (MAGIC, CAP–MBL, ACE–ENA). During the project, the PI maintained regular collaborations with Department of Energy Scientists at Brookhaven and Argonne National Laboratories. The project funded two students, who both obtained their M.S. in Atmospheric Science and are currently in the Atmospheric Science Ph.D. program. Both were lead authors on research papers funded by this grant.

54 ENVIRONMENTAL SCIENCES↗

Reduced basis approximations of parameterized dynamical partial differential equations via neural networks

Projection-based reduced order models are effective at approximating parameter-dependent differential equations that are parametrically separable. When parametric separability is not satisfied, which occurs in both linear and nonlinear problems, projection-based methods fail to adequately reduce the computational complexity. Devising alternative reduced order models is crucial for obtaining efficient and accurate approximations to expensive high-fidelity models. In this work, we develop a timestepping procedure for dynamical parameter-dependent problems, in which a neural-network is trained to propagate the coefficients of a reduced basis expansion. This results in an online stage with a computational cost independent of the size of the underlying problem. Here, we demonstrate our method on several parabolic partial differential equations, including a problem that is not parametrically separable.

97 MATHEMATICS AND COMPUTING↗

A North Sea in Situ Evaluation of the Fitch Wind Farm Parameterization Within the Mellor-Yamada-Nakanishi-Niino and 3D Planetary Boundary Layer Schemes

Wind resource assessments and wind power forecasts that account for wind farm wakes are sensitive to the choice of planetary boundary layer (PBL) scheme. This work compares the one-dimensional Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme with a three-dimensional PBL (3DPBL) scheme, evaluating predictions made with both schemes against two sets of North Sea in situ observations of wind farm wakes. The optimal PBL scheme varies based on the observations (FINO1 tower vs. aircraft), the quantity of interest (wind speed vs. turbulence kinetic energy [TKE]), and the error metric (bias, centered root mean square error [cRMSE], R2, and earth mover's distance [EMD]). Whereas 3DPBL wind speeds outperform MYNN wind speeds with respect to the cRMSE at the FINO1 site located at a single point within the turbine rotor layer, 3DPBL TKE bias is larger than MYNN TKE bias when compared to aircraft observations taken 100 m above a wind farm. Wind speeds in the aircraft region are ambiguous with regard to which PBL scheme is optimal. Aircraft MYNN wind speeds outperform 3DPBL wind speeds with respect to R2 and cRMSE but underperform with respect to bias and EMD. Future evaluations across broader temporal and spatial scales may offer further insight into model differences.

17 WIND ENERGY↗

Influence of explicit Phaeocystis parameterizations on the global distribution of marine dimethyl sulfide

Dimethyl sulfide (DMS) is a biogenic organosulfur compound which contributes strongly to marine aerosol mass and the determination of cloud condensation nuclei over the remote oceans. Since uncertainties in DMS flux to the atmosphere lead to large variations in climate forcing, the global DMS distribution has been the subject of increasingly complex dynamic simulations. DMS concentrations are directly controlled by marine ecosystems. Phaeocystis is a major DMS producer but is often omitted from global reduced sulfur mechanisms. Here we incorporate this phytoplankton group into the marine ecosystem-biogeochemical module of the Community Earth System Model. To examine its role in the ocean sulfur cycle, an earlier DMS model has been enhanced to include new knowledge gained over the last few years. Results from the baseline run show that simulated Phaeocystis biomass generally agrees with observations, with high concentrations near the Antarctic continent and between 50° and 60° north. Given the new explicit Phaeocystis representation, the DMS distribution shows significant improvements, especially regarding the amplitude and location of high-latitude peaks. The simulated global mean surface DMS value is 2.26 nM, comparable to an estimate of 2.34 nM from the latest climatology extrapolated based on observations. The total oceanic DMS source to the atmosphere is 20.4 Tg S/yr, on the low side of previous estimates. Comparisons with and without Phaeocystis show that the group dominates DMS distributions in temperate and cold waters, contributing 13% of the global flux. The proportion may increase as sea ice declines and should be considered in climate projections.

54 ENVIRONMENTAL SCIENCES↗

A decision tree algorithm for investigation of model biases related to dynamical cores and physical parameterizations: CESM/CAM EVALUATION BY DECISION TREES

An object-based evaluation method using a pattern recognition algorithm (i.e., classification trees) is applied to the simulated orographic precipitation for idealized experimental setups using the National Center of Atmospheric Research (NCAR) Community Atmosphere Model (CAM) with the finite volume (FV) and the Eulerian spectral transform dynamical cores with varying resolutions. Daily simulations were analyzed and three different types of precipitation features were identified by the classification tree algorithm. The statistical characteristics of these features (i.e., maximum value, mean value, and variance) were calculated to quantify the difference between the dynamical cores and changing resolutions. Even with the simple and smooth topography in the idealized setups, complexity in the precipitation fields simulated by the models develops quickly. The classification tree algorithm using objective thresholding successfully detected different types of precipitation features even as the complexity of the precipitation field increased. The results show that the complexity and the bias introduced in small-scale phenomena due to the spectral transform method of CAM Eulerian spectral dynamical core is prominent, and is an important reason for its dissimilarity from the FV dynamical core. The resolvable scales, both in horizontal and vertical dimensions, have significant effect on the simulation of precipitation. The results of this study also suggest that an efficient and informative study about the biases produced by GCMs should involve daily (or even hourly) output (rather than monthly mean) analysis over local scales.

54 ENVIRONMENTAL SCIENCES↗

Using precipitation, vertical root distribution, and satellite-retrieved vegetation information to parameterize water stress in a Penman-Monteith approach to evapotranspiration modeling under Mediterranean climate

Here, recent studies have shown that global Penman-Monteith equation based (PM-based) models poorly simulate water stress when estimating evapotranspiration (ET) in areas having a Mediterranean climate (AMC). In this study, we propose a novel approach using precipitation, vertical root distribution (VRD), and satellite-retrieved vegetation information to simulate water stress in a PM-based model (RS-WBPM) to address this issue. A multilayer water balance module is employed to simulate the soil water stress factor (SWSF) of multiple soil layers at different depths. The water stress factor (WSF) for surface evapotranspiration is determined by VRD information and SWSF in each layer. Additionally, four older PM-based models (PMOV) are evaluated at 27 flux sites in AMC. Results show that PMOV fails to estimate the magnitude or capture the variation of ET in summer at most sites, whereas RS-WBPM is successful. The daily ET resulting from RS-WBPM incorporating recommended VI (NDVI for shrub and EVI for other biomes) agrees well with observations, with R 2 = 0.60 ( RMSE = 18.72 W m -2 ) for all 27 sites and R 2 =0.62 ( RMSE 5 18.21 W m -2 ) for 25 nonagricultural sites. However, combined results from the optimum older PM-based models at specific sites show R 2 values of only 0.50 ( RMSE 5 20.74 W m -2 ) for all 27 sites. RS-WBPM is also found to outperform other ET models that also incorporate a soil water balance module. As all inputs of RS-WBPM are globally available, the results from RS-WBPM are encouraging and imply the potential of its implementation on a regional and global scale.

58 GEOSCIENCES↗

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗

Parameterized complex dielectric functions of CuIn 1−x Ga x Se 2 : applications in optical characterization of compositional non‐uniformities and depth profiles in materials and solar cells

Abstract In‐situ spectroscopic ellipsometry (SE) was employed to extract the complex dielectric functions ε = ε 1 + iε 2 over the spectral range of 0.75–6.5 eV for a set of polycrystalline CuIn 1− x Ga x Se 2 (CIGS) thin films with different alloy compositions x = [Ga]/{[In] + [Ga]}. For highest possible accuracy in ε for each CIGS thin film, specialized SE procedures were adopted including (i) deposition to a thickness of ~600 Å on smooth native oxide covered crystal silicon wafers, which minimizes the surface roughness on the film and thus the required corrections in data analysis, and (ii) measurement in‐situ , which minimizes ambient contamination and oxidation of the film surface. Assuming an analytical form for each of the ε spectra for these CIGS films, oscillator parameters were obtained in best fits, and these parameters were fit in turn to polynomials in x . With the resulting database of polynomial coefficients, the ε spectra for any composition of CIGS can be generated from the single parameter, x . In addition to enabling accurate contactless determination of bulk and surface roughness layer thicknesses of CIGS films by high speed multichannel SE, the database enables characterization of the composition and its profile with depth into these films, and even how the depth profile varies spatially within the plane of the films. In this study, depth profile parameters were found to correlate spatially with solar cell performance parameters. As a result, SE provides the capability of contactless compositional analysis of production‐scale CIGS photovoltaic modules at high speed. Copyright © 2016 John Wiley & Sons, Ltd.

Aryal, Puruswottam↗

Estimating and parameterizing mismatch power loss in bifacial photovoltaic systems

Nonuniform irradiance on the rear side of bifacial photovoltaic (PV) systems can cause additional mismatch loss, which may not be appropriately captured in PV energy production estimates and software. We evaluated several scenarios including horizontal single-axis tracking (HSAT) over natural ground-cover and rooftop-mounted systems over high albedo reflective roofs. We found mismatch losses of up to 1.5% annual loss for very close-mounted (0.15 m) rooftop systems, but losses for HSAT systems and high-ground-clearance rooftop systems were lower (<0.5%). A simplified empirical relationship was found that links the spatial variation of irradiance (specifically, the mean absolute difference of irradiance) to the resulting mismatch loss factor, with an R 2 better than .99. Furthermore, this approximate relationship was experimentally validated on mismatched PV modules, and it provides a basis for rapidly estimating bifacial mismatch loss factors for use in hourly PV performance simulations such as PVSyst or SAM. Additional parameters investigated include (a) climate dependence and module orientation, which were not found to have a strong impact on bifacial mismatch losses and (b) PV module fill factor and bifaciality ratio, which did have a strong linear impact on mismatch losses.

14 SOLAR ENERGY↗

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗

Role of Microphysical Parameterizations with Droplet Relative Dispersion in IAP AGCM 4.1

In previous studies we see that accurate descriptions of the cloud droplet effective radius (Re) and the autoconversion process of cloud droplets to raindrops (Au) can effectively improve simulated clouds and surface precipitation, and reduce the uncertainty of aerosol indirect effects in global climate models (GCMs). In this paper, we implement cloud microphysical schemes including two-moment Au and R e considering relative dispersion of the cloud droplet size distribution into version 4.1 of the Institute of Atmospheric Physics atmospheric GCM (IAP AGCM 4.1), which is the atmospheric component of the Chinese Academy of Sciences-Earth System model (CAS-ESM 1.0). An analysis of the effects of different schemes shows that the newly implemented schemes can improve both the simulated shortwave (SWCF) and longwave cloud radiative forcings (LWCF) as compared to the standard scheme in IAP AGCM 4.1. The new schemes also effectively enhance the large-scale precipitation, especially over low latitudes, although the influences of total precipitation are insignificant for different schemes. Further studies show that similar results can be found with the Community Atmosphere Model 5.1 (CAM5.1).

54 ENVIRONMENTAL SCIENCES↗

Parameterization of temperature and spectral distortions in future CMB experiments

CMB spectral distortions are induced by Compton collisions with electrons. We review the various schemes to characterize the anisotropic CMB with a non-Planckian spectrum. We advocate using logarithmically averaged temperature moments as the preferred language to describe these spectral distortions, both for theoretical modeling and observations. Numerical modeling is simpler, the moments are frame-independent, and in terms of scattering the mode truncation is exact.

79 ASTRONOMY AND ASTROPHYSICS↗

An Improved Discrete Least-Squares/Reduced-Basis Method for Parameterized Elliptic PDEs

It is shown that the computational efficiency of the discrete least-squares (DLS) approximation of solutions of stochastic elliptic PDEs is improved by incorporating a reduced-basis method into the DLS framework. In particular, we consider stochastic elliptic PDEs with an affine dependence on the random variable. The goal is to recover the entire solution map from the parameter space to the finite element space. To this end, first, a reduced-basis solution using a weak greedy algorithm is constructed, then a DLS approximation is determined by evaluating the reduced-basis approximation instead of the full finite element approximation. The main advantage of the new approach is that one only need apply the DLS operator to the coefficients of the reduced-basis expansion, resulting in huge savings in both the storage of the DLS coefficients and the online cost of evaluating the DLS approximation. In addition, the recently developed quasi-optimal polynomial space is also adopted in the new approach, resulting in superior convergence rates for a wider class of problems than previous analyzed. Numerical experiments are provided that illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗

Numerical Errors in Ice Microphysics Parameterizations and their Effects on Simulated Regional Climate

The major characteristics of ice microphysics in Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) bulk-type cloud microphysics originate from the diagnosed ice number concentration, which is a function of the cloud-ice mixing ratio. In this study, we correct numerical errors in ice microphysics processes of the WDM6, in which the cloud-ice shape is assumed as single bullets and examine the impact on regional climate simulations. Here, by rederiving the relationships between cloud microphysics characteristics, including the one linking the cloud-ice mixing ratio and number concentration, we remove numerical errors intrinsic to the description of cloud-ice characteristics in the original WDM6 microphysics scheme. The revised WDM6 is tested using a WRF framework for regional climate simulations over the East Asian region. We find that our correction to the WDM6 improves the model’s performance in capturing the observed distribution of the monsoon rain band. A reduction in cloud ice is significant in the revised WDM6, which strengthens the Western North Pacific High. By conducting the additional sensitivity experiment in which the characteristics of cloud-ice shape are revised as the one for the column type, our study also finds out that the impacts of the existing numerical errors on the simulated monsoon is as large as the ones of the changes in cloud-ice shape.

54 ENVIRONMENTAL SCIENCES↗