Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditional neutrality”

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 73 records · Page 4

Molecular Design of Functional Polymers for Silica Scale Inhibition

Silica polymerization, which involves the condensation reaction of silicic acid, is a fundamental process with wide-ranging implications in biological systems, material synthesis, and scale formation. The formation of a silica-based scale poses significant technological challenges to energy-efficient operations in various industrial processes, including heat exchangers and water treatment membranes. Despite the common strategy of applying functional polymers for inhibiting silica polymerization, the underlying mechanisms of inhibition remain elusive. In this study, we synthesized a series of nitrogen-containing polymers as silica inhibitors and elucidated the role of their molecular structures in stabilizing silicic acids. Polymers with both charged amine and uncharged amide groups in their backbones exhibit superior inhibition performance, retaining up to 430 ppm of reactive silica intact for 8 h under neutral pH conditions. In contrast, monomers of these amine/amide-containing polymers as well as polymers containing only amine or amide functionalities present insignificant inhibition. Molecular dynamics simulations reveal strong binding between the deprotonated silicic acid and a polymer when the amine groups in the polymer are protonated. Notably, an extended chain conformation of the polymer is crucial to prevent proximity between the interacting monomeric silica species, thereby facilitating effective silica inhibition. Furthermore, the hydrophobic nature of alkyl segments in polymer chains disrupts the hydration shell around the polymer, resulting in enhanced binding with ionized silicic acid precursors compared to monomers. Finally, our findings provide novel mechanistic insights into the stabilization of silicic acids with functional polymers, highlighting the molecular design principles of effective inhibitors for silica polymerization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Bicarbonate, Calcium, and pH on the Reactivity of As(V) and U(VI) Mixtures

Natural or anthropogenic processes can increase the concentration of uranium (U) and arsenic (As) above the maximum contaminant levels in water sources. Bicarbonate and calcium (Ca) can have major impacts on U speciation and can affect the reactivity between U and As. We therefore investigated the reactivity of aqueous U and As mixtures with bicarbonate and Ca for acidic and neutral pH conditions. In experiments performed with 1 mM U and As mixtures, 10 mM Ca, and without added bicarbonate (pCO 2 = 3.5), aqueous U decreased to <0.25 mM at pH 3 and 7. Aqueous As decreased the most at pH 3 (~0.125 mM). Experiments initiated with 0.005 mM As and U showed similar trends. X-ray spectroscopy (i.e., XAS and EDX) and diffraction indicated that U-As-Ca- and U-Ca-bearing solids resemble uranospinite [Ca(UO 2 ) 2 (AsO 4 ) 2 ·10H 2 O] and becquerelite [Ca(UO 2 ) 6 O 4 (OH) 6 ·8(H 2 O)]. These findings suggest that U-As-Ca-bearing solids formed in mixed solutions are stable at pH 3. However, the dissolution of U-As-Ca and U-Ca-bearing solids at pH 7 was observed in reactors containing 10 mM bicarbonate and Ca, suggesting a kinetic reaction of aqueous uranyl-calcium-carbonate complexation. Our study provides new insights regarding U and As mobilization for risk assessment and remediation strategies.

54 ENVIRONMENTAL SCIENCES↗

Understanding the Effect of Single Atom Cationic Defect Sites in an Al 2 O 3 (012) Surface on Altering Selenate and Sulfate Adsorption: An Ab Initio Study

Adsorption is a promising under-the-sink selenate remediation technique for distributed water systems. Recently it was shown that adsorption induced water network rearrangement control adsorption energetics on the α-Al 2 O 3 (012) surface. Here, we aim to elucidate the relative importance of the water network effects and surface cation identity on controlling selenate and sulfate adsorption energy using density functional theory calculations. Density functional theory (DFT) calculations predicted the adsorption energies of selenate and sulfate on nine transition metal cations (Sc–Cu) and two alkali metal cations (Ga and In) in the α-Al 2 O 3 (012) surface under simulated acidic and neutral pH conditions. We find that the water network effects had a larger impact on the adsorption energy than the cationic identity. However, cation identity secondarily controlled adsorption. Most cations decreased the adsorption energy, weakening the overall performance, the larger Sc and In cations enabled inner-sphere adsorption in acidic conditions because they relaxed outward from the surface, providing more space for adsorption. Additionally, only Ti induced Se selectivity over S by reducing the adsorbing selenate to selenite but not reducing the sulfate. Altogether, this study indicates that tuning water network structure will likely have a larger impact than tuning cation–selenate interactions for increasing adsorbate effectiveness.

Adsorption↗

