15 Matching Results

Search Results

Advanced search parameters have been applied.

Final Report on the “New Directions in the Variational Multiscale Formulation of Large Eddy Simulation of Turbulence”

Description: In the report we present a summary of the new models and algorithms developed by the PI and the students supported by this grant. These developments are described in detail in ten peer-reviewed journal articles that acknowledge support from this grant.
Date: July 16, 2013
Creator: Oberai, Assad A
Partner: UNT Libraries Government Documents Department

Milestone Report #2: Direct Evaporator Leak and Flammability Analysis Modifications and Optimization of the Organic Rankine Cycle to Improve the Recovery of Waste Heat

Description: The direct evaporator is a simplified heat exchange system for an Organic Rankine Cycle (ORC) that generates electricity from a gas turbine exhaust stream. Typically, the heat of the exhaust stream is transferred indirectly to the ORC by means of an intermediate thermal oil loop. In this project, the goal is to design a direct evaporator where the working fluid is evaporated in the exhaust gas heat exchanger. By eliminating one of the heat exchangers and the intermediate oil loop, the overall ORC system cost can be reduced by approximately 15%. However, placing a heat exchanger operating with a flammable hydrocarbon working fluid directly in the hot exhaust gas stream presents potential safety risks. The purpose of the analyses presented in this report is to assess the flammability of the selected working fluid in the hot exhaust gas stream stemming from a potential leak in the evaporator. Ignition delay time for cyclopentane at temperatures and pressure corresponding to direct evaporator operation was obtained for several equivalence ratios. Results of a computational fluid dynamic analysis of a pinhole leak scenario are given.
Date: September 1, 2013
Creator: Guillen, Donna Post
Partner: UNT Libraries Government Documents Department

Final Report on the “ New Directions in the Variational Multiscale Formulation of Large Eddy Simulation of Turbulence”

Description: The work performed as part of this grant may be described in the context of developing numerical methods for solving an abstract, nonlinear variational problem for which the straightforward Galerkin approximation on a space of reasonable dimension is inaccurate. Target applications include flows with turbulence and systems with shocks. With this in mind, an alternate finite dimensional problem is proposed in which a model term (with unknown parameters) is added to the Galerkin approximation to improve its performance. We have developed new formulations for this model term especially for turbulent flows and for systems with shocks. We have also developed a dynamic approach for evaluating the parameters that appear in this model term. We have conducted several numerical studies to confirm the usefulness of the proposed approach.
Date: July 13, 2007
Creator: Oberai, Assad A.
Partner: UNT Libraries Government Documents Department

Wind Turbine Modeling for Computational Fluid Dynamics: December 2010 - December 2012

Description: With the shortage of fossil fuel and the increasing environmental awareness, wind energy is becoming more and more important. As the market for wind energy grows, wind turbines and wind farms are becoming larger. Current utility-scale turbines extend a significant distance into the atmospheric boundary layer. Therefore, the interaction between the atmospheric boundary layer and the turbines and their wakes needs to be better understood. The turbulent wakes of upstream turbines affect the flow field of the turbines behind them, decreasing power production and increasing mechanical loading. With a better understanding of this type of flow, wind farm developers could plan better-performing, less maintenance-intensive wind farms. Simulating this flow using computational fluid dynamics is one important way to gain a better understanding of wind farm flows. In this study, we compare the performance of actuator disc and actuator line models in producing wind turbine wakes and the wake-turbine interaction between multiple turbines. We also examine parameters that affect the performance of these models, such as grid resolution, the use of a tip-loss correction, and the way in which the turbine force is projected onto the flow field.
Date: July 1, 2013
Creator: Tossas, L. A. M. & Leonardi, S.
Partner: UNT Libraries Government Documents Department

Adding Complex Terrain and Stable Atmospheric Condition Capability to the Simulator for On/Offshore Wind Farm Applications (SOWFA) (Presentation)

