Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wake dynamics”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

kynema-sgf [SWR-20-85]

Kynema-SGF (formerly AMR-Wind), wherein SGF stands for structured-grid fluid dynamics, is a massively parallel, block-structured adaptive-mesh, incompressible flow solver. The codebase was initiated in 2019 from incflo. The solver is built on top of the AMReX library. AMReX library provides the mesh data structures, mesh adaptivity, as well as the linear solvers used for solving the governing equations. Kynema-SGF is actively developed and maintained by a dedicated multi-institutional team from Lawrence Berkeley National Laboratory, National Laboratory of the Rockies, and Sandia National Laboratories. The primary applications for Kynema-SGF are: performing large-eddy simulations (LES) of atmospheric boundary layer (ABL) flows, simulating wind farm turbine-wake interactions using actuator disk or actuator line models for turbines, and as a background solver when coupled with a near-body solver (e.g., Kynema-UGF) with overset methodology to perform blade-resolved simulations of multiple wind turbines within a wind farm. For offshore applications, the ability to model the air-sea interaction effects and its impact on the ABL characteristics is another focus for the code development effort. As with other codes in the Kynema ecosystem, Kynema-SGF shares the following objectives: *an open, well-documented implementation of the state-of-the-art computational models for modeling wind farm flow physics at various fidelities that are backed by a comprehensive verification and validation (V&V) process; *be capable of performing the highest-fidelity simulations of flow fields within wind farms; and *be able to leverage the high-performance leadership class computing facilities available at DOE national laboratories.

Ananthan, Shreyas↗
DOE OSTI2021

Network Based Estimation of Wind Farm Power and Velocity Data Under Changing Wind Direction

This paper describes an estimation algorithm for velocity and power output signals in a wind farm under changing wind direction. A graph-theoretic definition describes the wind farm as a collection of nodes (turbines) and time-varying weighted edges (inter-turbine wake propagation) that change as a function of incoming wind direction. The velocity at each turbine is determined through a discrete input-output model. Changes in wind direction serve as the input and the output is defined in terms of a time-varying weighted adjacency matrix that depends on the time-delay of information propagation between turbines. These delays, which are defined in terms of the advection speed of the wind and the distance between the turbines, capture the delayed effect of wind direction changes on the inter-connectivity of the graph as the wind conditions at the farm inlet propagate through the turbine array. An event-based update framework is employed to capture time-dependent topology changes due to shifts in wind direction. Simulation results for dynamically changing wind inlet directions to a circular wind farm are compared to predictions from both the static and dynamic versions of the FLOw Redirection and Induction in Steady State (FLORIS) model. The approach is shown to enable real-time tracking of dynamic changes to wind farm power output within a framework that can be easily integrated into real-time, horizon-based, control strategies that typically do not account for wind direction changes.

distributed↗
DOE OSTI2022

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗
DOE OSTI2022

Understanding the role of flow dynamics in thermoacoustic combustion instability

Thermoacoustic combustion instability is one of the most challenging operational issues in several high-performance, low-emissions combustion technologies, including gas turbines, aircraft engines, rockets, and industrial boilers. Driven by the coupling between combustor acoustics and flame heat release rate fluctuations, thermoacoustic combustion instability can lead to reduced operability, increased emissions, and, in the most extreme cases, catastrophic failure of combustor components. The feedback loop between acoustics and combustion is often facilitated by fluid mechanic oscillations, referred to as “velocity coupling,” whereby acoustic oscillations drive flow fluctuations, which in turn create fluctuations in the flame. The character of these fluid mechanic oscillations is highly dependent on the structure of the flow field and the receptivity of the flow to external excitation. Combustor flow fields use features like fluid recirculation and shear to enhance flame holding and reduce emissions, but these are also the same features that can make the flow receptive to acoustic excitation or even drive self-excited oscillations. Here, in this paper, we discuss the basics of thermoacoustic instability with a focus on the role of hydrodynamic oscillations in typical combustor flows. To facilitate this discussion, we explore the hydrodynamic instability characteristics of several key combustor unit flows (wakes, swirling jets, etc.) and show how the hydrodynamic stability of a flow is an important consideration in determining a combustor’s propensity for thermoacoustic oscillations. Several examples of coupling between hydrodynamics and thermoacoustics are discussed to illustrate this important link. The paper concludes by discussing the potential for designing flow fields that are thermoacoustic instability resistant, either through a reduction in the receptivity of the flow or through nonlinear coupling mechanisms by which self-excited flow instabilities can suppress velocity-coupled combustion oscillations.

