scholarly journals Adjoint Optimization of Wind Plant Layouts

2016 ◽  
Author(s):  
Ryan N. King ◽  
Katherine Dykes ◽  
Peter Graf ◽  
Peter E. Hamlington

Abstract. Using adjoint optimization and three-dimensional Reynolds-averaged Navier Stokes (RANS) simulations, we present a new gradient-based approach for optimally siting wind turbines within utility-scale wind plants. By solving the adjoint equations of the flow model, the gradients needed for optimization are found at a cost that is independent of the number of control variables, thereby permitting optimization of large wind plants with many turbine locations. Moreover, compared to the common approach of superimposing prescribed wake deficits onto linearized flow models, the computational efficiency of the adjoint approach allows the use of higher-fidelity RANS flow models which can capture nonlinear turbulent flow physics within a wind plant. The RANS flow model is implemented in the Python finite element package FEniCS and the derivation of the adjoint equations is automated within the dolfin-adjoint framework. Gradient-based optimization of wind turbine locations is demonstrated on idealized test cases that reveal new optimization heuristics such as rotational symmetry, local speedups, and nonlinear wake curvature effects. Layout optimization is also demonstrated on more complex wind rose shapes, including a full annual energy production (AEP) layout optimization over 36 inflow directions and 5 windspeed bins.

2017 ◽  
Vol 2 (1) ◽  
pp. 115-131 ◽  
Author(s):  
Ryan N. King ◽  
Katherine Dykes ◽  
Peter Graf ◽  
Peter E. Hamlington

Abstract. Using adjoint optimization and three-dimensional steady-state Reynolds-averaged Navier–Stokes (RANS) simulations, we present a new gradient-based approach for optimally siting wind turbines within utility-scale wind plants. By solving the adjoint equations of the flow model, the gradients needed for optimization are found at a cost that is independent of the number of control variables, thereby permitting optimization of large wind plants with many turbine locations. Moreover, compared to the common approach of superimposing prescribed wake deficits onto linearized flow models, the computational efficiency of the adjoint approach allows the use of higher-fidelity RANS flow models which can capture nonlinear turbulent flow physics within a wind plant. The steady-state RANS flow model is implemented in the Python finite-element package FEniCS and the derivation and solution of the discrete adjoint equations are automated within the dolfin-adjoint framework. Gradient-based optimization of wind turbine locations is demonstrated for idealized test cases that reveal new optimization heuristics such as rotational symmetry, local speedups, and nonlinear wake curvature effects. Layout optimization is also demonstrated on more complex wind rose shapes, including a full annual energy production (AEP) layout optimization over 36 inflow directions and 5 wind speed bins.


2021 ◽  
Author(s):  
Michael LoCascio ◽  
Christopher Bay ◽  
Majid Bastankhah ◽  
Garrett Barter ◽  
Paul Fleming ◽  
...  

Abstract. Annual energy production (AEP) is often the objective function in wind plant layout optimization studies. The conventional method to compute AEP for a wind farm is to first evaluate power production for each wind direction and speed using either computational fluid dynamics simulations or engineering wake models. The AEP is then calculated by weighted-averaging (based on the wind rose at the wind farm site) the power produced across all wind directions. We propose a novel formulation for time-averaged wake velocity that incorporates an analytical integral of a wake deficit model across every wind direction. This approach computes the average flow field more efficiently, and layout optimization is an obvious application to exploit this benefit. The clear advantage of this new approach is that the layout optimization produces solutions with comparable AEP performance yet is completed about 700 times faster. The analytical integral and the use of a Fourier expansion to express the wind speed and wind direction frequency create a more smooth solution space for the gradient-based optimizer to excel compared with the discrete nature of the existing weighted-averaging power calculation.


Author(s):  
Benjamin Walther ◽  
Siva Nadarajah

This paper develops the discrete adjoint equations for a turbomachinery RANS solver and proposes a framework for fully-automatic gradient-based constrained aerodynamic shape optimization in a multistage turbomachinery environment. The systematic approach for the development of the discrete adjoint solver is discussed. Special emphasis is put on the development of the turbomachinery specific features of the adjoint solver, i.e. on the derivation of flow-consistent adjoint inlet/outlet boundary conditions and, to allow for a concurrent rotor/stator optimization and stage coupling, on the development of an exact adjoint counterpart to the non-reflective, conservative mixing-plane formulation used in the flow solver. The adjoint solver is validated by comparing its sensitivities with finite difference gradients obtained from the flow solver. A sequential quadratic programming algorithm is utilized to determine an improved blade shape based on the objective function gradient provided by the adjoint solution. The functionality of the proposed optimization method is demonstrated by the redesign of a single-stage transonic compressor. The objective is to maximize the isentropic efficiency while constraining the mass flow rate and the total pressure ratio.