Description: This presentation describes changes made to NREL's OpenFOAM-based wind plant aerodynamics solver so that it can compute the stably stratified atmospheric boundary layer and flow over terrain. Background about the flow solver, the Simulator for Off/Onshore Wind Farm Applications (SOWFA) is given, followed by details of the stable stratification/complex terrain modifications to SOWFA, along with some preliminary results calculations of a stable atmospheric boundary layer and flow over a simple set of hills.
Date: June 1, 2013
Creator: Churchfield, M. J.
Partner: UNT Libraries Government Documents Department

Large-Eddy Simulation of Wind-Plant Aerodynamics: Preprint

Description: In this work, we present results of a large-eddy simulation of the 48 multi-megawatt turbines composing the Lillgrund wind plant. Turbulent inflow wind is created by performing an atmospheric boundary layer precursor simulation and turbines are modeled using a rotating, variable-speed actuator line representation. The motivation for this work is that few others have done wind plant large-eddy simulations with a substantial number of turbines, and the methods for carrying out the simulations are varied. We wish to draw upon the strengths of the existing simulations and our growing atmospheric large-eddy simulation capability to create a sound methodology for performing this type of simulation. We have used the OpenFOAM CFD toolbox to create our solver.
Date: January 1, 2012
Creator: Churchfield, M. J.; Lee, S.; Moriarty, P. J.; Martinez, L. A.; Leonardi, S.; Vijayakumar, G. et al.
Partner: UNT Libraries Government Documents Department

Long-term Observations of the Convective Boundary Layer Using Insect Radar Returns at the SGP ARM Climate Research Facility

Description: A long-term study of the turbulent structure of the convective boundary layer (CBL) at the U.S. Department of Energy Atmospheric Radiation Measurement Program (ARM) Southern Great Plains (SGP) Climate Research Facility is presented. Doppler velocity measurements from insects occupying the lowest 2 km of the boundary layer during summer months are used to map the vertical velocity component in the CBL. The observations cover four summer periods (2004-08) and are classified into cloudy and clear boundary layer conditions. Profiles of vertical velocity variance, skewness, and mass flux are estimated to study the daytime evolution of the convective boundary layer during these conditions. A conditional sampling method is applied to the original Doppler velocity dataset to extract coherent vertical velocity structures and to examine plume dimension and contribution to the turbulent transport. Overall, the derived turbulent statistics are consistent with previous aircraft and lidar observations. The observations provide unique insight into the daytime evolution of the convective boundary layer and the role of increased cloudiness in the turbulent budget of the subcloud layer. Coherent structures (plumes-thermals) are found to be responsible for more than 80% of the total turbulent transport resolved by the cloud radar system. The extended dataset is suitable for evaluating boundary layer parameterizations and testing large-eddy simulations (LESs) for a variety of surface and cloud conditions.
Date: August 20, 2009
Creator: Chandra, A S; Kollias, P; Giangrande, S E & Klein, S A
Partner: UNT Libraries Government Documents Department

Mesoscale and Large-Eddy Simulations for Wind Energy

Description: Operational wind power forecasting, turbine micrositing, and turbine design require high-resolution simulations of atmospheric flow over complex terrain. The use of both Reynolds-Averaged Navier Stokes (RANS) and large-eddy (LES) simulations is explored for wind energy applications using the Weather Research and Forecasting (WRF) model. To adequately resolve terrain and turbulence in the atmospheric boundary layer, grid nesting is used to refine the grid from mesoscale to finer scales. This paper examines the performance of the grid nesting configuration, turbulence closures, and resolution (up to as fine as 100 m horizontal spacing) for simulations of synoptically and locally driven wind ramping events at a West Coast North American wind farm. Interestingly, little improvement is found when using higher resolution simulations or better resolved turbulence closures in comparison to observation data available for this particular site. This is true for week-long simulations as well, where finer resolution runs show only small changes in the distribution of wind speeds or turbulence intensities. It appears that the relatively simple topography of this site is adequately resolved by all model grids (even as coarse as 2.7 km) so that all resolutions are able to model the physics at similar accuracy. The accuracy of the results is shown in this paper to be more dependent on the parameterization of the land-surface characteristics such as soil moisture rather than on grid resolution.
Date: February 22, 2011
Creator: Marjanovic, N
Partner: UNT Libraries Government Documents Department