Ion heating in the PISCES-RF liquid-cooled high-power, steady-state, helicon plasma device

Radio frequency (RF) driven helicon plasma sources are commonly used for their ability to produce high-density argon plasmas ( n > 10 19 m -3 ) at relatively moderate powers (typical RF power < 2 kW). Typical electron temperatures are <10 eV and typical ion temperatures are <0.6 eV. A newly designed helicon antenna assembly (with concentric, double-layered, fully liquid-cooled RF-transparent windows) operates in steady-state at RF powers up to 10 kW. We report on the dependence of argon plasma density, electron temperature and ion temperature on RF power. At 10 kW, ion temperatures >2 eV in argon plasmas are measured with laser induced fluorescence, which is consistent with a simple volume averaged 0D power balance model. 1D Monte Carlo simulations of the neutral density profile for these plasma conditions show strong neutral depletion near the core and predict neutral temperatures well above room temperatures. The plasmas created in this high-power helicon source (when light ions are employed) are ideally suited for fusion divertor plasma-material interaction studies and negative ion production for neutral beams.

Physics↗

QTL × environment interactions underlie ionome divergence in switchgrass

Abstract Ionomics measures elemental concentrations in biological organisms and provides a snapshot of physiology under different conditions. In this study, we evaluate genetic variation of the ionome in outbred, perennial switchgrass in three environments across the species’ native range, and explore patterns of genotype-by-environment interactions. We grew 725 clonally replicated genotypes of a large full sib family from a four-way linkage mapping population, created from deeply diverged upland and lowland switchgrass ecotypes, at three common gardens. Concentrations of 18 mineral elements were determined in whole post-anthesis tillers using ion coupled plasma mass spectrometry (ICP-MS). These measurements were used to identify quantitative trait loci (QTL) with and without QTL-by-environment interactions (QTLxE) using a multi-environment QTL mapping approach. We found that element concentrations varied significantly both within and between switchgrass ecotypes, and GxE was present at both the trait and QTL level. Concentrations of 14 of the 18 elements were under some genetic control, and 77 QTL were detected for these elements. Seventy-four percent of QTL colocalized multiple elements, half of QTL exhibited significant QTLxE, and roughly equal numbers of QTL had significant differences in magnitude and sign of their effects across environments. The switchgrass ionome is under moderate genetic control and by loci with highly variable effects across environments.

59 BASIC BIOLOGICAL SCIENCES↗

Scaling Behavior of a Turbulent Kinetic Energy Closure Scheme for the Stably Stratified Atmosphere: A Steady-State Analysis

Abstract We present a turbulent kinetic energy (TKE) closure scheme for the stably stratified atmosphere in which the mixing lengths for momentum and heat are not parameterized in the same manner. The key difference is that, while the mixing length for heat tends toward the stability independent mixing length for momentum in neutrally stratified conditions, it tends toward one based on the Brunt–Väisälä time scale and square root of the TKE in the limit of large stability. This enables a unique steady-state solution for TKE to be obtained, which we demonstrate would otherwise be impossible if the mixing lengths were the same. Despite the model’s relative simplicity, it is shown to perform reasonably well with observational data from the 1999 Cooperative Atmosphere–Surface Exchange Study (CASES-99) using commonly employed model constants. Analyzing the scaling behavior of the nondimensional velocity and potential temperature gradients, or of the stability (correction) functions, reveals that for large stability the present model scales in the same manner as the first-order operational scheme of Viterbo et al. Alternatively, it appears as a blend of two cases of the TKE closure scheme of Baas et al. Critically, because a unique steady-state TKE can be obtained, the present model avoids the nonphysical behavior identified in one of the cases of Baas et al.

Meteorology & Atmospheric Sciences↗

Influence of the Indian Ocean Dipole on the large-scale circulation in South America

The influence of each phase of the Indian Ocean Dipole (IOD) on the largescale circulation in South America is investigated using rainfall observations, fully-coupled, large-ensemble, historical simulations (LENS), and forced experiments using the coupled model’s atmospheric component. IOD events often occur when El Niño Southern Oscillation (ENSO), the largest source of interannual variability of precipitation in South America, is active. To distinguish from effects of ENSO, only cases during neutral ENSO conditions are analyzed in LENS and observations. During the positive IOD polarity, a perturbation in the localWalker circulation leads to increased convection over equatorial South America, resulting in wet anomalies in the Amazon basin. This signal is the opposite of what is expected during El Niño events. Tropical convection anomalies in the Indian Ocean also force an extratropical Rossby wave train that reaches subtropical South America. During positive IOD, the moisture flux from the Amazon to central and southeastern Brazil weakens, resulting in a drying of the area associated with the South Atlantic Convergence Zone. Meanwhile, the South Atlantic Subtropical High strengthens, contributing to a drying in southeastern Brazil. During negative IOD, the induced wave train from the Indian Ocean leads to increased moisture transport to the La Plata basin, leading to wet anomalies in the region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating Implementations of the Immersed Boundary Method in the Weather Research and Forecasting Model