1999 ◽  
Vol 1 ◽  
pp. S86-S86
Author(s):  
R DESIMONE ◽  
G GLOMBITZA ◽  
C VAHL ◽  
H MEINZER ◽  
S HAGL

2019 ◽  
Vol 63 (5) ◽  
pp. 50401-1-50401-7 ◽  
Author(s):  
Jing Chen ◽  
Jie Liao ◽  
Huanqiang Zeng ◽  
Canhui Cai ◽  
Kai-Kuang Ma

Abstract For a robust three-dimensional video transmission through error prone channels, an efficient multiple description coding for multi-view video based on the correlation of spatial polyphase transformed subsequences (CSPT_MDC_MVC) is proposed in this article. The input multi-view video sequence is first separated into four subsequences by spatial polyphase transform and then grouped into two descriptions. With the correlation of macroblocks in corresponding subsequence positions, these subsequences should not be coded in completely the same way. In each description, one subsequence is directly coded by the Joint Multi-view Video Coding (JMVC) encoder and the other subsequence is classified into four sets. According to the classification, the indirectly coding subsequence selectively employed the prediction mode and the prediction vector of the counter directly coding subsequence, which reduces the bitrate consumption and the coding complexity of multiple description coding for multi-view video. On the decoder side, the gradient-based directional interpolation is employed to improve the side reconstructed quality. The effectiveness and robustness of the proposed algorithm is verified by experiments in the JMVC coding platform.


1986 ◽  
Vol 51 (11) ◽  
pp. 2489-2501
Author(s):  
Benitto Mayrhofer ◽  
Jana Mayrhoferová ◽  
Lubomír Neužil ◽  
Jaroslav Nývlt

A model is derived for a multi-stage crystallization with cross-current flows of the solution and the crystals being purified. The purity of the product is compared with that achieved in the countercurrent arrangement. A suitable function has been set up which allows the cross-current and countercurrent flow models to be compared and reduces substantially the labour of computation for the countercurrent arrangement. Using the recrystallization of KAl(SO4)2.12 H2O as an example, it is shown that, when the cross-current and countercurrent processes are operated at the same output, the countercurrent arrangement is more advantageous because its solvent consumption is lower.


Author(s):  
Susanne Charlotta Åberg ◽  
Annika Katarina Åberg ◽  
Kirsti Korkka-Niemi

AbstractGreater complexity in three-dimensional (3D) model structures yields more plausible groundwater recharge/discharge patterns, especially in groundwater/surface-water interactions. The construction of a 3D hydrostratigraphic model prior to flow modelling is beneficial when the hydraulic conductivity of geological units varies considerably. A workflow for 3D hydrostratigraphic modelling with Leapfrog Geo and flow modelling with MODFLOW-NWT was developed. It was used to evaluate how the modelling results for groundwater flow and recharge/discharge patterns differ when using simple or more complex hydrostratigraphic models. The workflow was applied to a study site consisting of complex Quaternary sediments underlain by fractured and weathered crystalline bedrock. Increasing the hydrostratigraphic detail appeared to improve the fit between the observed and simulated water table, and created more plausible groundwater flow patterns. Interlayered zones of low and high conductivity disperse the recharge/discharge patterns, increasing the vertical flow component. Groundwater flow was predominantly horizontal in models in which Quaternary sediments and bedrock were simplified as one layer per unit. It appears to be important to define the interlayered low-conductivity units, which can limit groundwater infiltration and also affect groundwater discharge patterns. Explicit modelling with Leapfrog Geo was found to be effective but time-consuming in the generation of scattered and thin-layered strata.


2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
András L. Szabó ◽  
Bitan Roy

Abstract We compute the effects of strong Hubbardlike local electronic interactions on three-dimensional four-component massless Dirac fermions, which in a noninteracting system possess a microscopic global U(1) ⊗ SU(2) chiral symmetry. A concrete lattice realization of such chiral Dirac excitations is presented, and the role of electron-electron interactions is studied by performing a field theoretic renormalization group (RG) analysis, controlled by a small parameter ϵ with ϵ = d−1, about the lower-critical one spatial dimension. Besides the noninteracting Gaussian fixed point, the system supports four quantum critical and four bicritical points at nonvanishing interaction couplings ∼ ϵ. Even though the chiral symmetry is absent in the interacting model, it gets restored (either partially or fully) at various RG fixed points as emergent phenomena. A representative cut of the global phase diagram displays a confluence of scalar and pseudoscalar excitonic and superconducting (such as the s-wave and p-wave) mass ordered phases, manifesting restoration of (a) chiral U(1) symmetry between two excitonic masses for repulsive interactions and (b) pseudospin SU(2) symmetry between scalar or pseudoscalar excitonic and superconducting masses for attractive interactions. Finally, we perturbatively study the effects of weak rotational symmetry breaking on the stability of various RG fixed points.


Sign in / Sign up

Export Citation Format

Share Document