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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Integration of nitrogen dynamics into the Noah-MP land surface model v1.1 for climate and environmental predictions

Climate and terrestrial biosphere models consider nitrogen an important factor in limiting plant carbon uptake, while operational environmental models view nitrogen as the leading pollutant causing eutrophication in water bodies. The community Noah land surface model with multi-parameterization options (Noah-MP) is unique in that it is the next-generation land surface model for the Weather Research and Forecasting meteorological model and for the operational weather/climate models in the National Centers for Environmental Prediction. Here in this study, we add a capability to Noah-MP to simulate nitrogen dynamics by coupling the Fixation and Uptake of Nitrogen (FUN) plant model and the Soil and Water Assessment Tool (SWAT) soil nitrogen dynamics. This model development incorporates FUN's state-of-the-art concept of carbon cost theory and SWAT's strength in representing the impacts of agricultural management on the nitrogen cycle. Parameterizations for direct root and mycorrhizal-associated nitrogen uptake, leaf retranslocation, and symbiotic biological nitrogen fixation are employed from FUN, while parameterizations for nitrogen mineralization, nitrification, immobilization, volatilization, atmospheric deposition, and leaching are based on SWAT. The coupled model is then evaluated at the Kellogg Biological Station – a Long Term Ecological Research site within the US Corn Belt. Results show that the model performs well in capturing the major nitrogen state/flux variables (e.g., soil nitrate and nitrate leaching). Furthermore, the addition of nitrogen dynamics improves the modeling of net primary productivity and evapotranspiration. The model improvement is expected to advance the capability of Noah-MP to simultaneously predict weather and water quality in fully coupled Earth system models.

58 GEOSCIENCES↗
NASA NTRS2023

The Sensitivity of the Equatorial Pacific ODZ to Particulate Organic Matter Remineralization in A Climate Model Under Pre-Industrial Conditions

Marine oxygen plays a fundamental role in regulating the transfer of organic carbon and nutrients to their dissolved inorganic forms, serving as the terminal electron acceptor for heterotrophic respiration. Oxygen can become limiting to these processes in coastal and open-ocean oxygen deficient zones (ODZs). The maintenance of ODZs depends on the balance between physical processes such as ventilation and biogeochemical processes such as remineralization. These two processes act in opposing directions on ODZs, with ventilation responsible for transporting oxygen from the surface, where it is near saturation with the partial pressure of O2 in the atmosphere, to deeper oxygen-depleted layers, and remineralization acting to consume oxygen through biogeochemical processes such as microbial degradation throughout the water column. Remineralization is represented in all CMIP6 models, but the actual parameterizations, as well as their magnitude, are widely varying. In this study, we examine the sensitivity to remineralization of the equatorial Pacific ODZ O2 ol/kg) in a model; the NASA GISS Ocean Biogeochemical Model (NOBM-G) embedded in the NASA-GISS coupled atmosphere-ocean model. We find that increasing the remineralization rate shoals the ODZ onset depth (i.e. the regionally averaged upper bound of the ODZ), and decreases ODZ volume and regional-averaged thickness. Changes in biological consumption and vertical convergence of oxygen are identified as the primary contributors to changes in ODZ onset depth. ODZ thickness is mainly influenced by the shoaling of the bottom boundary of the ODZ, which decreases with maximum remineralization rate due to decreasing oxygen consumption and increasing horizontal oxygen convergence in deep waters. On the other hand, vertically-averaged ODZ area has a more complex, non-monotonic relationship with maximum remineralization rate. While these findings suggest an important role for remineralization in determining the shape of the ODZ in our model, the relative importance of remineralization vs. other physical parameterizations remains to be established. While our results reflect the structure of our ocean biogeochemical model, the relationship of remineralization to ODZ characteristics in other models should be examined to better inform model inter-comparisons.

Oxygen-deficient zone↗
DOE OSTI2025

Ship‐Based Lidar Evaluation of Southern Ocean Low Clouds in the Storm‐Resolving General Circulation Model ICON and the ERA5 and MERRA‐2 Reanalyses