The terrain-following coordinate system used by many atmospheric models can cause numerical instabilities due to discretization errors as resolved terrain slopes increase and the grid becomes highly skewed. The immersed boundary (IB) method, which does not require the grid to conform to the terrain, has been shown to alleviate these errors, and has been used successfully for high-resolution atmospheric simulations over steep terrain, including vertical building surfaces. Since many previous applications of IB methods to atmospheric models have used very fine grid resolution (5 m or less), the present study seeks to evaluate IB method performance over a range of grid resolutions and aspect ratios. Two classes of IB algorithms, velocity reconstruction and shear stress reconstruction, are tested within the common framework of the Weather Research and Forecasting (WRF) Model. Performance is evaluated in two test cases, one with flat terrain and the other with the topography of Askervein Hill, both under neutrally stratified conditions. WRF-IB results are compared to similarity theory, observations, and native WRF results. Despite sensitivity to the location at which the IB intersects the model grid, the velocity reconstruction IB method shows consistent performance when used with a hybrid RANS/LES surface scheme. The shear stress reconstruction IB method is not sensitive to the grid intersection, but is less consistent and near-surface velocity errors can occur at coarse resolutions. This study represents an initial investigation of IB method variability across grid resolutions in WRF. Future work will focus on improving IB method performance at intermediate to coarse resolutions.

54 ENVIRONMENTAL SCIENCES↗

Mesoscale fractal whey protein particles derived from microscale linear-shaped protein assemblies (Part 1): Manufacturing method and particle characteristics

Whey protein isolates (WPI) are widely used in processed foods for their versatile functional properties. Modifying the structural properties of proteins by assembling them into mesoscale or microscale particles may improve their functionality and broaden their applications. This study aims to manufacture and characterize mesoscale whey protein particles (WPP) derived from WPI. Two types of WPP, WPP1 (0.05 mL/min) and WPP2 (0.25 mL/min), were prepared through a multistep approach involving liquid antisolvent (LAS) precipitation, heat treatment, and microfluidization. Liquid antisolvent precipitation was performed by injecting a 20% (wt/vol) WPI dispersion (pH 7) into an ethanol-glycerol mixture (75:25, vol/vol) under laminar flow, followed by heat treatment at 80°C for 20 min as a particle hardening step. This process produced stable fiber- and ribbon-shaped whey protein assemblies (WPA), which served as precursors to WPP. Subsequent microfluidization (150 MPa, 6 passages) reduced the size of WPA, yielding mesoscale WPP with irregular morphologies and a more uniform size distribution, as revealed by microscopy and dynamic light scattering. ζ-Potential and fluorescence labeling indicated higher surface charge and surface hydrophobicity of WPP compared with untreated WPI. The WPP showed internal mass fractal and surface fractal structures at larger length scales, analyzed using small-angle X-ray scattering. Fourier transform infrared spectroscopy demonstrated an increased fraction of intermolecular β-sheets in WPP, suggesting that hydrogen bonding contributed to their formation. Gel electrophoresis confirmed that disulfide bonds served as the primary cross-links stabilizing the WPP structure. Furthermore, turbidity measurements showed that WPP exhibited superior colloidal phase stability compared with untreated WPI and maintained high colloidal stability under both acidic and neutral pH conditions.

Antisolvent precipitation↗

On the Injection Scale of the Turbulence in the Partially Ionized Very Local Interstellar Medium

The cascade of magnetohydrodynamic (MHD) turbulence is subject to ion–neutral collisional damping and neutral viscous damping in the partially ionized interstellar medium. By examining the damping effects in the warm and partially ionized local interstellar medium, we find that the interstellar turbulence is damped by neutral viscosity at ~261 au and cannot account for the turbulent magnetic fluctuations detected by Voyager 1 and 2. The MHD turbulence measured by Voyager in the very local interstellar medium (VLISM) should be locally injected in the regime where ions are decoupled from neutrals for its cascade to survive the damping effects. With the imposed ion–neutral decoupling condition and the strong turbulence condition for the observed Kolmogorov magnetic energy spectrum, we find that the turbulence in the VLISM is sub-Alfvénic, and its largest possible injection scale is ~194 au.

79 ASTRONOMY AND ASTROPHYSICS↗

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗

