SEARCH · Engineering Papers

Results for “CFD simulations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

DOE OSTI2025

Dynamic separation of gases using microsieves

Separation of light weight molecules, such as nitrogen, argon, and oxygen, from heavier compounds can have significant impacts on energy capture, environmental monitoring, or isotopic applications. Large-scale gas separation techniques, like gas centrifugation and membrane mitigation, can be problematic as they impart tremendous energy and induce high mechanical stress onto the instrumentation. Microsieves, also known as micronozzles or microfunnels, are developed to create physical barriers to separate specific isotopes and gases. Separation is achieved using a converging and diverging micronozzle to impose supersonic gas flow around a curved wall, and it has been used for the separation of heavy actinide isotopes in low weight gas as well as separation of low weight gas compositions of nitrogen and argon back in 1900s. However, systematic reviews of this unique technology are lacking. The application of the Laval style nozzle, which has a converging/diverging entrance fundamental to the micronozzle, is included in this review due to its importance in industrial applications in uranium (U) isotope refinement. Using advanced computational fluid dynamic (CFD) simulations, the extent of gas separation can be modelled. Herein, we first examine the literature and survey recent advances on fabrication techniques for creating curved micronozzles, methods and separation principles used to design devices. Furthermore, we then follow with highlights of CFD simulations applied to evaluate the separation effects using microsieves. Finally, identification of the gap and recommendation for future development and applications are suggested for using intrinsic molecular features and fluidic dynamics in formulating separation strategies.

30 Microfluidics↗

DOE OSTI2021

EGR Dilution and Fuel Property Effects on High-Efficiency Spark-Ignition Flames

Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, new engine designs have moved towards increased tumble ratios and stroke-to-bore ratios. Increased tumble ratios correlate positively with increases in turbulent kinetic energy and improved fuel and residual gas mixing, all of which favor faster and more efficient combustion. Longer stroke-to-bore ratios allow higher geometric compression ratios and use of late intake valve closing to control peak compression pressures and temperatures. The addition of dilution to improve efficiency is limited by the resulting increase in combustion instabilities manifested by cycle-to-cycle variability. A number of effects - preferential diffusion, turbulence-combustion interactions, stochastic flow patterns, laminar-turbulent flame kernel transitions, and relative length and velocity scales between flame and turbulence - are believed to be responsible for the increase in cycle-to-cycle variations, where their contributions are likely interlinked. Several studies have shown the influence of stochastic flow characteristics on the nature of combustion instabilities, such as velocity patterns on flame kernel formation and cycle-to-cycle variations in residual gas. However, few have focused on the specific effects of fuel properties. The objective of this work is to contrast the effects of dilution on propane stoichiometric combustion against gasoline. Dilution tolerance experiments were conducted in a purpose-built high stroke-to-bore ratio single cylinder engine with both gasoline and LPG. Three-dimensional full cycle computational fluid dynamics (CFD) simulations employing a level-set combustion approach and Reynolds averaged Navier-Stokes (RANS) turbulence modeling was used to qualitatively assess the changes in length and velocity scales for turbulence and the flame. The experimental results showed that LPG can tolerate higher exhaust gas recirculation (EGR) dilution under a variety of conditions. Analysis of CFD simulations showed that propane flames are likely less sensitive to influences from the flow field due to less thickening of the flame and higher effective flame speeds.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Microstructural and infiltration properties of woven preforms during chemical vapor infiltration

Interface-resolved direct numerical simulations (DNSs) of chemical vapor infiltration (CVI) have been performed over a range of furnace-operating conditions (Thiele moduli) and for practical woven preform geometries. A level-set method is used to resolve the geometry of the initial preform at tow scale. The interface between the vapor and solid phase is then evolved in time through the entire CVI densification cycle, fully resolving the time-varying topology between the two phases. In contrast to previous level-set methods for CVI simulation, the physical reaction and diffusion processes govern the level-set movement in the current approach. The surface deposition kinetics is described by the usual one-step model. In this paper, the DNS data are used to study the evolving porosity, surface-to-volume ratio, and flow infiltration properties (permeability and effective diffusivities). Additionally, comparisons are made to popularly-assumed structure functions and the standard, Kozeny–Carmen porous media model commonly employed in modeled CFD simulations of CVI. The virtual DNS experiments reveal a Thiele modulus and preform geometry (fabric layup) dependence which the existing microstructural and infiltration models are not able to describe throughout the entire densification process. The DNS-based, woven geometry-specific correlations can be applied directly to mean-field, furnace-scale CFD simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2024