Global storm resolving models (GSRMs) represent the next generation of global climate models. One of them is a 5-km Icosahedral Nonhydrostatic Weather and Climate Model (ICON). Its high resolution means that parameterizations of convection and clouds, including subgrid-scale clouds, are omitted, relying on explicit simulation but necessarily utilizing microphysics and turbulence parameterizations. Standard-resolution (10–100 km) models, which use convection and cloud parameterizations, have substantial cloud biases over the Southern Ocean (SO), adversely affecting radiation and sea surface temperature. The SO is dominated by low clouds, which cannot be observed accurately from space due to overlapping clouds, attenuation, and ground clutter. We evaluated SO clouds in ICON and the ERA5 and MERRA-2 reanalyzes using approximately 2400 days of lidar observations and 2300 radiosonde profiles from 31 voyages and a Macquarie Island station during 2010–2021, compared to the model and reanalyzes using a ground-based lidar simulator. We found that ICON and the reanalyzes underestimate the total cloud fraction by about 10% and 20%, respectively. ICON and ERA5 overestimate the cloud occurrence peak at about 500 m, associated with underestimated lower tropospheric stability and overestimated lifting condensation level. The reanalyzes strongly underestimate fog and very low-level clouds, and MERRA-2 underestimates cloud occurrence at almost all heights. Outgoing shortwave radiation is overestimated in MERRA-2, implying a “too few, too bright” cloud problem. SO cloud and fog biases are a substantial issue in the analyzed model and reanalyzes and result in shortwave and longwave radiation biases.