Atmospheric and Wake Turbulence Impacts on Wind Turbine Fatigue Loading: Preprint

Description: Large-eddy simulations of atmospheric boundary layers under various stability and surface roughness conditions are performed to investigate the turbulence impact on wind turbines. In particular, the aeroelastic responses of the turbines are studied to characterize the fatigue loading of the turbulence present in the boundary layer and in the wake of the turbines. Two utility-scale 5 MW turbines that are separated by seven rotor diameters are placed in a 3 km by 3 km by 1 km domain. They are subjected to atmospheric turbulent boundary layer flow and data is collected on the structural response of the turbine components. The surface roughness was found to increase the fatigue loads while the atmospheric instability had a small influence. Furthermore, the downstream turbines yielded higher fatigue loads indicating that the turbulent wakes generated from the upstream turbines have significant impact.
Date: December 1, 2011
Creator: Lee, S.; Churchfield, M.; Moriarty, P.; Jonkman, J. & Michalakes, J.
Partner: UNT Libraries Government Documents Department

CFD validation in OECD/NEA t-junction benchmark.

Description: When streams of rapidly moving flow merge in a T-junction, the potential arises for large oscillations at the scale of the diameter, D, with a period scaling as O(D/U), where U is the characteristic flow velocity. If the streams are of different temperatures, the oscillations result in experimental fluctuations (thermal striping) at the pipe wall in the outlet branch that can accelerate thermal-mechanical fatigue and ultimately cause pipe failure. The importance of this phenomenon has prompted the nuclear energy modeling and simulation community to establish a benchmark to test the ability of computational fluid dynamics (CFD) codes to predict thermal striping. The benchmark is based on thermal and velocity data measured in an experiment designed specifically for this purpose. Thermal striping is intrinsically unsteady and hence not accessible to steady state simulation approaches such as steady state Reynolds-averaged Navier-Stokes (RANS) models.1 Consequently, one must consider either unsteady RANS or large eddy simulation (LES). This report compares the results for three LES codes: Nek5000, developed at Argonne National Laboratory (USA), and Cabaret and Conv3D, developed at the Moscow Institute of Nuclear Energy Safety at (IBRAE) in Russia. Nek5000 is based on the spectral element method (SEM), which is a high-order weighted residual technique that combines the geometric flexibility of the finite element method (FEM) with the tensor-product efficiencies of spectral methods. Cabaret is a 'compact accurately boundary-adjusting high-resolution technique' for fluid dynamics simulation. The method is second-order accurate on nonuniform grids in space and time, and has a small dispersion error and computational stencil defined within one space-time cell. The scheme is equipped with a conservative nonlinear correction procedure based on the maximum principle. CONV3D is based on the immersed boundary method and is validated on a wide set of the experimental and benchmark data. The numerical scheme has a very ...
Date: August 23, 2011
Creator: Obabko, A. V.; Fischer, P. F.; Tautges, T. J.; Karabasov, S.; Goloviznin, V. M.; Zaytsev, M. A. et al.
Partner: UNT Libraries Government Documents Department

Large-Eddy Simulation Study of Wake Propagation and Power Production in an Array of Tidal-Current Turbines: Preprint