CFD and systems engineering to minimize membrane-based carbon capture costs

In this study, we propose a superstructure for membrane-based carbon capture, addressing these challenges through the formulation of an optimization problem. The case study considers a flue gas with a flowrate of 1000 mol/s and a 30 % CO2 (molar basis), which is representative of steel plants. Reduced models, derived from rigorous CFD simulations of membrane modules, are utilized to enhance the efficiency of the optimization process. The reduced model includes specific input variables such as inlet flow rate, retentate and permeate pressures, and inlet CO2 concentration; while the CO2 recovery and purity in the permeate side are considered as output variables. In this way, we exploit the information from the CFD simulations in the superstructure, aiding in the selection of the optimal configuration for the multi-stage membrane process. The results underscore the efficiency of a multi-stage membrane design, showing capture costs around 45 $/t-CO2.

Pedrozo, Hector A.↗

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 OSTI2023

EXPERIMENTAL AND COMPUTATIONAL STUDY OF DIRECT CONTACT MEMBRANE DISTILLATION

Experimental and computational studies of direct contact membrane distillation (DCMD) are conducted. The permeate flux is measured in a module containing a flat sheet membrane. Polytetrafluoroethylene (PTFE) hydrophobic membrane (Membrane Solutions, FPB045A16) with a nominal pore size of 0.45 μm, the thickness of 179~239 μm, the porosity of 80% is used in the experiments. The feed solution is a binary mixture of water and NaCl. The inlet feed temperature is fixed at 70°C, and the feed Reynolds number is varied between 500 and 2000. Two inlet feed concentrations, 0 and 35g/L, are considered. The permeate flow rate was kept constant at Re of 1300 and temperature of 20°C. Three-dimensional steady-state computational fluid dynamics (CFD) simulations are performed using the laminar model. The Navier’s Stokes, energy, and mass transport equations in each channel coupled with flux boundary conditions imposed at the membrane surface are solved. The predicted module-averaged flux is compared to measured data for all operating conditions. Measured and predicted flux agrees well; validating both experiments and model, Temperature polarization in the feed and permeate channel and concentration polarization in the feed channel are characterized by validated CFD simulations.

Oztekin, Alparslan↗

DOE OSTI2025

Design Overview of a High-Pressure Helium Flow Visualization Apparatus for Blanket Cooling Studies

Cooling of the fusion blanket first wall remains a significant challenge given the adverse conditions of heat and particle flux encountered near the plasma. Helium emerges as an attractive cooling candidate because of its chemical and neutronic inertness and separability from hydrogenic species (e.g. tritium). Because of the low thermal mass of helium, optimization of these coolant channels is warranted to provide high heat transfer performance at low pumping costs. Increasingly, computational fluid dynamics (CFD) simulations are employed to model and optimize these flow channels, and accompanying experimental data are needed to validate the predictions of these models. To provide the aforementioned experimental data, a high-pressure helium flow visualization upgrade has been designed for the Helium Flow Loop Experiment facility. This apparatus was built to American Society of Mechanical Engineers boiler and pressure vessel standards to withstand operating pressure of 4 MPa and mated to high-pressure glass windows. Seedless flow visualization is performed via high-speed background oriented schlieren (BOS), with image correlation used for time-resolved two-dimensional velocimetry at frequencies in excess of 60 kHz. Rectangular flow channel test articles are additively manufactured via laser powder bed fusion and installed into this visualization apparatus, with one-sided heating supplied by resistive heaters. In conclusion, the chosen test geometries were informed by prior CFD simulations, and the helium flow structures observed via BOS (detachment, recirculation, etc.) will be used for the validation of these accompanying models, in support of the design and optimization of blanket cooling channel configurations.