Kuma, Peter [Stockholm Univ. (Sweden); Univ. of Ca↗
DOE OSTI2021

Modelica-based system modeling for studying control-related faults in chiller plants and boiler plants serving large office buildings

System modeling is critical when studying operation faults in chiller plants and boiler plants. However, current fault models have difficulties faithfully representing the operation of chiller plants and boiler plants under the effects of those faults, especially the control-related ones. In this study, we present a systematic method to develop high-fidelity models for approximating the behaviors of chiller plants and boiler plants under faulty conditions. Compared to existing ones, the resulting fault models have two advantages: first, they better characterize the dynamic patterns in the system operation. In those models, control architecture and control logic are faithfully implemented. Thus, they can be used to study control-related faults, such as incorrect staging control due to sensor bias and mistuned feedback control. Second, they are readily extensible and can support large-scale investigations to explore different faulty conditions/scenarios. Those models are established in a hierarchical structure while modules in each layer can be redeclared and parameterized at upper layers. In such a case, modifications to the models can be realized through model modifiers and the process can be easily streamlined with scripts. We applied the proposed models in a comprehensive fault impact evaluation of the 13 control-related faults of chiller and boiler plants. In this evaluation, the proposed plant model is coupled with the EnergyPlus thermal load model to study the impact of various faulty scenarios. Based on the evaluation results, we identified the faults that have the most significant impacts on the operation of the chiller and boiler plants, respectively. We also found that the relationship between the impacts of the studied faults and the severity level of the faults are highly non-linear

Huang, Sen↗
DOE OSTI2022

Parameterization of cloud transmittance for expeditious assessment and forecasting of all-sky DNI

Radiative transfer models require vast computing resources to solve cloud transmittance and reflectance from the radiative transfer equation. As a result, models offering precise simulation in operations often acquire individual cloud transmittance or reflectance from a lookup table precomputed for practicable scenarios. To further expedite the computation of global horizontal irradiance and to reduce the storage requirements, the Fast All-sky Radiation Model for Solar applications (FARMS) parameterized the lookup table using elementary functions with specified coefficients. This study extends FARMS direct normal irradiance (DNI) computation by utilizing hyperbolic tangent functions and various polynomial functions to parameterize the cloud transmittance for scattered solar radiation in the circumsolar region. The parameterization is implemented in FARMS with DNI (FARMS-DNI) and accounts for the circumsolar radiation when assessing or forecasting DNI. The evaluation, with long-term observations at the National Renewable Energy Laboratory's, Solar Radiation Research Laboratory, and the Atmospheric Radiation Measurement, Southern Great Plains, Central Facility, shows that the parameterized DNIs are virtually identical with those computed by coupling FARMS-DNI to a lookup table of cloud transmittance. This parameterization has diverse applications in radiative transfer models and numerical weather prediction models used to assess or forecast direct solar radiation.

14 SOLAR ENERGY↗
DOE OSTI2014

Probing non-standard gravity with the growth index: a background independent analysis

Measurements of the growth index of linear matter density fluctuations γ(z) provide a clue as to whether Einstein's field equations encompass gravity also on large cosmic scales, those where the expansion of the universe accelerates. We show that the information encoded in this function can be satisfactorily parameterized using a small set of coefficients γ{sub i}, in such a way that the true scaling of the growth index is recovered to better than 1% in most dark energy and dark gravity models. We find that the likelihood of current data, given this formalism and the Λ Cold Dark Matter (ΛCDM) expansion model of Planck, is maximal for γ{sub 0} = 0.74{sup +0.44}{sub −0.41} and γ{sub 1} = 0.01{sup +0.46}{sub −0.46}, a measurement compatible with the ΛCDM predictions (γ{sub 0} = 0.545, γ{sub 1} = −0.007). In addition, data tend to favor models predicting slightly less growth of structures than the Planck ΛCDM scenario. The main aim of the paper is to provide a prescription for routinely calculating, in an analytic way, the amplitude of the growth indices γ{sub i} in relevant cosmological scenarios, and to show that these parameters naturally define a space where predictions of alternative theories of gravity can be compared against growth data in a manner which is independent from the expansion history of the cosmological background. As the standard Ω-plane provides a tool to identify different expansion histories H(t) and their relation to various cosmological models, the γ-plane can thus be used to locate different growth rate histories f(t) and their relation to alternatives model of gravity. As a result, we find that the Dvali-Gabadadze-Porrati gravity model is rejected with a 95% confidence level. By simulating future data sets, such as those that a Euclid-like mission will provide, we also show how to tell apart ΛCDM predictions from those of more extreme possibilities, such as smooth dark energy models, clustering quintessence or parameterized post-Friedmann cosmological models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Response of cumulus clouds to large-scale forcing and cumulus parameterization

Using the cumulus ensemble model (Soong and Ogura, 1980), statistical properties of cumulus clouds that occur in response to the imposed large scale forcing were investigated both in the tropical and midlatitude situtations. A strong drying process was found to occur in the boundary layer in association with deep convection in a midlatitude case. The Arakawa-Schubert (1974) cumulus parameterization scheme was tested semi-prognostically against both the model result and the observation in a tropical rain event that occurred in GATE. The Cloud heating and drying effects predicted by the Arakawa-Schubert scheme were found to agree well both with the observation and the model result. However, it was also found that the Arakawa-Schuber scheme underestimates both condensation and evaporation rates substantially. An inclusion of the downdraft effect, as formulated by Johnson (1976), appears to alleviate this deficiency. This downdraft effect may be important in predicting the behavior of the boundary accurately.

Ogura, Y.↗
DOE OSTI2024

Understanding Processes Controlling the Temporal and Spatial Variations of PBL Structures Over the ARM SGP Site

The surface heat, moisture, and momentum fluxes are transferred to the atmosphere above through the planetary boundary layer (PBL), where vertical mixing due to turbulent eddies of different sizes plays critical roles. Therefore, reliably representing PBL processes in numerical models is critical for weather, climate, and air quality prediction. Currently, there are over ten PBL schemes that are selectable within the advanced research version of the Weather Research and Forecasting (WRF) model, indicative of the challenges in capturing the impacts of turbulence within the PBL in models. Further improvements in PBL parameterizations are needed for both weather and climate models, as emphasized in many recent national reports, but require an advanced understanding of the underlying boundary layer processes from observations. This project takes advantage of Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) investments in the atmospheric boundary layer observations and Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) simulations to characterize PBL structures, understand key physical processes controlling the mixed layer development, and to evaluate PBL parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Surrogate construction via weight parameterization of residual neural networks

