MC

Matt Churchfield

info

Please Note

9 records found

Journal article (2026) - Mehtab Ahmed Khan, Matthew J. Churchfield, Simon J. Watson
This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fri), and the aspect ratios of the ABL and wind farm (H̃i Sh). Simulation configurations that fully resolve AGWs-capturing at least one wavelength both horizontally and vertically-yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fri, multiple convergence–divergence zones develop under sub-critical conditions (Fri<1.0), while for super-critical conditions (Fri>1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃i: shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fri mitigates blockage, and increasing Sh mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fri, improving primarily with larger H̃i and Sh. ...
Journal article (2025) - M.A. Khan, D.J.N. Allaerts, S.J. Watson, Matthew Churchfield
Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, that contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches heavily relies on ad-hoc processes. For instance, the widely used Rayleigh damping method requires ad-hoc tuning to produce a setup which may be only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct an LES study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics.

In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This in turn will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths, (the representative wavelengths of the entire wave spectrum), to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency.

Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upwards propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow-outflow boundaries is of equal importance to the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time. ...
Journal article (2025) - M.A. Khan, S.J. Watson, D.J.N. Allaerts, Matthew Churchfield
Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (Hi), and inversion layer Froude number (Fri) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation. ...
Journal article (2025) - Luis A. Martinez-Tossas, Dries Allaerts, Emmanuel Branlard, Matthew J. Churchfield
The filtered lifting line theory presents a continuous form of the inviscid momentum equations of flow over a lifting device, such as a wing or rotor blade, using body forces without mathematical singularities. This theory is also consistent with an actuator line representation of a lifting device. In this work, we present a reformulation of the equations in terms of the local flow angle along the line, which allows solving the stand-alone equations using multivariate root-finding algorithms. This approach can be used to obtain a fast, computationally inexpensive solution of the loading distribution along a wing without the need to perform computational fluid dynamic simulations. We study the requirements in terms of resolution in the spanwise direction and establish the criteria for spacing and minimum amount of points required along the blade to obtain converged solutions. The solutions are compared to results from large-eddy simulations, and we observed excellent agreement with less than a percent difference in quantities along the blade between the methods. ...
Journal article (2024) - Majid Bastankhah, M. Becker, Matthew Churchfield, Caroline Draxl, Jay Prakash Goit, M.A. Khan, Luis A Martínez-Tossas, J.W. van Wingerden, S.J. Watson, More Authors...
Journal article (2023) - Sue Ellen Haupt, Branko Kosović, Larry K. Berg, Colleen M. Kaul, Matthew Churchfield, Jeffrey Mirocha, D.J.N. Allaerts, Thomas Brummet, Shannon Davis, More Authors...
The Mesoscale to Microscale Coupling team, part of the U.S. Department of Energy Atmosphere to Electrons (A2e) initiative, has studied various important challenges related to coupling mesoscale models to microscale models for the use case of wind energy development and operation. Several coupling methods and techniques for generating turbulence at the microscale that is subgrid to the mesoscale have been evaluated for a variety of cases. Case studies included flat-terrain, complex-terrain, and offshore environments. Methods were developed to bridge the terra incognita, which scales from about 100ĝ€¯m through the depth of the boundary layer. The team used wind-relevant metrics and archived code, case information, and assessment tools and is making those widely available. Lessons learned and discerned best practices are described in the context of the cases studied for the purpose of enabling further deployment of wind energy. ...
Journal article (2023) - Dries Allaerts, Eliot Quon, Matt Churchfield
Reproducing realistic date- and site-specific unsteady wind conditions in large-eddy simulations is becoming increasingly useful in wind energy. How to run a large-eddy simulation to match observed conditions, however, remains an open research question. One approach that has received considerable attention is mesoscale-to-microscale coupling, in which information about the mesoscale weather, most commonly acquired from a mesoscale numerical weather model, is passed on to a microscale model. In this paper, we demonstrate how the recently developed profile-assimilation technique, a form of mesoscale-to-microscale coupling, can be used to drive large-eddy simulations solely based on observed mean-flow profiles at a single location, bypassing the need for auxiliary mesoscale simulations. The new approach is evaluated for a diurnal cycle at the Scaled Wind Farm Technology site. Observed mean-flow profiles from the ground up to a height of 2 km are reconstructed by aggregating measurements from multiple instruments, and gaps in the data are infilled with natural neighbor interpolation. We perform nine simulations using various forcing approaches to deal with data limitations. The results show that it is indeed possible to drive microscale large-eddy simulation with observations using the profile-assimilation technique, notwithstanding large gaps in virtual potential temperature measurements. However, profile assimilation with vertical smoothing of the error between the desired and actual profiles is required. Without that smoothing, the microscale simulations develop unrealistically high turbulence levels under many situations. Finally, we show that simulated mesoscale data can account for missing observations, although care is needed as both data sources are not necessarily compatible. ...
Journal article (2021) - Caroline Draxl, Dries Allaerts, Eliot Quon, Matt Churchfield
To simulate the airflow through a wind farm across a wide range of atmospheric conditions, microscale models (e.g., large-eddy simulation, LES, models) have to be coupled with mesoscale models, because microscale models lack the atmospheric physical processes to represent time-varying local forcing. Here we couple mesoscale model outputs to a LES solver by applying mesoscale momentum- and temperature-budget components from the Weather Research and Forecasting model to the governing equations of the Simulator fOr Wind Farm Applications model. We test whether averaging the budget components affects the LES results with regard to quantities of interest to wind energy. Our study focuses on flat terrain during a quiescent diurnal cycle. The simulation results are compared with observations from a 200-m tall meteorological tower and a wind-profiling radar, by analyzing time series, profiles, rotor-averaged quantities, and spectra. While results show that averaging reduces the spatio-temporal variability of the mesoscale momentum-budget components, when coupled with the LES model, the mesoscale bias (in comparison with observations of wind speed and direction, and potential temperature) is not reduced. In contrast, the LES technique can correct for shear and veer. In both cases, however, averaging the budget components shows no significant impact on the mean flow quantities in the microscale and is not necessary when coupling mesocale budget components to the LES model. ...
Journal article (2020) - Dries Allaerts, Eliot Quon, Caroline Draxl, Matthew Churchfield
Mesoscale-to-microscale coupling (MMC) aims to address the limited scope of traditional large-eddy simulations by driving the microscale flow with information concerning large-scale weather patterns provided by mesoscale models. We present a new offline MMC technique for horizontally homogeneous microscale flow conditions, in which internal forcing terms are computed based on mesoscale time–height profiles of mean-flow quantities. The advantage of such an approach is that it can be used to drive a microscale simulation with either mesoscale or observational data, and that it does not rely on specific terms in the mesoscale budget equations, which are typically not part of the default output of a mesoscale solver. The performance of the proposed profile assimilation technique is assessed based on the simulation of a typical diurnal cycle over the Scaled Wind Farm Technology site in west Texas. Results indicate that simple data assimilation techniques lead to unphysically high levels of shear and turbulence caused by the algorithm’s inability to cope with inaccuracies in the mesoscale time–height profiles. Modifying the algorithm to account for vertical coherence in the mesoscale source terms gives the microscale solver a greater ability to correct the provided mesoscale time–height profiles, leading to improved predictions of shear and turbulence statistics. The resulting turbulence statistics are in good agreement with meteorological tower observations and simulation results obtained with state-of-the-art coupling techniques using mesoscale budget components. ...