42 ENGINEERING↗
DOE OSTI2023

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An Investigation of the Effects of Volume Fraction on Drag Coefficient of Non-Spherical Particles Using PR-DNS

Prediction of the drag coefficient is required in gas-particle multiphase flow modeling and simulation. Experimental data and correlations on the fixed-bed system of spherical particles with high volume fractions for various possible arrangements are available in the literature. However, the effect of volume fraction on the drag coefficient of non-spherical particles is not well studied. In solving the momentum equation, the volume fraction plays a vital role in determining the flow resistances. In this paper, we study the impact of volume fraction in the range of 0.069 to 0.65 on the drag coefficient using the computational fluid dynamics (CFD) simulation of air for Reynold number in the range of 10 to 10000 using particle resolved direct numerical solution (PR-DNS). Regular non-spherical particles such as a cube, tetrahedron, and spheroids are used in this study since their single particle’s drag coefficient data are available in the literature for comparison. For this work, the simulations are carried out in the Ansys Fluent using polyhedral mesh, which consumes significantly less computational time and power. The study showed the sphericity and volume fraction have significant impact on the bed pressure drop and average drag coefficient of the particles in the bed especially in high Reynolds number regime. The bed of the spheroid experiences the lowest drag being the most streamlined particle, and the particles with the edges result in a large drag coefficient due to flow separation at the discontinuity. The vector plots verify this behavior where large wake regions are observed behind the tetrahedron particle.

Mahyawansi, Pratik↗
DOE OSTI2022

Mechanisms of Low-Level Jet Formation in the U.S. Mid-Atlantic Offshore

Low level jets (LLJs) in the atmosphere exhibit a local windspeed maximum inthe boundary layer and are commonly observed both over land and in coastal environments. Because LLJs present strong positive shear beneath the maximum (or "jet nose") and negative shear above the nose, they pose a challenge to future offshore wind technology through their impacts on turbine performance and wakes. Summertime LLJs in the Great Plains have been attributed to frictional decoupling triggering an inertial oscillation, whilecoastal LLJs in California or the North Sea have additional driving mechanisms such as baroclinic forcing. By comparison, LLJs in the coastal US Mid-Atlantic have received less attention, with conflicting evidence about the mechanism that triggers these high shear events. Motivated by future wind energy development in the region, this work elucidates the atmospheric mechanisms of LLJ formation in the US Mid-Atlantic. Using observational data from two floating LiDAR buoys deployed bythe New York State Energy Research and Development Authority in the NY Bight, we identify four sustained LLJ events whose windspeeds and directions are well predicted by a WRF simulation (Weather Research and Forecasting Model). From analysis of the WRF data, we find that the four case studies provide evidence of concurrent inertial oscillation, reduced vertical mixing, and baroclinicity. In addition, we use simple atmospheric dynamics and a single-column-model to compare the relative contribution of each of these mechanisms to triggering the LLJ. By improving understanding and predictability of coastal mid-Atlantic LLJs, this work reduces uncertaintiesof wind energy deployment in the region, aiding the US's transition toward renewable energy.

atmospheric boundary layers↗
DOE OSTI2020

Downward jetting of a dynamic Leidenfrost drop

Jetting is a universal phenomenon frequently observed in nature and industries, for instance in rain drop impact, inkjet printing, spray cooling, fuel atomization, etc. In drop impact on a superheated surface, we observe the formation of a vapor cavity beneath the dynamic Leidenfrost drop and a consecutive downward ejection of a jet into the cavity, using ultrafast x-ray phase contrast imaging. We reveal that the cavity is induced mostly by the retraction of the drop and the jetting is caused by the convergence of capillary waves along the liquid-cavity interface. We find a jetting criterion based on the viscous damping of capillary waves: [OhWe 2 ] ≤ 66 ± 10. These results can provide important insight that leads to understanding and modeling of jets in nature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Sparse thermal data for cellular automata modeling of grain structure in additive manufacturing