Surrogate model development is a critical step for uncertainty quantification or other sample-intensive tasks for complex computational models. Here, in this work, we develop a multi-output surrogate form using a class of neural networks (NNs) that employ shortcut connections, namely Residual NNs (ResNets). ResNets are known to regularize the surrogate learning problem and improve the efficiency and accuracy of the resulting surrogate. Inspired by the continuous, Neural ODE analogy, we augment ResNets with weight parameterization strategy with respect to ResNet depth. Weight-parameterized ResNets regularize the NN surrogate learning problem and allow better generalization with a drastically reduced number of learnable parameters. We demonstrate that weight-parameterized ResNets are more accurate and efficient than conventional feed-forward multi-layer perceptron networks. We also compare various options for parameterization of the weights as functions of ResNet depth. We demonstrate the results on both synthetic examples and a large scale earth system model of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Dependence of Convective Cloud Microphysical Properties on Environmental Conditions during the TRACER and ESCAPE Field Campaigns: A Synergistic Approach of Observations, Machine Learning and Parcel Models

The sensitivity of convective clouds to aerosols and their interactions with environment, combined with limited observational constraints in parameterizations, introduces significant uncertainties in atmospheric models. Here, this study investigates the dependence of convective cloud microphysical properties on environmental conditions using a synergistic approach that combines unique observations from the TRACER and ESCAPE field campaigns, machine learning techniques, and parcel model simulations with a super-droplet microphysics scheme. A random forest algorithm identifies in-situ vertical velocity (w), temperature (T), and surface fine-mode aerosol mass concentration as the three most important environmental conditions influencing cloud properties including liquid water content (LWC), number concentration for particles with D max < 50 μm (N c ,<50), 50 μm ≤ D max ≤ 3000 μm (N c,50–3000 ), and droplet effective diameter (D e ). Results show that LWC, N c,<50 , and N c,50–3000 significantly increase with w in updrafts. Across w bins, as T decreases, LWC, D e , and N c,50–3000 increase, while N c,<50 decreases, which are closely linked to the distance above cloud bases. Warmer cloud bases yield higher LWC, greater N c,50–3000 , and smaller N c,<50 , while polluted environments produce greater N c,<50 . Parcel model simulations successfully replicate these observed dependencies. The simulation results indicate that warmer cloud bases enhance condensation generating larger droplets, and differences in droplet sizes are then amplified through collision-coalescence, resulting in a greater N c,50–3000 . Polluted conditions result in a greater N c,<50 primarily due to enhanced cloud condensation nuclei activation despite increased collision-coalescence rates compared to pristine conditions. This study provides observed quantitative patterns characterizing cloud microphysical properties as a function of key environmental parameters, offering valuable constraints for improving physics parameterizations and numerical models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2007

Empirical Storm-Time Correction to the International Reference Ionosphere Model E-Region Electron and Ion Density Parameterizations Using Observations from TIMED/SABER

The response of the ionospheric E-region to solar-geomagnetic storms can be characterized using observations of infrared 4.3 micrometers emission. In particular, we utilize nighttime TIMED/SABER measurements of broadband 4.3 micrometers limb emission and derive a new data product, the NO+(v) volume emission rate, which is our primary observation-based quantity for developing an empirical storm-time correction the IRI E-region electron density. In this paper we describe our E-region proxy and outline our strategy for developing the empirical storm model. In our initial studies, we analyzed a six day storm period during the Halloween 2003 event. The results of this analysis are promising and suggest that the ap-index is a viable candidate to use as a magnetic driver for our model.

Mertens, Christoper J.↗
DOE OSTI2022

Inferring relative dose-dependent color center populations in proton irradiated thoria single crystals using optical spectroscopy