Helium flow↗

DOE OSTI2023

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

Preliminary Design of Reactor Cavity Cooling System for a Horizontal Compact HTGR

The Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed MIGHTR in 3 years and support its commercialization as a safe and low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. This report documents the preliminary design study of the RCCS for the HC-HTGR. It includes the establishment of the design requirements, a high-level design study by initial scoping calculations, and preliminary performance calculations of the HC-HTGR RCCS design. Design requirements for the HC-HTGR RCCS have been established to guide preliminary design activities and scoping performance calculations. Initial scoping calculations including estimation of the water inventory, estimation of HVAC thermal capability, and a parametric study on loop dimensions by standalone RCCS analysis. Based on scoping calculation results, a set of baseline dimensions of the HC-HTGR RCCS was derived. A water panel modeling approach was investigated to explore various potential design options for the water panel under consideration for the HC-HTGR RCCS using RELAP5-3D. A test case study was performed to assess the prediction capability of two modeling approaches. The results were compared with CFD simulations conducted in constant RPV temperature and heat flux boundary conditions. It confirms the capability of the RELAP5-3D modeling approach to include all important heat transfer mechanisms expected in the HC-HTGR RCCS operation conditions. Then, a reference RELAP5-3D model for the 1/8 th of a compartment of the preliminary design of the HC-HTGR RCCS was developed. A preliminary performance analysis was conducted to evaluate single-phase natural circulation performance with different top tank temperature values and panel conduction performance in various operation conditions. From single-phase natural circulation performance analysis, the system operation mode was investigated in normal operating and limiting design conditions. It showed operation mode in a subcooled state with a proper top tank water cooling system. Parasitic heat loss by both internal air flow and RCCS was estimated, showing it satisfies maintaining below target maximum heat loss of the HC-HTGR RCCS. From the panel conduction performance analysis, two candidate materials for the riser tube such as carbon steel and stainless steel were compared in the thermal performance of HC-HTGR RCCS. From a single water panel test compared with CFD simulation results, it was confirmed that the current capability of the RELAP5-3D modeling approach for the water panel predicts the thermal conduction of two different materials of the water panel. Then, system-level thermal performance analysis was performed for 1/8 th of the compartment of the preliminary HC-HTGR RCCS design. It was first observed that the current preliminary HC-HTGR RCCS design had minimal impact on the overall thermal performance of the water panel by changing pipe material from carbon steel to stainless steel. From Argonne’s effort on the ongoing water-based NSTF testing program, several considerations other than the thermal performance point of view were addressed to be considered in selecting pipe materials.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗

DOE OSTI2024

Leading CFD-DEM simulation software for mastering particles

Overview of advances in DEM research and application of radiative modeling advancements in Aspherix.

DEM↗

DOE OSTI2022

CFD-DEM Simulations of Reacting Fluidized Beds with MFIX-Exa

SIAM Conference on Parallel Processing for Scientific Computing, Virtual, February 23-26, 2022

Porcu, Roberto↗

DOE OSTI2022

PIV and CFD-DEM simulation of gas jet fluidization of a 1mm ceramic particle bed

Video recording of presentation

Fullmer, William D.↗

DOE OSTI2020

Code-to-code-to-experiment validation of LES-ALM wind farm simulators

The aim of this work is to present a detailed code-to-code comparison of two Large-Eddy Simulation (LES) solvers. Corresponding experimental measurements are used as a reference to validate the quality of the CFD simulations. The comparison highlights the effects of solver order on the solutions, and it tries to answer the question of whether a high order solver is necessary to capture the main characteristics of a wind farm. Both solvers were used on different grids to study their convergence behavior. While both solvers show a good match with experimental measurements, it appears that the low order solver is more accurate and substantially cheaper in terms of computational cost.

17 WIND ENERGY↗

DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

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 project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. 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. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗

DOE OSTI2020

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

17 WIND ENERGY↗

DOE OSTI2021

Demonstrate moving-grid multi-turbine simulations primarily run on GPUs and propose improvements for successful KPP-2

17 WIND ENERGY↗