This article reviews the state-of-the-art numerical calculation of wind turbine wake aerodynamics. Different computational fluid dynamics techniques for modeling the rotor and the wake are discussed. Regarding rotor modeling, recent advances in the generalized actuator approach and the direct model are discussed, as far as it attributes to the wake description. For the wake, the focus is on the different turbulence models that are employed to study wake effects on downstream turbines.

This paper examines a class of explicit finite-difference advection schemes derived along the method of lines. An important application field is large-scale atmospheric transport. The paper therefore focuses on the demand of positivity. For the spatial discretization, attention is confined to conservative schemes using five points per direction. The fourth-order central scheme and the family of K-schemes, comprising the second-order central, the second-order upwind, and the third-order upwind biased, are studied. Positivity is enforced through flux limiting. It is concluded that the limited third-order upwind discretization is the best candidate from the four examined. For the time integration attention is confined to a number of explicit Runge-Kutta methods of orders two up to four. With regard to the demand of positivity, these integration methods turn out to behave almost equally and no best method could be identified. '~' 1995 Academic Press, Inc. l. INTRODUCTIONThe subject of this paper is the numerical solution of the partial differential equation for linear advection of a scalar quantity w in an arbitrary velocity field u, given byLinear advection is an important (classical) problem in computational fluid dynamics and has been the subject of numerous investigations. The central theme is how to approximate the advection term \7 · (uw), such that the resulting errors in both phase and amplitude are minimized and the computational cost is still affordable. An important application we have in mind concerns atmospheric transport of chemical species. Then w represents a concentration or density and u a wind field. In addition to the usual accuracy and efficiency requirements, here the main consideration is that the transported concentrations must remain positive, because in actual applications also chemical reactions are modeled for which positivity is a prerequisite for avoiding non-physical chemical instabilities. We emphasize *The research reported helongs to the projects EUSMOG and CIRK which are carried out in cooperation with the Air Laboratory of the RIVM-Tbe Dutch National Institute of Public Health and Environmental Protection. The RIVM is acknowledged for financial support. 35that the demand of positivity is important and that it severely restricts the choice of method, as it is essentially equivalent to the demand of avoiding numerical under-and overshoots in regions of strong variation.The research objective of this paper is to ex.amine a class of positive, finite-difference advection schemes which we consider promising for atmospheric transport applications and to select from this class the best possible candidate. We hereby follow the method-of-lines approach which means that the spatial discretization and temporal integration are considered separately.For the spatial discretization we confine ourselves to stencils using five points per (spatial) direction. We consider this a good starting point since a 5-point stencil is computationally attractive for the following reasons. First, a 5-point stenc...

a b s t r a c tA new formulation of Kapila's five-equation model for inviscid, non-heat-conducting, compressible two-fluid flow is derived, together with an appropriate numerical method. The new formulation uses flow equations based on conservation laws and exchange laws only. The two fluids exchange momentum and energy, for which exchange terms are derived from physical laws. All equations are written as a single system of equations in integral form. No equation is used to describe the topology of the two-fluid flow. Relations for the Riemann invariants of the governing equations are derived, and used in the construction of an Osher-type approximate Riemann solver. A consistent finite-volume discretization of the exchange terms is proposed. The exchange terms have distinct contributions in the cell interior and at the cell faces. For the exchange-term evaluation at the cell faces, the same Riemann solver as used for the flux evaluation is exploited. Numerical results are presented for two-fluid shock-tube and shock-bubble-interaction problems, the former also for a two-fluid mixture case. All results show good resemblance with reference results.

No abstract