In this work, we have utilized photoluminescence spectroscopy and optical ellipsometry to characterize the dose-dependence of the photoluminescence emission intensity and changes in optical absorption of thoria single crystals subject to irradiation with energetic protons at room- and elevated-temperatures. The photoluminescence peaks and the optical absorption bands are attributed to creation of new electronic states emerging from defects resulting from displacement damage. These bands are most likely associated with electrons trapped at the oxygen vacancy sites similar to color centers formed in other irradiated oxides and halides. Our experimental observations are supported by a standard density functional theory calculation of the electronic structure in pristine and oxygen vacancy-bearing thoria crystals. The dose-dependence of the intensity of the photoluminescence peaks is used to parameterize a rate theory model that estimates the concentration of color centers in the irradiated crystals. This parameterization provides optimized migration barrier parameters for oxygen interstitials and vacancies that simultaneously capture the optical response of the crystals irradiated at room- and elevated-temperature. These optical spectroscopy techniques offer a promising pathway to characterize the population of color centers formed at the sites of oxygen anion vacancies, particularly in irradiated nuclear fuels, where atomic-level defects cannot be readily imaged using electron microscopy. When combined with other direct and indirect characterization tools, our approach can provide new insight into defect formation and accumulation in energy materials over single atomic to extended length scales.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Improvement in the Modeled Representation of North American Monsoon Precipitation Using a Modified Kain–Fritsch Convective Parameterization Scheme

A commonly noted problem in the simulation of warm season convection in the North American monsoon region has been the inability of atmospheric models at the meso-β scales (10 s to 100 s of kilometers) to simulate organized convection, principally mesoscale convective systems. With the use of convective parameterization, high precipitation biases in model simulations are typically observed over the peaks of mountain ranges. To address this issue, the Kain–Fritsch (KF) cumulus parameterization scheme has been modified with new diagnostic equations to compute the updraft velocity, the convective available potential energy closure assumption, and the convective trigger function. The scheme has been adapted for use in the Weather Research and Forecasting (WRF). A numerical weather prediction-type simulation is conducted for the North American Monsoon Experiment Intensive Observing Period 2 and a regional climate simulation is performed, by dynamically downscaling. In both of these applications, there are notable improvements in the WRF model-simulated precipitation due to the better representation of organized, propagating convection. The use of the modified KF scheme for atmospheric model simulations may provide a more computationally economical alternative to improve the representation of organized convection, as compared to convective-permitting simulations at the kilometer scale or a super-parameterization approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Prognostic residual mean flow in an ocean general circulation model and its relation to prognostic Eulerian mean flow

Recent work has shown that taking the thickness-weighted average (TWA) of the Boussinesq equations in buoyancy coordinates results in exact equations governing the prognostic residual mean flow where eddy–mean flow interactions appear in the horizontal momentum equations as the divergence of the Eliassen–Palm flux tensor (EPFT). It has been proposed that, given the mathematical tractability of the TWA equations, the physical interpretation of the EPFT, and its relation to potential vorticity fluxes, the TWA is an appropriate framework for modeling ocean circulation with parameterized eddies. The authors test the feasibility of this proposition and investigate the connections between the TWA framework and the conventional framework used in models, where Eulerian mean flow prognostic variables are solved for. Using the TWA framework as a starting point, this study explores the well-known connections between vertical transfer of horizontal momentum by eddy form drag and eddy overturning by the bolus velocity, used by Greatbatch and Lamb and Gent and McWilliams to parameterize eddies. After implementing the TWA framework in an ocean general circulation model, we verify our analysis by comparing the flows in an idealized Southern Ocean configuration simulated using the TWA and conventional frameworks with the same mesoscale eddy parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Vertical Structure and Ice Production Processes of Shallow Convective Postfrontal Clouds over the Southern Ocean in MARCUS. Part II: Modeling Study

Abstract Part I of this series presented a detailed overview of postfrontal mixed-phase clouds observed during the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign. In Part II, we focus on a multiday (23–26 February 2018) case with the aim of understanding ice production as well as model sensitivity to ice process parameterizations using the Weather Research and Forecasting (WRF) Model. The control simulation with the Predicted Particle Properties (P3) microphysics scheme underestimates the ice content and overestimates the supercooled liquid water content, contrary to the bias common in global climate models. The simulations targeted at ice production processes show negligible sensitivity to cloud droplet number concentrations. Further, neither increasing ice nuclei particle (INP) concentrations to an unrealistic level nor adjusting it to MARCUS field estimations alone guarantees more ice production in the model. However, the simulated clouds are found to be highly sensitive to the implementation of immersion freezing, the thresholding of condensation/deposition freezing initiation, and the rime splintering process. By increasing immersion freezing of cloud droplets, relaxing thresholds for condensation/deposition freezing, or removing rime splintering thresholds, the model significantly improves its performance in producing ice. The relaxation of the immersion freezing temperature threshold to the observed cloud-top temperature suggests an in-cloud seeder–feeder mechanism. The results of this work call for an increase in observations of INP, especially over the remote Southern Ocean and at relatively high temperatures, and measurements of ice particle size distributions to better constrain ice nucleating processes in models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Lower length scale model for palladium attack of silicon carbide in TRISO fuel