The distribution of neutral hydrogen in the colour–magnitude plane of galaxies

ABSTRACT We present the conditional H i (neutral hydrogen) mass function (HIMF) conditioned on observed optical properties, Mr (r-band absolute magnitude), and Cur (u − r colour), for a sample of 7709 galaxies from Arecibo Legacy Fast ALFA (40 per cent data release – α.40) which overlaps with a common volume in SDSS DR7. Based on the conditional HIMF, we find that the luminous red, luminous blue, and faint blue populations dominate the total HIMF at the high-mass end, knee, and the low-mass end, respectively. We use the conditional HIMF to derive the underlying distribution function of ΩH i (H i density parameter), p(ΩH i), in the colour–magnitude plane of galaxies. The distribution, p(ΩH i), peaks in the blue cloud at $\mathit{ M_{{r}}}^{\text{max}}=-19.25, \mathit{ C_{{ur}}}^{\text{max}}=1.44$ but is skewed. It has a long tail towards faint blue galaxies, and luminous red galaxies. We argue that p(ΩH i) can be used to reveal the underlying relation between cold gas, stellar mass, and the star formation rate in an unbiased way, that is, the derived relation does not suffer from survey or sample selection.

Dutta, Saili (ORCID:0000000288585845)↗

Rapid Crystallization of Zeolites with Controllable Defects: Disentangling Fluoride Concentration and pH Using NH 4 F

Zeolite synthesis is typically conducted either under basic conditions or in neutral fluoride media using hydrofluoric acid (HF). While basic (OH − ) conditions generally result in faster zeolite crystallization, they can also increase the likelihood of framework defects and crystal intergrowths. In contrast, synthesis in neutral fluoride media tends to produce fewer defects because fluoride balances positive charges from structure-directing agents. However, this method often requires significantly longer crystallization times and involves the handling of dangerous HF. In the present study, we pursue the best of both synthesis conditions, rapid syntheses with controllable defect concentrations, by disentangling of mineralizing agent and chargebalancing agent using ammonium fluoride (NH 4 F) as an alternative to HF. We have investigated the use of NH 4 F in the syntheses of siliceous and aluminum-containing zeolite A (LTA, small pore), ZSM-5 (MFI, medium pore), and siliceous Beta(*BEA, large pore). The crystallization times of all four zeolites decreased substantially with an increasing NH 4 F concentration. Crystallization times were reduced from 24 to 4 h (Si-LTA), 96 to 36 h (Al-LTA), 240 to 6 h (ZSM-5), and 24 to 3 h (Si-*BEA). Additionally, increasing the NH 4 F concentration in the synthesis mixtures decreases the defect densities of siliceous zeolites. Raman spectroscopy, along with 29 Si MAS NMR, 19 F MAS NMR, 13 C MAS NMR, and fluorine elemental analysis of Si-LTA samples confirms that the reduction in charged defects (Si−O − ) is due to the higher incorporation of F − within the double four-membered ring (D4R) present in the LTA samples. We show that the accelerated crystallization is due to the role of F − in enhancing the silica mineralization rate (formation of silicon hexafluoride species) and stabilizing D4Rs under basic conditions. As a result, combining basic and fluoride-mediated synthesis could therefore be advantageous for faster zeolite production and improved control over structural properties for a wide variety of zeolite structures.

Anions↗

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗

Insights into oxygen diffusion in rare earth disilicate environmental barrier coatings

Environmental barrier coatings (EBC) are crucial for the use of SiC-based ceramic matrix composites in high-temperature combustion environments, yet knowledge of oxygen diffusion in these coatings is limited. This study investigates oxygen diffusion dynamics in the β-RE 2 Si 2 O 7 system to minimize oxygen penetration in rare earth disilicates. We analyze defect formation energy under varying oxygen conditions, identifying key diffusion mechanisms. In oxygen-rich environments, the most favorable neutral interstitial oxygen diffuses along the [110] direction. In oxygen-poor conditions, neutral oxygen vacancies rotate around Y and Si atoms, exhibiting a diffusivity of 6.59×10 −22 m 2 /s at 1500 K for β-Yb 2 Si 2 O 7 . Under intermediate oxygen levels, charged interstitial oxygen diffuses via concerted interstitialcy along the [001] direction with a diffusivity of 6.21×10 −17 m 2 /s. Additionally, alloying rare earth Y with Er and Yb increases diffusion barriers, contributing to improved EBC performance in extreme environments. The insights gained provides valuable guidance for designing robust coatings tailored to withstand extreme operational environments.

36 MATERIALS SCIENCE↗