Grain growth in the wake of the melt pool formed during alloy-based additive manufacturing (AM) is complex and multifaceted, depending on parameters governing heat transport, fluid flow, and solidification itself. Cellular automata (CA) models have proven effective in providing computationally efficient and physically sound predictions of grain structure for several AM problems, but their efficiency is tied to the performance of heat transport models. CA models use only a small portion of the problem's temperature data (near the moving melt pool boundary), and much of the CA calculations do not affect the final result due to re-melting of material. Coupling of and communication between heat transport and solidification models, and eliminating operations irrelevant towards final grain structure prediction, will be necessary for using these methods for efficient simulation of large parts. In this work, we introduce a procedure of decoupling the CA from temperature field simulation, using files of relevant temperature data written by the heat transport model. This approach is validated against the standard coupling approach using data obtained through the computational fluid dynamics software OpenFOAM. Negligible differences are seen in grain size, volume, and texture distributions for multilayer simulation of test problems, while the quantity of temperature data for these test problems was reduced by four orders of magnitude (from 100s of GB to 10s of MB) and the code performance sped up by a factor of around 50. Furthermore, variability in microstructure as a function of cell size, substrate, time step, and nucleation parameters is studied, and it is found that cell sizes less than or equal to 1.67 μm and sufficiently small time steps yield Αstatistically equivalent microstructures. Finally, a potential use case for this CA approach—the layer-wise convergence in grain structure starting from extremes in initial grain size—is examined. This approach's ability to simulate expected trends in nucleation and epitaxial grain growth for large regions of microstructure, simulated independently of heat transport models themselves, should prove useful for investigation of various microstructure uncertainties and prediction of part-scale experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Reference floating wind array designs for three representative regions

This work presents the systematic development of three open-source reference floating wind array designs. The designs are tailored to representative site conditions for three regions of the United States: Humboldt Bay off the coast of California, the Gulf of Maine, and the Gulf of America. We adopted existing reference designs for the individual 15 MW turbines, semisubmersible floating platforms, substations, mooring systems, and power cables – integrating and adapting them as needed for each location. We adapted existing dynamic cable designs to use larger conductor sizes to meet the arrays' power transmission requirements, and we set up redundant mooring systems for each substation. The layout of each array is a uniform-grid design optimized to approximately minimize the levelized cost of energy (LCOE) within a square lease area while satisfying spatial constraints. These constraints ensure adequate clearances between adjacent turbines and between underwater components during the layout optimization to prevent clashing and ensure that all components reside within the lease boundaries. Substations are included to allow accounting for intra-array cable costs. They are placed within the uniform grid to maintain the navigability of the arrays. For each feasible layout considered, annual energy production and cable routing costs are calculated and updated in the LCOE objective function. After the optimization, we refined the cable routing with a mix of algorithmic and manual methods to ensure that the cables avoid mooring system components and approach the substation with adequate clearances. We confirmed the suitability of each reference array's layout by comparing the wake losses at each wind heading angle to the wind rose, observing that the optimized layouts largely avoid wake losses in the predominant wind directions. These reference arrays provide open-source baseline designs to enable future research and innovation of floating wind technology at the array scale.

17 WIND ENERGY↗
DOE OSTI2024

Vortex bursting and associated twist dynamics on helical vortex tubes and vortex rings

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to$0.0625$, as well as vortex rings with a large radius-to-core size ratio$10$, both at a circulation-based Reynolds number$5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Mechanics↗
DOE OSTI2023

Progress in the validation of rotor aerodynamic codes using field data

Within the framework of the fourth phase of the International Energy Agency (IEA) Wind Task 29, a large comparison exercise between measurements and aeroelastic simulations has been carried out featuring three simulation cases in axial, sheared and yawed inflow conditions. Results were obtained from more than 19 simulation tools originating from 12 institutes, ranging in fidelity from blade element momentum (BEM) to computational fluid dynamics (CFDs) and compared to state-of-the-art field measurements from the 2 MW DanAero turbine. More than 15 different variable types ranging from lifting-line variables to blade surface pressures, loads and velocities have been compared for the different conditions, resulting in over 250 comparison plots. The result is a unique insight into the current status and accuracy of rotor aerodynamic modeling. For axial flow conditions, a good agreement was found between the various code types, where a dedicated grid sensitivity study was necessary for the CFD simulations. However, compared to wind tunnel experiments on rotors featuring controlled conditions, it remains a challenge to achieve good agreement of absolute levels between simulations and measurements in the field. For sheared inflow conditions, uncertainties due to rotational and unsteady effects on airfoil data result in the CFD predictions standing out above the codes that need input of sectional airfoil data. However, it was demonstrated that using CFD-synthesized airfoil data is an effective means to bypass this shortcoming. For yawed flow conditions, it was observed that modeling of the skewed wake effect is still problematic for BEM codes where CFD and free vortex wake codes inherently model the underlying physics correctly. The next step is a comparison in turbulent inflow conditions, which is featured in IEA Wind Task 47. Doing this analysis in cooperation under the auspices of the IEA Wind Technology Collaboration Program (TCP) has led to many mutual benefits for the participants. The large size of the consortium brought ample manpower for the analysis where the learning process by combining several complementary experiences and modeling techniques gave valuable insights that could not be found when the analysis is carried out individually.