TRistructural ISOtropic (TRISO) particle fuels rely on silicon carbide (SiC) as the primary barrier for metallic fission product (FP) release. Palladium (Pd) generated by fission degrades the SiC layer, resulting in the formation of lamellar layers of palladium silicides (PdxSi) and carbon (C) perpendicular to the direction of attack. The Pd attack has been hypothesized to be responsible for failure of the SiC layer and enhance FP release. To better understand and quantify Pd attack of SiC in TRISO particles, a multiscale, mechanistic model of Pd transport is being developed by the NEAMS program. Previous work provided an initial hypothesis for modeling lamellar microstructure formation in SiC due to Pd attack using a phase-field model. The work described in this report builds on the previous model by using molecular dynamics (MD) simulations to parameterize the phase-field model kinetics, and build a reduced order model in BISON using the improved mesoscale Pd penetration model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Evaluation of Preindustrial to Present-day Black Carbon and its Albedo Forcing from Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP)

As a part of the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), we evaluate the historical black carbon (BC) aerosols simulated by 8 ACCMIP models against the observations including 12 ice core records, a long-term surface mass concentrations and recent Arctic BC snowpack measurements. We also estimate BC albedo forcing by performing additional simulations using the NCAR Community Land and Sea-Ice model 4 with prescribed meteorology from 1996-2000, which includes the SNICAR BC-snow model. We evaluated the vertical profile of BC snow concentrations from these offline simulations to using recent BC snowpack measurements. Despite using the same BC emissions, global BC burden differs by approximately a factor of 3 among models due to the differences in aerosol removal parameterizations and simulated meteorology among models; 34 Gg to 103 Gg in 1850 and 82 Gg to 315 Gg in 2000. However,models agree well on 2.5~3 times increase in the global BC burden from preindustrial to present-day, which matches with the 2.5 times increase in BC emissions. We find a large model diversity at both NH and SH high latitude regions for BC burden and at SH high latitude regions for deposition fluxes. The ACCMIP simulations match the observed BC mass concentrations well in Europe and North America except at Jungfrauch and Ispra. However, the models fail to capture the Arctic BC seasonality due tosevere underestimations during winter and spring. Compared to recent snowpack measurements, the simulated vertically resolved BC snow concentrations are, on average, within a factor of 2-3 of observations except for Greenland and Arctic Ocean. However, model and observation differ widely due to missing interannual variations in emissions and possibly due to the choice of the prescribed meteorology period (i.e., 1996-2000).

ACCMIP↗
DOE OSTI2020

Simulated wind farm wake sensitivity to configuration choices in the Weather Research and Forecasting model version 3.8.1

Wakes from wind farms can extend over 50 km downwind in stably stratified conditions. These wakes can undermine power production at downwind turbines, adversely impacting revenue. As such, wind farm wake impacts must be considered in wind resource assessments, especially in regions of dense wind farm development. The open-source Weather Research and Forecasting (WRF) numerical weather prediction model includes a wind farm parameterization to estimate wind farm wake effects, but model configuration choices can influence the resulting predictions of wind farm wakes. These choices include vertical resolution, horizontal resolution, and whether or not to include the addition of turbulent kinetic energy generated by the rotating wind turbines. Despite the sensitivity to model configuration, no clear guidance currently exists for these options. Here we compare simulated wind farm wakes produced by varying model configurations with meteorological observations near a land-based wind farm in flat terrain over several diurnal cycles. A WRF configuration comprised of horizontal resolutions of 3 km or 1 km paired with a vertical resolution of 10 m provides the most accurate representation of wind farm wake effects, such as the correct surface warming and elevated wind speed deficit. The inclusion of turbine-generated turbulence is also critical to produce accurate surface warming and should not be omitted.

17 WIND ENERGY↗