Description: This paper presents our initial work in performing large-eddy simulations of tidal turbine array flows. First, a horizontally-periodic precursor simulation is performed to create turbulent flow data. Then that data is used as inflow into a tidal turbine array two rows deep and infinitely wide. The turbines are modeled using rotating actuator lines, and the finite-volume method is used to solve the governing equations. In studying the wakes created by the turbines, we observed that the vertical shear of the inflow combined with wake rotation causes lateral wake asymmetry. Also, various turbine configurations are simulated, and the total power production relative to isolated turbines is examined. Staggering consecutive rows of turbines in the simulated configurations allows the greatest efficiency using the least downstream row spacing. Counter-rotating consecutive downstream turbines in a non-staggered array shows a small benefit. This work has identified areas for improvement, such as the use of a larger precursor domain to better capture elongated turbulent structures, the inclusion of salinity and temperature equations to account for density stratification and its effect on turbulence, improved wall shear stress modelling, and the examination of more array configurations.
Date: July 1, 2012
Creator: Churchfield, M. J.; Li, Y. & Moriarty, P. J.
Partner: UNT Libraries Government Documents Department

Recovery Discontinuous Galerkin Jacobian-free Newton-Krylov Method for all-speed flows

Description: There is an increasing interest to develop the next generation simulation tools for the advanced nuclear energy systems. These tools will utilize the state-of-art numerical algorithms and computer science technology in order to maximize the predictive capability, support advanced reactor designs, reduce uncertainty and increase safety margins. In analyzing nuclear energy systems, we are interested in compressible low-Mach number, high heat flux flows with a wide range of Re, Ra, and Pr numbers. Under these conditions, the focus is placed on turbulent heat transfer, in contrast to other industries whose main interest is in capturing turbulent mixing. Our objective is to develop singlepoint turbulence closure models for large-scale engineering CFD code, using Direct Numerical Simulation (DNS) or Large Eddy Simulation (LES) tools, requireing very accurate and efficient numerical algorithms. The focus of this work is placed on fully-implicit, high-order spatiotemporal discretization based on the discontinuous Galerkin method solving the conservative form of the compressible Navier-Stokes equations. The method utilizes a local reconstruction procedure derived from weak formulation of the problem, which is inspired by the recovery diffusion flux algorithm of van Leer and Nomura [?] and by the piecewise parabolic reconstruction [?] in the finite volume method. The developed methodology is integrated into the Jacobianfree Newton-Krylov framework [?] to allow a fully-implicit solution of the problem.
Date: July 1, 2008
Creator: Park, HyeongKae; Nourgaliev, Robert; Mousseau, Vincent & Knoll, Dana
Partner: UNT Libraries Government Documents Department

Topical report : CFD analysis for the applicability of the natural convection shutdown heat removal test facility (NSTF) for the simulation of the VHTR RCCS.

Description: The Very High Temperature gas cooled reactor (VHTR) is one of the GEN IV reactor concepts that have been proposed for thermochemical hydrogen production and other process-heat applications like coal gasification. The United States Department of Energy has selected the VHTR for further research and development, aiming to demonstrate emissions-free electricity and hydrogen production at a future time. One of the major safety advantages of the VHTR is the potential for passive decay heat removal by natural circulation of air in a Reactor Cavity Cooling System (RCCS). The air-side of the RCCS is very similar to the Reactor Vessel Auxiliary Cooling System (RVACS) that has been proposed for the PRISM reactor design. The design and safety analysis of the RVACS have been based on extensive analytical and experimental work performed at ANL. The Natural Convection Shutdown Heat Removal Test Facility (NSTF) at ANL that simulates at full scale the air-side of the RVACS was built to provide experimental support for the design and analysis of the PRISM RVACS system. The objective of this work is to demonstrate that the NSTF facility can be used to generate RCCS experimental data: to validate CFD and systems codes for the analysis of the RCCS; and to support the design and safety analysis of the RCCS. At this time no reference design is available for the NGNP. The General Atomics (GA) gas turbine - modular helium reactor (GT-MHR) has been used in many analyses as a starting reference design. In the GT-MHR the reactor outlet temperature is 850 C, while the target outlet reactor temperature in VHTR is 1000 C. VHTR scoping studies with a reactor outlet temperature of 1000 C have been performed at GA and INEL. Although the reactor outlet temperature in the VHTR is significantly higher than in the GT-MHR, the ...
Date: May 16, 2007
Creator: Tzanos, C. P. (Nuclear Engineering Division)
Partner: UNT Libraries Government Documents Department