17 WIND ENERGY↗
DOE OSTI2025

AMR-Wind: A Performance-Portable, High-Fidelity Flow Solver for Wind Farm Simulations

We present AMR-Wind, a verified and validated high-fidelity computational-fluid-dynamics code for wind farm flows. AMR-Wind is a block-structured, adaptive-mesh, incompressible-flow solver that enables predictive simulations of the atmospheric boundary layer and wind plants. It is a highly scalable code designed for parallel high-performance computing with a specific focus on performance portability for current and future computing architectures, including graphical processing units (GPUs). In this paper, we detail the governing equations, the numerical methods, and the turbine models. Establishing a foundation for the correctness of the code, we present the results of formal verification and validation. The verification studies, which include a novel actuator line test case, indicate that AMR-Wind is spatially and temporally second-order accurate. The validation studies demonstrate that the key physics capabilities implemented in the code, including actuator disk models, actuator line models, turbulence models, and large eddy simulation (LES) models for atmospheric boundary layers, perform well in comparison to reference data from established computational tools and theory. We conclude with a demonstration simulation of a 12-turbine wind farm operating in a turbulent atmospheric boundary layer, detailing computational performance and realistic wake interactions.

17 WIND ENERGY↗
DOE OSTI2023

Analysis of Narrow Bipolar Events Using Mode Decomposition Methods

Multi-resolution analysis methods can reveal the underlying physical dynamics of nonstationary signals, such as those from lightning. In this paper we demonstrate the application of two multi-resolution analysis methods: Ensemble Empirical Mode Decomposition (EEMD) and Variational Mode Decomposition (VMD) in a comparative way in the analysis of electric field change waveforms from lightning. EEMD and VMD decompose signals into a set of Intrinsic Mode Functions (IMFs). The IMFs can be combined using distance and divergence metrics to obtain noise reduction or to obtain new waveforms that isolate the physical processes of interest while removing irrelevant components of the original signal. We apply the EEMD and VMD methods to the observations of three close Narrow Bipolar Events (NBEs) that were reported by Rison et al. The ΔE observations reveal the occurrence of complex oscillatory processes after the main NBE sferic. We show that both EEMD and VMD are able to isolate the oscillations from the main NBE, with VMD being more effective of the two methods since it requires the least user supervision. The oscillations are found to begin at the end of the NBEs' downward fast positive breakdown and appear to be produced by a half-wavelength standing wave within a weakly-conducting resonant ionization cavity left behind in the wake of the streamer-based NBE event. In conclusion, additional analysis shows that one of the NBEs was likely initiated by an energetic cosmic ray shower, and also corrects a misinterpretation in the literature that fast breakdown is an artifact of NBE-like events in interferometer observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2021

The effects of laser polarization and wavelength on injection dynamics of a laser wakefield accelerator

Here, we investigate the effects of laser polarization and wavelength on electron injection dynamics in a laser wakefield accelerator. During the ionization process, electrons gain residual momentum and kinetic energy via above threshold ionization, which has a strong dependence on laser polarization. A circularly polarized laser pulse results in a much higher residual momentum and kinetic energy gain for the ionized electrons compared with the linearly polarized case. This residual momentum results in particle injection because of the sensitivity of particle trapping to the initial conditions and enhanced the total injected beam charge in both experiments and particle-in-cell simulations. Due to the strong correlation of above threshold ionization with laser wavelength, in this work we extended the investigation to long wavelength (up to 20 μm) drive pulses using particle-in-cell simulations. Owing to the gain in kinetic energy, it may be expected that the charge trapped would consistently increase for circular polarization with increasing laser wavelength, but this was not observed. Instead, there are oscillations with wavelength in the relative trapped charge between linear and circular polarization cases, which arise because of ionization and heating effects on the plasma. Our studies highlight the complex interplay between several different physical effects, including injection regimes—above threshold ionization assisted injection, wave-breaking injection by carrier-envelope-phase effects and ionization injection—ionization gradient induced laser pulse evolution, and thermal modifications to the wake structure that need considering when extrapolating laser wakefield acceleration to different wavelength regimes.

Ma, Y. (ORCID:0000000239194856)↗