Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A large-scale numerical computation of five wind farms was performed as a part of the American WAKE experimeNt (AWAKEN). This high-fidelity computation used the ExaWind/AMR-Wind LES solver to simulate a 100 km × 100 km domain containing 541 turbines under unstable atmospheric conditions matching previous measurements. The turbines were represented by Joukowski and OpenFAST coupled actuator disk models. Results of this qualitative comparison illustrate the interactions of wind farms with large-scale ABL structures in the flow, as well as the extent of downstream wake penetration in the flow and blockage effects around wind farms.

Abstract not provided.

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, capturing the thin boundary layers, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Nalu-Wind is part of the ExaWind code suite.

Turbine generator power from simulations using Actuator Line Models and Actuator Disk Models with a Filtered Lifting Line Correction are compared to field data of a V27 turbine. Preliminary results of the wake characteristics are also presented. Turbine quantities of interest from traditional ALM and ADM with the Gaussian kernel (ϵ) set at the optimum value for matching power production and that resolve the kernel at all mesh sizes are also presented. The atmospheric boundary layer is simulated using Nalu-Wind, a Large Eddy Simulation code which is part of the ExaWind code suite. The effect of mesh resolution on quantities of interest is also examined.

The Jet Propulsion Laboratory has a keen interest in exploring icy moons in the solar system, particularly Jupiter's Europa. Successful exploration of the moon's surface includes planetary protection initiatives to prevent the introduction of viable organisms from Earth to Europa. To that end, the Europa lander requires a Terminal Sterilization Subsystem (TSS) to rid the lander of viable organisms that would potentially contaminate the moon's environment. Sandia National Laboratories has been developing a TSS architecture, relying heavily on computational models to support TSS development. Sandia's TSS design approach involves using energetic material to thermally sterilize lander components at the end of the mission. A hierarchical modeling approach was used for system development and analysis, where simplified systems were constructed to perform empirical tests for evaluating energetic material formulation development and assist in developing computational models with multiple tiers of physics fidelity. Computational models have been developed using multiple Sandia-native computational tools. Three experimental systems and corresponding computational models have been developed: Tube, Sub-Box Small, and Sub-Box Large systems. This paper presents an explanation of the application context of the TSS along with an overview description of a small portion of the TSS development from a modeling and simulation perspective, specifically highlighting verification, validation, and uncertainty quantification (VVUQ) aspects of the modeling and simulation work. Multiple VVUQ approaches were implemented during TSS development, including solution verification, calibration, uncertainty quantification, global sensitivity analysis, and validation. This paper is not intended to express the design results or parameter values used to model the TSS but to communicate the approaches used and how the results of the VVUQ efforts were used and interpreted to assist system development.

The complexity and associated uncertainties involved with atmospheric-turbine-wake interactions produce challenges for accurate wind farm predictions of generator power and other important quantities of interest (QoIs), even with state-of-the-art high-fidelity atmospheric and turbine models. A comprehensive computational study was undertaken with consideration of simulation methodology, parameter selection, and mesh refinement on atmospheric, turbine, and wake QoIs to identify capability gaps in the validation process. For neutral atmospheric boundary layer conditions, the massively parallel large eddy simulation (LES) code Nalu-Wind was used to produce high-fidelity computations for experimental validation using high-quality meteorological, turbine, and wake measurement data collected at the Department of Energy/Sandia National Laboratories Scaled Wind Farm Technology (SWiFT) facility located at Texas Tech University's National Wind Institute. The wake analysis showed the simulated lidar model implemented in Nalu-Wind was successful at capturing wake profile trends observed in the experimental lidar data.

We study wind turbine wakes of rotors operating at high thrust coefficients (CT > 24/25) using large-eddy simulations with a rotating actuator disk model. Wind turbine wakes at high thrust coefficients are different from wakes at low thrust coefficients. Wakes behave differently at high thrust, with increased turbulence and faster recovery. Lower induction in the wake is achieved because wakes in high-thrust conditions recover much faster than in normal operating conditions. This enhanced recovery is possible thanks to the turbulence generated in the near wake. We explore the mechanism behind this behavior and propose a simple model to reproduce it. We also propose a Gaussian fit for the wakes under high-thrust conditions and use it use it to initialize an Ainslie type model within the FAST.Farm framework.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

Abstract not provided.

Milestone Description: Enhance Nalu-Wind's actuator disc model through hardening, documenting, stress-testing, verifying, and validating. Existing workflows will be improved by reducing the data output stream, and by making the analysis capabilities more modular and generally better. These model capabilities are needed by other A2e areas, namely Wake Dynamics, AWAKEN, and VV&UQ.

We investigate the spatial organization and temporal dynamics of large-scale, coherent structures in turbulent Rayleigh-Bénard convection via direct numerical simulation of a 6.3 aspect-ratio cylinder with Rayleigh and Prandtl numbers of and, respectively. Fourier modal decomposition is performed to investigate the structural organization of the coherent turbulent motions by analysing the length scales, time scales and the underlying dynamical processes that are ultimately responsible for the large-scale structure formation and evolution. We observe a high level of rotational symmetry in the large-scale structure in this study and that the structure is well described by the first four azimuthal Fourier modes. Two different large-scale organizations are observed during the duration of the simulation and these patterns are dominated spatially and energetically by azimuthal Fourier modes with frequencies of 2 and 3. Studies of the transition between these two large-scale patterns, radial and vertical variations in the azimuthal energy spectra, as well as the spatial and modal variations in the system's correlation time are conducted. Rotational dynamics are observed for individual Fourier modes and the global structure with strong similarities to the dynamics that have been reported for unit aspect-ratio domains in prior works. It is shown that the large-scale structures have very long correlation time scales, on the order of hundreds to thousands of free-fall time units, and that they are the primary source for a horizontal inhomogeneity within the system that can be observed during a finite, but a very long-time simulation or experiment.

Abstract not provided.