Effects of pH on an IDP conformational ensemble explored by molecular dynamics simulation

The conformational ensemble of intrinsically disordered proteins, such as α-synuclein, are responsible for their function and malfunction. Misfolding of α-synuclein can lead to neurodegenerative diseases, and the ability to study their conformations and those of other intrinsically disordered proteins under varying physiological conditions can be crucial to understanding and preventing pathologies. In contrast to well-folded peptides, a consensus feature of IDPs is their low hydropathy and high charge, which makes their conformations sensitive to pH perturbation. We examine a prominent member of this subset of IDPs, α-synuclein, using a divide-and-conquer scheme that provides enhanced sampling of IDP structural ensembles. We constructed conformational ensembles of α-synuclein under neutral (pH ~ 7) and low (pH ~ 3) pH conditions and compared our results with available information obtained from smFRET, SAXS, and NMR studies. Specifically, α-synuclein has been found to in a more compact state at low pH conditions and the structural changes observed are consistent with those from experiments. We also characterize the conformational and dynamic differences between these ensembles and discussed the implication on promoting pathogenic fibril formation. We find that under low pH conditions, neutralization of negatively charged residues leads to compaction of the C-terminal portion of α-synuclein while internal reorganization allows α-synuclein to maintain its overall end-to-end distance. Here, we also observe different levels of intra-protein interaction between three regions of α-synuclein at varying pH and a shift towards more hydrophilic interactions with decreasing pH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

RANS simulations of neutral atmospheric boundary layer flow over complex terrain with comparisons to field measurements

Micro-scale Reynolds-averaged Navier-Stokes (RANS) simulations of the neutral atmospheric boundary layer (ABL) over complex terrain and a comparison of the results with conditionally averaged met-tower data are presented. A robust conditional sampling procedure for the meteorological tower (met-tower) data to identify near-neutral conditions based on a criterion for the turbulence intensity is developed. The conditionally averaged wind data on 14 met-towers are used for the model validation. The ABL flow simulations are conducted over complex terrain which includes a prominent hill using the OpenFOAM-based simulator for on/offshore wind farm applications (SOWFA) with the k-ϵ and the SST k-ω turbulence models. The discretization of the production term in the transport equation for the turbulent kinetic energy (TKE) is modified to greatly reduce the commonly observed nonphysical near surface TKE peak. The driving inflow is generated through an iterative approach using a precursor method to reproduce the measured wind statistics at the reference tower. Both of the RANS models are able to capture the flow behavior windward of the hill. The SST k-ω model predicts more intense flow separation than the k-ϵ model downstream of the steepest sections of the hill. The wind statistics predicted at the location of the met-towers by both of the RANS models are fairly consistent. Overall, the comparisons of the direction, mean, and standard deviation of the wind between the simulations and the tower data show reasonable agreement except for the differences of the mean wind speeds at four met-towers located closer to the main ridge of the hill in a region of strong terrain variations.

17 WIND ENERGY↗

Observations of wind farm wake recovery at an operating wind farm

Abstract. The interplay of momentum surrounding wind farms significantly influences wake recovery, affecting the speed at which wakes return to their freestream velocities. Under stable atmospheric conditions, wind farm wakes can extend over considerable distances, leading to sustained vertical momentum flux downstream, with variations observed throughout the diurnal cycle. Particularly in regions such as the US Great Plains, stable conditions can induce low-level jets (LLJs), impacting wind farm performance and power output. This study examines the implications of wake recovery using long-term observations of vertical momentum flux profiles across diverse atmospheric conditions. In these observations, several key findings were observed, such as (a) LLJ heights being altered downstream of a wind farm, especially when the LLJs are below 250 m above ground level; (b) a notable impact of LLJ height on wake recovery being observed using momentum flux profiles at upwind and downwind locations, wherein LLJs between 250 and 500 m above ground level resulted in larger momentum transfer within the wake (i.e., smaller velocity deficit) compared to LLJs below 250 m above ground level; (c) the largest momentum flux variability being observed during stable atmospheric conditions, with non-negligible variability observed during neutral and unstable atmospheric conditions; (d) detection of wake effects almost always being observed throughout the atmospheric boundary layer height; and finally (e) enhancement of wake recovery being observed in the presence of propagating gravity waves. These insights deepen our understanding of the intricate dynamics governing wake recovery in wind farms, advancing efforts to model and predict their behavior across varying atmospheric contexts. In addition, the performance of large-eddy-simulation-based semi-empirical internal boundary layer height model estimates incorporating real-world atmospheric and turbine inputs was evaluated using observations during LLJ conditions.

17 WIND ENERGY↗