Stochastic models for turbulent reacting flows
Kerstein, A.
1993-12-01
The goal of this program is to develop and apply stochastic models of various processes occurring within turbulent reacting flows in order to identify the fundamental mechanisms governing these flows, to support experimental studies of these flows, and to further the development of comprehensive turbulent reacting flow models.
Approximate Model for Turbulent Stagnation Point Flow.
Dechant, Lawrence
2016-01-01
Here we derive an approximate turbulent self-similar model for a class of favorable pressure gradient wedge-like flows, focusing on the stagnation point limit. While the self-similar model provides a useful gross flow field estimate this approach must be combined with a near wall model is to determine skin friction and by Reynolds analogy the heat transfer coefficient. The combined approach is developed in detail for the stagnation point flow problem where turbulent skin friction and Nusselt number results are obtained. Comparison to the classical Van Driest (1958) result suggests overall reasonable agreement. Though the model is only valid near the stagnation region of cylinders and spheres it nonetheless provides a reasonable model for overall cylinder and sphere heat transfer. The enhancement effect of free stream turbulence upon the laminar flow is used to derive a similar expression which is valid for turbulent flow. Examination of free stream enhanced laminar flow suggests that the rather than enhancement of a laminar flow behavior free stream disturbance results in early transition to turbulent stagnation point behavior. Excellent agreement is shown between enhanced laminar flow and turbulent flow behavior for high levels, e.g. 5% of free stream turbulence. Finally the blunt body turbulent stagnation results are shown to provide realistic heat transfer results for turbulent jet impingement problems.
Eiamsa-ard, Smith [Department of Mechanical Engineering, Faculty of Engineering, Mahanakorn University of Technology, Bangkok 10530 (Thailand); Seemawute, Panida [Department of Civil Engineering, Faculty of Engineering, Mahanakorn University of Technology, Bangkok 10530 (Thailand); Wongcharee, Khwanchit [Department of Chemical Engineering, Faculty of Engineering, Mahanakorn University of Technology, Bangkok 10530 (Thailand)
2010-09-15
Effects of peripherally-cut twisted tape insert on heat transfer, friction loss and thermal performance factor characteristics in a round tube were investigated. Nine different peripherally-cut twisted tapes with constant twist ratio (y/W = 3.0) and different three tape depth ratios (DR = d/W = 0.11, 0.22 and 0.33), each with three different tape width ratios (WR = w/W = 0.11, 0.22 and 0.33) were tested. Besides, one typical twisted tape was also tested for comparison. The measurement of heat transfer rate was conducted under uniform heat flux condition while that of friction factor was performed under isothermal condition. Tests were performed with Reynolds number in a range from 1000 to 20,000, using water as a working fluid. The experimental results revealed that both heat transfer rate and friction factor in the tube equipped with the peripherally-cut twisted tapes were significantly higher than those in the tube fitted with the typical twisted tape and plain tube, especially in the laminar flow regime. The higher turbulence intensity of fluid in the vicinity of the tube wall generated by the peripherally-cut twisted tape compared to that induced by the typical twisted tape is referred as the main reason for achieved results. The obtained results also demonstrated that as the depth ratio increased and width ratio decreased, the heat transfer enhancement increased. Over the range investigated, the peripherally-cut twisted tape enhanced heat transfer rates in term of Nusselt numbers up to 2.6 times (turbulent regime) and 12.8 times (laminar regime) of that in the plain tube. These corresponded to the maximum performance factors of 1.29 (turbulent regime) and 4.88 (laminar regime). (author)
ASCR Workshop on Turbulent Flow Simulations at the Exascale:...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
ASCR Workshop on Turbulent Flow Simulations at the Exascale: Opportunities and Challenges ASCR Workshop on Turbulent Flow Simulations at the Exascale: Opportunities and Challenges...
Direct numerical simulation of turbulent reacting flows
Chen, J.H.
1993-12-01
The development of turbulent combustion models that reflect some of the most important characteristics of turbulent reacting flows requires knowledge about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between the turbulence and the chemistry is so strong in certain regimes that is is very difficult to isolate the role played by one individual phenomenon. Direct numerical simulation (DNS) is an extremely useful tool to study in detail the turbulence-chemistry interactions in certain well defined regimes. Globally, non-premixed flames are controlled by two limiting cases: the fast chemistry limit, where the turbulent fluctuations. In between these two limits, finite-rate chemical effects are important and the turbulence interacts strongly with the chemical processes. This regime is important because industrial burners operate in regimes in which, locally the flame undergoes extinction, or is at least in some nonequilibrium condition. Furthermore, these nonequilibrium conditions strongly influence the production of pollutants. To quantify the finite-rate chemistry effect, direct numerical simulations are performed to study the interaction between an initially laminar non-premixed flame and a three-dimensional field of homogeneous isotropic decaying turbulence. Emphasis is placed on the dynamics of extinction and on transient effects on the fine scale mixing process. Differential molecular diffusion among species is also examined with this approach, both for nonreacting and reacting situations. To address the problem of large-scale mixing and to examine the effects of mean shear, efforts are underway to perform large eddy simulations of round three-dimensional jets.
Optical monitor for observing turbulent flow
Albrecht, Georg F.; Moore, Thomas R.
1992-01-01
The present invention provides an apparatus and method for non-invasively monitoring turbulent fluid flows including anisotropic flows. The present invention uses an optical technique to filter out the rays travelling in a straight line, while transmitting rays with turbulence induced fluctuations in time. The output is two dimensional, and can provide data regarding the spectral intensity distribution, or a view of the turbulence in real time. The optical monitor of the present invention comprises a laser that produces a coherent output beam that is directed through a fluid flow, which phase-modulates the beam. The beam is applied to a temporal filter that filters out the rays in the beam that are straight, while substantially transmitting the fluctuating, turbulence-induced rays. The temporal filter includes a lens and a photorefractive crystal such as BaTiO.sub.3 that is positioned in the converging section of the beam near the focal plane. An imaging system is used to observe the filtered beam. The imaging system may take a photograph, or it may include a real time camera that is connected to a computer. The present invention may be used for many purposes including research and design in aeronautics, hydrodynamics, and combustion.
Statistical theory of turbulent incompressible multimaterial flow
Kashiwa, B.
1987-10-01
Interpenetrating motion of incompressible materials is considered. ''Turbulence'' is defined as any deviation from the mean motion. Accordingly a nominally stationary fluid will exhibit turbulent fluctuations due to a single, slowly moving sphere. Mean conservation equations for interpenetrating materials in arbitrary proportions are derived using an ensemble averaging procedure, beginning with the exact equations of motion. The result is a set of conservation equations for the mean mass, momentum and fluctuational kinetic energy of each material. The equation system is at first unclosed due to integral terms involving unknown one-point and two-point probability distribution functions. In the mean momentum equation, the unclosed terms are clearly identified as representing two physical processes. One is transport of momentum by multimaterial Reynolds stresses, and the other is momentum exchange due to pressure fluctuations and viscous stress at material interfaces. Closure is approached by combining careful examination of multipoint statistical correlations with the traditional physical technique of kappa-epsilon modeling for single-material turbulence. This involves representing the multimaterial Reynolds stress for each material as a turbulent viscosity times the rate of strain based on the mean velocity of that material. The multimaterial turbulent viscosity is related to the fluctuational kinetic energy kappa, and the rate of fluctuational energy dissipation epsilon, for each material. Hence a set of kappa and epsilon equations must be solved, together with mean mass and momentum conservation equations, for each material. Both kappa and the turbulent viscosities enter into the momentum exchange force. The theory is applied to (a) calculation of the drag force on a sphere fixed in a uniform flow, (b) calculation of the settling rate in a suspension and (c) calculation of velocity profiles in the pneumatic transport of solid particles in a pipe.
Adaptive LES Methodology for Turbulent Flow Simulations
Oleg V. Vasilyev
2008-06-12
Although turbulent flows are common in the world around us, a solution to the fundamental equations that govern turbulence still eludes the scientific community. Turbulence has often been called one of the last unsolved problem in classical physics, yet it is clear that the need to accurately predict the effect of turbulent flows impacts virtually every field of science and engineering. As an example, a critical step in making modern computational tools useful in designing aircraft is to be able to accurately predict the lift, drag, and other aerodynamic characteristics in numerical simulations in a reasonable amount of time. Simulations that take months to years to complete are much less useful to the design cycle. Much work has been done toward this goal (Lee-Rausch et al. 2003, Jameson 2003) and as cost effective accurate tools for simulating turbulent flows evolve, we will all benefit from new scientific and engineering breakthroughs. The problem of simulating high Reynolds number (Re) turbulent flows of engineering and scientific interest would have been solved with the advent of Direct Numerical Simulation (DNS) techniques if unlimited computing power, memory, and time could be applied to each particular problem. Yet, given the current and near future computational resources that exist and a reasonable limit on the amount of time an engineer or scientist can wait for a result, the DNS technique will not be useful for more than 'unit' problems for the foreseeable future (Moin & Kim 1997, Jimenez & Moin 1991). The high computational cost for the DNS of three dimensional turbulent flows results from the fact that they have eddies of significant energy in a range of scales from the characteristic length scale of the flow all the way down to the Kolmogorov length scale. The actual cost of doing a three dimensional DNS scales as Re{sup 9/4} due to the large disparity in scales that need to be fully resolved. State-of-the-art DNS calculations of isotropic turbulence
Dynamic Multiscale Averaging (DMA) of Turbulent Flow
Richard W. Johnson
2012-09-01
A new approach called dynamic multiscale averaging (DMA) for computing the effects of turbulent flow is described. The new method encompasses multiple applications of temporal and spatial averaging, that is, multiscale operations. Initially, a direct numerical simulation (DNS) is performed for a relatively short time; it is envisioned that this short time should be long enough to capture several fluctuating time periods of the smallest scales. The flow field variables are subject to running time averaging during the DNS. After the relatively short time, the time-averaged variables are volume averaged onto a coarser grid. Both time and volume averaging of the describing equations generate correlations in the averaged equations. These correlations are computed from the flow field and added as source terms to the computation on the next coarser mesh. They represent coupling between the two adjacent scales. Since they are computed directly from first principles, there is no modeling involved. However, there is approximation involved in the coupling correlations as the flow field has been computed for only a relatively short time. After the time and spatial averaging operations are applied at a given stage, new computations are performed on the next coarser mesh using a larger time step. The process continues until the coarsest scale needed is reached. New correlations are created for each averaging procedure. The number of averaging operations needed is expected to be problem dependent. The new DMA approach is applied to a relatively low Reynolds number flow in a square duct segment. Time-averaged stream-wise velocity and vorticity contours from the DMA approach appear to be very similar to a full DNS for a similar flow reported in the literature. Expected symmetry for the final results is produced for the DMA method. The results obtained indicate that DMA holds significant potential in being able to accurately compute turbulent flow without modeling for practical
Compound cooling flow turbulator for turbine component
Lee, Ching-Pang; Jiang, Nan; Marra, John J; Rudolph, Ronald J
2014-11-25
Multi-scale turbulation features, including first turbulators (46, 48) on a cooling surface (44), and smaller turbulators (52, 54, 58, 62) on the first turbulators. The first turbulators may be formed between larger turbulators (50). The first turbulators may be alternating ridges (46) and valleys (48). The smaller turbulators may be concave surface features such as dimples (62) and grooves (54), and/or convex surface features such as bumps (58) and smaller ridges (52). An embodiment with convex turbulators (52, 58) in the valleys (48) and concave turbulators (54, 62) on the ridges (46) increases the cooling surface area, reduces boundary layer separation, avoids coolant shadowing and stagnation, and reduces component mass.
Kinetic Theory of Turbulent Multiphase Flow | The Ames Laboratory
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Kinetic Theory of Turbulent Multiphase Flow FWP/Project Description: Project Leader(s): Rodney Fox It is proposed to further the present understanding of turbulent gas-solid fluidized-bed reactors from the conceptual standpoint of kinetic theory and turbulence modeling. The primary purpose is to provide a theoretical underpinning for the construction of computer codes to better understand and predict multiphase flow behavior in polydisperse gas-solid fluidized-bed reactors and risers. In
Subgrid models for mass and thermal diffusion in turbulent mixing
Sharp, David H; Lim, Hyunkyung; Li, Xiao - Lin; Gilmm, James G
2008-01-01
We are concerned with the chaotic flow fields of turbulent mixing. Chaotic flow is found in an extreme form in multiply shocked Richtmyer-Meshkov unstable flows. The goal of a converged simulation for this problem is twofold: to obtain converged solutions for macro solution features, such as the trajectories of the principal shock waves, mixing zone edges, and mean densities and velocities within each phase, and also for such micro solution features as the joint probability distributions of the temperature and species concentration. We introduce parameterized subgrid models of mass and thermal diffusion, to define large eddy simulations (LES) that replicate the micro features observed in the direct numerical simulation (DNS). The Schmidt numbers and Prandtl numbers are chosen to represent typical liquid, gas and plasma parameter values. Our main result is to explore the variation of the Schmidt, Prandtl and Reynolds numbers by three orders of magnitude, and the mesh by a factor of 8 per linear dimension (up to 3200 cells per dimension), to allow exploration of both DNS and LES regimes and verification of the simulations for both macro and micro observables. We find mesh convergence for key properties describing the molecular level of mixing, including chemical reaction rates between the distinct fluid species. We find results nearly independent of Reynolds number for Re 300, 6000, 600K . Methodologically, the results are also new. In common with the shock capturing community, we allow and maintain sharp solution gradients, and we enhance these gradients through use of front tracking. In common with the turbulence modeling community, we include subgrid scale models with no adjustable parameters for LES. To the authors' knowledge, these two methodologies have not been previously combined. In contrast to both of these methodologies, our use of Front Tracking, with DNS or LES resolution of the momentum equation at or near the Kolmogorov scale, but without resolving the
Quantitative imaging of turbulent and reacting flows
Paul, P.H.
1993-12-01
Quantitative digital imaging, using planar laser light scattering techniques is being developed for the analysis of turbulent and reacting flows. Quantitative image data, implying both a direct relation to flowfield variables as well as sufficient signal and spatial dynamic range, can be readily processed to yield two-dimensional distributions of flowfield scalars and in turn two-dimensional images of gradients and turbulence scales. Much of the development of imaging techniques to date has concentrated on understanding the requisite molecular spectroscopy and collision dynamics to be able to determine how flowfield variable information is encoded into the measured signal. From this standpoint the image is seen as a collection of single point measurements. The present effort aims at realizing necessary improvements in signal and spatial dynamic range, signal-to-noise ratio and spatial resolution in the imaging system as well as developing excitation/detection strategies which provide for a quantitative measure of particular flowfield scalars. The standard camera used for the study is an intensified CCD array operated in a conventional video format. The design of the system was based on detailed modeling of signal and image transfer properties of fast UV imaging lenses, image intensifiers and CCD detector arrays. While this system is suitable for direct scalar imaging, derived quantities (e.g. temperature or velocity images) require an exceptionally wide dynamic range imaging detector. To apply these diagnostics to reacting flows also requires a very fast shuttered camera. The authors have developed and successfully tested a new type of gated low-light level detector. This system relies on fast switching of proximity focused image-diode which is direct fiber-optic coupled to a cooled CCD array. Tests on this new detector show significant improvements in detection limit, dynamic range and spatial resolution as compared to microchannel plate intensified arrays.
Gilmore, Mark A.
2013-06-27
Final Report for grant DE-FG02-06ER54898. The dynamics and generation of intermittent plasma turbulent structures, widely known as "blobs" have been studied in the presence of sheared plasma flows in a controlled laboratory experiment.
Turbulence-chemistry interactions in reacting flows
Barlow, R.S.; Carter, C.D.
1993-12-01
Interactions between turbulence and chemistry in nonpremixed flames are investigated through multiscalar measurements. Simultaneous point measurements of major species, NO, OH, temperature, and mixture fraction are obtained by combining spontaneous Raman scattering, Rayleigh scattering, and laser-induced fluorescence (LIF). NO and OH fluorescence signals are converted to quantitative concentrations by applying shot-to-shot corrections for local variations of the Boltzmann fraction and collisional quenching rate. These measurements of instantaneous thermochemical states in turbulent flames provide insights into the fundamental nature of turbulence-chemistry interactions. The measurements also constitute a unique data base for evaluation and refinement of turbulent combustion models. Experimental work during the past year has focused on three areas: (1) investigation of the effects of differential molecular diffusion in turbulent combustion: (2) experiments on the effects of Halon CF{sub 3}Br, a fire retardant, on the structure of turbulent flames of CH{sub 4} and CO/H{sub 2}/N{sub 2}; and (3) experiments on NO formation in turbulent hydrogen jet flames.
Direct Numerical Simulations of High Reynolds Number Turbulent Channel Flow
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
| Argonne Leadership Computing Facility Visualization of the spanwise vorticity in a turbulent channel. S. Hoyas and O. Flores while they were at Universidad Politecnica de Madrid Direct Numerical Simulations of High Reynolds Number Turbulent Channel Flow PI Name: Robert Moser PI Email: rmoser@ices.utexas.edu Institution: University of Texas Allocation Program: INCITE Allocation Hours at ALCF: 175 Million Year: 2013 Research Domain: Engineering Approximately 28% of U.S. energy resources are
Simulations of Turbulent Flows with Strong Shocks and Density Variations
Zhong, Xiaolin
2012-12-13
In this report, we present the research efforts made by our group at UCLA in the SciDAC project Simulations of turbulent flows with strong shocks and density variations. We use shock-fitting methodologies as an alternative to shock-capturing schemes for the problems where a well defined shock is present. In past five years, we have focused on development of high-order shock-fitting Navier-Stokes solvers for perfect gas flow and thermochemical non-equilibrium flow and simulation of shock-turbulence interaction physics for very strong shocks. Such simulation has not been possible before because the limitation of conventional shock capturing methods. The limitation of shock Mach number is removed by using our high-order shock-fitting scheme. With the help of DOE and TeraGrid/XSEDE super computing resources, we have obtained new results which show new trends of turbulence statistics behind the shock which were not known before. Moreover, we are also developing tools to consider multi-species non-equilibrium flows. The main results are in three areas: (1) development of high-order shock-fitting scheme for perfect gas flow, (2) Direct Numerical Simulation (DNS) of interaction of realistic turbulence with moderate to very strong shocks using super computing resources, and (3) development and implementation of models for computation of mutli-species non-quilibrium flows with shock-fitting codes.
A generic model for transport in turbulent shear flows
Newton, Andrew P. L.; Kim, Eun-Jin [Department of Applied Mathematics, University of Sheffield, Sheffield, S3 7RH (United Kingdom)
2011-05-15
Turbulence regulation by large-scale shear flows is crucial for a predictive modeling of transport in plasma. In this paper the suppression of turbulent transport by large-scale flows is studied numerically by measuring the turbulent diffusion D{sub t} and scalar amplitude
Lynn, Jacob W.; Quataert, Eliot; Chandran, Benjamin D. G.; Parrish, Ian J.
2014-08-10
We use analytic estimates and numerical simulations of test particles interacting with magnetohydrodynamic (MHD) turbulence to show that subsonic MHD turbulence produces efficient second-order Fermi acceleration of relativistic particles. This acceleration is not well described by standard quasi-linear theory but is a consequence of resonance broadening of wave-particle interactions in MHD turbulence. We provide momentum diffusion coefficients that can be used for astrophysical and heliospheric applications and discuss the implications of our results for accretion flows onto black holes. In particular, we show that particle acceleration by subsonic turbulence in radiatively inefficient accretion flows can produce a non-thermal tail in the electron distribution function that is likely important for modeling and interpreting the emission from low-luminosity systems such as Sgr A* and M87.
On the laminar to turbulent flow transition in diabatic helically coiled pipe flow
Cioncolini, Andrea; Santini, Lorenzo [Department of Nuclear Engineering, Politecnico di Milano, via Ponzio 34/3, 20133 Milano (Italy)
2006-07-15
Recently the authors experimentally investigated the turbulence emergence process in adiabatic coiled pipe flow. The results of such an investigation compared favorably with existing experimental evidence and revealed as well some new and striking features of the turbulence emergence process in coiled pipes that were not observed in previous research. The objective of the present investigation is to confirm such findings with diabatic flow through coiled pipes. (author)
Large-eddy simulation of turbulent circular jet flows
Jones, S. C.; Sotiropoulos, F.; Sale, M. J.
2002-07-01
This report presents a numerical method for carrying out large-eddy simulations (LES) of turbulent free shear flows and an application of a method to simulate the flow generated by a nozzle discharging into a stagnant reservoir. The objective of the study was to elucidate the complex features of the instantaneous flow field to help interpret the results of recent biological experiments in which live fish were exposed to the jet shear zone. The fish-jet experiments were conducted at the Pacific Northwest National Laboratory (PNNL) under the auspices of the U.S. Department of Energys Advanced Hydropower Turbine Systems program. The experiments were designed to establish critical thresholds of shear and turbulence-induced loads to guide the development of innovative, fish-friendly hydropower turbine designs.
Convection Heat Transfer in Three-Dimensional Turbulent Separated/Reattached Flow
Bassem F. Armaly
2007-10-31
The measurements and the simulation of convective heat transfer in separated flow have been a challenge to researchers for many years. Measurements have been limited to two-dimensional flow and simulations failed to predict accurately turbulent heat transfer in the separated and reattached flow region (prediction are higher than measurements by more than 50%). A coordinated experimental and numerical effort has been initiated under this grant for examining the momentum and thermal transport in three-dimensional separated and reattached flow in an effort to provide new measurements that can be used for benchmarking and for improving the simulation capabilities of 3-D convection in separated/reattached flow regime. High-resolution and non-invasive measurements techniques are developed and employed in this study to quantify the magnitude and the behavior of the three velocity components and the resulting convective heat transfer. In addition, simulation capabilities are developed and employed for improving the simulation of 3-D convective separated/reattached flow. Such basic measurements and simulation capabilities are needed for improving the design and performance evaluation of complex (3-D) heat exchanging equipment. Three-dimensional (3-D) convective air flow adjacent to backward-facing step in rectangular channel is selected for the experimental component of this study. This geometry is simple but it exhibits all the complexities that appear in any other separated/reattached flow, thus making the results generated in this study applicable to any other separated and reattached flow. Boundary conditions, inflow, outflow, and wall thermal treatment in this geometry can be well measured and controlled. The geometry can be constructed with optical access for non-intrusive measurements of the flow and thermal fields. A three-component laser Doppler velocimeter (LDV) is employed to measure simultaneously the three-velocity components and their turbulent fluctuations
Chemical preconcentrator with integral thermal flow sensor
Manginell, Ronald P.; Frye-Mason, Gregory C.
2003-01-01
A chemical preconcentrator with integral thermal flow sensor can be used to accurately measure fluid flow rate in a microanalytical system. The thermal flow sensor can be operated in either constant temperature or constant power mode and variants thereof. The chemical preconcentrator with integral thermal flow sensor can be fabricated with the same MEMS technology as the rest of the microanlaytical system. Because of its low heat capacity, low-loss, and small size, the chemical preconcentrator with integral thermal flow sensor is fast and efficient enough to be used in battery-powered, portable microanalytical systems.
Aerosol deposition in bends with turbulent flow
McFarland, A.R.; Gong, H.; Wente, W.B.
1997-08-01
The losses of aerosol particles in bends were determined numerically for a broad range of design and operational conditions. Experimental data were used to check the validity of the numerical model, where the latter employs a commercially available computational fluid dynamics code for characterizing the fluid flow field and Lagrangian particle tracking technique for characterizing aerosol losses. Physical experiments have been conducted to examine the effect of curvature ratio and distortion of the cross section of bends. If it curvature ratio ({delta} = R/a) is greater than about 4, it has little effect on deposition, which is in contrast with the recommendation given in ANSI N13.1-1969 for a minimum curvature ratio of 10. Also, experimental results show that if the tube cross section is flattened by 25% or less, the flattening also has little effect on deposition. Results of numerical tests have been used to develop a correlation of aerosol penetration through a bend as a function of Stokes number (Stk), curvature ratio ({delta}) and the bend angle ({theta}). 17 refs., 10 figs., 2 tabs.
Flow topologies and turbulence scales in a jet-in-cross-flow
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Oefelein, Joseph C.; Ruiz, Anthony M.; Lacaze, Guilhem
2015-04-03
This study presents a detailed analysis of the flow topologies and turbulence scales in the jet-in-cross-flow experiment of [Su and Mungal JFM 2004]. The analysis is performed using the Large Eddy Simulation (LES) technique with a highly resolved grid and time-step and well controlled boundary conditions. This enables quantitative agreement with the first and second moments of turbulence statistics measured in the experiment. LES is used to perform the analysis since experimental measurements of time-resolved 3D fields are still in their infancy and because sampling periods are generally limited with direct numerical simulation. A major focal point is the comprehensivemoreĀ Ā» characterization of the turbulence scales and their evolution. Time-resolved probes are used with long sampling periods to obtain maps of the integral scales, Taylor microscales, and turbulent kinetic energy spectra. Scalar-fluctuation scales are also quantified. In the near-field, coherent structures are clearly identified, both in physical and spectral space. Along the jet centerline, turbulence scales grow according to a classical one-third power law. However, the derived maps of turbulence scales reveal strong inhomogeneities in the flow. From the modeling perspective, these insights are useful to design optimized grids and improve numerical predictions in similar configurations.Ā«Ā less
Flow topologies and turbulence scales in a jet-in-cross-flow
Oefelein, Joseph C.; Ruiz, Anthony M.; Lacaze, Guilhem
2015-04-03
This study presents a detailed analysis of the flow topologies and turbulence scales in the jet-in-cross-flow experiment of [Su and Mungal JFM 2004]. The analysis is performed using the Large Eddy Simulation (LES) technique with a highly resolved grid and time-step and well controlled boundary conditions. This enables quantitative agreement with the first and second moments of turbulence statistics measured in the experiment. LES is used to perform the analysis since experimental measurements of time-resolved 3D fields are still in their infancy and because sampling periods are generally limited with direct numerical simulation. A major focal point is the comprehensive characterization of the turbulence scales and their evolution. Time-resolved probes are used with long sampling periods to obtain maps of the integral scales, Taylor microscales, and turbulent kinetic energy spectra. Scalar-fluctuation scales are also quantified. In the near-field, coherent structures are clearly identified, both in physical and spectral space. Along the jet centerline, turbulence scales grow according to a classical one-third power law. However, the derived maps of turbulence scales reveal strong inhomogeneities in the flow. From the modeling perspective, these insights are useful to design optimized grids and improve numerical predictions in similar configurations.
Coherent structures in ion temperature gradient turbulence-zonal flow
Singh, Rameswar; Singh, R.; Kaw, P.; Gürcan, Ö. D.; Diamond, P. H.
2014-10-15
Nonlinear stationary structure formation in the coupled ion temperature gradient (ITG)-zonal flow system is investigated. The ITG turbulence is described by a wave-kinetic equation for the action density of the ITG mode, and the longer scale zonal mode is described by a dynamic equation for the m?=?n?=?0 component of the potential. Two populations of trapped and untrapped drift wave trajectories are shown to exist in a moving frame of reference. This novel effect leads to the formation of nonlinear stationary structures. It is shown that the ITG turbulence can self-consistently sustain coherent, radially propagating modulation envelope structures such as solitons, shocks, and nonlinear wave trains.
Turbulence structure in free-surface channel flows
Rashidi, M.; Banerjee, S.
1988-09-01
A turbulence structure in horizontal liquid streams bounded by a free surface and a wall has been investigated using 10--25 ..mu..m oxygen bubbles as tracers. High speed video movies indicate that the dominant flow structure is caused by the periodic ejection of intensely turbulent fluid with low streamwise momentum from the wall region into the relatively quiescent bulk fluid which it displaces and mixes with slowly. The motion of these bursts is constrained by the free interface. Between bursts and the interface a high speed region with a steep velocity gradient develops as a consequence. This in turn causes progress of the burst fluid toward the interface to slow down and eventually to turn back toward the wall, giving rise to characteristic rolling structures, which rotate clockwise if the flow is viewed as going from left to right. To complement the video studies, quantitative data were obtained by analyzing bubble streak lines generated by photography of optically chopped flashes. These data show that in the vicinity of the interface the velocity fluctuations normal to it are damped whereas those parallel to it are enhanced. Analysis of conditional samples of the data indicate that fluid with relatively low streamwise momentum tends to move toward the interface while fluid with high momentum moves away giving rise to rotating structures that roll along with the flow in agreement with the video studies. A high degree of correlation between ejection events near the wall and the fluid motion near the interface also confirm that the bursts extend across the flow stream. This has important implications for surface renewal theories of turbulent transport at fluid--fluid interfaces.
Prediction of turbulent buoyant flow using an RNG {kappa}-{epsilon} model
Gan, G.
1998-02-06
Buoyant flows occur in various engineering practices such as heating, ventilation, and air-conditioning of buildings. This phenomenon is particularly important in rooms with displacement ventilation, where supply air velocities are generally very low (< 0.2 m/s) so that the predominant indoor airflow is largely due to thermal buoyancy created by internal heat sources such as occupants and equipment. This type of ventilation system has been shown to be an effective means to remove excess heat and achieve good indoor air quality. Here, numerical predictions were carried out for turbulent natural convection in two tall air cavities. The standard and RNG {kappa}-{epsilon} turbulence models were used for the predictions. The predicted results were compared with experimental data from the literature, and good agreement between prediction and measurement was obtained. Improved prediction was achieved using the RNG {kappa}-{epsilon} model in comparison with the standard {kappa}-{epsilon} model. The principal parameters for the improvement were investigated.
Ando, Yasutaka; Tobe, Shogo [Ashikaga Institute of Technology, 268-1 Omae, Ashikaga, Tochigi 326-8558 (Japan); Tahara, Hirokazu [Osaka Institute of Technology, 5-16-1 Omiya, Asahi-Ku, Osaka 535-8585 (Japan)
2008-02-21
In this study, to provide continuous plasma atmosphere on the substrate surface in the case of atmospheric thermal plasma CVD, TiO{sub 2} film deposition by thermal plasma CVD using laminar plasma jet was carried out. For comparison, the film deposition using turbulence plasma jet was conducted as well. Consequently, transition of the plasma jet from laminar to turbulent occurred on the condition of over 3.5 1/min in Ar working gas flow rate and the plasma jet became turbulent on the condition of over 10 1/min. In the case of the turbulent plasma jet use, anatase rich titanium oxide film could be obtained though plasma jet could not contact with the surface of the substrate continuously even on the condition that feedstock material was injected into the plasma jet. On the other hand,, in the case of laminar gas flow rate, the plasma jet could contact with the substrate continuously without melt down of the substrate during film deposition. Besides, titanium oxide film could be obtained even in the case of the laminar plasma jet use. From these results, this technique was thought to have high potential for atmospheric thermal plasma CVD.
Rosa, B.; Parishani, H.; Ayala, O.; Wang, L.-P.
2015-01-15
In this paper, we study systematically the effects of forcing time scale in the large-scale stochastic forcing scheme of Eswaran and Pope [āAn examination of forcing in direct numerical simulations of turbulence,ā Comput. Fluids 16, 257 (1988)] on the simulated flow structures and statistics of forced turbulence. Using direct numerical simulations, we find that the forcing time scale affects the flow dissipation rate and flow Reynolds number. Other flow statistics can be predicted using the altered flow dissipation rate and flow Reynolds number, except when the forcing time scale is made unrealistically large to yield a Taylor microscale flow Reynolds number of 30 and less. We then study the effects of forcing time scale on the kinematic collision statistics of inertial particles. We show that the radial distribution function and the radial relative velocity may depend on the forcing time scale when it becomes comparable to the eddy turnover time. This dependence, however, can be largely explained in terms of altered flow Reynolds number and the changing range of flow length scales present in the turbulent flow. We argue that removing this dependence is important when studying the Reynolds number dependence of the turbulent collision statistics. The results are also compared to those based on a deterministic forcing scheme to better understand the role of large-scale forcing, relative to that of the small-scale turbulence, on turbulent collision of inertial particles. To further elucidate the correlation between the altered flow structures and dynamics of inertial particles, a conditional analysis has been performed, showing that the regions of higher collision rate of inertial particles are well correlated with the regions of lower vorticity. Regions of higher concentration of pairs at contact are found to be highly correlated with the region of high energy dissipation rate.
Assessment of a hybrid finite element and finite volume code for turbulent incompressible flows
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Xia, Yidong; Wang, Chuanjin; Luo, Hong; Christon, Mark; Bakosi, Jozsef
2015-12-15
Hydra-TH is a hybrid finite-element/finite-volume incompressible/low-Mach flow simulation code based on the Hydra multiphysics toolkit being developed and used for thermal-hydraulics applications. In the present work, a suite of verification and validation (V&V) test problems for Hydra-TH was defined to meet the design requirements of the Consortium for Advanced Simulation of Light Water Reactors (CASL). The intent for this test problem suite is to provide baseline comparison data that demonstrates the performance of the Hydra-TH solution methods. The simulation problems vary in complexity from laminar to turbulent flows. A set of RANS and LES turbulence models were used in themoreĀ Ā» simulation of four classical test problems. Numerical results obtained by Hydra-TH agreed well with either the available analytical solution or experimental data, indicating the verified and validated implementation of these turbulence models in Hydra-TH. Where possible, we have attempted some form of solution verification to identify sensitivities in the solution methods, and to suggest best practices when using the Hydra-TH code.Ā«Ā less
Assessment of a hybrid finite element and finite volume code for turbulent incompressible flows
Xia, Yidong; Wang, Chuanjin; Luo, Hong; Christon, Mark; Bakosi, Jozsef
2015-12-15
Hydra-TH is a hybrid finite-element/finite-volume incompressible/low-Mach flow simulation code based on the Hydra multiphysics toolkit being developed and used for thermal-hydraulics applications. In the present work, a suite of verification and validation (V&V) test problems for Hydra-TH was defined to meet the design requirements of the Consortium for Advanced Simulation of Light Water Reactors (CASL). The intent for this test problem suite is to provide baseline comparison data that demonstrates the performance of the Hydra-TH solution methods. The simulation problems vary in complexity from laminar to turbulent flows. A set of RANS and LES turbulence models were used in the simulation of four classical test problems. Numerical results obtained by Hydra-TH agreed well with either the available analytical solution or experimental data, indicating the verified and validated implementation of these turbulence models in Hydra-TH. Where possible, we have attempted some form of solution verification to identify sensitivities in the solution methods, and to suggest best practices when using the Hydra-TH code.
Large-eddy simulations of turbulent flow for grid-to-rod fretting in nuclear reactors
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Bakosi, J.; Christon, M. A.; Lowrie, R. B.; Pritchett-Sheats, L. A.; Nourgaliev, R. R.
2013-07-12
The grid-to-rod fretting (GTRF) problem in pressurized water reactors is a flow-induced vibration problem that results in wear and failure of the fuel rods in nuclear assemblies. In order to understand the fluid dynamics of GTRF and to build an archival database of turbulence statistics for various configurations, implicit large-eddy simulations of time-dependent single-phase turbulent flow have been performed in 3 Ć 3 and 5 Ć 5 rod bundles with a single grid spacer. To assess the computational mesh and resolution requirements, a method for quantitative assessment of unstructured meshes with no-slip walls is described. The calculations have been carriedmoreĀ Ā» out using Hydra-TH, a thermal-hydraulics code developed at Los Alamos for the Consortium for Advanced Simulation of Light water reactors, a United States Department of Energy Innovation Hub. Hydra-TH uses a second-order implicit incremental projection method to solve the singlephase incompressible Navier-Stokes equations. The simulations explicitly resolve the large scale motions of the turbulent flow field using first principles and rely on a monotonicity-preserving numerical technique to represent the unresolved scales. Each series of simulations for the 3 Ć 3 and 5 Ć 5 rod-bundle geometries is an analysis of the flow field statistics combined with a mesh-refinement study and validation with available experimental data. Our primary focus is the time history and statistics of the forces loading the fuel rods. These hydrodynamic forces are believed to be the key player resulting in rod vibration and GTRF wear, one of the leading causes for leaking nuclear fuel which costs power utilities millions of dollars in preventive measures. As a result, we demonstrate that implicit large-eddy simulation of rod-bundle flows is a viable way to calculate the excitation forces for the GTRF problem.Ā«Ā less
DNS/LES of Complex Turbulent Flows | Argonne Leadership Computing Facility
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Autoignition of a turbulent hydrogen jet injected into hot oxygen Autoignition of a turbulent hydrogen jet injected into hot oxygen. The green isocontours illustrate the turbulent structures in the fuel, the blue isocontours show the HO2 radical (an important precursor to autoignition), and the orange isosurfaces illustrate temperature. Credit: Rajapadiyan Asaithambi and Krishnan Mahesh, University of Minnesota DNS/LES of Complex Turbulent Flows PI Name: Krishnan Mahesh PI Email: mahesh@umn.edu
Organized Oscillations of Initially-Turbulent Flow Past a Cavity
J.C. Lin; D. Rockwell
2002-09-17
Flow past an open cavity is known to give rise to self-sustained oscillations in a wide variety of configurations, including slotted-wall, wind and water tunnels, slotted flumes, bellows-type pipe geometries, high-head gates and gate slots, aircraft components and internal piping systems. These cavity-type oscillations are the origin of coherent and broadband sources of noise and, if the structure is sufficiently flexible, flow-induced vibration as well. Moreover, depending upon the state of the cavity oscillation, substantial alterations of the mean drag may be induced. In the following, the state of knowledge of flow past cavities, based primarily on laminar inflow conditions, is described within a framework based on the flow physics. Then, the major unresolved issues for this class of flows will be delineated. Self-excited cavity oscillations have generic features, which are assessed in detail in the reviews of Rockwell and Naudascher, Rockwell, Howe and Rockwell. These features, which are illustrated in the schematic of Figure 1, are: (i) interaction of a vorticity concentration(s) with the downstream corner, (ii) upstream influence from this corner interaction to the sensitive region of the shear layer formed from the upstream corner of the cavity; (iii) conversion of the upstream influence arriving at this location to a fluctuation in the separating shear layer; and (iv) amplification of this fluctuation in the shear layer as it develops in the streamwise direction. In view of the fact that inflow shear-layer in the present investigation is fully turbulent, item (iv) is of particular interest. It is generally recognized, at least for laminar conditions at separation from the leading-corner of the cavity, that the disturbance growth in the shear layer can be described using concepts of linearized, inviscid stability theory, as shown by Rockwell, Sarohia, and Knisely and Rockwell. As demonstrated by Knisely and Rockwell, on the basis of experiments interpreted
Flow and Thermal Behavior of an EGS Reservoir - Geothermal Code...
Office of Scientific and Technical Information (OSTI)
Conference: Flow and Thermal Behavior of an EGS Reservoir - Geothermal Code Comparison Study Citation Details In-Document Search Title: Flow and Thermal Behavior of an EGS ...
Transition from thermal to turbulent equilibrium with a resulting electromagnetic spectrum
Ziebell, L. F.; Yoon, P. H.; School of Space Research, Kyung Hee University, Yongin, Gyeonggi 446-701 ; Gaelzer, R.; Instituto de FĆsica e MatemĆ”tica, UFPel, Pelotas, RS ; Pavan, J.
2014-01-15
A recent paper [Ziebell et al., Phys. Plasmas 21, 010701 (2014)] discusses a new type of radiation emission process for plasmas in a state of quasi-equilibrium between the particles and enhanced Langmuir turbulence. Such a system may be an example of the so-called āturbulent quasi-equilibrium.ā In the present paper, it is shown on the basis of electromagnetic weak turbulence theory that an initial thermal equilibrium state (i.e., only electrostatic fluctuations and Maxwellian particle distributions) transitions toward the turbulent quasi-equilibrium state with enhanced electromagnetic radiation spectrum, thus demonstrating that the turbulent quasi-equilibrium discussed in the above paper correctly describes the weakly turbulent plasma dynamically interacting with electromagnetic fluctuations, while maintaining a dynamical steady-state in the average sense.
Turbulence and bias-induced flows in simple magnetized toroidal plasmas
Li, B.; Rogers, B. N.; Ricci, P.; Gentle, K. W.; Bhattacharjee, A.
2011-05-15
Turbulence and bias-induced flows in simple magnetized toroidal plasmas are explored with global three-dimensional fluid simulations, focusing on the parameters of the Helimak experiment. The simulations show that plasma turbulence and transport in the regime of interest are dominated by the ideal interchange instability. The application of a bias voltage alters the structure of the plasma potential, resulting in the equilibrium sheared flows.These bias-induced vertical flows located in the gradient region appear to reduce the radial extent of turbulent structures,and thereby lower the radial plasma transport on the low field side.
Laminar and turbulent nozzle-jet flows and their acoustic near-field
Bühler, Stefan; Obrist, Dominik; Kleiser, Leonhard
2014-08-15
We investigate numerically the effects of nozzle-exit flow conditions on the jet-flow development and the near-field sound at a diameter-based Reynolds number of Re{sub D} = 18?100 and Mach number Ma = 0.9. Our computational setup features the inclusion of a cylindrical nozzle which allows to establish a physical nozzle-exit flow and therefore well-defined initial jet-flow conditions. Within the nozzle, the flow is modeled by a potential flow core and a laminar, transitional, or developing turbulent boundary layer. The goal is to document and to compare the effects of the different jet inflows on the jet flow development and the sound radiation. For laminar and transitional boundary layers, transition to turbulence in the jet shear layer is governed by the development of Kelvin-Helmholtz instabilities. With the turbulent nozzle boundary layer, the jet flow development is characterized by a rapid changeover to a turbulent free shear layer within about one nozzle diameter. Sound pressure levels are strongly enhanced for laminar and transitional exit conditions compared to the turbulent case. However, a frequency and frequency-wavenumber analysis of the near-field pressure indicates that the dominant sound radiation characteristics remain largely unaffected. By applying a recently developed scaling procedure, we obtain a close match of the scaled near-field sound spectra for all nozzle-exit turbulence levels and also a reasonable agreement with experimental far-field data.
Turbulent Flow Inside Pipes with Two-Dimensional Rib Roughness
Energy Science and Technology Software Center (OSTI)
1994-01-24
A commonly used internal enhancement for single-phase forced-convective turbulent flow applications is tranverse and/or near tranverse ribs. These enhanced surfaces consist of a uniform inside diameter with periodic and discrete disruption of ribs. Enhanced tubes of this type are made by an extrusion process and are used in some condensers and evaporators in refrigeration systems. Tubes of this type fall into an enhancement category called separation and reattachment that has been identified as one ofmoreĀ Ā» the most energy efficient. Lacking are prediction methods that are mechanistic based that can be used to calculate the heat-transfer coefficients and friction-factors for tubes with this enhancement type. This program calculates the Nusselt number and friction factor for enhanced tubes with tranverse, rectangular ribs with a spacing exceeding the reattachment length. The input quantities are the enhancement height, spacing, and the width. The Nusselt number and friction factor are calculated for a specific Reynolds number or for a range of Reynolds numbers. Users of the program are heat-exchanger designers, enhanced tubing suppliers, and research organizations or academia who are developing or validating prediction methods. The manufacturers of refrigeration heat exchangers and enhanced tube suppliers are potential users of this software.Ā«Ā less
Direct Simulation of Turbulent Swept Flow Over a Wire in a Channel |
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Argonne Leadership Computing Facility Turbulent Swept Flow Over a Wire in a Channel Authors: RANJAN, R., PANTANO, C., FISCHER, P. Turbulent swept flow over a cylindrical wire placed on a wall of a channel is investigated using direct numerical simulations. This geometry is a model of the flow through the wire-wrapped fuel pins, the heat exchanger, typical of many nuclear reactor designs. Mean flow along and across the wire axis is imposed, leading to the formation of separated flow regions.
ASCR Workshop on Turbulent Flow Simulations at the Exascale: Opportunities and Challenges
Broader source: Energy.gov [DOE]
The need for accurateĀ simulation of turbulent flows is evident across the US Department of Energy applied-scienceĀ andĀ engineering portfolio, including combustion, plasma physics, nuclear-reactor...
Direct Simulation of Fully Resolved Vaporizing Droplets in a Turbulent Flow
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
| Argonne Leadership Computing Facility Direct Simulation of Fully Resolved Vaporizing Droplets in a Turbulent Flow PI Name: Said Elghobashi PI Email: selghoba@uci.edu Institution: University of California, Irvine Allocation Program: INCITE Allocation Hours at ALCF: 20,000,000 Year: 2012 Research Domain: Engineering This project will enhance the understanding of liquid droplet vaporization and mixing processes in a turbulent flow. All liquid fuel combustion devices, mobile or stationary, use
Caughey, David
2010-10-08
A Symposium on Turbulence and Combustion was held at Cornell University on August 3-4, 2009. The overall goal of the Symposium was to promote future advances in the study of turbulence and combustion, through an unique forum intended to foster interactions between leading members of these two research communities. The Symposium program consisted of twelve invited lectures given by world-class experts in these fields, two poster sessions consisting of nearly 50 presentations, an open forum, and other informal activities designed to foster discussion. Topics covered in the lectures included turbulent dispersion, wall-bounded flows, mixing, finite-rate chemistry, and others, using experiment, modeling, and computations, and included perspectives from an international community of leading researchers from academia, national laboratories, and industry.
GPU accelerated flow solver for direct numerical simulation of turbulent flows
Salvadore, Francesco; Botti, Michela
2013-02-15
Graphical processing units (GPUs), characterized by significant computing performance, are nowadays very appealing for the solution of computationally demanding tasks in a wide variety of scientific applications. However, to run on GPUs, existing codes need to be ported and optimized, a procedure which is not yet standardized and may require non trivial efforts, even to high-performance computing specialists. In the present paper we accurately describe the porting to CUDA (Compute Unified Device Architecture) of a finite-difference compressible NavierāStokes solver, suitable for direct numerical simulation (DNS) of turbulent flows. Porting and validation processes are illustrated in detail, with emphasis on computational strategies and techniques that can be applied to overcome typical bottlenecks arising from the porting of common computational fluid dynamics solvers. We demonstrate that a careful optimization work is crucial to get the highest performance from GPU accelerators. The results show that the overall speedup of one NVIDIA Tesla S2070 GPU is approximately 22 compared with one AMD Opteron 2352 Barcelona chip and 11 compared with one Intel Xeon X5650 Westmere core. The potential of GPU devices in the simulation of unsteady three-dimensional turbulent flows is proved by performing a DNS of a spatially evolving compressible mixing layer.
Direct numerical simulation of turbulent flow in a rotating square duct
Dai, Yi-Jun; Huang, Wei-Xi Xu, Chun-Xiao; Cui, Gui-Xiang
2015-06-15
A fully developed turbulent flow in a rotating straight square duct is simulated by direct numerical simulations at Re{sub ?} = 300 and 0 ? Ro{sub ?} ? 40. The rotating axis is parallel to two opposite walls of the duct and normal to the main flow. Variations of the turbulence statistics with the rotation rate are presented, and a comparison with the rotating turbulent channel flow is discussed. Rich secondary flow patterns in the cross section are observed by varying the rotation rate. The appearance of a pair of additional vortices above the pressure wall is carefully examined, and the underlying mechanism is explained according to the budget analysis of the mean momentum equations.
Radial flow nuclear thermal rocket (RFNTR)
Leyse, Carl F.
1995-01-01
A radial flow nuclear thermal rocket fuel assembly includes a substantially conical fuel element having an inlet side and an outlet side. An annular channel is disposed in the element for receiving a nuclear propellant, and a second, conical, channel is disposed in the element for discharging the propellant. The first channel is located radially outward from the second channel, and separated from the second channel by an annular fuel bed volume. This fuel bed volume can include a packed bed of loose fuel beads confined by a cold porous inlet frit and a hot porous exit frit. The loose fuel beads include ZrC coated ZrC-UC beads. In this manner, nuclear propellant enters the fuel assembly axially into the first channel at the inlet side of the element, flows axially across the fuel bed volume, and is discharged from the assembly by flowing radially outward from the second channel at the outlet side of the element.
Radial flow nuclear thermal rocket (RFNTR)
Leyse, Carl F.
1995-11-07
A radial flow nuclear thermal rocket fuel assembly includes a substantially conical fuel element having an inlet side and an outlet side. An annular channel is disposed in the element for receiving a nuclear propellant, and a second, conical, channel is disposed in the element for discharging the propellant. The first channel is located radially outward from the second channel, and separated from the second channel by an annular fuel bed volume. This fuel bed volume can include a packed bed of loose fuel beads confined by a cold porous inlet frit and a hot porous exit frit. The loose fuel beads include ZrC coated ZrC-UC beads. In this manner, nuclear propellant enters the fuel assembly axially into the first channel at the inlet side of the element, flows axially across the fuel bed volume, and is discharged from the assembly by flowing radially outward from the second channel at the outlet side of the element.
Cioncolini, Andrea; Santini, Lorenzo [Department of Nuclear Engineering, Politecnico di Milano, via Ponzio 34/3, 20133 Milano (Italy)
2006-03-01
An experimental study was carried out to investigate the transition from laminar to turbulent flow in helically coiled pipes. Twelve coils have been tested, with ratios of coil diameter to tube diameter ranging from 6.9 to 369, and the interaction between turbulence emergence and coil curvature has been analyzed from direct observation of the experimental friction factor profiles. The experimental data compare favorably with existing results and reveal new features that apparently were not observed in previous research. (author)
Self-sustaining turbulence in a restricted nonlinear model of plane Couette flow
Thomas, Vaughan L.; Gayme, Dennice F.; Lieu, Binh K.; JovanoviÄ, Mihailo R.; Farrell, Brian F.; Ioannou, Petros J.
2014-10-15
This paper demonstrates the maintenance of self-sustaining turbulence in a restricted nonlinear (RNL) model of plane Couette flow. The RNL system is derived directly from the Navier-Stokes equations and permits higher resolution studies of the dynamical system associated with the stochastic structural stability theory (S3T) model, which is a second order approximation of the statistical state dynamics of the flow. The RNL model shares the dynamical restrictions of the S3T model but can be easily implemented by reducing a DNS code so that it retains only the RNL dynamics. Comparisons of turbulence arising from DNS and RNL simulations demonstrate that the RNL system supports self-sustaining turbulence with a mean flow as well as structural and dynamical features that are consistent with DNS. These results demonstrate that the simplified RNL system captures fundamental aspects of fully developed turbulence in wall-bounded shear flows and motivate use of the RNL/S3T framework for further study of wall-turbulence.
Hydrokinetic canal measurements: inflow velocity, wake flow velocity, and turbulence
Gunawan, Budi
2014-06-11
The dataset consist of acoustic Doppler current profiler (ADCP) velocity measurements in the wake of a 3-meter diameter vertical-axis hydrokinetic turbine deployed in Roza Canal, Yakima, WA, USA. A normalized hub-centerline wake velocity profile and two cross-section velocity contours, 10 meters and 20 meters downstream of the turbine, are presented. Mean velocities and turbulence data, measured using acoustic Doppler velocimeter (ADV) at 50 meters upstream of the turbine, are also presented. Canal dimensions and hydraulic properties, and turbine-related information are also included.
Turbulent Flow Effects on the Biological Performance of Hydro-Turbines
Richmond, Marshall C.; Romero Gomez, Pedro DJ
2014-08-25
The hydro-turbine industry uses Computational Fluid Dynamics (CFD) tools to predict the flow conditions as part of the design process for new and rehabilitated turbine units. Typically the hydraulic design process uses steady-state simulations based on Reynolds-Averaged Navier-Stokes (RANS) formulations for turbulence modeling because these methods are computationally efficient and work well to predict averaged hydraulic performance, e.g. power output, efficiency, etc. However, in view of the increasing emphasis on environmental concerns, such as fish passage, the consideration of the biological performance of hydro-turbines is also required in addition to hydraulic performance. This leads to the need to assess whether more realistic simulations of the turbine hydraulic environment -those that resolve unsteady turbulent eddies not captured in steady-state RANS computations- are needed to better predict the occurrence and extent of extreme flow conditions that could be important in the evaluation of fish injury and mortality risks. In the present work, we conduct unsteady, eddy-resolving CFD simulations on a Kaplan hydro-turbine at a normal operational discharge. The goal is to quantify the impact of turbulence conditions on both the hydraulic and biological performance of the unit. In order to achieve a high resolution of the incoming turbulent flow, Detached Eddy Simulation (DES) turbulence model is used. These transient simulations are compared to RANS simulations to evaluate whether extreme hydraulic conditions are better captured with advanced eddy-resolving turbulence modeling techniques. The transient simulations of key quantities such as pressure and hydraulic shear flow that arise near the various components (e.g. wicket gates, stay vanes, runner blades) are then further analyzed to evaluate their impact on the statistics for the lowest absolute pressure (nadir pressures) and for the frequency of collisions that are known to cause mortal injury in fish
Thermal shallow water models of geostrophic turbulence in Jovian atmospheres
Warneford, Emma S. Dellar, Paul J.
2014-01-15
Conventional shallow water theory successfully reproduces many key features of the Jovian atmosphere: a mixture of coherent vortices and stable, large-scale, zonal jets whose amplitude decreases with distance from the equator. However, both freely decaying and forced-dissipative simulations of the shallow water equations in Jovian parameter regimes invariably yield retrograde equatorial jets, while Jupiter itself has a strong prograde equatorial jet. Simulations by Scott and Polvani [Equatorial superrotation in shallow atmospheres, Geophys. Res. Lett. 35, L24202 (2008)] have produced prograde equatorial jets through the addition of a model for radiative relaxation in the shallow water height equation. However, their model does not conserve mass or momentum in the active layer, and produces mid-latitude jets much weaker than the equatorial jet. We present the thermal shallow water equations as an alternative model for Jovian atmospheres. These equations permit horizontal variations in the thermodynamic properties of the fluid within the active layer. We incorporate a radiative relaxation term in the separate temperature equation, leaving the mass and momentum conservation equations untouched. Simulations of this model in the Jovian regime yield a strong prograde equatorial jet, and larger amplitude mid-latitude jets than the Scott and Polvani model. For both models, the slope of the non-zonal energy spectra is consistent with the classic Kolmogorov scaling, and the slope of the zonal energy spectra is consistent with the much steeper spectrum observed for Jupiter. We also perform simulations of the thermal shallow water equations for Neptunian parameter values, with a radiative relaxation time scale calculated for the same 25 mbar pressure level we used for Jupiter. These Neptunian simulations reproduce the broad, retrograde equatorial jet and prograde mid-latitude jets seen in observations. The much longer radiative time scale for the colder planet Neptune
Flow and Thermal Behavior of an EGS Reservoir - Geothermal Code...
Office of Scientific and Technical Information (OSTI)
of an EGS Reservoir - Geothermal Code Comparison Study Citation Details In-Document Search Title: Flow and Thermal Behavior of an EGS Reservoir - Geothermal Code Comparison Study ...
Radial Spreading of Drift-Wave-Zonal-Flow Turbulence via Soliton Formation
Guo Zehua; Chen Liu; Zonca, Fulvio
2009-07-31
The self-consistent spatiotemporal evolution of a drift-wave (DW) radial envelope and a zonal-flow (ZF) amplitude is investigated in a slab model. The stationary solution of the coupled partial differential equations in a simple limit yields the formation of DW-ZF soliton structures, which propagate radially with speed depending on the envelope peak amplitude. Additional interesting physics, e.g., the generation, destruction, collision, and reflection of solitons, as well as turbulence bursting can also be observed due to the effects of linear growth or damping, dissipation, equilibrium nonuniformities and soliton dynamics. The propagation of soliton causes significant radial spreading of DW turbulence and therefore can affect transport scaling with the system size by broadening of the turbulent region. The correspondence of the present analysis with the description of DW-ZF interactions in toroidal geometry is also discussed.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Schilling, Oleg; Mueschke, Nicholas J.
2010-10-18
Data from a 1152X760X1280 direct numerical simulation (DNS) of a transitional Rayleigh-Taylor mixing layer modeled after a small Atwood number water channel experiment is used to comprehensively investigate the structure of mean and turbulent transport and mixing. The simulation had physical parameters and initial conditions approximating those in the experiment. The budgets of the mean vertical momentum, heavy-fluid mass fraction, turbulent kinetic energy, turbulent kinetic energy dissipation rate, heavy-fluid mass fraction variance, and heavy-fluid mass fraction variance dissipation rate equations are constructed using Reynolds averaging applied to the DNS data. The relative importance of mean and turbulent production, turbulent dissipationmore »and destruction, and turbulent transport are investigated as a function of Reynolds number and across the mixing layer to provide insight into the flow dynamics not presently available from experiments. The analysis of the budgets supports the assumption for small Atwood number, Rayleigh/Taylor driven flows that the principal transport mechanisms are buoyancy production, turbulent production, turbulent dissipation, and turbulent diffusion (shear and mean field production are negligible). As the Reynolds number increases, the turbulent production in the turbulent kinetic energy dissipation rate equation becomes the dominant production term, while the buoyancy production plateaus. Distinctions between momentum and scalar transport are also noted, where the turbulent kinetic energy and its dissipation rate both grow in time and are peaked near the center plane of the mixing layer, while the heavy-fluid mass fraction variance and its dissipation rate initially grow and then begin to decrease as mixing progresses and reduces density fluctuations. All terms in the transport equations generally grow or decay, with no qualitative change in their profile, except for the pressure flux contribution to the total turbulent kinetic
Isolation of Metals from Liquid Wastes: Reactive Scavenging in Turbulent Thermal Reactors
Jost O.L. Wendt; Alan R. Kerstein; Alexander Scheeline; Arne Pearlstein; William Linak
2003-08-06
The Overall project demonstrated that toxic metals (cesium Cs and strontium Sr) in aqueous and organic wastes can be isolated from the environment through reaction with kaolinite based sorbent substrates in high temperature reactor environments. In addition, a state-of-the art laser diagnostic tool to measure droplet characteristic in practical 'dirty' laboratory environments was developed, and was featured on the cover of a recent edition of the scientific journal ''applied Spectroscopy''. Furthermore, great strides have been made in developing a theoretical model that has the potential to allow prediction of the position and life history of every particle of waste in a high temperature, turbulent flow field, a very challenging problem involving as it does, the fundamentals of two phase turbulence and of particle drag physics.
Flow, Mixing and Combustion of Transient Turbulent Gaseous Jets in Confined
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Cylindrical Geometries | Argonne Leadership Computing Facility Flow, Mixing and Combustion of Transient Turbulent Gaseous Jets in Confined Cylindrical Geometries PI Name: Christos Frouzakis PI Email: frouzakis@lav.mavt.ethz.ch Institution: Swiss Federal Institute of Technology Zurich (ETHZ) Allocation Program: ESP Year: 2015 Research Domain: Engineering Tier 2 Code Development Project Numerical Methods/Algorithms Direct numerical simulations for this project will be based on the open source
The effect of diamagnetic flows on turbulent driven ion toroidal rotation
Lee, J. P.; Barnes, M.; Parra, F. I.; Belli, E. A.; Candy, J.
2014-05-15
Turbulent momentum redistribution determines the radial profile of rotation in a tokamak. The momentum transport driven by diamagnetic flow effects is an important piece of the radial momentum transport for sub-sonic rotation, which is often observed in experiments. In a non-rotating state, the diamagnetic flow and the E × B flow must cancel. The diamagnetic flow and the E × B flow have different effects on the turbulent momentum flux, and this difference in behavior induces intrinsic rotation. The momentum flux is evaluated using gyrokinetic equations that are corrected to higher order in the ratio of the poloidal Larmor radius to the minor radius, which requires evaluation of the diamagnetic corrections to Maxwellian equilibria. To study the momentum transport due to diamagnetic flow effects, three experimental observations of ion rotation are examined. First, a strong pressure gradient at the plasma edge is shown to result in a significant inward momentum transport due to the diamagnetic effect, which may explain the observed peaking of rotation in a high confinement mode. Second, the direction of momentum transport is shown to change as collisionality increases, which is qualitatively consistent with the observed reversal of intrinsic rotation by varying plasma density and current. Last, the dependence of the intrinsic momentum flux on the magnetic shear is found, and it may explain the observed rotation changes in the presence of lower hybrid current drive.
Simulations of Turbulent Flows with Strong Shocks and Density Variations: Final Report
Sanjiva Lele
2012-10-01
The target of this SciDAC Science Application was to develop a new capability based on high-order and high-resolution schemes to simulate shock-turbulence interactions and multi-material mixing in planar and spherical geometries, and to study Rayleigh-Taylor and Richtmyer-Meshkov turbulent mixing. These fundamental problems have direct application in high-speed engineering flows, such as inertial confinement fusion (ICF) capsule implosions and scramjet combustion, and also in the natural occurrence of supernovae explosions. Another component of this project was the development of subgrid-scale (SGS) models for large-eddy simulations of flows involving shock-turbulence interaction and multi-material mixing, that were to be validated with the DNS databases generated during the program. The numerical codes developed are designed for massively-parallel computer architectures, ensuring good scaling performance. Their algorithms were validated by means of a sequence of benchmark problems. The original multi-stage plan for this five-year project included the following milestones: 1) refinement of numerical algorithms for application to the shock-turbulence interaction problem and multi-material mixing (years 1-2); 2) direct numerical simulations (DNS) of canonical shock-turbulence interaction (years 2-3), targeted at improving our understanding of the physics behind the combined two phenomena and also at guiding the development of SGS models; 3) large-eddy simulations (LES) of shock-turbulence interaction (years 3-5), improving SGS models based on the DNS obtained in the previous phase; 4) DNS of planar/spherical RM multi-material mixing (years 3-5), also with the two-fold objective of gaining insight into the relevant physics of this instability and aiding in devising new modeling strategies for multi-material mixing; 5) LES of planar/spherical RM mixing (years 4-5), integrating the improved SGS and multi-material models developed in stages 3 and 5. This final report is
MODELING STRATEGIES FOR UNSTEADY TURBULENT FLOWS IN THE LOWER PLENUM OF THE VHTR
Richard W. Johnson
2006-09-01
Validation simulations are presented for turbulent flow in a staggered tube bank, geometry similar to that in the lower plenum of a block very high temperature reactor. Steady 2D RANS predictions are compared to unsteady 2D RANS results and experiment. The unsteady calculations account for the fact that nonturbulent fluctuations (due to vortex-shedding) are present in the flow. The unsteady computations are shown to predict the mean variables and the total shear stress quite well. Previous workers have presented results that indicated that 3D simulations were necessary to obtain reasonable results. Best practices are based on requirements for the ASME Journal of Fluids Engineering.
Internal (Annular) and Compressible External (Flat Plate) Turbulent Flow Heat Transfer Correlations.
Dechant, Lawrence; Smith, Justin
2016-01-01
Here we provide a discussion regarding the applicability of a family of traditional heat transfer correlation based models for several (unit level) heat transfer problems associated with flight heat transfer estimates and internal flow heat transfer associated with an experimental simulation design (Dobranich 2014). Variability between semi-empirical free-flight models suggests relative differences for heat transfer coefficients on the order of 10%, while the internal annular flow behavior is larger with differences on the order of 20%. We emphasize that these expressions are strictly valid only for the geometries they have been derived for e.g. the fully developed annular flow or simple external flow problems. Though, the application of flat plate skin friction estimate to cylindrical bodies is a traditional procedure to estimate skin friction and heat transfer, an over-prediction bias is often observed using these approximations for missile type bodies. As a correction for this over-estimate trend, we discuss a simple scaling reduction factor for flat plate turbulent skin friction and heat transfer solutions (correlations) applied to blunt bodies of revolution at zero angle of attack. The method estimates the ratio between axisymmetric and 2-d stagnation point heat transfer skin friction and Stanton number solution expressions for sub-turbulent Reynolds numbers %3C1x10 4 . This factor is assumed to also directly influence the flat plate results applied to the cylindrical portion of the flow and the flat plate correlations are modified by
Flow field computation of the NREL S809 airfoil using various turbulence models
Chang, Y.L.; Yang, S.L.; Arici, O. [Michigan Technological Univ., Houghton, MI (United States). Mechanical Engineering-Engineering Mechanics Dept.
1996-10-01
Performance comparison of three popular turbulence models, namely Baldwin-Lomas algebraic model, Chien`s Low-Reynolds-Number {kappa}-{epsilon} model, and Wilcox`s Low-Reynolds-Number {kappa}-{omega} model, is given. These models were applied to calculate the flow field around the National Renewable Energy Laboratory S809 airfoil using Total Variational Diminishing scheme. Numerical results of C{sub P}, C{sub L}, and C{sub D} are presented along with the Delft experimental data. It is shown that all three models perform well for attached flow, i.e., no flow separation at low angles of attack. However, at high angles of attack with flow separation, convergence characteristics show Wilcox`s model outperforms the other models. Results of this study will be used to guide the authors in their dynamic stall research.
Theory of Fine-scale Zonal Flow Generation From Trapped Electron Mode Turbulence
Lu Wang and T.S. Hahm
2009-06-11
Most existing zonal flow generation theory has been developed with a usual assumption of qr??” << 1 (qr is the radial wave number of zonal flow, and ??” is the ion poloidal gyrora- dius). However, recent nonlinear gyrokinetic simulations of trapped electron mode (TEM) turbulence exhibit a relatively short radial scale of the zonal flows with qr??” ~ 1 [Z. Lin et al., IAEA-CN/TH/P2-8 (2006); D. Ernst et al., Phys. Plasmas 16, 055906 (2009)]. This work reports an extension of zonal flow growth calculation to this short wavelength regime via the wave kinetics approach. A generalized expression for the polarization shielding for arbitrary radial wavelength [Lu Wang and T.S. Hahm, to appear in Phys. Plasmas (2009)] which extends the Rosenbluth-Hinton formula in the long wavelength limit is applied.
Response Relationship Between Juvenile Salmon and an Autonomous Sensor in Turbulent Flows
Richmond, Marshall C.; Deng, Zhiqun; McKinstry, Craig A.; Mueller, Robert P.; Carlson, Thomas J.; Dauble, Dennis D.
2009-01-27
Juvenile fall chinook salmon (Oncorhynchus tshawythscha) and an autonomous sensor device (Sensor Fish) were exposed to turbulent shear flows in order to determine how hydraulic conditions effected fish injury response. Studies were designed to establish correlation metrics between Sensor Fish device measurements and live fish injuries by conducting concurrent releases in a range of turbulent shear flows. Comparisons were made for two exposure scenarios. In the fast-fish-to-slow-water scenario, test fish were carried by the fast-moving water of a submerged turbulent jet and exposed into the standing water of a flume. In the slow-fish-to-fast-water scenario, test fish were introduced into a turbulent jet from standing water through an introduction tube placed just outside the edge of the jet. Motion-tracking analysis was performed on high-speed, high-resolution digital videos of all the releases at water jet velocities ranging from 3 to 22.9 m Ā· s^{-1}. Velocities of the Sensor Fish were very similar to those of live fish, but maximum accelerations of live fish were larger than those by Sensor Fish for all the nozzle velocities of both cenarios. A 10% probability of major injury threshold was found to occur at sensor fish accelerations of 513 and 260 (m Ā· s^{-2}) for the fast-fish-to-slow-water and slow-fish-to-fast-water scenarios, respectively. The findings provide a linkage between laboratory experiments of fish injury, field survival studies, and numerical modeling.
Bodo, G.; Rossi, P.; Cattaneo, F.; Mignone, A.
2012-12-20
We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero mean magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. The calculations begin from an isothermal state spanning three scale heights above and below the mid-plane. After a long transient the layers settle to a stationary state in which thermal losses out of the boundaries are balanced by dissipative heating. We identify two regimes. The first is a conductive regime in which the heat is transported mostly by conduction and the density decreases with height. In the limit of large thermal diffusivity this regime resembles the more familiar isothermal case. The second is the convective regime, observed at smaller values of the thermal diffusivity, in which the layer becomes unstable to overturning motions, the heat is carried mostly by advection, and the density becomes nearly constant throughout the layer. In this latter constant-density regime we observe evidence for large-scale dynamo action leading to a substantial increase in transport efficiency relative to the conductive case.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Wu, Xiaohua; Moin, Parviz; Adrian, Ronald J.; Baltzer, Jon R.
2015-06-15
We report that the precise dynamics of breakdown in pipe transition is a century-old unresolved problem in fluid mechanics. We demonstrate that the abruptness and mysteriousness attributed to the Osborne Reynolds pipe transition can be partially resolved with a spatially developing direct simulation that carries weakly but finitely perturbed laminar inflow through gradual rather than abrupt transition arriving at the fully developed turbulent state. Our results with this approach show during transition the energy norms of such inlet perturbations grow exponentially rather than algebraically with axial distance. When inlet disturbance is located in the core region, helical vortex filaments evolvemoreĀ Ā» into large-scale reverse hairpin vortices. The interaction of these reverse hairpins among themselves or with the near-wall flow when they descend to the surface from the core produces small-scale hairpin packets, which leads to breakdown. When inlet disturbance is near the wall, certain quasi-spanwise structure is stretched into a Lambda vortex, and develops into a large-scale hairpin vortex. Small-scale hairpin packets emerge near the tip region of the large-scale hairpin vortex, and subsequently grow into a turbulent spot, which is itself a local concentration of small-scale hairpin vortices. This vortex dynamics is broadly analogous to that in the boundary layer bypass transition and in the secondary instability and breakdown stage of natural transition, suggesting the possibility of a partial unification. Under parabolic base flow the friction factor overshoots Moodyās correlation. Plug base flow requires stronger inlet disturbance for transition. Finally, accuracy of the results is demonstrated by comparing with analytical solutions before breakdown, and with fully developed turbulence measurements after the completion of transition.Ā«Ā less
Experimental investigation of 3-D turbulent free shear flow past propellers and windmills
Kotb, M.A.
1984-01-01
An experimental investigation of the flowfield region near a 0.49 m (1.615 ft), 3 bladed, horizontal axis rotor was conducted in the VPI 2 m x 2 m (6 ft x 6 ft) wind tunnel. Two different modes of operations were studied - propeller and windmill. For each case, tests were run with a uniform approach flow as a base line condition, and the main test series was run with a variable mesh wire grid upstream to produce an approach flow with an essentially linear velocity gradient. The results are compared to elucidate the effects of the non-uniform approach flow. Several types of measurements are reported. First are gross quantities such as overall thrust and power. The second type of measurements are mean (in the turbulence sense) quantities obtained with a five port yawhead tube. All three components of mean velocity and static pressure were obtained. The third type of measurements were made with an x-wire anemometer and an r.m.s. meter. These measurements yield all components of the turbulence intensities and stresses at a point averaged over many passes of the rotor blades.
Uribe, A. L. [University of Chicago, Chicago, IL 60637 (United States); Klahr, H.; Henning, Th., E-mail: uribe@oddjob.uchicago.edu [Max-Planck-Institut fuer Astronomie, Heidelberg (Germany)
2013-06-01
We have performed three-dimensional magnetohydrodynamical simulations of stellar accretion disks, using the PLUTO code, and studied the accretion of gas onto a Jupiter-mass planet and the structure of the circumplanetary gas flow after opening a gap in the disk. We compare our results with simulations of laminar, yet viscous disks with different levels of an {alpha}-type viscosity. In all cases, we find that the accretion flow across the surface of the Hill sphere of the planet is not spherically or azimuthally symmetric, and is predominantly restricted to the mid-plane region of the disk. Even in the turbulent case, we find no significant vertical flow of mass into the Hill sphere. The outer parts of the circumplanetary disk are shown to rotate significantly below Keplerian speed, independent of viscosity, while the circumplanetary disk density (therefore the angular momentum) increases with viscosity. For a simulation of a magnetized turbulent disk, where the global averaged alpha stress is {alpha}{sub MHD} = 10{sup -3}, we find the accretion rate onto the planet to be M-dot {approx}2 Multiplication-Sign 10{sup -6}M{sub J} yr{sup -1} for a gap surface density of 12 g cm{sup -2}. This is about a third of the accretion rate obtained in a laminar viscous simulation with equivalent {alpha} parameter.
WENO schemes on arbitrary unstructured meshes for laminar, transitional and turbulent flows
Tsoutsanis, Panagiotis, E-mail: panagiotis.tsoutsanis@cranfield.ac.uk; Antoniadis, Antonios Foivos, E-mail: a.f.antoniadis@cranfield.ac.uk; Drikakis, Dimitris, E-mail: d.drikakis@cranfield.ac.uk
2014-01-01
This paper presents the development and implementation of weighted-essentially-non-oscillatory (WENO) schemes for viscous flows on arbitrary unstructured grids. WENO schemes up to fifth-order accurate have been implemented in conjunction with hybrid and non-hybrid unstructured grids. The schemes are investigated with reference to numerical and experimental results for the TaylorGreen vortex, as well as for laminar and turbulent flows around a sphere, and the turbulent shock-wave boundary layer interaction flow problem. The results show that the accuracy of the schemes depends on the arbitrariness of shape and orientation of the unstructured mesh elements, as well as the compactness of directional stencils. The WENO schemes provide a more accurate numerical framework compared to second-order and third-order total variation diminishing (TVD) methods, however, the fifth-order version of the schemes is computationally too expensive to make the schemes practically usable. On the other hand, the third-order variant offers an excellent numerical framework in terms of accuracy and computational cost compared to the fifth-order WENO and second-order TVD schemes. Parallelisation of the CFD code (henceforth labelled as UCNS3D), where the schemes have been implemented, shows that the present methods offer very good scalable performance.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Schilling, Oleg; Mueschke, Nicholas J.
2010-10-18
Data from a 1152X760X1280 direct numerical simulation (DNS) of a transitional Rayleigh-Taylor mixing layer modeled after a small Atwood number water channel experiment is used to comprehensively investigate the structure of mean and turbulent transport and mixing. The simulation had physical parameters and initial conditions approximating those in the experiment. The budgets of the mean vertical momentum, heavy-fluid mass fraction, turbulent kinetic energy, turbulent kinetic energy dissipation rate, heavy-fluid mass fraction variance, and heavy-fluid mass fraction variance dissipation rate equations are constructed using Reynolds averaging applied to the DNS data. The relative importance of mean and turbulent production, turbulent dissipationmoreĀ Ā» and destruction, and turbulent transport are investigated as a function of Reynolds number and across the mixing layer to provide insight into the flow dynamics not presently available from experiments. The analysis of the budgets supports the assumption for small Atwood number, Rayleigh/Taylor driven flows that the principal transport mechanisms are buoyancy production, turbulent production, turbulent dissipation, and turbulent diffusion (shear and mean field production are negligible). As the Reynolds number increases, the turbulent production in the turbulent kinetic energy dissipation rate equation becomes the dominant production term, while the buoyancy production plateaus. Distinctions between momentum and scalar transport are also noted, where the turbulent kinetic energy and its dissipation rate both grow in time and are peaked near the center plane of the mixing layer, while the heavy-fluid mass fraction variance and its dissipation rate initially grow and then begin to decrease as mixing progresses and reduces density fluctuations. All terms in the transport equations generally grow or decay, with no qualitative change in their profile, except for the pressure flux contribution to the total turbulent kinetic
Dye laser amplifier including a low turbulence, stagnation-free dye flow configuration
Davin, James
1992-01-01
A large (high flow rate) dye laser amplifier in which a continuous replenished supply of dye is excited by a first light beam, specifically a copper vapor laser beam, in order to amplify the intensity of a second different light beam, specifically a dye beam, passing through the dye is disclosed herein. This amplifier includes a dye cell defining a dye chamber through which a continuous stream of dye is caused to pass at a flow rate of for example 30 gallons/minute, a specifically designed support vessel for containing the dye cell and a screen device for insuring that the dye stream passes into the dye cell in a substantially turbulent free, stagnation-free manner.
Dye laser amplifier including a low turbulence, stagnation-free dye flow configuration
Davin, J.
1992-12-01
A large (high flow rate) dye laser amplifier in which a continuous replenished supply of dye is excited by a first light beam, specifically a copper vapor laser beam, in order to amplify the intensity of a second different light beam, specifically a dye beam, passing through the dye is disclosed herein. This amplifier includes a dye cell defining a dye chamber through which a continuous stream of dye is caused to pass at a flow rate of for example 30 gallons/minute, a specifically designed support vessel for containing the dye cell and a screen device for insuring that the dye stream passes into the dye cell in a substantially turbulent free, stagnation-free manner. 9 figs.
Sippola, Mark R.; Nazaroff, William W.
2002-06-01
This report reviews published experimental and theoretical investigations of particle deposition from turbulent flows and considers the applicability of this body of work to the specific case of particle deposition from flows in the ducts of heating, ventilating and air conditioning (HVAC) systems. Particle deposition can detrimentally affect the performance of HVAC systems and it influences the exposure of building occupants to a variety of air pollutants. The first section of this report describes the types of HVAC systems under consideration and discusses the components, materials and operating parameters commonly found in these systems. The second section reviews published experimental investigations of particle deposition rates from turbulent flows and considers the ramifications of the experimental evidence with respect to HVAC ducts. The third section considers the structure of turbulent airflows in ventilation ducts with a particular emphasis on turbulence investigations that have been used as a basis for particle deposition models. The final section reviews published literature on predicting particle deposition rates from turbulent flows.
Evaluating the Efficacy of Wavelet Configurations on Turbulent-Flow Data
Li, Shaomeng; Gruchalla, Kenny; Potter, Kristin; Clyne, John; Childs, Hank
2015-10-25
I/O is increasingly becoming a significant constraint for simulation codes and visualization tools on modern supercomputers. Data compression is an attractive workaround, and, in particular, wavelets provide a promising solution. However, wavelets can be applied in multiple configurations, and the variations in configuration impact accuracy, storage cost, and execution time. While the variation in these factors over wavelet configurations have been explored in image processing, they are not well understood for visualization and analysis of scientific data. To illuminate this issue, we evaluate multiple wavelet configurations on turbulent-flow data. Our approach is to repeat established analysis routines on uncompressed and lossy-compressed versions of a data set, and then quantitatively compare their outcomes. Our findings show that accuracy varies greatly based on wavelet configuration, while storage cost and execution time vary less. Overall, our study provides new insights for simulation analysts and visualization experts, who need to make tradeoffs between accuracy, storage cost, and execution time.
Experiments measuring particle deposition from fully developed turbulent flow in ventilation ducts
Sippola, Mark R.; Nazaroff, William W.
2003-08-01
Particle deposition in ventilation ducts influences particle exposures of building occupants and may lead to a variety of indoor air quality concerns. Experiments have been performed in a laboratory to study the effects of particle size and air speed on deposition rates of particles from turbulent air flows in galvanized steel and internally insulated ducts with hydraulic diameters of 15.2 cm. The duct systems were constructed of materials typically found in commercial heating, ventilating and air conditioning (HVAC) systems. In the steel duct system, experiments with nominal particle sizes of 1, 3, 5, 9 and 16 {micro}m were conducted at each of three nominal air speeds: 2.2, 5.3 and 9.0 m/s. In the insulated duct system, deposition rates of particles with nominal sizes of 1, 3, 5, 8 and 13 {micro}m were measured at nominal air speeds of 2.2, 5.3 and 8.8 m/s. Fluorescent techniques were used to directly measure the deposition velocities of monodisperse fluorescent particles to duct surfaces (floor, wall and ceiling) at two straight duct sections where the turbulent flow profile was fully developed. In steel ducts, deposition rates were higher to the duct floor than to the wall, which were, in turn, greater than to the ceiling. In insulated ducts, deposition was nearly the same to the duct floor, wall and ceiling for a given particle size and air speed. Deposition to duct walls and ceilings was greatly enhanced in insulated ducts compared to steel ducts. Deposition velocities to each of the three duct surface orientations in both systems were found to increase with increasing particle size or air velocity over the ranges studied. Deposition rates measured in the current experiments were in general agreement with the limited observations of similar systems by previous researchers.
Optimization of a Two-Fluid Hydrodynamic Model of Churn-Turbulent Flow
Donna Post Guillen
2009-07-01
A hydrodynamic model of two-phase, churn-turbulent flows is being developed using the computational multiphase fluid dynamics (CMFD) code, NPHASE-CMFD. The numerical solutions obtained by this model are compared with experimental data obtained at the TOPFLOW facility of the Institute of Safety Research at the Forschungszentrum Dresden-Rossendorf. The TOPFLOW data is a high quality experimental database of upward, co-current air-water flows in a vertical pipe suitable for validation of computational fluid dynamics (CFD) codes. A five-field CMFD model was developed for the continuous liquid phase and four bubble size groups using mechanistic closure models for the ensemble-averaged Navier-Stokes equations. Mechanistic models for the drag and non-drag interfacial forces are implemented to include the governing physics to describe the hydrodynamic forces controlling the gas distribution. The closure models provide the functional form of the interfacial forces, with user defined coefficients to adjust the force magnitude. An optimization strategy was devised for these coefficients using commercial design optimization software. This paper demonstrates an approach to optimizing CMFD model parameters using a design optimization approach. Computed radial void fraction profiles predicted by the NPHASE-CMFD code are compared to experimental data for four bubble size groups.
Quenching and anisotropy of hydromagnetic turbulent transport
Karak, Bidya Binay; Brandenburg, Axel; Rheinhardt, Matthias; KĆ¤pylĆ¤, Petri J.; KĆ¤pylĆ¤, Maarit J.
2014-11-01
Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the Ī± effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Roberts flow, isotropically forced non-helical turbulence, and rotating thermal convection. We see significant quenching only when the mean magnetic field is larger than the equipartition value of the turbulence. Expressing the magnetic field in terms of the equipartition value of the quenched flows, we obtain for the quenching exponents of the turbulent magnetic diffusivity about 1.3, 1.1, and 1.3 for Roberts flow, forced turbulence, and convection, respectively. However, when the magnetic field is expressed in terms of the equipartition value of the unquenched flows, these quenching exponents become about 4, 1.5, and 2.3, respectively. For the Ī± effect, the exponent is about 1.3 for the Roberts flow and 2 for convection in the first case, but 4 and 3, respectively, in the second. In convection, the quenching of turbulent pumping follows the same power law as turbulent diffusion, while for the coefficient describing the Ī©ĆJ effect nearly the same quenching exponent is obtained as for Ī±. For forced turbulence, turbulent diffusion proportional to the second derivative along the mean magnetic field is quenched much less, especially for larger values of the magnetic Reynolds number. However, we find that in corresponding axisymmetric mean-field dynamos with dominant toroidal field the quenched diffusion coefficients are the same for the poloidal and toroidal field constituents.
Structure of hydrogen-rich transverse jets in a vitiated turbulent flow
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Lyra, Sgouria; Wilde, Benjamin; Kolla, Hemanth; Seitzman, Jerry M.; Lieuwen, Timothy C.; Chen, Jacqueline H.
2014-11-24
Our paper reports the results of a joint experimental and numerical study of the flow characteristics and flame structure of a hydrogen rich jet injected normal to a turbulent, vitiated crossflow of lean methane combustion products. Simultaneous high-speed stereoscopic PIV and OH PLIF measurements were obtained and analyzed alongside three-dimensional direct numerical simulations of inert and reacting JICF with detailed H2/COH2/CO chemistry. Both the experiment and the simulation reveal that, contrary to most previous studies of reacting JICF stabilized in low-to-moderate temperature air crossflow, the present conditions lead to a burner-attached flame that initiates uniformly around the burner edge. SignificantmoreĀ Ā» asymmetry is observed, however, between the reaction zones located on the windward and leeward sides of the jet, due to the substantially different scalar dissipation rates. The windward reaction zone is much thinner in the near field, while also exhibiting significantly higher local and global heat release than the much broader reaction zone found on the leeward side of the jet. The unsteady dynamics of the windward shear layer, which largely control the important jet/crossflow mixing processes in that region, are explored in order to elucidate the important flow stability implications arising in the inert and reacting JICF. The paper concludes with an analysis of the ignition, flame characteristics, and global structure of the burner-attached flame. FurthermoreChemical explosive mode analysis (CEMA) shows that the entire windward shear layer, and a large region on the leeward side of the jet, are highly explosive prior to ignition and are dominated by non-premixed flame structures after ignition. The predominantly mixing limited nature of the flow after ignition is examined by computing the Takeno flame index, which shows that ~70% of the heat release occurs in non-premixed regions.Ā«Ā less
Structure of hydrogen-rich transverse jets in a vitiated turbulent flow
Lyra, Sgouria; Wilde, Benjamin; Kolla, Hemanth; Seitzman, Jerry M.; Lieuwen, Timothy C.; Chen, Jacqueline H.
2014-11-24
Our paper reports the results of a joint experimental and numerical study of the flow characteristics and flame structure of a hydrogen rich jet injected normal to a turbulent, vitiated crossflow of lean methane combustion products. Simultaneous high-speed stereoscopic PIV and OH PLIF measurements were obtained and analyzed alongside three-dimensional direct numerical simulations of inert and reacting JICF with detailed H2/COH2/CO chemistry. Both the experiment and the simulation reveal that, contrary to most previous studies of reacting JICF stabilized in low-to-moderate temperature air crossflow, the present conditions lead to a burner-attached flame that initiates uniformly around the burner edge. Significant asymmetry is observed, however, between the reaction zones located on the windward and leeward sides of the jet, due to the substantially different scalar dissipation rates. The windward reaction zone is much thinner in the near field, while also exhibiting significantly higher local and global heat release than the much broader reaction zone found on the leeward side of the jet. The unsteady dynamics of the windward shear layer, which largely control the important jet/crossflow mixing processes in that region, are explored in order to elucidate the important flow stability implications arising in the inert and reacting JICF. The paper concludes with an analysis of the ignition, flame characteristics, and global structure of the burner-attached flame. FurthermoreChemical explosive mode analysis (CEMA) shows that the entire windward shear layer, and a large region on the leeward side of the jet, are highly explosive prior to ignition and are dominated by non-premixed flame structures after ignition. The predominantly mixing limited nature of the flow after ignition is examined by computing the Takeno flame index, which shows that ~70% of the heat release occurs in non-premixed regions.
Schuster, Eugenio
2014-05-02
The strong coupling between the different physical variables involved in the plasma transport phenomenon and the high complexity of its dynamics call for a model-based, multivariable approach to profile control where those predictive models could be exploited. The overall objective of this project has been to extend the existing body of work by investigating numerically and experimentally active control of unstable fluctuations, including fully developed turbulence and the associated cross-field particle transport, via manipulation of flow profiles in a magnetized laboratory plasma device. Fluctuations and particle transport can be monitored by an array of electrostatic probes, and Ex#2;B flow profiles can be controlled via a set of biased concentric ring electrodes that terminate the plasma column. The goals of the proposed research have been threefold: i- to develop a predictive code to simulate plasma transport in the linear HELCAT (HELicon-CAThode) plasma device at the University of New Mexico (UNM), where the experimental component of the proposed research has been carried out; ii- to establish the feasibility of using advanced model-based control algorithms to control cross-field turbulence-driven particle transport through appropriate manipulation of radial plasma flow profiles, iii- to investigate the fundamental nonlinear dynamics of turbulence and transport physics. Lehigh University (LU), including Prof. Eugenio Schuster and one full-time graduate student, has been primarily responsible for control-oriented modeling and model-based control design. Undergraduate students have also participated in this project through the National Science Foundation Research Experience for Undergraduate (REU) program. The main goal of the LU Plasma Control Group has been to study the feasibility of controlling turbulence-driven transport by shaping the radial poloidal flow profile (i.e., by controlling flow shear) via biased concentric ring electrodes.
Lee, S. N.; Tak, N. I.; Kim, M. H.; Noh, J. M.
2012-07-01
The existing study of Spall et al. shows that only {nu}{sup 2}-f turbulence model well matches with the experimental data of Shehata and McEligot which were obtained under strongly heated gas flows. Significant over-predictions in those literatures were observed in the convective heat transfer with the other famous turbulence models such as the k-{epsilon} and k-{omega} models. In spite of such good evidence about the performance of the{nu}{sup 2}-f model, the application of the {nu}{sup 2}-f model to the thermo-fluid analysis of a prismatic core is very rare. In this paper, therefore, the convective heat transfer of the coolant flow in a prismatic core has been investigated using the {nu}{sup 2}-f model. Computational fluid dynamics (CFD) calculations have been carried out for the typical unit cell geometry of a prismatic fuel column with typical operating conditions of prismatic designs. The tested Reynolds numbers of the coolant flow are 10,000, 20,000, 30,000 and 50,000. The predicted Nusselt numbers with the {nu}{sup 2}-f model are compared with the results by the other turbulence models (k-{epsilon} and SST) as well as the empirical correlations. (authors)
Single-Phase, Turbulent Heat-Transfer Friction-Factor Data Base Flow Enhanced Tb
Energy Science and Technology Software Center (OSTI)
1994-01-21
Heat-exchanger designers need to know what type of performance improvement can be obtained before they will consider enhanced tubes. In particular, they need access to the heat-transfer coefficients and friction-factor values of enhanced tube types that are commercially available. To compile these data from the numerous publications and reports in the open literature is a formidable task that can discourage the designer from using them. A computer program that contains a comprehensive data base withmoreĀ Ā» a search feature would be a handy tool for the designer to obtain an estimate of the performance improvement that can be obtained with a particular enhanced tube geometry. In addition, it would be a valuable tool for researchers who are developing and/or validating new prediction methods. This computer program can be used to obtain friction-factor and/or heat-transfer data for a broad range of internally enhanced tube geometries with forced-convective turbulent flow. The program has search features; that is the user can select data for tubes with a particular enhancement geometry range or data obtained from a particular source or publication. The friction factor data base contains nearly 5,000 points and the heat-transfer data base contains more than 4,700 points. About 360 different tube geometries are included from the 36 different sources. Data for tubes with similar geometries and the same and/or different types can be easily extracted with the sort feature of this data base and compared. Users of the program are heat-exchanger designers, enhanced tubing suppliers, and research organizations or academia who are developing or validating prediction methods.Ā«Ā less
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Minamoto, Yuki; Kolla, Hemanth; Grout, Ray W.; Gruber, Andrea; Chen, Jacqueline H.
2015-07-24
Here, three-dimensional direct numerical simulation results of a transverse syngas fuel jet in turbulent cross-flow of air are analyzed to study the influence of varying volume fractions of CO relative to H2 in the fuel composition on the near field flame stabilization. The mean flame stabilizes at a similar location for CO-lean and CO-rich cases despite the trend suggested by their laminar flame speed, which is higher for the CO-lean condition. To identify local mixtures having favorable mixture conditions for flame stabilization, explosive zones are defined using a chemical explosive mode timescale. The explosive zones related to flame stabilization aremoreĀ Ā» located in relatively low velocity regions. The explosive zones are characterized by excess hydrogen transported solely by differential diffusion, in the absence of intense turbulent mixing or scalar dissipation rate. The conditional averages show that differential diffusion is negatively correlated with turbulent mixing. Moreover, the local turbulent Reynolds number is insufficient to estimate the magnitude of the differential diffusion effect. Alternatively, the Karlovitz number provides a better indicator of the importance of differential diffusion. A comparison of the variations of differential diffusion, turbulent mixing, heat release rate and probability of encountering explosive zones demonstrates that differential diffusion predominantly plays an important role for mixture preparation and initiation of chemical reactions, closely followed by intense chemical reactions sustained by sufficient downstream turbulent mixing. The mechanism by which differential diffusion contributes to mixture preparation is investigated using the Takeno Flame Index. The mean Flame Index, based on the combined fuel species, shows that the overall extent of premixing is not intense in the upstream regions. However, the Flame Index computed based on individual contribution of H2 or CO species reveals that hydrogen
Minamoto, Yuki; Kolla, Hemanth; Grout, Ray W.; Gruber, Andrea; Chen, Jacqueline H.
2015-07-24
Here, three-dimensional direct numerical simulation results of a transverse syngas fuel jet in turbulent cross-flow of air are analyzed to study the influence of varying volume fractions of CO relative to H_{2} in the fuel composition on the near field flame stabilization. The mean flame stabilizes at a similar location for CO-lean and CO-rich cases despite the trend suggested by their laminar flame speed, which is higher for the CO-lean condition. To identify local mixtures having favorable mixture conditions for flame stabilization, explosive zones are defined using a chemical explosive mode timescale. The explosive zones related to flame stabilization are located in relatively low velocity regions. The explosive zones are characterized by excess hydrogen transported solely by differential diffusion, in the absence of intense turbulent mixing or scalar dissipation rate. The conditional averages show that differential diffusion is negatively correlated with turbulent mixing. Moreover, the local turbulent Reynolds number is insufficient to estimate the magnitude of the differential diffusion effect. Alternatively, the Karlovitz number provides a better indicator of the importance of differential diffusion. A comparison of the variations of differential diffusion, turbulent mixing, heat release rate and probability of encountering explosive zones demonstrates that differential diffusion predominantly plays an important role for mixture preparation and initiation of chemical reactions, closely followed by intense chemical reactions sustained by sufficient downstream turbulent mixing. The mechanism by which differential diffusion contributes to mixture preparation is investigated using the Takeno Flame Index. The mean Flame Index, based on the combined fuel species, shows that the overall extent of premixing is not intense in the upstream regions. However, the Flame Index computed based on individual contribution of H_{2} or CO species reveals that
Enhanced thermal and gas flow performance in a three-way catalytic...
Office of Energy Efficiency and Renewable Energy (EERE) Indexed Site
Enhanced thermal and gas flow performance in a three-way catalytic converter through use of insulation within the ceramic monolith Emissions performance comparison of conventional ...
Contribution to the numerical study of turbulence in high intensity discharge lamps
Kaziz, S.; Ben Ahmed, R.; Helali, H.; Gazzah, H.; Charrada, K. [Unite d'Etude des Milieux Ionises et Reactifs, IPEIM, 5019 route de Kairouan Monastir (Tunisia)
2011-07-15
We present in this paper a comparison between results obtained with a laminar and turbulent models for high-pressure mercury arc. The two models are based on the resolution of bidimensional time-dependent equations by a semi-implicit finite-element code. The numerical computation of turbulent model is solved with large eddy simulation model; this approach takes into account the various scales of turbulence by a filtering method on each scale. The results show the quantitative influence of turbulence on the flow fields and also the difference between laminar and turbulent effects on the dynamic thermal behaviour and on the characteristics of the discharge.
Comparative study of the vorticity field in turbulent flows: Theory, experiments, computations
Levich, E.
1991-09-01
The Goal of the project was to understand the role of topology vortex lines in general and the helicity invariant (inviscid) in particular for turbulent dynamics. The project consisted of three main ingredients: theoretical, numerical and experimental. The achievements and failures of the above are separately reported in this paper.
Levich, E.
1991-09-01
The Goal of the project was to understand the role of topology vortex lines in general and the helicity invariant (inviscid) in particular for turbulent dynamics. The project consisted of three main ingredients: theoretical, numerical and experimental. The achievements and failures of the above are separately reported in this paper.
Arima, T.; Sonoda, T.; Shirotori, M.; Tamura, A.; Kikuchi, K.
1999-01-01
The authors have developed a computer simulation code for three-dimensional viscous flow in turbomachinery based on the time-averaged compressible Navier-Stokes equations and a low-Reynolds-number {kappa}-{epsilon} turbulence model. It is described in detail in this paper. The code is used to compute the flow fields for two types of rotor (a transonic fan NASA Rotor 67 and a transonic axial compressor NASA rotor 37), and numerical results are compared to experimental data based on aerodynamic probe and laser anemometer measurements. In the case of Rotor 67, calculated and experimental results are compared under the design speed to validate the code. The calculated results show good agreement with the experimental data, such as the rotor performance map and the spanwise distribution of total pressure, total temperature, and flow angle downstream of the rotor. In the case of Rotor 37, detailed comparisons between the numerical results and the experimental data are made under the design speed condition to assess the overall quality of the numerical solution. Furthermore, comparisons under the part-speed condition are used to investigate a flow field without passage shock. The results are well predicted qualitatively. However, considerable quantitative discrepancies remain in predicting the flow near the tip. In order to assess the predictive capabilities of the developed code, computed flow structures are presented with the experimental data for each rotor and the cause of the discrepancies is discussed.
Punjabi, Sangeeta B.; Sahasrabudhe, S. N.; Das, A. K.; Joshi, N. K.; Mangalvedekar, H. A.; Kothari, D. C.
2014-01-15
This paper provides 2D comparative study of results obtained using laminar and turbulent flow model for RF (radio frequency) Inductively Coupled Plasma (ICP) torch. The study was done for the RF-ICP torch operating at 50?kW DC power and 3?MHz frequency located at BARC. The numerical modeling for this RF-ICP torch is done using ANSYS software with the developed User Defined Function. A comparative study is done between laminar and turbulent flow model to investigate how temperature and flow fields change when using different operating conditions such as (a) swirl and no swirl velocity for sheath gas flow rate, (b) variation in sheath gas flow rate, and (c) variation in plasma gas flow rate. These studies will be useful for different material processing applications.
Heat Transfer and Friction-Factor Methods Turbulent Flow Inside Pipes 3d Rough
Energy Science and Technology Software Center (OSTI)
1994-01-21
Three-dimensional roughened internally enhanced tubes have been shown to be one of the most energy efficient for turbulent, forced convection applications. However, there is only one prediction method presented in the open literature and that is restricted to three-dimensional sand-grain roughness. Other roughness types are being proposed: hemispherical sectors, truncated cones, and full and truncated pyramids. There are no validated heat-transfer and friction-factor prediction methods for these different roughness shapes that can be used inmoreĀ Ā» the transition and fully rough region. This program calculates the Nusselt number and friction factor values, for a broad range of three-dimensional roughness types such as hemispherical sectors, truncated cones, and full and truncated pyramids. Users of this program are heat-exchangers designers, enhanced tubing suppliers, and research organizations or academia who are developing or validating prediction methods.Ā«Ā less
Talbot, L.; Cheng, R.K.
1993-12-01
Turbulent combustion is the dominant process in heat and power generating systems. Its most significant aspect is to enhance the burning rate and volumetric power density. Turbulent mixing, however, also influences the chemical rates and has a direct effect on the formation of pollutants, flame ignition and extinction. Therefore, research and development of modern combustion systems for power generation, waste incineration and material synthesis must rely on a fundamental understanding of the physical effect of turbulence on combustion to develop theoretical models that can be used as design tools. The overall objective of this program is to investigate, primarily experimentally, the interaction and coupling between turbulence and combustion. These processes are complex and are characterized by scalar and velocity fluctuations with time and length scales spanning several orders of magnitude. They are also influenced by the so-called {open_quotes}field{close_quotes} effects associated with the characteristics of the flow and burner geometries. The authors` approach is to gain a fundamental understanding by investigating idealized laboratory flames. Laboratory flames are amenable to detailed interrogation by laser diagnostics and their flow geometries are chosen to simplify numerical modeling and simulations and to facilitate comparison between experiments and theory.
Heat Flow in VC-2A and VC-2B, and Constraints on the Thermal...
to: navigation, search OpenEI Reference LibraryAdd to library Conference Paper: Heat Flow in VC-2A and VC-2B, and Constraints on the Thermal Regime of the Valles Caldera, New...
Electromagnetic Transport From Microtearing Mode Turbulence
Guttenfelder, W; Kaye, S M; Nevins, W M; Wang, E; Bell, R E; Hammett, G W; LeBlanc, B P; Mikkelsen, D R
2011-03-23
This Letter presents non-linear gyrokinetic simulations of microtearing mode turbulence. The simulations include collisional and electromagnetic effects and use experimental parameters from a high beta discharge in the National Spherical Torus Experiment (NSTX). The predicted electron thermal transport is comparable to that given by experimental analysis, and it is dominated by the electromagnetic contribution of electrons free streaming along the resulting stochastic magnetic field line trajectories. Experimental values of flow shear can significantly reduce the predicted transport.
Sudden Viscous Dissipation of Compressing Turbulence
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Davidovits, Seth; Fisch, Nathaniel J.
2016-03-11
Here we report compression of turbulent plasma can amplify the turbulent kinetic energy, if the compression is fast compared to the viscous dissipation time of the turbulent eddies. A sudden viscous dissipation mechanism is demonstrated, whereby this amplified turbulent kinetic energy is rapidly converted into thermal energy, suggesting a new paradigm for fast ignition inertial fusion.
Internal thermal coupling in direct-flow coaxial vacuum tube collectors
Glembin, J.; Rockendorf, G.; Scheuren, J.
2010-07-15
This investigation covers the impact of low flow rates on the efficiency of coaxial vacuum tube collectors. Measurements show an efficiency reduction of 10% if reducing the flow rate from 78 kg/m{sup 2} h to 31 kg/m{sup 2} h for a collector group with 60 parallel vacuum tubes with a coaxial flow conduit at one-sided connection. For a more profound understanding a model of the coaxial tube was developed which defines the main energy fluxes including the internal thermal coupling. The tube simulations show a non-linear temperature profile along the tube with the maximum temperature in the outer pipe. Due to heat transfer to the entering flow this maximum is not located at the fluid outlet. The non-linearity increases with decreasing flow rates. The experimentally determined flow distribution allows simulating the measured collector array. The simulation results confirm the efficiency decrease at low flow rates. The flow distribution has a further impact on efficiency reduction, but even at an ideal uniform flow, a considerable efficiency reduction at low flow rates is to be expected. As a consequence, low flow rates should be prevented for coaxial tube collectors, thus restricting the possible operation conditions. The effect of constructional modifications like diameter or material variations is presented. Finally the additional impact of a coaxial manifold design is discussed. (author)
Flow distribution analysis on the cooling tube network of ITER thermal shield
Nam, Kwanwoo; Chung, Wooho; Noh, Chang Hyun; Kang, Dong Kwon; Kang, Kyoung-O; Ahn, Hee Jae; Lee, Hyeon Gon
2014-01-29
Thermal shield (TS) is to be installed between the vacuum vessel or the cryostat and the magnets in ITER tokamak to reduce the thermal radiation load to the magnets operating at 4.2K. The TS is cooled by pressurized helium gas at the inlet temperature of 80K. The cooling tube is welded on the TS panel surface and the composed flow network of the TS cooling tubes is complex. The flow rate in each panel should be matched to the thermal design value for effective radiation shielding. This paper presents one dimensional analysis on the flow distribution of cooling tube network for the ITER TS. The hydraulic cooling tube network is modeled by an electrical analogy. Only the cooling tube on the TS surface and its connecting pipe from the manifold are considered in the analysis model. Considering the frictional factor and the local loss in the cooling tube, the hydraulic resistance is expressed as a linear function with respect to mass flow rate. Sub-circuits in the TS are analyzed separately because each circuit is controlled by its own control valve independently. It is found that flow rates in some panels are insufficient compared with the design values. In order to improve the flow distribution, two kinds of design modifications are proposed. The first one is to connect the tubes of the adjacent panels. This will increase the resistance of the tube on the panel where the flow rate is excessive. The other design suggestion is that an orifice is installed at the exit of tube routing where the flow rate is to be reduced. The analysis for the design suggestions shows that the flow mal-distribution is improved significantly.
Tidal Flow Turbulence Measurements
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
What)can)we)measure?) * Acous+c)Doppler)Current)Profilers)(ADCP)) - Temporal)noise) - Spa+al)aliasing)) * Acous+c)Doppler)Velocimeters)(ADVs)) * Other)op+ons) Northwest)Na+onal)Mar...
Thermal/chemical degradation of ceramic cross-flow filter materials
Alvin, M.A.; Lane, J.E.; Lippert, T.E.
1989-11-01
This report summarizes the 14-month, Phase 1 effort conducted by Westinghouse on the Thermal/Chemical Degradation of Ceramic Cross-Flow Filter Materials program. In Phase 1 expected filter process conditions were identified for a fixed-bed, fluid-bed, and entrained-bed gasification, direct coal fired turbine, and pressurized fluidized-bed combustion system. Ceramic cross-flow filter materials were also selected, procured, and subjected to chemical and physical characterization. The stability of each of the ceramic cross-flow materials was assessed in terms of potential reactions or phase change as a result of process temperature, and effluent gas compositions containing alkali and fines. In addition chemical and physical characterization was conducted on cross-flow filters that were exposed to the METC fluid-bed gasifier and the New York University pressurized fluidized-bed combustor. Long-term high temperature degradation mechanisms were proposed for each ceramic cross-flow material at process operating conditions. An experimental bench-scale test program is recommended to be conducted in Phase 2, generating data that support the proposed cross-flow filter material thermal/chemical degradation mechanisms. Papers on the individual subtasks have been processed separately for inclusion on the data base.
Effects Of Thermal Exchange On Material Flow During Steel Thixoextrusion Process
Becker, Eric; Gu Guochao; Langlois, Laurent; Bigot, Regis; Pesci, Raphael
2011-01-17
Semisolid processing is an innovative technology for near net-shape production of components, where the metallic alloys are processed in the semisolid state. Taking advantage of the thixotropic behavior of alloys in the semisolid state, significant progress has been made in semisolid processing. However, the consequences of such behavior on the flow during thixoforming are still not completely understood. To explore and better understand the influence of the different parameters on material flow during thixoextrusion process, thixoextrusion experiments were performed using the low carbon steel C38. The billet was partially melted at high solid fraction. Effects of various process parameters including the initial billet temperature, the temperature of die, the punch speed during process and the presence of a Ceraspray layer at the interface of tool and billet were investigated through experiments and simulation. After analyzing the results thus obtained, it was identified that the aforementioned parameters mainly affect thermal exchanges between die and part. The Ceraspray layer not only plays a lubricant role, but also acts as a thermal barrier at the interface of tool and billet. Furthermore, the thermal effects can affect the material flow which is composed of various distinct zones.
An experimental investigation of thermal mixing and combustion in supersonic flows
Srikrishnan, A.R.; Kurian, J.; Sriramulu, V.
1996-12-01
A radially lobed nozzle (petal nozzle) is being increasingly recognized as a potential candidate for promoting mixing in compressible flows. An experimental investigation has been conducted to study its effectiveness in improving thermal mixing and combustion in supersonic flow. A hot gas jet issuing supersonically from a lobed nozzle mixes with a cold supersonic jet in a circular mixing tube. The two jets issue coaxially. A detailed survey of the flow field inside the mixing duct reveals that nearly complete thermal mixing (as exemplified by the nearly uniform temperature distribution) could be achieved in a short distance when a lobed nozzle is employed. The results also indicate the presence of large-scale vortices in the flow field downstream of the lobed nozzle. Having thus created a field in which mixing is good, supersonic combustion was then attempted. Kerosene was introduced into the hot stream issuing from the lobed nozzle and it burned mainly in the mixing tube, which served as a supersonic combustor. Resulting temperature and pressure rises were measured and the supersonic combustion efficiency was found to be of the order of 60%. The performance of a conventional conical nozzle was found to be much inferior to that of the petal nozzle under identical conditions.
Advances in compressible turbulent mixing
Dannevik, W.P.; Buckingham, A.C.; Leith, C.E.
1992-01-01
This volume includes some recent additions to original material prepared for the Princeton International Workshop on the Physics of Compressible Turbulent Mixing, held in 1988. Workshop participants were asked to emphasize the physics of the compressible mixing process rather than measurement techniques or computational methods. Actual experimental results and their meaning were given precedence over discussions of new diagnostic developments. Theoretical interpretations and understanding were stressed rather than the exposition of new analytical model developments or advances in numerical procedures. By design, compressibility influences on turbulent mixing were discussed--almost exclusively--from the perspective of supersonic flow field studies. The papers are arranged in three topical categories: Foundations, Vortical Domination, and Strongly Coupled Compressibility. The Foundations category is a collection of seminal studies that connect current study in compressible turbulent mixing with compressible, high-speed turbulent flow research that almost vanished about two decades ago. A number of contributions are included on flow instability initiation, evolution, and transition between the states of unstable flow onset through those descriptive of fully developed turbulence. The Vortical Domination category includes theoretical and experimental studies of coherent structures, vortex pairing, vortex-dynamics-influenced pressure focusing. In the Strongly Coupled Compressibility category the organizers included the high-speed turbulent flow investigations in which the interaction of shock waves could be considered an important source for production of new turbulence or for the enhancement of pre-existing turbulence. Individual papers are processed separately.
Thermal characteristics of air flow cooling in the lithium ion batteries experimental chamber
Lukhanin A.; Rohatgi U.; Belyaev, A.; Fedorchenko, D.; Khazhmuradov, M.; Lukhanin, O; Rudychev, I.
2012-07-08
A battery pack prototype has been designed and built to evaluate various air cooling concepts for the thermal management of Li-ion batteries. The heat generation from the Li-Ion batteries was simulated with electrical heat generation devices with the same dimensions as the Li-Ion battery (200 mm x 150 mm x 12 mm). Each battery simulator generates up to 15W of heat. There are 20 temperature probes placed uniformly on the surface of the battery simulator, which can measure temperatures in the range from -40 C to +120 C. The prototype for the pack has up to 100 battery simulators and temperature probes are recorder using a PC based DAQ system. We can measure the average surface temperature of the simulator, temperature distribution on each surface and temperature distributions in the pack. The pack which holds the battery simulators is built as a crate, with adjustable gap (varies from 2mm to 5mm) between the simulators for air flow channel studies. The total system flow rate and the inlet flow temperature are controlled during the test. The cooling channel with various heat transfer enhancing devices can be installed between the simulators to investigate the cooling performance. The prototype was designed to configure the number of cooling channels from one to hundred Li-ion battery simulators. The pack is thermally isolated which prevents heat transfer from the pack to the surroundings. The flow device can provide the air flow rate in the gap of up to 5m/s velocity and air temperature in the range from -30 C to +50 C. Test results are compared with computational modeling of the test configurations. The present test set up will be used for future tests for developing and validating new cooling concepts such as surface conditions or heat pipes.
The high-energy-density counterpropagating shear experiment and turbulent self-heating
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Doss, F. W.; Fincke, J. R.; Loomis, E. N.; Welser-Sherrill, L.; Flippo, K. A.
2013-12-06
The counterpropagating shear experiment has previously demonstrated the ability to create regions of shockdriven shear, balanced symmetrically in pressure and experiencing minimal net drift. This allows for the creation of a high-Mach-number high-energy-density shear environment. New data from the counterpropagating shear campaign is presented, and both hydrocode modeling and theoretical analysis in the context of a Reynolds-averaged-Navier-Stokes model suggest turbulent dissipation of energy from the supersonic flow bounding the layer is a significant driver in its expansion. A theoretical minimum shear flow Mach number threshold is suggested for substantial thermal-turbulence coupling.
Evidence for Radial Flow of Thermal Dileptons in High-Energy Nuclear Collisions
Arnaldi, R.; Colla, A.; Cortese, P.; Ferretti, A.; Oppedisano, C.; Scomparin, E.; Banicz, K.; Damjanovic, S.; Castor, J.; Devaux, A.; Fargeix, J.; Force, P.; Manso, F.; Chaurand, B.; Cicalo, C.; Falco, A. de; Floris, M.; Masoni, A.; Puddu, G.; Serci, S.
2008-01-18
The NA60 experiment at the CERN SPS has studied low-mass dimuon production in 158A GeV In-In collisions. An excess of pairs above the known meson decays has been reported before. We now present precision results on the associated transverse momentum spectra. The slope parameter T{sub eff} extracted from the spectra rises with dimuon mass up to the {rho}, followed by a sudden decline above. While the initial rise is consistent with the expectations for radial flow of a hadronic decay source, the decline signals a transition to an emission source with much smaller flow. This may well represent the first direct evidence for thermal radiation of partonic origin in nuclear collisions.
Thermally induced flow oscillation in vertical two-phase natural circulation loop
Lee, Sang Yong; Ishii, Mamoru
1988-01-01
In order to study the two-phase natural circulation during a small break loss of coolant accident in LWR, simulation experiments have been performed using Freon-113 boiling and condensation loop. In quasi-steady state, the flow became relatively stabilized and certain regular patterns of flow oscillations were detected with ranges of periods in 8-/approximately/35 seconds and 2.5-/approximately/4 minutes. In order to find out the nature of these oscillations, one-dimensional field equations for the single-phase (liquid) and two-phase region were set up, and these field equations were integrated along the loop. The homogeneous flow model was used for the two-phase region. Then the characteristic equation was derived using perturbation method. Thermal non-equilibrium and compressibility of each phase were not considered in the present analysis. The characteristic equation derived can be used to obtain the stability criteria. A simplified approach showed that the short-period oscillation were the manometer oscillation. The longer period oscillations were the density wave oscillation which had the period of oscillations close to the residence time of a fluid around the loop.
Turbulent electron transport in edge pedestal by electron temperature gradient turbulence
Singh, R.; Institute for Plasma Research, Bhat Gandhinagar, Gujarat 2382 428 ; Jhang, Hogun; Diamond, P. H.; CMTFO and CASS, University of California, San Diego 92093-0424, California
2013-11-15
We present a model for turbulent electron thermal transport at the edge pedestal in high (H)-mode plasmas based on electron temperature gradient (ETG) turbulence. A quasi-linear analysis of electrostatic toroidal ETG modes shows that both turbulent electron thermal diffusivity and hyper-resistivity exhibits the Ohkawa scaling in which the radial correlation length of turbulence becomes the order of electron skin depth. Combination of the Ohkawa scales and the plasma current dependence results in a novel confinement scaling inside the pedestal region. It is also shown that ETG turbulence induces a thermoelectric pinch, which may accelerate the density pedestal formation.
Simulation of microtearing turbulence in national spherical torus experiment
Guttenfelder, W.; Kaye, S. M.; Bell, R. E.; Hammett, G. W.; LeBlanc, B. P.; Mikkelsen, D. R.; Ren, Y.; Candy, J.; Nevins, W. M.; Wang, E.; Zhang, J.; Crocker, N. A.; Yuh, H.
2012-05-15
Thermal energy confinement times in National Spherical Torus Experiment (NSTX) dimensionless parameter scans increase with decreasing collisionality. While ion thermal transport is neoclassical, the source of anomalous electron thermal transport in these discharges remains unclear, leading to considerable uncertainty when extrapolating to future spherical tokamak (ST) devices at much lower collisionality. Linear gyrokinetic simulations find microtearing modes to be unstable in high collisionality discharges. First non-linear gyrokinetic simulations of microtearing turbulence in NSTX show they can yield experimental levels of transport. Magnetic flutter is responsible for almost all the transport ({approx}98%), perturbed field line trajectories are globally stochastic, and a test particle stochastic transport model agrees to within 25% of the simulated transport. Most significantly, microtearing transport is predicted to increase with electron collisionality, consistent with the observed NSTX confinement scaling. While this suggests microtearing modes may be the source of electron thermal transport, the predictions are also very sensitive to electron temperature gradient, indicating the scaling of the instability threshold is important. In addition, microtearing turbulence is susceptible to suppression via sheared E Multiplication-Sign B flows as experimental values of E Multiplication-Sign B shear (comparable to the linear growth rates) dramatically reduce the transport below experimental values. Refinements in numerical resolution and physics model assumptions are expected to minimize the apparent discrepancy. In cases where the predicted transport is strong, calculations suggest that a proposed polarimetry diagnostic may be sensitive to the magnetic perturbations associated with the unique structure of microtearing turbulence.
The interaction of high-speed turbulence with flames: Turbulent flame speed
Poludnenko, A.Y.; Oran, E.S. [Laboratory for Computational Physics and Fluid Dynamics, Naval Research Laboratory, Washington, DC 20375 (United States)
2011-02-15
Direct numerical simulations of the interaction of a premixed flame with driven, subsonic, homogeneous, isotropic, Kolmogorov-type turbulence in an unconfined system are used to study the mechanisms determining the turbulent flame speed, S{sub T}, in the thin reaction zone regime. High intensity turbulence is considered with the r.m.s. velocity 35 times the laminar flame speed, S{sub L}, resulting in the Damkoehler number Da=0.05. The simulations were performed with Athena-RFX, a massively parallel, fully compressible, high-order, dimensionally unsplit, reactive-flow code. A simplified reaction-diffusion model, based on the one-step Arrhenius kinetics, represents a stoichiometric H{sub 2}-air mixture under the assumption of the Lewis number Le=1. Global properties and the internal structure of the flame were analyzed in an earlier paper, which showed that this system represents turbulent combustion in the thin reaction zone regime. This paper demonstrates that: (1) The flame brush has a complex internal structure, in which the isosurfaces of higher fuel mass fractions are folded on progressively smaller scales. (2) Global properties of the turbulent flame are best represented by the structure of the region of peak reaction rate, which defines the flame surface. (3) In the thin reaction zone regime, S{sub T} is predominantly determined by the increase of the flame surface area, A{sub T}, caused by turbulence. (4) The observed increase of S{sub T} relative to S{sub L} exceeds the corresponding increase of A{sub T} relative to the surface area of the planar laminar flame, on average, by {approx}14%, varying from only a few percent to as high as {approx}30%. (5) This exaggerated response is the result of tight flame packing by turbulence, which causes frequent flame collisions and formation of regions of high flame curvature >or similar 1/{delta}{sub L}, or ''cusps,'' where {delta}{sub L} is the thermal width of the laminar flame. (6) The local flame speed in the cusps
Monte, F. de; Galli, G.; Marcotullio, F.
1996-12-31
A closed form-expression for the effectiveness of the heat exchangers and regenerator of a Stirling cycle machine is given. This result may be used in a simple way in order to evaluate their effect on the machine performance. The proposed method, indeed, allows the actual cycle gas temperatures in the heater and cooler to be obtained readily, once the geometry of the heater, cooler and regenerator is known and some quantities characterizing the engine dynamics (strokes, frequency and phase angle of the moving elements) and its heat-exchange processes (inlet temperatures of the heating and cooling fluids, and their volumetric flow rates) are measured. Thus, an immediate indication about the effectiveness of the heat exchangers and regenerator as well as about the machine thermal efficiency may be obtained. The availability of a closed-form expression for the heater, regenerator and cooler effectiveness is useful especially for those engines, like the free-piston Stirling engines, whose design requires the application of analytically based optimization criteria.
Pressure atomizer having multiple orifices and turbulent generation feature
VanBrocklin, Paul G.; Geiger, Gail E.; Moran, Donald James; Fournier, Stephane
2002-01-01
A pressure atomizer includes a silicon plate having a top surface and a bottom surface. A portion of the top surface defines a turbulent chamber. The turbulent chamber is peripherally bounded by the top surface of the plate. The turbulent chamber is recessed a predetermined depth relative to the top surface. The silicon plate further defines at least one flow orifice. Each flow orifice extends from the bottom surface of the silicon plate to intersect with and open into the turbulent chamber. Each flow orifice is in fluid communication with the turbulent chamber.
First-order particle acceleration in magnetically driven flows
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Beresnyak, Andrey; Li, Hui
2016-03-02
In this study, we demonstrate that particles are regularly accelerated while experiencing curvature drift in flows driven by magnetic tension. Some examples of such flows include spontaneous turbulent reconnection and decaying magnetohydrodynamic turbulence, where a magnetic field relaxes to a lower-energy configuration and transfers part of its energy to kinetic motions of the fluid. We show that this energy transfer, which normally causes turbulent cascade and heating of the fluid, also results in a first-order acceleration of non-thermal particles. Since it is generic, this acceleration mechanism is likely to play a role in the production of non-thermal particle distribution inmoreĀ Ā» magnetically dominant environments such as the solar chromosphere, pulsar magnetospheres, jets from supermassive black holes, and Ī³-ray bursts.Ā«Ā less
Anh Bui; Nam Dinh; Brian Williams
2013-09-01
In addition to validation data plan, development of advanced techniques for calibration and validation of complex multiscale, multiphysics nuclear reactor simulation codes are a main objective of the CASL VUQ plan. Advanced modeling of LWR systems normally involves a range of physico-chemical models describing multiple interacting phenomena, such as thermal hydraulics, reactor physics, coolant chemistry, etc., which occur over a wide range of spatial and temporal scales. To a large extent, the accuracy of (and uncertainty in) overall model predictions is determined by the correctness of various sub-models, which are not conservation-laws based, but empirically derived from measurement data. Such sub-models normally require extensive calibration before the models can be applied to analysis of real reactor problems. This work demonstrates a case study of calibration of a common model of subcooled flow boiling, which is an important multiscale, multiphysics phenomenon in LWR thermal hydraulics. The calibration process is based on a new strategy of model-data integration, in which, all sub-models are simultaneously analyzed and calibrated using multiple sets of data of different types. Specifically, both data on large-scale distributions of void fraction and fluid temperature and data on small-scale physics of wall evaporation were simultaneously used in this works calibration. In a departure from traditional (or common-sense) practice of tuning/calibrating complex models, a modern calibration technique based on statistical modeling and Bayesian inference was employed, which allowed simultaneous calibration of multiple sub-models (and related parameters) using different datasets. Quality of data (relevancy, scalability, and uncertainty) could be taken into consideration in the calibration process. This work presents a step forward in the development and realization of the CIPS Validation Data Plan at the Consortium for Advanced Simulation of LWRs to enable
Can we characterize turbulence in premixed flames?
Lipatnikov, A.N. [Department of Applied Mechanics, Chalmers University of Technology, Gothenburg, 412 96 (Sweden)
2009-06-15
Modeling of premixed turbulent combustion involves averaging reaction rates in turbulent flows. The focus of most approaches to resolving this problem has been placed on determining the dependence of the mean rate w of product creation on the laminar flame speed S{sub L}, the rms turbulence velocity u', etc. The goal of the present work is to draw attention to another issue: May the input quantity u{sup '} for a model of w= w(u'/S{sub L},..) be considered to be known? The point is that heat release substantially affects turbulence and, hence, turbulence characteristics in premixed flames should be modeled. However, standard moment methods for numerically simulating turbulent flows do not allow us to evaluate the true turbulence characteristics in a flame. For instance, the Reynolds stresses in premixed flames are affected not only by turbulence itself, but also by velocity jump across flamelets. A common way to resolving this problem consists of considering the Reynolds stresses conditioned on unburned (or burned) mixture to be the true turbulence characteristics. In the present paper, this widely accepted but never proved hypothesis is put into question, first, by considering simple model constant-density problems (flame motion in an oscillating one-dimensional laminar flow; flame stabilized in a periodic shear, one-dimensional, laminar flow; turbulent mixing). In all the cases, the magnitude of velocity fluctuations, calculated using the conditioned Reynolds stresses, is affected by the intermittency of reactants and products and, hence, is not the true rms velocity. Second, the above claim is further supported by comparing balance equations for the mean and conditioned Reynolds stresses. The conditioned Reynolds stresses do not characterize the true turbulence in flames, because conditional averaging cuts off flow regions characterized by either high or low velocities. (author)
Cyclone separator having boundary layer turbulence control
Krishna, Coimbatore R.; Milau, Julius S.
1985-01-01
A cyclone separator including boundary layer turbulence control that is operable to prevent undue build-up of particulate material at selected critical areas on the separator walls, by selectively varying the fluid pressure at those areas to maintain the momentum of the vortex, thereby preventing particulate material from inducing turbulence in the boundary layer of the vortical fluid flow through the separator.
Aqueous Solution Vessel Thermal Model Development II
Buechler, Cynthia Eileen
2015-10-28
The work presented in this report is a continuation of the work described in the May 2015 report, āAqueous Solution Vessel Thermal Model Developmentā. This computational fluid dynamics (CFD) model aims to predict the temperature and bubble volume fraction in an aqueous solution of uranium. These values affect the reactivity of the fissile solution, so it is important to be able to calculate them and determine their effects on the reaction. Part A of this report describes some of the parameter comparisons performed on the CFD model using Fluent. Part B describes the coupling of the Fluent model with a Monte-Carlo N-Particle (MCNP) neutron transport model. The fuel tank geometry is the same as it was in the May 2015 report, annular with a thickness-to-height ratio of 0.16. An accelerator-driven neutron source provides the excitation for the reaction, and internal and external water cooling channels remove the heat. The model used in this work incorporates the Eulerian multiphase model with lift, wall lubrication, turbulent dispersion and turbulence interaction. The buoyancy-driven flow is modeled using the Boussinesq approximation, and the flow turbulence is determined using the k-Ļ Shear-Stress-Transport (SST) model. The dispersed turbulence multiphase model is employed to capture the multiphase turbulence effects.
Lagrangian-Averaged Scale-Dependent subfilter turbulence model
Energy Science and Technology Software Center (OSTI)
2011-03-01
LASD are Fortran 90 modules that compute the stresses and scalar fluxes arising from unrelolved scales of turbulence, required for large-eddy eimulations of fluid flows.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
SjĆ¶berg, Ylva; Coon, Ethan; K. Sannel, A. Britta; Pannetier, Romain; Harp, Dylan; Frampton, Andrew; Painter, Scott L.; Lyon, Steve W.
2016-02-05
Modeling and observation of ground temperature dynamics are the main tools for understanding current permafrost thermal regimes and projecting future thaw. Until recently, most studies on permafrost have focused on vertical ground heat fluxes. Groundwater can transport heat in both lateral and vertical directions but its influence on ground temperatures at local scales in permafrost environments is not well understood. In this paper, we combine field observations from a subarctic fen in the sporadic permafrost zone with numerical simulations of coupled water and thermal fluxes. At the Tavvavuoma study site in northern Sweden, ground temperature profiles and groundwater levels weremoreĀ Ā» observed in boreholes. These observations were used to set up one- and two-dimensional simulations down to 2 m depth across a gradient of permafrost conditions within and surrounding the fen. Two-dimensional scenarios representing the fen under various hydraulic gradients were developed to quantify the influence of groundwater flow on ground temperature. Our observations suggest that lateral groundwater flow significantly affects ground temperatures. This is corroborated by modeling results that show seasonal ground ice melts 1 month earlier when a lateral groundwater flux is present. Further, although the thermal regime may be dominated by vertically conducted heat fluxes during most of the year, isolated high groundwater flow rate events such as the spring freshet are potentially important for ground temperatures. Finally, as sporadic permafrost environments often contain substantial portions of unfrozen ground with active groundwater flow paths, knowledge of this heat transport mechanism is important for understanding permafrost dynamics in these environments.Ā«Ā less
ENERGY EFFICIENT THERMAL MANAGEMENT FOR NATURAL GAS ENGINE AFTERTREATMENT VIA ACTIVE FLOW CONTROL
David K. Irick; Ke Nguyen
2004-04-01
The project is focused on the development of an energy efficient aftertreatment system capable of reducing NOx and methane by 90% from lean-burn natural gas engines by applying active exhaust flow control. Compared to conventional passive flow-through reactors, the proposed scheme cuts supplemental energy by 50%-70%. The system consists of a Lean NOx Trap (LNT) system and an oxidation catalyst. Through alternating flow control, a major amount of engine exhaust flows through a large portion of the LNT system in the absorption mode, while a small amount of exhaust goes through a small portion of the LNT system in the regeneration or desulfurization mode. By periodically reversing the exhaust gas flow through the oxidation catalyst, a higher temperature profile is maintained in the catalyst bed resulting in greater efficiency of the oxidation catalyst at lower exhaust temperatures. The project involves conceptual design, theoretical analysis, computer simulation, prototype fabrication, and empirical studies. This report details the progress during the first twelve months of the project. The primary activities have been to develop the bench flow reactor system, develop the computer simulation and modeling of the reverse-flow oxidation catalyst, install the engine into the test cell, and begin design of the LNT system.
Energy Efficient Thermal Management for Natural Gas Engine Aftertreatment via Active Flow Control
David K. Irick; Ke Nguyen; Vitacheslav Naoumov; Doug Ferguson
2006-04-01
The project is focused on the development of an energy efficient aftertreatment system capable of reducing NOx and methane by 90% from lean-burn natural gas engines by applying active exhaust flow control. Compared to conventional passive flow-through reactors, the proposed scheme cuts supplemental energy by 50%-70%. The system consists of a Lean NOx Trap (LNT) system and an oxidation catalyst. Through alternating flow control, a major amount of engine exhaust flows through a large portion of the LNT system in the absorption mode, while a small amount of exhaust goes through a small portion of the LNT system in the regeneration or desulfurization mode. By periodically reversing the exhaust gas flow through the oxidation catalyst, a higher temperature profile is maintained in the catalyst bed resulting in greater efficiency of the oxidation catalyst at lower exhaust temperatures. The project involves conceptual design, theoretical analysis, computer simulation, prototype fabrication, and empirical studies. This report details the progress during the first twelve months of the project. The primary activities have been to develop the bench flow reactor system, develop the computer simulation and modeling of the reverse-flow oxidation catalyst, install the engine into the test cell, and begin design of the LNT system.
Energy Efficient Thermal Management for Natural Gas Engine Aftertreatment via Active Flow Control
David K. Irick; Ke Nguyen; Vitacheslav Naoumov; Doug Ferguson
2005-04-01
The project is focused on the development of an energy efficient aftertreatment system capable of reducing NOx and methane by 90% from lean-burn natural gas engines by applying active exhaust flow control. Compared to conventional passive flow-through reactors, the proposed scheme cuts supplemental energy by 50%-70%. The system consists of a Lean NOx Trap (LNT) system and an oxidation catalyst. Through alternating flow control, a major amount of engine exhaust flows through a large portion of the LNT system in the absorption mode, while a small amount of exhaust goes through a small portion of the LNT system in the regeneration or desulfurization mode. By periodically reversing the exhaust gas flow through the oxidation catalyst, a higher temperature profile is maintained in the catalyst bed resulting in greater efficiency of the oxidation catalyst at lower exhaust temperatures. The project involves conceptual design, theoretical analysis, computer simulation, prototype fabrication, and empirical studies. This report details the progress during the first twelve months of the project. The primary activities have been to develop the bench flow reactor system, develop the computer simulation and modeling of the reverse-flow oxidation catalyst, install the engine into the test cell, and begin design of the LNT system.
Cleveland, J.C.
1988-01-01
The analyses reported here have been conducted at the Oak Ridge National Laboratory (ORNL) for the US Nuclear Regulatory Commission's (NRC's) Division of Regulatory Applications of the Office of Nuclear Regulatory Research. The short-term thermal response of the Modular High-Temperature Gas-Cooled Reactor (MHTGR) is analyzed for a range of flow and reactivity transients. These include loss of forced circulation (LOFC) without scram, moisture ingress, spurious withdrawal of a control rod group, hypothetical large and rapid positive reactivity insertion, and a rapid core cooling event. The coupled heat transfer-neutron kinetics model is also described.
D. D. Blackwell; K. W. Wisian; M. C. Richards; J. L. Steele
2000-04-01
Several activities related to geothermal resources in the western United States are described in this report. A database of geothermal site-specific thermal gradient and heat flow results from individual exploration wells in the western US has been assembled. Extensive temperature gradient and heat flow exploration data from the active exploration of the 1970's and 1980's were collected, compiled, and synthesized, emphasizing previously unavailable company data. Examples of the use and applications of the database are described. The database and results are available on the world wide web. In this report numerical models are used to establish basic qualitative relationships between structure, heat input, and permeability distribution, and the resulting geothermal system. A series of steady state, two-dimensional numerical models evaluate the effect of permeability and structural variations on an idealized, generic Basin and Range geothermal system and the results are described.
Yang, Q. Q. Zhong, F. C. E-mail: fczhong@dhu.edu.cn; Jia, M. N.; Xu, G. S. E-mail: fczhong@dhu.edu.cn; Wang, L.; Wang, H. Q.; Chen, R.; Yan, N.; Liu, S. C.; Chen, L.; Li, Y. L.; Liu, J. B.
2015-06-15
The power fall-off width in the H-mode scrape-off layer (SOL) in tokamaks shows a strong inverse dependence on the plasma current, which was noticed by both previous multi-machine scaling work [T. Eich et al., Nucl. Fusion 53, 093031 (2013)] and more recent work [L. Wang et al., Nucl. Fusion 54, 114002 (2014)] on the Experimental Advanced Superconducting Tokamak. To understand the underlying physics, probe measurements of three H-mode discharges with different plasma currents have been studied in this work. The results suggest that a higher plasma current is accompanied by a stronger EĆB shear and a shorter radial correlation length of turbulence in the SOL, thus resulting in a narrower power fall-off width. A simple model has also been applied to demonstrate the suppression effect of EĆB shear on turbulence in the SOL and shows relatively good agreement with the experimental observations.
Imaging Fluid Flow in Geothermal Wells Using Distributed Thermal Perturbation Sensing
Broader source: Energy.gov [DOE]
Project objective: A New Geothermal Well Imaging Tool. 1.To develop a robust and easily deployable DTPS for monitoring in geothermal wells; and 2. Develop the associated analysis methodology for flow imaging; andĀwhen possible by wellbore conditionsĀto determine in situthermal conductivity and basal heat flux.
Supercritical carbon dioxide tubular flow under temporally varying thermal boundary condition
Son, H. M.; Halimi, B.; Suh, K. Y.
2012-07-01
During transient operation of fusion power plants the amount of thermal energy transferred from plasma to surrounding blanket modules will be varied over time, and will affect behavior of the working fluid inside the blanket and power conversion system where the coolant is in a supercritical state. Transient behavior of the power is in pulsed state in tokamak. The Optimized Supercritical Cycle Operation (OSCO) loop is constructed to investigate the thermohydraulic characteristics of the supercritical fluid under temporally varying thermal boundary condition. In this study the tube outer wall temperature data are measured for abrupt change in thermal power as a preliminary power transient test. The OSCO test conditions are selected to include the erratic behavior of the supercritical fluid under pseudo-critical condition during transient. In order to incorporate the delayed response of utilized thermocouples, a time constant is applied to adjust the obtained results. Along with the experimental study, computational fluid dynamic software is used to perform detailed analysis over the test section geometry. The preliminary test results are presented for comparison against the available correlations from the literature. (authors)
Uzdensky, Dmitri A.
2013-10-01
In this paper, we consider two outstanding intertwined problems in modern high-energy astrophysics: (1) the vertical-thermal structure of an optically thick accretion disk heated by the dissipation of magnetohydrodynamic turbulence driven by the magnetorotational instability (MRI), and (2) determining the fraction of the accretion power released in the corona above the disk. For simplicity, we consider a gas-pressure-dominated disk and assume a constant opacity. We argue that the local turbulent dissipation rate due to the disruption of the MRI channel flows by secondary parasitic instabilities should be uniform across most of the disk, almost up to the disk photosphere. We then obtain a self-consistent analytical solution for the vertical thermal structure of the disk, governed by the balance between the heating by MRI turbulence and the cooling by radiative diffusion. Next, we argue that the coronal power fraction is determined by the competition between the Parker instability, viewed as a parasitic instability feeding off of MRI channel flows, and other parasitic instabilities. We show that the Parker instability inevitably becomes important near the disk surface, leading to a certain lower limit on the coronal power. While most of the analysis in this paper focuses on the case of a disk threaded by an externally imposed vertical magnetic field, we also discuss the zero net flux case, in which the magnetic field is produced by the MRI dynamo itself, and show that most of our arguments and conclusions should be valid in this case as well.
Vaidyaraman, S.; Lackey, W.J.; Agrawal, P.K.; Freeman, G.B.; Langman, M.D.
1995-10-01
Carbon fiber-carbon matrix composites were fabricated using the forced flow-thermal gradient chemical vapor infiltration (FCVI) process. Preforms were prepared by stacking 40 layers of plain weave carbon cloth in a graphite holder. The preforms were infiltrated using propylene, propane, and methane. The present work showed that the FCVI process is well suited for fabricating carbon-carbon composites; without optimization of the process, the authors have achieved uniform and thorough densification. Composites with porosities as low as 7% were fabricated in 8--12 h. The highest deposition rate obtained in the present study was {approximately}3 {micro}m/h which is more than an order of magnitude faster than the typical value of 0.1--0.25 {micro}m/h for the isothermal process. It was also found that the use of propylene and propane as reagents resulted in faster infiltration compared to methane.
Magnetic and thermal energy flow during disruptions in DIII-D
Hyatt, A.W.; Lee, R.L.; Humphreys, D.A.; Kellman, A.G.; Taylor, P.L.; Cuthbertson, J.W.; Lasnier, C.J.
1996-07-01
The authors present results from disruption experiments where they measure magnetic energy flow across a closed surface surrounding the plasma using a Poynting flux analysis to measure the electromagnetic power, bolometers to measure radiation power and IR scanners to measure radiation and particle heat conduction to the divertor. The initial and final stored energies within the volume are found using the full equilibrium reconstruction code EFIT. From this analysis they calculate an energy balance and find that they can account for all energy deposited on the first wall and the divertor to within about 10%.
Flux-driven simulations of turbulence collapse
Park, G. Y.; Kim, S. S.; Jhang, Hogun; Rhee, T.; Diamond, P. H.; Xu, X. Q.
2015-03-15
Using three-dimensional nonlinear simulations of tokamak turbulence, we show that an edge transport barrier (ETB) forms naturally once input power exceeds a threshold value. Profiles, turbulence-driven flows, and neoclassical coefficients are evolved self-consistently. A slow power ramp-up simulation shows that ETB transition is triggered by the turbulence-driven flows via an intermediate phase which involves coherent oscillation of turbulence intensity and EĆB flow shear. A novel observation of the evolution is that the turbulence collapses and the ETB transition begins when R{sub T}ā>ā1 at tā=āt{sub R} (R{sub T}: normalized Reynolds power), while the conventional transition criterion (Ļ{sub EĆB}>Ī³{sub lin} where Ļ{sub EĆB} denotes mean flow shear) is satisfied only after tā=āt{sub C} (ā>t{sub R}), when the mean flow shear grows due to positive feedback.
Gyrokinetic Simulation of Global Turbulent Transport Properties in Tokamak Experiments
Wang, W.X.; Lin, Z.; Tang, W.M.; Lee, W.W.; Ethier, S.; Lewandowski, J.L.V.; Rewoldt, G.; Hahm, T.S.; Manickam, J.
2006-01-01
A general geometry gyro-kinetic model for particle simulation of plasma turbulence in tokamak experiments is described. It incorporates the comprehensive influence of noncircular cross section, realistic plasma profiles, plasma rotation, neoclassical (equilibrium) electric fields, and Coulomb collisions. An interesting result of global turbulence development in a shaped tokamak plasma is presented with regard to nonlinear turbulence spreading into the linearly stable region. The mutual interaction between turbulence and zonal flows in collisionless plasmas is studied with a focus on identifying possible nonlinear saturation mechanisms for zonal flows. A bursting temporal behavior with a period longer than the geodesic acoustic oscillation period is observed even in a collisionless system. Our simulation results suggest that the zonal flows can drive turbulence. However, this process is too weak to be an effective zonal flow saturation mechanism.
An h-adaptive finite element method for turbulent heat transfer
Carriington, David B [Los Alamos National Laboratory
2009-01-01
A two-equation turbulence closure model (k-{omega}) using an h-adaptive grid technique and finite element method (FEM) has been developed to simulate low Mach flow and heat transfer. These flows are applicable to many flows in engineering and environmental sciences. Of particular interest in the engineering modeling areas are: combustion, solidification, and heat exchanger design. Flows for indoor air quality modeling and atmospheric pollution transport are typical types of environmental flows modeled with this method. The numerical method is based on a hybrid finite element model using an equal-order projection process. The model includes thermal and species transport, localized mesh refinement (h-adaptive) and Petrov-Galerkin weighting for the stabilizing the advection. This work develops the continuum model of a two-equation turbulence closure method. The fractional step solution method is stated along with the h-adaptive grid method (Carrington and Pepper, 2002). Solutions are presented for 2d flow over a backward-facing step.
Reaction and diffusion in turbulent combustion
Pope, S.B.
1993-12-01
The motivation for this project is the need to obtain a better quantitative understanding of the technologically-important phenomenon of turbulent combustion. In nearly all applications in which fuel is burned-for example, fossil-fuel power plants, furnaces, gas-turbines and internal-combustion engines-the combustion takes place in a turbulent flow. Designers continually demand more quantitative information about this phenomenon-in the form of turbulent combustion models-so that they can design equipment with increased efficiency and decreased environmental impact. For some time the PI has been developing a class of turbulent combustion models known as PDF methods. These methods have the important virtue that both convection and reaction can be treated without turbulence-modelling assumptions. However, a mixing model is required to account for the effects of molecular diffusion. Currently, the available mixing models are known to have some significant defects. The major motivation of the project is to seek a better understanding of molecular diffusion in turbulent reactive flows, and hence to develop a better mixing model.
A Multi-Year Plan for Enhancing Turbulence Modeling in Hydra-TH Revised and Updated Version 2.0
Smith, Thomas M.; Berndt, Markus; Baglietto, Emilio; Magolan, Ben
2015-10-01
The purpose of this report is to document a multi-year plan for enhancing turbulence modeling in Hydra-TH for the Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Hydra-TH is being developed to the meet the high- fidelity, high-Reynolds number CFD based thermal hydraulic simulation needs of the program. This work is being conducted within the thermal hydraulics methods (THM) focus area. This report is an extension of THM CASL milestone L3:THM.CFD.P10.02 [33] (March, 2015) and picks up where it left off. It will also serve to meet the requirements of CASL THM level three milestone, L3:THM.CFD.P11.04, scheduled for completion September 30, 2015. The objectives of this plan will be met by: maturation of recently added turbulence models, strategic design/development of new models and systematic and rigorous testing of existing and new models and model extensions. While multi-phase turbulent flow simulations are important to the program, only single-phase modeling will be considered in this report. Large Eddy Simulation (LES) is also an important modeling methodology. However, at least in the first year, the focus is on steady-state Reynolds Averaged Navier-Stokes (RANS) turbulence modeling.
Maeyama, S. Nakata, M.; Miyato, N.; Yagi, M.; Ishizawa, A.; Watanabe, T.-H.; Idomura, Y.
2014-05-15
Electromagnetic turbulence driven by kinetic ballooning modes (KBMs) in high-? plasma is investigated based on the local gyrokinetic model. Analysis of turbulent fluxes, norms, and phases of fluctuations shows that KBM turbulence gives narrower spectra and smaller phase factors than those in ion-temperature-gradient (ITG)-driven turbulence. This leads to the smaller transport fluxes in KBM turbulence than those in ITG turbulence even when they have similar linear growth rates. From the analysis of the entropy balance relation, it is found that the entropy transfer from ions to electrons through the field-particle interactions mainly drives electron perturbations, which creates radial twisted modes by rapid parallel motions of electrons in a sheared magnetic geometry. The nonlinear coupling between the dominant unstable mode and its twisted modes is important for the saturation of KBM turbulence, in contrast to the importance of zonal flow shearing in ITG turbulence. The coupling depends on the flux-tube domain with the one-poloidal-turn parallel length and on the torus periodicity constraint.
Large Eddy Simulation of Stable Boundary Layer Turbulent Processes in Complex Terrain
Eric D. Skyllingstad
2005-01-26
Research was performed using a turbulence boundary layer model to study the behavior of cold, dense flows in regions of complex terrain. Results show that flows develop a balance between turbulent entrainment of warm ambient air and dense, cold air created by surface cooling. Flow depth and strength is a function of downslope distance, slope angle and angle changes, and the ambient air temperature.
Large eddy simulation of mixing between hot and cold sodium flows - comparison with experiments
Simoneau, J.P.; Noe, H.; Menant, B.
1995-09-01
The large eddy simulation is becoming a potential powerful tool for the calculation of turbulent flows. In nuclear liquid metal cooled fast reactors, the knowledge of the turbulence characteristics is of great interest for the prediction and the analysis of thermal stripping phenomena. The objective of this paper is to give a contribution in the evaluation of the large eddy simulation technique is an individual case. The problem chosen is the case of the mixing between hot and cold sodium flows. The computations are compared with available sodium tests. This study shows acceptable qualitative results but the simple model used is not able to predict the turbulence characteristics. More complex models including larger domains around the fluctuating zone and fluctuating boundary conditions could be necessary. Validation works are continuing.
A Two-length Scale Turbulence Model for Single-phase Multi-fluid Mixing
Schwarzkopf, J. D.; Livescu, D.; Baltzer, J. R.; Gore, R. A.; Ristorcelli, J. R.
2015-09-08
A two-length scale, second moment turbulence model (Reynolds averaged Navier-Stokes, RANS) is proposed to capture a wide variety of single-phase flows, spanning from incompressible flows with single fluids and mixtures of different density fluids (variable density flows) to flows over shock waves. The two-length scale model was developed to address an inconsistency present in the single-length scale models, e.g. the inability to match both variable density homogeneous Rayleigh-Taylor turbulence and Rayleigh-Taylor induced turbulence, as well as the inability to match both homogeneous shear and free shear flows. The two-length scale model focuses on separating the decay and transport length scales, as the two physical processes are generally different in inhomogeneous turbulence. This allows reasonable comparisons with statistics and spreading rates over such a wide range of turbulent flows using a common set of model coefficients. The specific canonical flows considered for calibrating the model include homogeneous shear, single-phase incompressible shear driven turbulence, variable density homogeneous Rayleigh-Taylor turbulence, Rayleigh-Taylor induced turbulence, and shocked isotropic turbulence. The second moment model shows to compare reasonably well with direct numerical simulations (DNS), experiments, and theory in most cases. The model was then applied to variable density shear layer and shock tube data and shows to be in reasonable agreement with DNS and experiments. Additionally, the importance of using DNS to calibrate and assess RANS type turbulence models is highlighted.
Greif, Ralph (University of California, Berkeley, CA); Evans, Gregory Herbert; Kearney, Sean Patrick (Sandia National Laboratories, Albuquerque, NM); Laskowski, Gregory Michael
2006-02-01
Heat transfer to and from a circular cylinder in a cross-flow of water at low Reynolds number was studied both experimentally and numerically. The experiments were carried out in a high aspect ratio water channel. The test section inflow temperature and velocity, channel lower surface temperature and cylinder surface temperature were controlled to yield either laminar or turbulent flow for a desired Richardson number. When the lower surface was unheated, the temperatures of the lower surface and water upstream of the cylinder were maintained approximately equal and the flow was laminar. When the lower surface was heated, turbulence intensities as high as 20% were measured several cylinder diameters upstream of the cylinder due to turbulent thermal plumes produced by heating the lower surface. Variable property, two-dimensional simulations were undertaken using a variant of the u{sup 2}-f turbulence model with buoyancy production of turbulence accounted for by a simple gradient diffusion model. Predicted and measured heat flux distributions around the cylinder are compared for values of the Richardson number, Gr{sub d}/Re{sub d}{sup 2} from 0.3 to 9.3. For laminar flow, the predicted and measured heat flux results agreed to within the experimental uncertainty. When the lower surface was heated, and the flow was turbulent, there was qualitative agreement between predicted and measured heat flux distributions around the cylinder. However the predicted spatially averaged Nusselt number was from 37% to 53% larger than the measured spatially averaged Nusselt number. Additionally, spatially averaged Nusselt numbers are compared to correlations in the literature for mixed convection heat transfer to/from cylinders in cross-flow. The results presented here are larger than the correlation values. This is believed to be due to the effects of buoyancy-induced turbulence resulting from heating the lower surface and the proximity of the cylinder to that surface.
The deterministic chaos and random noise in turbulent jet
Yao, Tian-Liang; Liu, Hai-Feng Xu, Jian-Liang; Li, Wei-Feng
2014-06-01
A turbulent flow is usually treated as a superposition of coherent structure and incoherent turbulence. In this paper, the largest Lyapunov exponent and the random noise in the near field of round jet and plane jet are estimated with our previously proposed method of chaotic time series analysis [T. L. Yao, et al., Chaos 22, 033102 (2012)]. The results show that the largest Lyapunov exponents of the round jet and plane jet are in direct proportion to the reciprocal of the integral time scale of turbulence, which is in accordance with the results of the dimensional analysis, and the proportionality coefficients are equal. In addition, the random noise of the round jet and plane jet has the same linear relation with the Kolmogorov velocity scale of turbulence. As a result, the random noise may well be from the incoherent disturbance in turbulence, and the coherent structure in turbulence may well follow the rule of chaotic motion.
Turbulent natural and mixed convection along a vertical plate
Abu-Mulaweh, H.I.; Armaly, B.F.; Chen, T.S.; Zhao, J.Z.
1997-07-01
Measurements of turbulent boundary-layer air flow in natural and mixed convection adjacent to an isothermal vertical flat plate are reported. Laser-Doppler velocimeter and cold wire anemometer were used, respectively, to measure simultaneously the mean turbulent velocity and temperature distributions were measured for a temperature difference, {Delta}T, of 30 C between the heated wall and the free stream air at a fixed location x = 3 m (with a corresponding Grashof number Gr{sub x} = 8.55 x 10{sup 10}), and for a range of free stream velocities 0 m/s {le} U{sub {infinity} } {le} 0.41 m/s. The effect of small free stream velocity on the turbulent natural convection is examined. These results reveal that the introduction of small free stream velocity on turbulent natural convection flow suppresses turbulence and decreases the heat transfer rate from the heated wall.
Radiosonde measurements of turbulence
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Turbulence detection on aerial platforms using orientation sensors R. Giles Harrison, Robin J. Hogan, George W. Rogers, Alyssa M. Heath and Keri A. Nicoll Department of Meteorology University of Reading, UK r.g.harrison@reading.ac.uk 2 "Bumpiness" (or turbulence), still remains an aircraft hazard, even at cruising altitudes How can in-situ direct position and orientation sensing can be used to detect turbulence affecting aerial platforms? * Investigated using balloon platforms,
Chatterjee, Rupa; Srivastava, Dinesh K.
2009-02-15
We calculate the elliptic flow of thermal photons from Au+Au collisions at RHIC energies for a range of values for the formation time {tau}{sub 0} but a fixed entropy (or particle rapidity density). The results are found to be quite sensitive to {tau}{sub 0}. The value of v{sub 2} for photons decreases as {tau}{sub 0} decreases and admits a larger contribution from the quark gluon plasma phase, which has a smaller v{sub 2}. The elliptic flow coefficient for hadrons, however, is only marginally dependent on {tau}{sub 0}.
Sabli, Nordin; Talib, Zainal Abidin; Yunus, Wan Mahmood Mat; Zainal, Zulkarnain; Hilal, Hikmat S.; Fujii, Masatoshi
2014-03-05
This work describes a new technique to enhance photoresponse of metal chalcogenide-based semiconductor film electrodes deposited by thermal vacuum evaporation under argon gas flow from synthesized Cu{sub 2}SnSe{sub 3} sources. SnSe formation with Cu-doped was obtained under higher argon gas flow rate (V{sub A} = 25 cm{sup 3}/min). Higher value of photoresponse was observed for films deposited under V{sub A} = 25 cm{sup 3}/min which was 9.1%. This finding indicates that Cu atoms inside the SnSe film were important to increase carrier concentrations that promote higher photoresponse.
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
... Furthermore, combined spontaneous Raman spectroscopy, Rayleigh scatter, and OHCO-laser-induced fluorescence measurements have been used to quantify in situ turbulent flame mixture ...
NO concentration imaging in turbulent nonpremixed flames
Schefer, R.W.
1993-12-01
The importance of NO as a pollutant species is well known. An understanding of the formation characteristics of NO in turbulent hydrocarbon flames is important to both the desired reduction of pollutant emissions and the validation of proposed models for turbulent reacting flows. Of particular interest is the relationship between NO formation and the local flame zone, in which the fuel is oxidized and primary heat release occurs. Planar imaging of NO provides the multipoint statistics needed to relate NO formation to the both the flame zone and the local turbulence characteristics. Planar imaging of NO has been demonstrated in turbulent flames where NO was seeded into the flow at high concentrations (2000 ppm) to determine the gas temperature distribution. The NO concentrations in these experiments were significantly higher than those expected in typical hydrocarbon-air flames, which require a much lower detectability limit for NO measurements. An imaging technique based on laser-induced fluorescence with sufficient sensitivity to study the NO formation mechanism in the stabilization region of turbulent lifted-jet methane flames.
Gyrokinetic Particle Simulation of Turbulent Transport in Burning Plasmas (GPS - TTBP) Final Report
Chame, Jacqueline
2011-05-27
The goal of this project is the development of the Gyrokinetic Toroidal Code (GTC) Framework and its applications to problems related to the physics of turbulence and turbulent transport in tokamaks,. The project involves physics studies, code development, noise effect mitigation, supporting computer science efforts, diagnostics and advanced visualizations, verification and validation. Its main scientific themes are mesoscale dynamics and non-locality effects on transport, the physics of secondary structures such as zonal flows, and strongly coherent wave-particle interaction phenomena at magnetic precession resonances. Special emphasis is placed on the implications of these themes for rho-star and current scalings and for the turbulent transport of momentum. GTC-TTBP also explores applications to electron thermal transport, particle transport; ITB formation and cross-cuts such as edge-core coupling, interaction of energetic particles with turbulence and neoclassical tearing mode trigger dynamics. Code development focuses on major initiatives in the development of full-f formulations and the capacity to simulate flux-driven transport. In addition to the full-f -formulation, the project includes the development of numerical collision models and methods for coarse graining in phase space. Verification is pursued by linear stability study comparisons with the FULL and HD7 codes and by benchmarking with the GKV, GYSELA and other gyrokinetic simulation codes. Validation of gyrokinetic models of ion and electron thermal transport is pursed by systematic stressing comparisons with fluctuation and transport data from the DIII-D and NSTX tokamaks. The physics and code development research programs are supported by complementary efforts in computer sciences, high performance computing, and data management.
Bursting frequency prediction in turbulent boundary layers
LIOU,WILLIAM W.; FANG,YICHUNG
2000-02-01
The frequencies of the bursting events associated with the streamwise coherent structures of spatially developing incompressible turbulent boundary layers were predicted using global numerical solution of the Orr-Sommerfeld and the vertical vorticity equations of hydrodynamic stability problems. The structures were modeled as wavelike disturbances associated with the turbulent mean flow. The global method developed here involves the use of second and fourth order accurate finite difference formula for the differential equations as well as the boundary conditions. An automated prediction tool, BURFIT, was developed. The predicted resonance frequencies were found to agree very well with previous results using a local shooting technique and measured data.
Sridharan, Kumar; Anderson, Mark; Allen, Todd; Corradini, Michael
2012-01-30
The goal of this NERI project was to perform research on high temperature fluoride and chloride molten salts towards the long-term goal of using these salts for transferring process heat from high temperature nuclear reactor to operation of hydrogen production and chemical plants. Specifically, the research focuses on corrosion of materials in molten salts, which continues to be one of the most significant challenges in molten salts systems. Based on the earlier work performed at ORNL on salt properties for heat transfer applications, a eutectic fluoride salt FLiNaK (46.5% LiF-11.5%NaF-42.0%KF, mol.%) and a eutectic chloride salt (32%MgCl2-68%KCl, mole %) were selected for this study. Several high temperature candidate Fe-Ni-Cr and Ni-Cr alloys: Hastelloy-N, Hastelloy-X, Haynes-230, Inconel-617, and Incoloy-800H, were exposed to molten FLiNaK with the goal of understanding corrosion mechanisms and ranking these alloys for their suitability for molten fluoride salt heat exchanger and thermal storage applications. The tests were performed at 850C for 500 h in sealed graphite crucibles under an argon cover gas. Corrosion was noted to occur predominantly from dealloying of Cr from the alloys, an effect that was particularly pronounced at the grain boundaries Alloy weight-loss due to molten fluoride salt exposure correlated with the initial Cr-content of the alloys, and was consistent with the Cr-content measured in the salts after corrosion tests. The alloys weight-loss was also found to correlate to the concentration of carbon present for the nominally 20% Cr containing alloys, due to the formation of chromium carbide phases at the grain boundaries. Experiments involving molten salt exposures of Incoloy-800H in Incoloy-800H crucibles under an argon cover gas showed a significantly lower corrosion for this alloy than when tested in a graphite crucible. Graphite significantly accelerated alloy corrosion due to the reduction of Cr from solution by graphite and formation
KOHURT, P. , KALINICHENKO, S.D.; KROSHILIN, A.E.; KROSHILIN, V.E.; SMIRNOV, A.V.
2006-06-04
In this paper we present recent results of the application of the thermal-hydraulic code BAGIRA to the analysis of complex two-phase flows in nuclear power plants primary loops. In particular, we performed benchmark numerical simulation of an integral LOCA experiment performed on a test facility modeling the primary circuit of VVER-1000. In addition, we have also analyzed the flow patterns in the VVER-1000 steam generator vessel for stationary and transient operation regimes. For both of these experiments we have compared the numerical results with measured data. Finally, we demonstrate the capabilities of BAGIRA by modeling a hypothetical severe accident for a VVER-1000 type nuclear reactor. The numerical analysis, which modeled all stages of the hypothetical severe accident up to the complete ablation of the reactor cavity bottom, shows the importance of multi-dimensional flow effects.
Wang, Y.; Pan, J; Mao, Y; Li, Z; Li, L; Hsiao, B
2010-01-01
The present Article reports the relationships between molecular orientation, formation, and spatial distribution of {gamma}-crystals in metallocene-made isotactic polypropylene (m-iPP) samples prepared by two industrial processes: conventional injection molding (CIM) and oscillatory shear injection molding (OSIM), in which combined thermal and flow fields typically exist. In particular, spatial distributions of crystallinity, fraction of {gamma}-crystal (f{gamma}) with respect to {alpha}-crystal, and lamella-branched shish-kebab structure in the shaped samples were characterized by synchrotron two-dimensional (2D) wide-angle X-ray diffraction (WAXD) and small-angle X-ray scattering (SAXS) techniques. The results showed that the crystallinity in any given region of OSIM samples was always higher than that of CIM samples. The value of f{gamma} increased monotonously from skin to core in CIM samples, whereas the corresponding f{gamma} increased nonmonotonically in OSIM samples. The spatial distribution of {gamma}-crystal in OSIM samples can be explained by the epitaxial arrangement between {gamma}- and {alpha}-crystal in a lamella-branched shish-kebab structure. In the proposed model, the parent lamellae of {alpha}-crystal provide secondary nucleation sites for daughter lamellae of {alpha}-crystal and {gamma}-crystal, and the different content of parent lamellae results in varying amounts of {gamma}-crystal. In OSIM samples, the smallest parent-daughter ratio ([R] = 1.38) in the core region led to the lowest fraction of {gamma}-crystal (0.57), but relatively higher {gamma}-crystal content (0.69) at 600 and 1200 {micro}m depth of the samples (corresponding to [R] of 4.5 and 5.8, respectively). This is consistent with the proposed model where more parent lamellae provide more nucleation sites for crystallization, thus resulting in higher content of {gamma}-crystal. The melting behavior of CIM and OSIM samples was studied by differential scanning calorimetery (DSC). The
Toward the Theory of Turbulence in Magnetized Plasmas
Boldyrev, Stanislav
2013-07-26
The goal of the project was to develop a theory of turbulence in magnetized plasmas at large scales, that is, scales larger than the characteristic plasma microscales (ion gyroscale, ion inertial scale, etc.). Collisions of counter-propagating Alfven packets govern the turbulent cascade of energy toward small scales. It has been established that such an energy cascade is intrinsically anisotropic, in that it predominantly supplies energy to the modes with mostly field-perpendicular wave numbers. The resulting energy spectrum of MHD turbulence, and the structure of the fluctuations were studied both analytically and numerically. A new parallel numerical code was developed for simulating reduced MHD equations driven by an external force. The numerical setting was proposed, where the spectral properties of the force could be varied in order to simulate either strong or weak turbulent regimes. It has been found both analytically and numerically that weak MHD turbulence spontaneously generates a condensate, that is, concentration of magnetic and kinetic energy at small k{sub {parallel}}. A related topic that was addressed in the project is turbulent dynamo action, that is, generation of magnetic field in a turbulent flow. We were specifically concentrated on the generation of large-scale magnetic field compared to the scales of the turbulent velocity field. We investigate magnetic field amplification in a turbulent velocity field with nonzero helicity, in the framework of the kinematic Kazantsev-Kraichnan model.
RELAXATION PROCESSES IN SOLAR WIND TURBULENCE
Servidio, S.; Carbone, V.; Gurgiolo, C.; Goldstein, M. L.
2014-07-10
Based on global conservation principles, magnetohydrodynamic (MHD) relaxation theory predicts the existence of several equilibria, such as the Taylor state or global dynamic alignment. These states are generally viewed as very long-time and large-scale equilibria, which emerge only after the termination of the turbulent cascade. As suggested by hydrodynamics and by recent MHD numerical simulations, relaxation processes can occur during the turbulent cascade that will manifest themselves as local patches of equilibrium-like configurations. Using multi-spacecraft analysis techniques in conjunction with Cluster data, we compute the current density and flow vorticity and for the first time demonstrate that these localized relaxation events are observed in the solar wind. Such events have important consequences for the statistics of plasma turbulence.
Core density turbulence in the HSX Stellarator
Deng, C. B.; Brower, D. L.; Anderson, D. T.; Anderson, F. S. B.; Briesemeister, Alexis R.; Likin, K. M.
2015-10-23
Broadband turbulent density fluctuations are explored in the helically symmetric stellarator experiment (HSX) by investigating changes related to plasma heating power and location. No fluctuation response is observed to occur with large changes in electron temperature and its gradient, thereby eliminating temperature gradient as a driving mechanism. Instead, measurements reveal that density turbulence varies inversely with electron density scale length. This response is consistent with density gradient drive as one might expect for trapped electron mode (TEM) turbulence. In general, the plasma stored energy and particle confinement are higher for discharges with reduced fluctuations in the plasma core. When the density fluctuation amplitude is reduced, increased plasma rotation is also evident suggesting a role is being played by intrinsic plasma flow.
Energy Science and Technology Software Center (OSTI)
2008-01-25
BOUT is a parallelized 3D nonlocal electromagnetic turbulence code. The principal calculations are the boundary plasma turbulence in a realistic magnetic geometry. BOUT uses fluid Braginskii equations for plasma vorticity, density, electron and ion temperature and Parallel mementum. With sources added in the core-edge region and sinks in the scrape-off-layer (SOL), BOUT follows the self-consistent profile evolution together with turbulence. BOUT also includes coupling to a magnetohyfrodynamic equlibrium (EFIT package) and a two-dimensional hydrodynamic edgemoreĀ Ā» transport model (UEDGE package).Ā«Ā less
Observations of Edge Turbulence
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Edge Turbulence near the X-point of Alcator C-Mod APS-2007 (1) J.L. Terry, S.J. Zweben*, B. LaBombard, I. Cziegler, O. Grulke + , D.P. Stotler* MIT - Plasma Science and Fusion...
Characterization of Fuego for laminar and turbulent natural convection heat transfer.
Francis, Nicholas Donald, Jr. (,; .)
2005-08-01
A computational fluid dynamics (CFD) analysis is conducted for internal natural convection heat transfer using the low Mach number code Fuego. The flow conditions under investigation are primarily laminar, transitional, or low-intensity level turbulent flows. In the case of turbulent boundary layers at low-level turbulence or transitional Reynolds numbers, the use of standard wall functions no longer applies, in general, for wall-bounded flows. One must integrate all the way to the wall in order to account for gradients in the dependent variables in the viscous sublayer. Fuego provides two turbulence models in which resolution of the near-wall region is appropriate. These models are the v2-f turbulence model and a Launder-Sharma, low-Reynolds number turbulence model. Two standard geometries are considered: the annulus formed between horizontal concentric cylinders and a square enclosure. Each geometry emphasizes wall shear flow and complexities associated with turbulent or near turbulent boundary layers in contact with a motionless core fluid. Overall, the Fuego simulations for both laminar and turbulent flows compared well to measured data, for both geometries under investigation, and to a widely accepted commercial CFD code (FLUENT).
Scaled Experimental Modeling of VHTR Plenum Flows
ICONE 15
2007-04-01
Abstract The Very High Temperature Reactor (VHTR) is the leading candidate for the Next Generation Nuclear Power (NGNP) Project in the U.S. which has the goal of demonstrating the production of emissions free electricity and hydrogen by 2015. Various scaled heated gas and water flow facilities were investigated for modeling VHTR upper and lower plenum flows during the decay heat portion of a pressurized conduction-cooldown scenario and for modeling thermal mixing and stratification (thermal striping) in the lower plenum during normal operation. It was concluded, based on phenomena scaling and instrumentation and other practical considerations, that a heated water flow scale model facility is preferable to a heated gas flow facility and to unheated facilities which use fluids with ranges of density to simulate the density effect of heating. For a heated water flow lower plenum model, both the Richardson numbers and Reynolds numbers may be approximately matched for conduction-cooldown natural circulation conditions. Thermal mixing during normal operation may be simulated but at lower, but still fully turbulent, Reynolds numbers than in the prototype. Natural circulation flows in the upper plenum may also be simulated in a separate heated water flow facility that uses the same plumbing as the lower plenum model. However, Reynolds number scaling distortions will occur at matching Richardson numbers due primarily to the necessity of using a reduced number of channels connected to the plenum than in the prototype (which has approximately 11,000 core channels connected to the upper plenum) in an otherwise geometrically scaled model. Experiments conducted in either or both facilities will meet the objectives of providing benchmark data for the validation of codes proposed for NGNP designs and safety studies, as well as providing a better understanding of the complex flow phenomena in the plenums.
The effect of solids concentration on self-induced turbulence
Kenning, V.; Crowe, C.T.
1994-12-31
A model to predict the turbulence intensity due to the solid particles in a simple flow has been developed. The flow is one in which all the turbulence is due to the presence of the particles. The model accounts for energy input through the loss of potential energy by the solid particles and energy loss by turbulent dissipation. The predictions are compared for various solids concentrations and particle sizes. The peak turbulence intensity is seen to b3e reached faster for higher solid concentrations. The peak is also higher for higher concentrations. In all cases, a peak value of turbulent intensity is reached if the supply of particles is maintained. The case in which the supply of particles are no longer available to supply the turbulence with energy. When normalized by the peak value, the turbulence was seen to decay more rapidly for higher concentrations of solid particles in the present model. An experimental study will be conducted to compare with the current model.
Numerical Simulations of Small Non-spherical Particles in Turbulence |
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Argonne Leadership Computing Facility Numerical Simulations of Small Non-spherical Particles in Turbulence Event Sponsor: Mathematics and Computer Science Division LANS Seminar Start Date: Aug 31 2016 - 3:00pm Building/Room: Building 240/Room 1404-1405 Location: Argonne National Laboratory Speaker(s): Nimish Pujara Speaker(s) Title: UC Berkeley Motivated by the ubiquity of natural particles in turbulent flows in the natural environment as well as in many industrial processes, we investigate
Direct Numerical Simulation of Compressible, Turbulent Flow ...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
The computational mesh for this direct numerical simulation was over 33 billion cells, and was run on up to 102,400 cores under a DoD HPCMP Frontier Project. Nicholas Bisek and ...
Feizollahi, F.; Quapp, W.J.; Hempill, H.G.; Groffie, F.J.
1994-07-01
This document accompanies a full report which describes the testing and evaluation of ten different methods for incinerating mixed low-level radioactive wastes. It consists of flowsheets and diagrams of a rotary kiln, pyrolysis methods, a plasma furnace, a fixed hearth, and thermal desorption methods.
James E. O'Brien; Piyush Sabharwall; SuJong Yoon
2001-09-01
A new high-temperature multi-fluid, multi-loop test facility for advanced nuclear applications is under development at the Idaho National Laboratory. The facility will include three flow loops: high-temperature helium, molten salt, and steam/water. Molten salts have been identified as excellent candidate heat transport fluids for primary or secondary coolant loops, supporting advanced high temperature and small modular reactors (SMRs). Details of some of the design aspects and challenges of this facility, which is currently in the conceptual design phase, are discussed. A preliminary design configuration will be presented, with the required characteristics of the various components. The loop will utilize advanced high-temperature compact printed-circuit heat exchangers (PCHEs) operating at prototypic intermediate heat exchanger (IHX) conditions. The initial configuration will include a high-temperature (750Ā°C), high-pressure (7 MPa) helium loop thermally integrated with a molten fluoride salt (KF-ZrF4) flow loop operating at low pressure (0.2 MPa) at a temperature of ~450Ā°C. Experiment design challenges include identification of suitable materials and components that will withstand the required loop operating conditions. Corrosion and high temperature creep behavior are major considerations. The facility will include a thermal energy storage capability designed to support scaled process heat delivery for a variety of hybrid energy systems and grid stabilization strategies. Experimental results obtained from this research will also provide important data for code ve
Particle dispersion in homogeneous turbulence using the one-dimensional turbulence model
Sun, Guangyuan; Lignell, David O.; Hewson, John C.; Gin, Craig R.
2014-10-09
Lagrangian particle dispersion is studied using the one-dimensional turbulence (ODT) model in homogeneous decaying turbulence configurations. The ODT model has been widely and successfully applied to a number of reacting and nonreacting flow configurations, but only limited application has been made to multiphase flows. We present a version of the particle implementation and interaction with the stochastic and instantaneous ODT eddy events. The model is characterized by comparison to experimental data of particle dispersion for a range of intrinsic particle time scales and body forces. Particle dispersion, velocity, and integral time scale results are presented. Moreover, the particle implementation introduces a single model parameter Ī² p , and sensitivity to this parameter and behavior of the model are discussed. Good agreement is found with experimental data and the ODT model is able to capture the particle inertial and trajectory crossing effects. Our results serve as a validation case of the multiphase implementations of ODT for extensions to other flow configurations.
Bell, Geoffrey C.; Feustel, Helmut E.; Dickerhoff, Darryl J.
2002-01-01
A fume hood is provided having an adequate level of safety while reducing the amount of air exhausted from the hood. A displacement flow fume hood works on the principal of a displacement flow which displaces the volume currently present in the hood using a push-pull system. The displacement flow includes a plurality of air supplies which provide fresh air, preferably having laminar flow, to the fume hood. The displacement flow fume hood also includes an air exhaust which pulls air from the work chamber in a minimally turbulent manner. As the displacement flow produces a substantially consistent and minimally turbulent flow in the hood, inconsistent flow patterns associated with contaminant escape from the hood are minimized. The displacement flow fume hood largely reduces the need to exhaust large amounts of air from the hood. It has been shown that exhaust air flow reductions of up to 70% are possible without a decrease in the hood's containment performance. The fume hood also includes a number of structural adaptations which facilitate consistent and minimally turbulent flow within a fume hood.
Rutland, Christopher J.
2009-04-26
The Terascale High-Fidelity Simulations of Turbulent Combustion (TSTC) project is a multi-university collaborative effort to develop a high-fidelity turbulent reacting flow simulation capability utilizing terascale, massively parallel computer technology. The main paradigm of the approach is direct numerical simulation (DNS) featuring the highest temporal and spatial accuracy, allowing quantitative observations of the fine-scale physics found in turbulent reacting flows as well as providing a useful tool for development of sub-models needed in device-level simulations. Under this component of the TSTC program the simulation code named S3D, developed and shared with coworkers at Sandia National Laboratories, has been enhanced with new numerical algorithms and physical models to provide predictive capabilities for turbulent liquid fuel spray dynamics. Major accomplishments include improved fundamental understanding of mixing and auto-ignition in multi-phase turbulent reactant mixtures and turbulent fuel injection spray jets.
Turbulence assessment at potential turbine sites
Daniels, A.
1996-12-31
As opposed to a fixed anemometer, the Tala kite is free to move in the air. The motion of the kite is not random, it moves with or against the speed gradient towards the center of passing turbulence events of higher or lower speeds thus allowing the kite to measure event maximum or minimum speed rather than the speed at some unknown distance from the event center like a fixed anemometer. This behavior is confirmed both by a theoretical aerodynamics analysis of the kite motion and by data from a field study where kite and hot film anemometer (HFA) events, defined by the rain flow count method, were compared with flap events on a rotating turbine blade. The HFAs simulated too few events lasting too long while the kites reproduced both the number of events and event periods remarkably close. It is concluded that the kite is the optimal tool for measuring turbulence at potential turbine sites. Kite turbulence can form the bases for economic return estimates and an example is given where less windy sites could be more economical than other more turbulent higher speed sites. 13 refs., 8 figs.
A spray-suppression model for turbulent combustion
DESJARDIN,PAUL E.; TIESZEN,SHELDON R.; GRITZO,LOUIS A.
2000-02-14
A spray-suppression model that captures the effects of liquid suppressant on a turbulent combusting flow is developed and applied to a turbulent diffusion flame with water spray suppression. The spray submodel is based on a stochastic separated flow approach that accounts for the transport and evaporation of liquid droplets. Flame extinguishment is accounted for by using a perfectly stirred reactor (PSR) submodel of turbulent combustion. PSR pre-calculations of flame extinction times are determined using CHEMKIN and are compared to local turbulent time scales of the flow to determine if local flame extinguishment has occurred. The PSR flame extinguishment and spray submodels are incorporated into Sandia's flow fire simulation code, VULCAN, and cases are run for the water spray suppression studies of McCaffrey for turbulent hydrogen-air jet diffusion flames. Predictions of flame temperature decrease and suppression efficiency are compared to experimental data as a function of water mass loading using three assumed values of drop sizes. The results show that the suppression efficiency is highly dependent on the initial droplet size for a given mass loading. A predicted optimal suppression efficiency was observed for the smallest class of droplets while the larger drops show increasing suppression efficiency with increasing mass loading for the range of mass loadings considered. Qualitative agreement to the experiment of suppression efficiency is encouraging, however quantitative agreement is limited due to the uncertainties in the boundary conditions of the experimental data for the water spray.
Implementation and Validation of the BHR Turbulence Model in the FLAG
Office of Scientific and Technical Information (OSTI)
Hydrocode (Technical Report) | SciTech Connect Technical Report: Implementation and Validation of the BHR Turbulence Model in the FLAG Hydrocode Citation Details In-Document Search Title: Implementation and Validation of the BHR Turbulence Model in the FLAG Hydrocode The BHR-2 turbulence model, developed at Los Alamos National Laboratory for variable density and compressible flows, is implemented in an Arbitrary Lagrangian-Eulerian hydrocode, FLAG. The BHR-2 formulation is discussed, with
Collisionless inter-species energy transfer and turbulent heating in drift wave turbulence
Zhao, L.; Diamond, P. H.
2012-08-15
We reconsider the classic problems of calculating 'turbulent heating' and collisionless inter-species transfer of energy in drift wave turbulence. These issues are of interest for low collisionality, electron heated plasmas, such as ITER, where collisionless energy transfer from electrons to ions is likely to be significant. From the wave Poynting theorem at steady state, a volume integral over an annulus r{sub 1}
A Two-length Scale Turbulence Model for Single-phase Multi-fluid Mixing
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Schwarzkopf, J. D.; Livescu, D.; Baltzer, J. R.; Gore, R. A.; Ristorcelli, J. R.
2015-09-08
A two-length scale, second moment turbulence model (Reynolds averaged Navier-Stokes, RANS) is proposed to capture a wide variety of single-phase flows, spanning from incompressible flows with single fluids and mixtures of different density fluids (variable density flows) to flows over shock waves. The two-length scale model was developed to address an inconsistency present in the single-length scale models, e.g. the inability to match both variable density homogeneous Rayleigh-Taylor turbulence and Rayleigh-Taylor induced turbulence, as well as the inability to match both homogeneous shear and free shear flows. The two-length scale model focuses on separating the decay and transport length scales,moreĀ Ā» as the two physical processes are generally different in inhomogeneous turbulence. This allows reasonable comparisons with statistics and spreading rates over such a wide range of turbulent flows using a common set of model coefficients. The specific canonical flows considered for calibrating the model include homogeneous shear, single-phase incompressible shear driven turbulence, variable density homogeneous Rayleigh-Taylor turbulence, Rayleigh-Taylor induced turbulence, and shocked isotropic turbulence. The second moment model shows to compare reasonably well with direct numerical simulations (DNS), experiments, and theory in most cases. The model was then applied to variable density shear layer and shock tube data and shows to be in reasonable agreement with DNS and experiments. Additionally, the importance of using DNS to calibrate and assess RANS type turbulence models is highlighted.Ā«Ā less
The interaction of high-speed turbulence with flames: Global properties and internal flame structure
Poludnenko, A.Y.; Oran, E.S. [Laboratory for Computational Physics and Fluid Dynamics, Naval Research Laboratory, Washington, DC 20375 (United States)
2010-05-15
We study the dynamics and properties of a turbulent flame, formed in the presence of subsonic, high-speed, homogeneous, isotropic Kolmogorov-type turbulence in an unconfined system. Direct numerical simulations are performed with Athena-RFX, a massively parallel, fully compressible, high-order, dimensionally unsplit, reactive flow code. A simplified reaction-diffusion model represents a stoichiometric H{sub 2}-air mixture. The system being modeled represents turbulent combustion with the Damkoehler number Da=0.05 and with the turbulent velocity at the energy injection scale 30 times larger than the laminar flame speed. The simulations show that flame interaction with high-speed turbulence forms a steadily propagating turbulent flame with a flame brush width approximately twice the energy injection scale and a speed four times the laminar flame speed. A method for reconstructing the internal flame structure is described and used to show that the turbulent flame consists of tightly folded flamelets. The reaction zone structure of these is virtually identical to that of the planar laminar flame, while the preheat zone is broadened by approximately a factor of two. Consequently, the system evolution represents turbulent combustion in the thin reaction zone regime. The turbulent cascade fails to penetrate the internal flame structure, and thus the action of small-scale turbulence is suppressed throughout most of the flame. Finally, our results suggest that for stoichiometric H{sub 2}-air mixtures, any substantial flame broadening by the action of turbulence cannot be expected in all subsonic regimes. (author)
Test particle study of ion transport in drift type turbulence
Vlad, M.; Spineanu, F.
2013-12-15
Ion transport regimes in drift type turbulence are determined in the frame of a realistic model for the turbulence spectrum based on numerical simulations. The model includes the drift of the potential with the effective diamagnetic velocity, turbulence anisotropy, and dominant waves. The effects of the zonal flow modes are also analyzed. A semi-analytical method that is able to describe trajectory stochastic trapping or eddying is used for obtaining the transport coefficients as function of the parameters of the turbulence. Analytical approximations of the transport coefficients are derived from the results. They show the transition from Bohm to gyro-Bohm scaling as plasma size increases in very good agreement with the numerical simulations.
Li, P.W.; Daisaka, H.; Kawaguchi, Y.; Yabe, A.; Hishida, K.; Maeda, M.
1999-07-01
The turbulent characteristics of a surfactant water solution in changing from drag-reducing flow to turbulent flow inside a two-dimensional smooth channel and in changing from turbulent flow to drag-reducing flow in the same channel with a mesh plug were investigated through LDV measurement in this study. The mesh plug was used to exert high shear stress to destroy micelle structures in the surfactant solution so that turbulence could be produced for better heat transfer. The two-component LDV system was installed on a movable platform, which could be moved streamwise of the flow to measure the two-dimensional velocity at different stations downstream from the mesh plug. The surfactant tested was Cetyltrimethyl ammonium chloride (C{sub 16}H{sub 33}N(CH{sub 3}){sub 3}Cl, abbreviated as CTAC). Local tap water was used as solvent and same weight concentration of sodium salicylate was used as the counter-ion material. The investigation of turbulent parameters for the drag-reducing flow with increasing Reynolds number showed that when the Reynolds number exceeded the drag-reducing region, the turbulent character was the same as that of water. The turbulent parameters of surfactant flow downstream the mesh plug showed that the high heat transfer region had the same turbulent intensity as that of water flow. As the critical Reynolds number was approached, it became easier to obtain such a turbulent region by mesh plug. In such cases, the mesh helped to create high wall shear stress and therefore to destroy the super-ordered structures of rod-like micelles for introducing turbulence. However, it was found that the turbulent intensities of the velocity gradually decreased to the same as those of drag-reducing flow downstream from the mesh because the mesh plug only produced a local high shear stress.
A turbulent transport network model in MULTIFLUX coupled with TOUGH2
Danko, G.; Bahrami, D.; Birkholzer, J.T.
2011-02-15
A new numerical method is described for the fully iterated, conjugate solution of two discrete submodels, involving (a) a transport network model for heat, moisture, and airflows in a high-permeability, air-filled cavity; and (b) a variably saturated fractured porous medium. The transport network submodel is an integrated-parameter, computational fluid dynamics solver, describing the thermal-hydrologic transport processes in the flow channel system of the cavity with laminar or turbulent flow and convective heat and mass transport, using MULTIFLUX. The porous medium submodel, using TOUGH2, is a solver for the heat and mass transport in the fractured rock mass. The new model solution extends the application fields of TOUGH2 by integrating it with turbulent flow and transport in a discrete flow network system. We present demonstrational results for a nuclear waste repository application at Yucca Mountain with the most realistic model assumptions and input parameters including the geometrical layout of the nuclear spent fuel and waste with variable heat load for the individual containers. The MULTIFLUX and TOUGH2 model elements are fully iterated, applying a programmed reprocessing of the Numerical Transport Code Functionalization model-element in an automated Outside Balance Iteration loop. The natural, convective airflow field and the heat and mass transport in a representative emplacement drift during postclosure are explicitly solved in the new model. The results demonstrate that the direction and magnitude of the air circulation patterns and all transport modes are strongly affected by the heat and moisture transport processes in the surrounding rock, justifying the need for a coupled, fully iterated model solution such as the one presented in the paper.
Stagnation Region Heat Transfer Augmentation at Very High Turbulence Levels
Ames, Forrest; Kingery, Joseph E.
2015-06-17
A database for stagnation region heat transfer has been extended to include heat transfer measurements acquired downstream from a new high intensity turbulence generator. This work was motivated by gas turbine industry heat transfer designers who deal with heat transfer environments with increasing Reynolds numbers and very high turbulence levels. The new mock aero-combustor turbulence generator produces turbulence levels which average 17.4%, which is 37% higher than the older turbulence generator. The increased level of turbulence is caused by the reduced contraction ratio from the liner to the exit. Heat transfer measurements were acquired on two large cylindrical leading edge test surfaces having a four to one range in leading edge diameter (40.64 cm and 10.16 cm). Gandvarapu and Ames [1] previously acquired heat transfer measurements for six turbulence conditions including three grid conditions, two lower turbulence aero-combustor conditions, and a low turbulence condition. The data are documented and tabulated for an eight to one range in Reynolds numbers for each test surface with Reynolds numbers ranging from 62,500 to 500,000 for the large leading edge and 15,625 to 125,000 for the smaller leading edge. The data show augmentation levels of up to 136% in the stagnation region for the large leading edge. This heat transfer rate is an increase over the previous aero-combustor turbulence generator which had augmentation levels up to 110%. Note, the rate of increase in heat transfer augmentation decreases for the large cylindrical leading edge inferring only a limited level of turbulence intensification in the stagnation region. The smaller cylindrical leading edge shows more consistency with earlier stagnation region heat transfer results correlated on the TRL (Turbulence, Reynolds number, Length scale) parameter. The downstream regions of both test surfaces continue to accelerate the flow but at a much lower rate than the leading edge. Bypass transition occurs
Simulations of drift resistive ballooning L-mode turbulence in the edge plasma of the DIII-D tokamak
Cohen, B. I.; Umansky, M. V.; Nevins, W. M.; Makowski, M. A.; Boedo, J. A.; Rudakov, D. L.; McKee, G. R.; Yan, Z.; Groebner, R. J.
2013-05-15
Results from simulations of electromagnetic drift-resistive ballooning turbulence for tokamak edge turbulence in realistic single-null geometry are reported. The calculations are undertaken with the BOUT three-dimensional fluid code that solves Braginskii-based fluid equations [X. Q. Xu and R. H. Cohen, Contrib. Plasma Phys. 36, 158 (1998)]. The simulation setup models L-mode edge plasma parameters in the actual magnetic geometry of the DIII-D tokamak [J. L. Luxon et al., Fusion Sci. Technol. 48, 807 (2002)]. The computations track the development of drift-resistive ballooning turbulence in the edge region to saturation. Fluctuation amplitudes, fluctuation spectra, and particle and thermal fluxes are compared to experimental data near the outer midplane from Langmuir probe and beam-emission-spectroscopy for a few well-characterized L-mode discharges in DIII-D. The simulations are comprised of a suite of runs in which the physics model is varied to include more fluid fields and physics terms. The simulations yield results for fluctuation amplitudes, correlation lengths, particle and energy fluxes, and diffusivities that agree with measurements within an order of magnitude and within factors of 2 or better for some of the data. The agreement of the simulations with the experimental measurements varies with respect to including more physics in the model equations within the suite of models investigated. The simulations show stabilizing effects of sheared E × B poloidal rotation (imposed zonal flow) and of lower edge electron temperature and density.
The selection of turbulence models for prediction of room airflow
Nielsen, P.V.
1998-10-01
The airflow in buildings involves a combination of many different flow elements. It is, therefore, difficult to find an adequate, all-round turbulence model covering all aspects. Consequently, it is appropriate and economical to choose turbulence models according to the situation that is to be predicted. This paper discusses the use of different turbulence models and their advantages in given situations. As an example, it is shown that a simple zero-equation model can be used for the prediction of special situations as flow with a low level of turbulence. A zero-equation model with compensation for room dimensions and velocity level also is discussed. A {kappa}-{epsilon} model expanded by damping functions is used to improve the prediction of the flow in a room ventilated by displacement ventilation. The damping functions especially take into account the turbulence level and the vertical temperature gradient. Low Reynolds number models (LNR models) are used to improve the prediction of evaporation-controlled emissions from building material, which is shown by an example. Finally, large eddy simulation (LES) of room airflow is discussed and demonstrated.
Developments in the Gyrofluid approach to tokamak turbulence simulations
Hammett, G.W.; Beer, M.A.; Dorland, W.; Smith, S.A.; Cowley, S.C.
1993-06-01
A status report is given on recent developments in the gyrofluid approach to simulating tokamak turbulence. ``Gyrofluid`` (or ``gyro-Landau fluid``) equations attempt to extend the range of validity of fluid equations to a more collisionless regime typical of tokamaks, by developing fluid models of important kinetic effects such as Landau-damping and gyro-orbit averaging. The fluid moments approach should converge if enough moments are kept, though this may require a large number of moments for some processes. Toroidal gyrofluid equations have been extended from 4 to 6 moments, and to include the {mu} {gradient}B magnetic mirroring force. An efficient field-line coordinate system for toroidal turbulence simulations (useful for both particle and fluid simulations) is presented. Nonlinear 3-D simulations of toroidal ITG-driven turbulence indicate that turbulence-generated sheared flows play. an important role in the development and saturation of the turbulence. There is a strong enhancement of the flows when the electrons are assumed adiabatic on each flux surface, which is partially offset by toroidal drift effects which reduce the flows.
PLASMA EMISSION BY WEAK TURBULENCE PROCESSES
Ziebell, L. F.; Gaelzer, R.; Yoon, P. H.; Pavan, J. E-mail: rudi.gaelzer@ufrgs.br E-mail: joel.pavan@ufpel.edu.br
2014-11-10
The plasma emission is the radiation mechanism responsible for solar type II and type III radio bursts. The first theory of plasma emission was put forth in the 1950s, but the rigorous demonstration of the process based upon first principles had been lacking. The present Letter reports the first complete numerical solution of electromagnetic weak turbulence equations. It is shown that the fundamental emission is dominant and unless the beam speed is substantially higher than the electron thermal speed, the harmonic emission is not likely to be generated. The present findings may be useful for validating reduced models and for interpreting particle-in-cell simulations.
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Turbulent Combustion - Sandia Energy Energy Search Icon Sandia Home Locations Contact Us Employee Locator Energy & Climate Secure & Sustainable Energy Future Stationary Power Energy Conversion Efficiency Solar Energy Wind Energy Water Power Supercritical CO2 Geothermal Natural Gas Safety, Security & Resilience of the Energy Infrastructure Energy Storage Nuclear Power & Engineering Grid Modernization Battery Testing Nuclear Energy Defense Waste Management Programs Advanced Nuclear
Identifying Turbulent Structures through Topological Segmentation
Bremer, Peer-Timo; Gruber, Andrea; Bennett, Janine C.; Gyulassy, Attila; Kolla, Hemanth; Chen, Jacqueline H.; Grout, Ray W.
2016-01-01
A new method of extracting vortical structures from a turbulent flow is proposed whereby topological segmentation of an indicator function scalar field is used to identify the regions of influence of the individual vortices. This addresses a long-standing challenge in vector field topological analysis: indicator functions commonly used produce a scalar field based on the local velocity vector field; reconstructing regions of influence for a particular structure requires selecting a threshold to define vortex extent. In practice, the same threshold is rarely meaningful throughout a given flow. By also considering the topology of the indicator field function, the characteristics of vortex strength and extent can be separated and the ambiguity in the choice of the threshold reduced. The proposed approach is able to identify several types of vortices observed in a jet in cross-flow configuration simultaneously where no single threshold value for a selection of common indicator functions appears able to identify all of these vortex types.
Convection causes enhanced magnetic turbulence in accretion disks in outburst
Hirose, Shigenobu; Blaes, Omer; Coleman, Matthew S. B.; Krolik, Julian H.; Sano, Takayoshi
2014-05-20
We present the results of local, vertically stratified, radiation magnetohydrodynamic (MHD) shearing box simulations of magneto-rotational instability (MRI) turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near 10{sup 4} K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio ?, possibly by generating vertical magnetic field that may seed the axisymmetric MRI, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the ?-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.
Development of one-equation transition/turbulence models
Edwards, J.R.; Roy, C.J.; Blottner, F.G.; Hassan, H.A.
2000-01-14
This paper reports on the development of a unified one-equation model for the prediction of transitional and turbulent flows. An eddy viscosity--transport equation for nonturbulent fluctuation growth based on that proposed by Warren and Hassan is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittency function based on the work of Dhawan and Narasimha. The model predicts both the onset and extent of transition. Low-speed test cases include transitional flow over a flat plate, a single element airfoil, and a multi-element airfoil in landing configuration. High-speed test cases include transitional Mach 3.5 flow over a 5{degree} cone and Mach 6 flow over a flared-cone configuration. Results are compared with experimental data, and the grid-dependence of selected predictions is analyzed.
Numerical calculation of two-phase turbulent jets
Saif, A.A.
1995-05-01
Two-phase turbulent round jets were numerically simulated using a multidimensional two-phase CFD code based on the two-fluid model. The turbulence phenomena were treated with the standard k-{epsilon} model. It was modified to take into account the additional dissipation of turbulent kinetic energy by the dispersed phase. Within the context of the two-fluid model it is more appropriate and physically justified to treat the diffusion by an interfacial force in the momentum equation. In this work, the diffusion force and the additional dissipation effect by the dispersed phase were modeled starting from the classical turbulent energy spectrum analysis. A cut-off frequency was proposed to decrease the dissipation effect by the dispersed phase when large size particles are introduced in the flow. The cut-off frequency combined with the bubble-induced turbulence effect allows for an increase in turbulence for large particles. Additional care was taken in choosing the right kind of experimental data from the literature so that a good separate effect test was possible for their models. The models predicted the experimental data very closely and they were general enough to predict extreme limit cases: water-bubble and air-droplet jets.
Turbulent Equipartition Theory of Toroidal Momentum Pinch
T.S. Hahm, P.H. Diamond, O.D. Gurcan, and G. Rewaldt
2008-01-31
The mode-independet part of magnetic curvature driven turbulent convective (TuroCo) pinch of the angular momentum density [Hahm et al., Phys. Plasmas 14,072302 (2007)] which was originally derived from the gyrokinetic equation, can be interpreted in terms of the turbulent equipartition (TEP) theory. It is shown that the previous results can be obtained from the local conservation of "magnetically weighted angular momentum density," nmi U|| R/B2, and its homogenization due to turbulent flows. It is also demonstrated that the magnetic curvature modification of the parallel acceleration in the nonlinear gyrokinetic equation in the laboratory frame, which was shown to be responsible for the TEP part of the TurCo pinch of angular momentum density in the previous work, is closely related to the Coriolis drift coupling to the perturbed electric field. In addition, the origin of the diffusive flux in the rotating frame is highlighted. Finally, it is illustratd that there should be a difference in scalings between the momentum pinch originated from inherently toroidal effects and that coming from other mechanisms which exist in a simpler geometry.
Fully Developed Turbulent Mixing in an Annular Sector
Lim, Hyun-Kyung; Zhou, Yijie; de Almeida, Valmor F; Glimm, James G
2014-01-01
We review recent progress on the characterization of turbulent mixing fluid flow and relate these ideas to high-speed, two-phase Couette flow with application to mixing in a centrifugal contactor. The general ideas are more broadly applicable and have been applied to the study of Rayleigh-Taylor and Richtmyer-Meshkov fluid mixing, combustion in the engine of a scram jet and the analysis of inertial confinement pellet simulations.
HIERARCHICAL STRUCTURE OF MAGNETOHYDRODYNAMIC TURBULENCE IN POSITION-POSITION-VELOCITY SPACE
Burkhart, Blakesley; Lazarian, A.; Goodman, Alyssa; Rosolowsky, Erik
2013-06-20
Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium (ISM) that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic position-position-velocity (PPV) emission cubes of isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the presence of self-gravity and the global sonic and Alfvenic Mach numbers. Simulations with higher Alfvenic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observe a strong dependency on the sonic and Alfvenic Mach numbers and self-gravity when we apply the statistical moments (i.e., mean, variance, skewness, kurtosis) to the leaf and node distribution of the dendrogram. Simulations with self-gravity, larger magnetic field and higher sonic Mach number have dendrogram distributions with higher statistical moments. Application of the dendrogram to three-dimensional density cubes, also known as position-position-position (PPP) cubes, reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening, and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions for obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP data.
Practical application of large eddy simulation to film cooling flow analysis on gas turbine airfoils
Takata, T.; Takeishi, K.; Kawata, Y.; Tsuge, A.
1999-07-01
Large eddy simulation (LES) using body-fitted coordinates is applied to solve film cooling flow on turbine blades. The turbulent model was tuned using the experimental flow field and adiabatic film cooling effectiveness measurements for a single row of holes on a flat plate surface. The results show the interaction between the main stream boundary layer and injected film cooling air generates kidney and horseshoe shaped vortices. Comparison of the temperature distribution between experimental results and present analysis has been conducted. The non-dimensional temperature distribution at x/d = 1 is dome style and quantitatively agrees with experimental results. LES was also applied to solve film cooling on a turbine airfoil. If LES was applied to solve whole flow field domain large CPU time would make the solution impractical. LES, using body-fitted coordinates, is applied to solve the non-isotropic film cooling flow near the turbine blade. The cascade flow domain, with a pitch equal to one film cooling hole spacing, is solved using {kappa}-{epsilon} model. By using such a hybrid numerical method, CPU time is reduced and numerical accuracy is insured. The analytical results show the interaction between the flow blowing through film cooling holes and mainstream on the suction and pressure surfaces of the turbine airfoil. They also show the fundamental structure of the film cooling air flow is governed by arch internal secondary flow and horseshoe vortices which have a similar structure to film cooling air flow blowing through a cooling hole on a flat plate. In the flow field, the effect of turbulent structure on curvature (relaminarization) and flow pattern, involving the interaction between main flow and the cooling jet, are clearly shown. Film cooling effectiveness on the blade surface is predicted from the results of the thermal field calculation and is compared with the test result.
Turbulent Flame Propagation Characteristics of High Hydrogen Content Fuels
Seitzman, Jerry; Lieuwen, Timothy
2014-09-30
evidence that the leading points model can provide useful predictions of turbulent flame speed over a wide range of operating conditions and flow geometries.
Sun, K.; Chenu, A.; Mikityuk, K.; Krepel, J.; Chawla, R.
2012-07-01
The core behaviour of a large (3600 MWth) sodium-cooled fast reactor (SFR) is investigated in this paper with the use of a coupled TRACE/PARCS model. The SFR neutron spectrum is characterized by several performance advantages, but also leads to one dominating neutronics drawback - a positive sodium void reactivity. This implies a positive reactivity effect when sodium coolant is removed from the core. In order to evaluate such feedback in terms of the dynamics, a representative unprotected loss-of-flow (ULOF) transient, i.e. flow run-down without SCRAM in which sodium boiling occurs, is analyzed. Although analysis of a single transient cannot allow general conclusions to be drawn, it does allow better understanding of the underlying physics and can lead to proposals for improving the core response during such an accident. The starting point of this study is the reference core design considered in the framework of the Collaborative Project on the European Sodium Fast Reactor (CP-ESFR). To reduce the void effect, the core has been modified by introducing an upper sodium plenum (along with a boron layer) and by reducing the core height-to-diameter ratio. For the ULOF considered, a sharp increase in core power results in melting of the fuel in the case of the reference core. In the modified core, a large dryout leads to melting of the clad. It seems that, for the hypothetical event considered, fuel failure cannot be avoided with just improvement of the neutronics design; therefore, thermal-hydraulics optimization has been considered. An innovative assembly design is proposed to prevent sodium vapour blocking the fuel channel. This results in preventing a downward propagation of the sodium boiling to the core center, thus limiting it to the upper region. Such a void map introduces a negative coolant density reactivity feedback, which dominates the total reactivity change. As a result, the power level and the fuel temperature are effectively reduced, and a large dryout
Hall MHD Stability and Turbulence in Magnetically Accelerated Plasmas
H. R. Strauss
2012-11-27
The object of the research was to develop theory and carry out simulations of the Z pinch and plasma opening switch (POS), and compare with experimental results. In the case of the Z pinch, there was experimental evidence of ion kinetic energy greatly in excess of the ion thermal energy. It was thought that this was perhaps due to fine scale turbulence. The simulations showed that the ion energy was predominantly laminar, not turbulent. Preliminary studies of a new Z pinch experiment with an axial magnetic field were carried out. The axial magnetic is relevant to magneto - inertial fusion. These studies indicate the axial magnetic field makes the Z pinch more turbulent. Results were also obtained on Hall magnetohydrodynamic instability of the POS.
Magnetized Turbulent Dynamo in Protogalaxies
Leonid Malyshkin; Russell M. Kulsrud
2002-01-28
The prevailing theory for the origin of cosmic magnetic fields is that they have been amplified to their present values by the turbulent dynamo inductive action in the protogalactic and galactic medium. Up to now, in calculation of the turbulent dynamo, it has been customary to assume that there is no back reaction of the magnetic field on the turbulence, as long as the magnetic energy is less than the turbulent kinetic energy. This assumption leads to the kinematic dynamo theory. However, the applicability of this theory to protogalaxies is rather limited. The reason is that in protogalaxies the temperature is very high, and the viscosity is dominated by magnetized ions. As the magnetic field strength grows in time, the ion cyclotron time becomes shorter than the ion collision time, and the plasma becomes strongly magnetized. As a result, the ion viscosity becomes the Braginskii viscosity. Thus, in protogalaxies the back reaction sets in much earlier, at field strengths much lower than those which correspond to field-turbulence energy equipartition, and the turbulent dynamo becomes what we call the magnetized turbulent dynamo. In this paper we lay the theoretical groundwork for the magnetized turbulent dynamo. In particular, we predict that the magnetic energy growth rate in the magnetized dynamo theory is up to ten times larger than that in the kinematic dynamo theory. We also briefly discuss how the Braginskii viscosity can aid the development of the inverse cascade of magnetic energy after the energy equipartition is reached.
Fichtl, G.H.
1983-09-01
When designing a wind energy converison system (WECS), it may be necessary to take into account the distribution of wind across the disc of rotation. The specific engineering applications include structural strength, fatigue, and control. This wind distribution consists of two parts, namely that associated with the mean wind profile and that associated with the turbulence velocity fluctuation field. The work reported herein is aimed at the latter, namely the distribution of turbulence velocity fluctuations across the WECS disk of rotation. A theory is developed for the two-time covariance matrix for turbulence velocity vector components for wind energy conversion system (WECS) design. The theory is developed for homogeneous and iotropic turbulance with the assumption that Taylor's hypothesis is valid. The Eulerian turbulence velocity vector field is expanded about the hub of the WECS. Formulae are developed for the turbulence velocity vector component covariance matrix following the WECS blade elements. It is shown that upon specification of the turbulence energy spectrum function and the WECS rotation rate, the two-point, two-time covariance matrix of the turbulent flow relative to the WECS bladed elements is determined. This covariance matrix is represented as the sum of nonstationary and stationary contributions. Generalized power spectral methods are used to obtain two-point, double frequency power spectral density functions for the turbulent flow following the blade elements. The Dryden turbulence model is used to demonstrate the theory. A discussion of linear system response analysis is provided to show how the double frequency turbulence spectra might be used to calculate response spectra of a WECS to turbulent flow. Finally the spectrum of the component of turbulence normal to the WECS disc of rotation, following the blade elements, is compared with experimental results.
UNDERSTANDING GALAXY OUTFLOWS AS THE PRODUCT OF UNSTABLE TURBULENT SUPPORT
Scannapieco, Evan
2013-02-01
The interstellar medium is a multiphase gas in which turbulent support is as important as thermal pressure. Sustaining this configuration requires both continuous turbulent stirring and continuous radiative cooling to match the decay of turbulent energy. While this equilibrium can persist for small turbulent velocities, if the one-dimensional velocity dispersion is larger than Almost-Equal-To 35 km s{sup -1}, the gas moves into an unstable regime that leads to rapid heating. I study the implications of this turbulent runaway, showing that it causes a hot gas outflow to form in all galaxies with a gas surface density above Almost-Equal-To 50 M{sub Sun} pc{sup -2}, corresponding to a star formation rate per unit area of Almost-Equal-To 0.1 M{sub Sun} yr{sup -1} kpc{sup -2}. For galaxies with v{sub esc} {approx}> 200 km s{sup -1}, the sonic point of this hot outflow should lie interior to the region containing cold gas and stars, while for galaxies with smaller escape velocities, the sonic point should lie outside this region. This leads to efficient cold cloud acceleration in higher mass galaxies, while in lower mass galaxies, clouds may be ejected by random turbulent motions rather than accelerated by the wind. Finally, I show that energy balance cannot be achieved at all for turbulent media above a surface density of Almost-Equal-To 10{sup 5} M{sub Sun} pc{sup -2}.
Study of natural circulation in a VHTR after a LOFA using different turbulence models
Yu-Hsin Tung; Yuh-Ming Ferng; Richard W. Johnson; Ching-Chang Chieng
2013-10-01
Natural convection currents in the core are anticipated in the event of the failure of the gas circulator in a prismatic gas-cooled very high temperature reactor (VHTR). The paths that the helium coolant takes in forming natural circulation loops and the effective heat transport are of interest. The heated flow in the reactor core is turbulent during normal operating conditions and at the beginning of the LOFA with forced convection, but the flow may significantly be slowed down after the event and laminarized with mixed convection. In the present study, the potential occurrence and effective heat transport of natural circulation are demonstrated using computational fluid dynamic (CFD) calculations with different turbulence models as well as laminar flow. Validations and recommendation on turbulence model selection are conducted. The study concludes that large loop natural convection is formed due to the enhanced turbulence levels by the buoyancy effect and the turbulent regime near the interface of upper plenum and flow channels increases the flow resistance for channel flows entering upper plenum and thus less heat can be removed from the core than the prediction by laminar flow assumption.
National Solar Thermal Test Facility
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
SunShot Grand Challenge: Regional Test Centers National Solar Thermal Test Facility HomeTag:National Solar Thermal Test Facility Molten Nitrate Salt Initial Flow Testing is a ...
ENSEMBLE SIMULATIONS OF PROTON HEATING IN THE SOLAR WIND VIA TURBULENCE AND ION CYCLOTRON RESONANCE
Cranmer, Steven R.
2014-07-01
Protons in the solar corona and heliosphere exhibit anisotropic velocity distributions, violation of magnetic moment conservation, and a general lack of thermal equilibrium with the other particle species. There is no agreement about the identity of the physical processes that energize non-Maxwellian protons in the solar wind, but a traditional favorite has been the dissipation of ion cyclotron resonant AlfvĆ©n waves. This paper presents kinetic models of how ion cyclotron waves heat protons on their journey from the corona to interplanetary space. It also derives a wide range of new solutions for the relevant dispersion relations, marginal stability boundaries, and nonresonant velocity-space diffusion rates. A phenomenological model containing both cyclotron damping and turbulent cascade is constructed to explain the suppression of proton heating at low alpha-proton differential flow speeds. These effects are implemented in a large-scale model of proton thermal evolution from the corona to 1 AU. A Monte Carlo ensemble of realistic wind speeds, densities, magnetic field strengths, and heating rates produces a filled region of parameter space (in a plane described by the parallel plasma beta and the proton temperature anisotropy ratio) similar to what is measured. The high-beta edges of this filled region are governed by plasma instabilities and strong heating rates. The low-beta edges correspond to weaker proton heating and a range of relative contributions from cyclotron resonance. On balance, the models are consistent with other studies that find only a small fraction of the turbulent power spectrum needs to consist of ion cyclotron waves.
THE VIOLATION OF THE TAYLOR HYPOTHESIS IN MEASUREMENTS OF SOLAR WIND TURBULENCE
Klein, K. G.; Howes, G. G.; TenBarge, J. M.
2014-08-01
Motivated by the upcoming Solar Orbiter and Solar Probe Plus missions, qualitative and quantitative predictions are made for the effects of the violation of the Taylor hypothesis on the magnetic energy frequency spectrum measured in the near-Sun environment. The synthetic spacecraft data method is used to predict observational signatures of the violation for critically balanced AlfvĆ©nic turbulence or parallel fast/whistler turbulence. The violation of the Taylor hypothesis can occur in the slow flow regime, leading to a shift of the entire spectrum to higher frequencies, or in the dispersive regime, in which the dissipation range spectrum flattens at high frequencies. It is found that AlfvĆ©nic turbulence will not significantly violate the Taylor hypothesis, but whistler turbulence will. The flattening of the frequency spectrum is therefore a key observational signature for fast/whistler turbulence.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Ren, Y.; Wang, W. X.; LeBlanc, B. P.; Guttenfelder, W.; Kaye, S. M.; Ethier, S.; Mazzucato, E.; Lee, K. C.; Domier, C. W.; Bell, R.; et al
2015-11-03
In this letter, we report the first observation of the fast response of electron-scale turbulence to auxiliary heating cessation in National Spherical Torus eXperiment [Ono et al., Nucl. Fusion 40, 557 (2000)]. The observation was made in a set of RF-heated L-mode plasmas with toroidal magnetic field of 0.55 T and plasma current of 300 kA. It is observed that electron-scale turbulence spectral power (measured with a high-k collective microwave scattering system) decreases significantly following fast cessation of RF heating that occurs in less than 200 Ī¼s. The large drop in the turbulence spectral power has a short time delaymoreĀ Ā» of about 1ā2 ms relative to the RF cessation and happens on a time scale of 0.5ā1 ms, much smaller than the energy confinement time of about 10 ms. Power balance analysis shows a factor of about 2 decrease in electron thermal diffusivity after the sudden drop of turbulence spectral power. Measured small changes in equilibrium profiles across the RF cessation are unlikely able to explain this sudden reduction in the measured turbulence and decrease in electron thermal transport, supported by local linear stability analysis and both local and global nonlinear gyrokinetic simulations. Furthermore, the observations imply that nonlocal flux-driven mechanism may be important for the observed turbulence and electron thermal transport.Ā«Ā less
Particle dispersion in homogeneous turbulence using the one-dimensional turbulence model
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sun, Guangyuan; Lignell, David O.; Hewson, John C.; Gin, Craig R.
2014-10-09
Lagrangian particle dispersion is studied using the one-dimensional turbulence (ODT) model in homogeneous decaying turbulence configurations. The ODT model has been widely and successfully applied to a number of reacting and nonreacting flow configurations, but only limited application has been made to multiphase flows. We present a version of the particle implementation and interaction with the stochastic and instantaneous ODT eddy events. The model is characterized by comparison to experimental data of particle dispersion for a range of intrinsic particle time scales and body forces. Particle dispersion, velocity, and integral time scale results are presented. Moreover, the particle implementation introducesmoreĀ Ā» a single model parameter Ī² p , and sensitivity to this parameter and behavior of the model are discussed. Good agreement is found with experimental data and the ODT model is able to capture the particle inertial and trajectory crossing effects. Our results serve as a validation case of the multiphase implementations of ODT for extensions to other flow configurations.Ā«Ā less
Development of a One-Equation Transition/Turbulence Model
EDWARDS,JACK R.; ROY,CHRISTOPHER J.; BLOTTNER,FREDERICK G.; HASSAN,HASSAN A.
2000-09-26
This paper reports on the development of a unified one-equation model for the prediction of transitional and turbulent flows. An eddy viscosity - transport equation for non-turbulent fluctuation growth based on that proposed by Warren and Hassan (Journal of Aircraft, Vol. 35, No. 5) is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittence function based on the work of Dhawan and Narasimha (Journal of Fluid Mechanics, Vol. 3, No. 4). The model predicts both the onset and extent of transition. Low-speed test cases include transitional flow over a flat plate, a single element airfoil, and a multi-element airfoil in landing configuration. High-speed test cases include transitional Mach 3.5 flow over a 5{degree} cone and Mach 6 flow over a flared-cone configuration. Results are compared with experimental data, and the spatial accuracy of selected predictions is analyzed.
Tokuhiro, A.T.; Kimura, N.; Nishimura, M.; Kobayashi, J.; Miyakoshi, H.
1999-07-01
The thermal-hydraulic mixing of three quasi-planar vertical water jets was experimentally and numerically investigated. The central jet was initially 5 C lower in temperature than the other two. The hydraulic diameter and average exit velocity-based Reynolds and Richardson numbers were, Re{sub D} = 2 x 10{sup 4}, Ri{sub D} = 0.002. Besides temperature measurements from a traversing array of 37 thermocouples, velocity measurements were made using laser and ultrasound Doppler velocimetries (LDV and UDV). In parallel the in-house code, CASCADE, featuring a {kappa}-{epsilon} turbulence model was used to simulate the experimental flow configuration. A comparison of the experimental and numerical results showed that code validation by LDV/UDV was possible and in particular that time-averaged field and frequency characteristics of transversely swaying jets, even under Reynolds averaging of the conservation equations, could be simulated. A representative comparison of the amplitude of oscillation is shown in Figure A-1 with an inset of the visualized flow and sample time-series of the temperature fluctuations at the position indicated. The difference in the predominant frequency, the numerically predicted {approximately}1.6 Hz versus the experimental {approximately}2.25 Hz, is attributed to the turbulence model that overestimate thus effective fluid viscosity. Development of an accurate numerical simulation is of relevance to the design of the liquid metal fast breeder reactor (LMFBR), where the lack of mixing of the cold sodium may initiate thermal striping; that is, poorly mixed hot and cold streams may thermally stress the components onto which they impinge. Turbulent mixing of jets is equally of general interest to environmental and material processing flows.
Boundary Layer Cloud Turbulence Characteristics
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Boundary Layer Cloud Turbulence Characteristics Virendra Ghate Bruce Albrecht Parameter Observational Readiness (/10) Modeling Need (/10) Cloud Boundaries 9 9 Cloud Fraction Variance Skewness Up/Downdraft coverage Dominant Freq. signal Dissipation rate ??? Observation-Modeling Interface
Edge turbulence and transport: Text and ATF modeling
Ritz, C.P.; Rhodes, T.L.; Lin, H.; Rowan, W.L.; Bengtson, R.; Wootton, A.J. . Fusion Research Center); Carreras, B.A.; Leboeuf, J.N.; Lee, D.K.; Harris, J.; Hidalgo, C.; Bell, J.D.; Holmes, J.A.; Isler, R.; Lynch, V.E.; Uckan, T. ); Diamond, P.H.; Ware, A.S. ); Thayer, D.R. (Science Applications Inter
1990-01-01
We present experimental results on edge turbulence and transport from the tokamak TEXT and the torsatron ATF. The measured electrostatic fluctuations can explain the edge transport of particles and energy. Certain drive (radiation) and stabilizing (velocity shear) terms are suggested by the results. The experimental fluctuation levels and spectral widths can be reproduced by considering the nonlinear evolution of the reduced MHD equations, incorporating a thermal drive from line radiation. In the tokamak limit (with toroidal electric field) the model corresponds to the resistivity gradient mode, while in the currentless torsatron or stellarator limit it corresponds to a thermally driven drift wave.
Terascale High-Fidelity Simulations of Turbulent Combustion with Detailed Chemistry
Im, Hong G [University of Michigan] [University of Michigan; Trouve, Arnaud [University of Maryland] [University of Maryland; Rutland, Christopher J [University of Wisconsin] [University of Wisconsin; Chen, Jacqueline H [Sandia National Laboratories] [Sandia National Laboratories
2012-08-13
The TSTC project is a multi-university collaborative effort to develop a high-fidelity turbulent reacting flow simulation capability utilizing terascale, massively parallel computer technology. The main paradigm of our approach is direct numerical simulation (DNS) featuring highest temporal and spatial accuracy, allowing quantitative observations of the fine-scale physics found in turbulent reacting flows as well as providing a useful tool for development of sub-models needed in device-level simulations. The code named S3D, developed and shared with Chen and coworkers at Sandia National Laboratories, has been enhanced with new numerical algorithms and physical models to provide predictive capabilities for spray dynamics, combustion, and pollutant formation processes in turbulent combustion. Major accomplishments include improved characteristic boundary conditions, fundamental studies of auto-ignition in turbulent stratified reactant mixtures, flame-wall interaction, and turbulent flame extinction by water spray. The overarching scientific issue in our recent investigations is to characterize criticality phenomena (ignition/extinction) in turbulent combustion, thereby developing unified criteria to identify ignition and extinction conditions. The computational development under TSTC has enabled the recent large-scale 3D turbulent combustion simulations conducted at Sandia National Laboratories.
Terascale High-Fidelity Simulations of Turbulent Combustion with Detailed Chemistry
Hong G. Im; Arnaud Trouve; Christopher J. Rutland; Jacqueline H. Chen
2009-02-02
The TSTC project is a multi-university collaborative effort to develop a high-fidelity turbulent reacting flow simulation capability utilizing terascale, massively parallel computer technology. The main paradigm of our approach is direct numerical simulation (DNS) featuring highest temporal and spatial accuracy, allowing quantitative observations of the fine-scale physics found in turbulent reacting flows as well as providing a useful tool for development of sub-models needed in device-level simulations. The code named S3D, developed and shared with Chen and coworkers at Sandia National Laboratories, has been enhanced with new numerical algorithms and physical models to provide predictive capabilities for spray dynamics, combustion, and pollutant formation processes in turbulent combustion. Major accomplishments include improved characteristic boundary conditions, fundamental studies of auto-ignition in turbulent stratified reactant mixtures, flame-wall interaction, and turbulent flame extinction by water spray. The overarching scientific issue in our recent investigations is to characterize criticality phenomena (ignition/extinction) in turbulent combustion, thereby developing unified criteria to identify ignition and extinction conditions. The computational development under TSTC has enabled the recent large-scale 3D turbulent combustion simulations conducted at Sandia National Laboratories.
None
2016-07-05
Thermal rectifiers using linear nanostructures as core thermal conductors have been fabricated. A high mass density material is added preferentially to one end of the nanostructures to produce an axially non-uniform mass distribution. The resulting nanoscale system conducts heat asymmetrically with greatest heat flow in the direction of decreasing mass density. Thermal rectification has been demonstrated for linear nanostructures that are electrical insulators, such as boron nitride nanotubes, and for nanostructures that are conductive, such as carbon nanotubes.
On the interaction between turbulence and a planar rarefaction
Johnson, Bryan M.
2014-04-01
The modeling of turbulence, whether it be numerical or analytical, is a difficult challenge. Turbulence is amenable to analysis with linear theory if it is subject to rapid distortions, i.e., motions occurring on a timescale that is short compared to the timescale for nonlinear interactions. Such an approach (referred to as rapid distortion theory) could prove useful for understanding aspects of astrophysical turbulence, which is often subject to rapid distortions, such as supernova explosions or the free-fall associated with gravitational instability. As a proof of principle, a particularly simple problem is considered here: the evolution of vorticity due to a planar rarefaction in an ideal gas. Analytical solutions are obtained for incompressive modes having a wave vector perpendicular to the distortion; as in the case of gradient-driven instabilities, these are the modes that couple most strongly to the mean flow. Vorticity can either grow or decay in the wake of a rarefaction front, and there are two competing effects that determine which outcome occurs: entropy fluctuations couple to the mean pressure gradient to produce vorticity via baroclinic effects, whereas vorticity is damped due to the conservation of angular momentum as the fluid expands. Whether vorticity grows or decays depends upon the ratio of entropic to vortical fluctuations at the location of the front; growth occurs if this ratio is of order unity or larger. In the limit of purely entropic fluctuations in the ambient fluid, a strong rarefaction generates vorticity with a turbulent Mach number on the order of the rms of the ambient entropy fluctuations. The analytical results are shown to compare well with results from two- and three-dimensional numerical simulations. Analytical solutions are also derived in the linear regime of Reynolds-averaged turbulence models. This highlights an inconsistency in standard turbulence models that prevents them from accurately capturing the physics of
The time evolution of turbulent parameters in reversed-field pinch plasmas
Titus, J. B.; Alexander, Brandon; Johnson, J. A. III
2013-04-28
Turbulence is abundant in fully ionized fusion plasmas, with unique turbulent characteristics in different phases of the discharge. Using Fourier and chaos-based techniques, a set of parameters have been developed to profile the time evolution of turbulence in high temperature, fusion plasmas, specifically in self-organized, reversed-field pinch plasma in the Madison Symmetric Torus. With constant density and plasma current, the turbulence profile is measured during ramp-up, magnetic reconnection, and increased confinement phases. During magnetic reconnection, a scan of plasma current is performed with a constant density. Analysis revealed that the energy associated with turbulence (turbulent energy) is found to increase when changes in magnetic energy occur and is correlated to edge ion temperatures. As the turbulent energy increases with increasing current, the rate at which this energy flow between scales (spectral index) and anti-persistence of the fluctuations increases (Hurst exponent). These turbulent parameters are then compared to the ramp-up phase and increased confinement regime.
Power-law wrinkling turbulence-flame interaction model for astrophysical flames
Jackson, Aaron P.; Townsley, Dean M.; Calder, Alan C.
2014-04-01
We extend a model for turbulence-flame interactions (TFI) to consider astrophysical flames with a particular focus on combustion in Type Ia supernovae. The inertial range of the turbulent cascade is nearly always under-resolved in simulations of astrophysical flows, requiring the use of a model in order to quantify the effects of subgrid-scale wrinkling of the flame surface. We provide implementation details to extend a well-tested TFI model to low-Prandtl number flames for use in the compressible hydrodynamics code FLASH. A local, instantaneous measure of the turbulent velocity is calibrated for FLASH and verification tests are performed. Particular care is taken to consider the relation between the subgrid rms turbulent velocity and the turbulent flame speed, especially for high-intensity turbulence where the turbulent flame speed is not expected to scale with the turbulent velocity. Finally, we explore the impact of different TFI models in full-star, three-dimensional simulations of Type Ia supernovae.
Garcia, L.; Carreras, B. A.; Llerena, I.; Calvo, I. [Departamento de Fisica, Universidad Carlos III de Madrid, 28911 Leganes, Madrid (Spain); Departament d'Algebra i Geometria, Facultat de Matematiques, Universitat de Barcelona, Barcelona (Spain); Laboratorio Nacional de Fusion, Asociacion EURATOM-CIEMAT, 28040 Madrid (Spain)
2009-10-15
For the resistive pressure-gradient-driven turbulence model, the transition from laminar regime to fully developed turbulence is not simple and goes through several phases. For low values of the plasma parameter {beta}, a single quasicoherent structure forms. As {beta} increases, several of these structures may emerge and in turn take the dominant role. Finally, at high {beta}, fully developed turbulence with a broad spectrum is established. A suitable characterization of this transition can be given in terms of topological properties of the flow. Here, we analyze these properties that provide an understanding of the turbulence-induced transport and give a measure of the breaking of the homogeneity of the turbulence. To this end, an approach is developed that allows discriminating between topological properties of plasma turbulence flows that are relevant to the transport dynamics and the ones that are not. This is done using computational homology tools and leads to a faster convergence of numerical results for a fixed level of resolution than previously presented in Phys. Rev. E 78, 066402 (2008)
STATIONARITY IN SOLAR WIND FLOWS
Perri, S.; Balogh, A. E-mail: a.balogh@imperial.ac.u
2010-05-01
By using single-point measurements in space physics it is possible to study a phenomenon only as a function of time. This means that we cannot have direct access to information about spatial variations of a measured quantity. However, the investigation of the properties of turbulence and of related phenomena in the solar wind widely makes use of an approximation frequently adopted in hydrodynamics under certain conditions, the so-called Taylor hypothesis; indeed, the solar wind flow has a bulk velocity along the radial direction which is much higher than the velocity of a single turbulent eddy embedded in the main flow. This implies that the time of evolution of the turbulent features is longer than the transit time of the flow through the spacecraft position, so that the turbulent field can be considered frozen into the solar wind flow. This assumption allows one to easily associate time variations with spatial variations and stationarity to homogeneity. We have investigated, applying criteria for weak stationarity to Ulysses magnetic field data in different solar wind regimes, at which timescale and under which conditions the hypothesis of stationarity, and then of homogeneity, of turbulence in the solar wind is well justified. We extend the conclusions of previous studies by Matthaeus and Goldstein to different parameter ranges in the solar wind. We conclude that the stationarity assumption in the inertial range of turbulence on timescales of 10 minutes to 1 day is reasonably satisfied in fast and uniform solar wind flows, but that in mixed, interacting fast, and slow solar wind streams the assumption is frequently only marginally valid.
Dissipation of turbulence in the wake of a wind turbine
Lundquist, J. K.; Bariteau, L.
2014-11-06
The wake of a wind turbine is characterized by increased turbulence and decreased wind speed. Turbines are generally deployed in large groups in wind farms, and so the behaviour of an individual wake as it merges with other wakes and propagates downwind is critical in assessing wind-farm power production. This evolution depends on the rate of turbulence dissipation in the wind-turbine wake, which has not been previously quantified in field-scale measurements. In situ measurements of winds and turbulence dissipation from the wake region of a multi-MW turbine were collected using a tethered lifting system (TLS) carrying a payload of high-rate turbulence probes. Ambient flow measurements were provided from sonic anemometers on a meteorological tower located near the turbine. Good agreement between the tower measurements and the TLS measurements was established for a case without a wind-turbine wake. When an operating wind turbine is located between the tower and the TLS so that the wake propagates to the TLS, the TLS measures dissipation rates one to two orders of magnitude higher in the wake than outside of the wake. These data, collected between two and three rotor diameters D downwind of the turbine, document the significant enhancement of turbulent kinetic energy dissipation rate within the wind-turbine wake. These wake measurements suggest that it may be useful to pursue modelling approaches that account for enhanced dissipation. Furthermore. comparisons of wake and non-wake dissipation rates to mean wind speed, wind-speed variance, and turbulence intensity are presented to facilitate the inclusion of these measurements in wake modelling schemes.
Dissipation of turbulence in the wake of a wind turbine
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Lundquist, J. K.; Bariteau, L.
2014-11-06
The wake of a wind turbine is characterized by increased turbulence and decreased wind speed. Turbines are generally deployed in large groups in wind farms, and so the behaviour of an individual wake as it merges with other wakes and propagates downwind is critical in assessing wind-farm power production. This evolution depends on the rate of turbulence dissipation in the wind-turbine wake, which has not been previously quantified in field-scale measurements. In situ measurements of winds and turbulence dissipation from the wake region of a multi-MW turbine were collected using a tethered lifting system (TLS) carrying a payload of high-ratemoreĀ Ā» turbulence probes. Ambient flow measurements were provided from sonic anemometers on a meteorological tower located near the turbine. Good agreement between the tower measurements and the TLS measurements was established for a case without a wind-turbine wake. When an operating wind turbine is located between the tower and the TLS so that the wake propagates to the TLS, the TLS measures dissipation rates one to two orders of magnitude higher in the wake than outside of the wake. These data, collected between two and three rotor diameters D downwind of the turbine, document the significant enhancement of turbulent kinetic energy dissipation rate within the wind-turbine wake. These wake measurements suggest that it may be useful to pursue modelling approaches that account for enhanced dissipation. Furthermore. comparisons of wake and non-wake dissipation rates to mean wind speed, wind-speed variance, and turbulence intensity are presented to facilitate the inclusion of these measurements in wake modelling schemes.Ā«Ā less
Large Eddy Simulations of Combustor Liner Flows | Argonne Leadership
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Computing Facility Turbulent flow going through a multi-hole coupon geometry Turbulent flow going through a multi-hole coupon geometry. Density iso surface colored by velocity magnitude. Joseph Insley, Argonne National Laboratory; Lee Shunn, Cascade Technologies; Anne Dord, GE Aviation; and Yonduck Sung, GE Global Research Large Eddy Simulations of Combustor Liner Flows PI Name: Anne Dord PI Email: anne.dord@ge.com Institution: GE Aviation Allocation Program: INCITE Allocation Hours at ALCF:
Fletcher, C.D.; Ghan, L.S.; Determan, J.C.; Nielsen, H.H. )
1994-04-01
A system model of the Advanced Neutron Source Reactor (ANSR) has been developed and used to perform conceptual safety analyses. To better represent thermal-hydraulic behavior in the unique geometry and conditions of the ANSR core, three specific changes in the RELAP5/MOD3 computer code were implemented: a turbulent forced-convection heat transfer correlation, a critical heat flux correlation, and an interfacial drag correlation. The system model includes representations of the ANSR core, heat exchanger coolant loops, and the pressurizing and letdown systems. Analyses of ANSR station blackout and loss-of-flow accident scenarios are described. The results show that the core can survive without exceeding the flow excursion or critical heat flux thermal limits defined for the conceptual safety analysis, if the proper mitigation options are provided.
Turbulence at Hydroelectric Power Plants and its Potential Effects on Fish.
Cada, Glenn F.; Odeh, Mufeed
2001-01-01
The fundamental influence of fluid dynamics on aquatic organisms is receiving increasing attention among aquatic ecologists. For example, the importance of turbulence to ocean plankton has long been a subject of investigation (Peters and Redondo 1997). More recently, studies have begun to emerge that explicitly consider the effects of shear and turbulence on freshwater invertebrates (Statzner et al. 1988; Hart et al. 1996) and fishes (Pavlov et al. 1994, 1995). Hydraulic shear stress and turbulence are interdependent natural fluid phenomena that are important to fish, and consequently it is important to develop an understanding of how fish sense, react to, and perhaps utilize these phenomena under normal river flows. The appropriate reaction to turbulence may promote movement of migratory fish or prevent displacement of resident fish. It has been suggested that one of the adverse effects of flow regulation by hydroelectric projects is the reduction of normal turbulence, particularly in the headwaters of reservoirs, which can lead to disorientation and slowing of migration (Williams et al. 1996; Coutant et al. 1997; Coutant 1998). On the other hand, greatly elevated levels of shear and turbulence may be injurious to fish; injuries can range from removal of the mucous layer on the body surface to descaling to torn opercula, popped eyes, and decapitation (Neitzel et al. 2000a,b). Damaging levels of fluid stress can occur in a variety of circumstances in both natural and man-made environments. This paper discusses the effects of shear stress and turbulence on fish, with an emphasis on potentially damaging levels in man-made environments. It defines these phenomena, describes studies that have been conducted to understand their effects, and identifies gaps in our knowledge. In particular, this report reviews the available information on the levels of turbulence that can occur within hydroelectric power plants, and the associated biological effects. The final section
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sankaran, Ramanan; Hawkes, Evatt R.; Yoo, Chun Sang; Chen, Jacqueline H.
2015-06-22
Direct numerical simulations of three-dimensional spatially-developing turbulent Bunsen flames were performed at three different turbulence intensities. We performed these simulations using a reduced methaneāair chemical mechanism which was specifically tailored for the lean premixed conditions simulated here. A planar-jet turbulent Bunsen flame configuration was used in which turbulent preheated methaneāair mixture at 0.7 equivalence ratio issued through a central jet and was surrounded by a hot laminar coflow of burned products. The turbulence characteristics at the jet inflow were selected such that combustion occured in the thin reaction zones (TRZ) regime. At the lowest turbulence intensity, the conditions fall onmoreĀ Ā» the boundary between the TRZ regime and the corrugated flamelet regime, and progressively moved further into the TRZ regime by increasing the turbulent intensity. The data from the three simulations was analyzed to understand the effect of turbulent stirring on the flame structure and thickness. Furthermore, statistical analysis of the data showed that the thermal preheat layer of the flame was thickened due to the action of turbulence, but the reaction zone was not significantly affected. A global and local analysis of the burning velocity of the flame was performed to compare the different flames. Detailed statistical averages of the flame speed were also obtained to study the spatial dependence of displacement speed and its correlation to strain rate and curvature.Ā«Ā less
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sankaran, Ramanan; Hawkes, Evatt R.; Yoo, Chun Sang; Chen, Jacqueline H.
2015-06-22
Direct numerical simulations of three-dimensional spatially-developing turbulent Bunsen flames were performed at three different turbulence intensities. We performed these simulations using a reduced methaneair chemical mechanism which was specifically tailored for the lean premixed conditions simulated here. A planar-jet turbulent Bunsen flame configuration was used in which turbulent preheated methaneair mixture at 0.7 equivalence ratio issued through a central jet and was surrounded by a hot laminar coflow of burned products. The turbulence characteristics at the jet inflow were selected such that combustion occured in the thin reaction zones (TRZ) regime. At the lowest turbulence intensity, the conditions fall onmore »the boundary between the TRZ regime and the corrugated flamelet regime, and progressively moved further into the TRZ regime by increasing the turbulent intensity. The data from the three simulations was analyzed to understand the effect of turbulent stirring on the flame structure and thickness. Furthermore, statistical analysis of the data showed that the thermal preheat layer of the flame was thickened due to the action of turbulence, but the reaction zone was not significantly affected. A global and local analysis of the burning velocity of the flame was performed to compare the different flames. Detailed statistical averages of the flame speed were also obtained to study the spatial dependence of displacement speed and its correlation to strain rate and curvature.« less
Sankaran, Ramanan; Hawkes, Evatt R.; Yoo, Chun Sang; Chen, Jacqueline H.
2015-06-22
Direct numerical simulations of three-dimensional spatially-developing turbulent Bunsen flames were performed at three different turbulence intensities. We performed these simulations using a reduced methaneāair chemical mechanism which was specifically tailored for the lean premixed conditions simulated here. A planar-jet turbulent Bunsen flame configuration was used in which turbulent preheated methaneāair mixture at 0.7 equivalence ratio issued through a central jet and was surrounded by a hot laminar coflow of burned products. The turbulence characteristics at the jet inflow were selected such that combustion occured in the thin reaction zones (TRZ) regime. At the lowest turbulence intensity, the conditions fall on the boundary between the TRZ regime and the corrugated flamelet regime, and progressively moved further into the TRZ regime by increasing the turbulent intensity. The data from the three simulations was analyzed to understand the effect of turbulent stirring on the flame structure and thickness. Furthermore, statistical analysis of the data showed that the thermal preheat layer of the flame was thickened due to the action of turbulence, but the reaction zone was not significantly affected. A global and local analysis of the burning velocity of the flame was performed to compare the different flames. Detailed statistical averages of the flame speed were also obtained to study the spatial dependence of displacement speed and its correlation to strain rate and curvature.
Magnetohydrodynamic turbulence: Observation and experiment
Brown, M. R.; Schaffner, D. A.; Weck, P. J.
2015-05-15
We provide a tutorial on the paradigms and tools of magnetohydrodynamic (MHD) turbulence. The principal paradigm is that of a turbulent cascade from large scales to small, resulting in power law behavior for the frequency power spectrum for magnetic fluctuations E{sub B}(f). We will describe five useful statistical tools for MHD turbulence in the time domain: the temporal autocorrelation function, the frequency power spectrum, the probability distribution function of temporal increments, the temporal structure function, and the permutation entropy. Each of these tools will be illustrated with an example taken from MHD fluctuations in the solar wind. A single dataset from the Wind satellite will be used to illustrate all five temporal statistical tools.
Supercomputers Capture Turbulence in the Solar Wind
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Supercomputers Capture Turbulence in the Solar Wind Supercomputers Capture Turbulence in the Solar Wind Berkeley Lab visualizations could help scientists forecast destructive space weather December 16, 2013 Linda Vu, +1 510 495 2402, lvu@lbl.gov eddies1.jpg This visualization zooms in on current sheets revealing the "cascade of turbulence" in the solar wind occurring down to electron scales. This is a phenomenon common in fluid dynamics-turbulent energy injected at large eddies is
Noise correction of turbulent spectra obtained from Acoustic Doppler Velocimeters
Durgesh, Vibhav; Thomson, Jim; Richmond, Marshall C.; Polagye, Brian
2014-03-02
Accurately estimated auto-spectral density functions are essential for characterization of turbulent flows, and they also have applications in computational fluid dynamics modeling, site and inflow characterization for hydrokinetic turbines, and inflow turbulence generation. The Acoustic Doppler Velocimeter (ADV) provides single-point temporally resolved data, that are used to characterize turbulent flows in rivers, seas, and oceans. However, ADV data are susceptible to contamination from various sources, including instrument noise, which is the intrinsic limit to the accuracy of acoustic velocity measurements. Due to the presence of instrument noise, the spectra obtained are altered at high frequencies. The focus of this study is to develop a robust and effective method for accurately estimating auto-spectral density functions from ADV data by reducing or removing the spectral contribution derived from instrument noise. For this purpose, the āNoise Auto-Correlationā (NAC) approach was developed, which exploits the correlation properties of instrument noise to identify and remove its contribution from spectra. The spectra estimated using the NAC approach exhibit increased fidelity and a slope of -5/3 in the inertial range, which is typically observed for turbulent flows. Finally, this study also compares the effectiveness of low-pass Gaussian filters in removing instrument noise with that of the NAC approach. For the data used in this study, both the NAC and Gaussian filter approaches are observed to be capable of removing instrument noise at higher frequencies from the spectra. However, the NAC results are closer to the expected frequency power of -5/3 in the inertial sub-range.
Numerical Analysis of Coolant Flow and Heat Transfer in ITER Diagnostic First Wall
Khodak, A.; Loesser, G.; Zhai, Y.; Udintsev, V.; Klabacha, J.; Wang, W.; Johnson, D.; Feder, R.
2015-07-24
We performed numerical simulations of the ITER Diagnostic First Wall (DFW) using ANSYS workbench. During operation DFW will include solid main body as well as liquid coolant. Thus thermal and hydraulic analysis of the DFW was performed using conjugated heat transfer approach, in which heat transfer was resolved in both solid and liquid parts, and simultaneously fluid dynamics analysis was performed only in the liquid part. This approach includes interface between solid and liquid part of the systemAnalysis was performed using ANSYS CFX software. CFX software allows solution of heat transfer equations in solid and liquid part, and solution of the flow equations in the liquid part. Coolant flow in the DFW was assumed turbulent and was resolved using Reynolds averaged Navier-Stokes equations with Shear Stress Transport turbulence model. Meshing was performed using CFX method available within ANSYS. The data cloud for thermal loading consisting of volumetric heating and surface heating was imported into CFX Volumetric heating source was generated using Attila software. Surface heating was obtained using radiation heat transfer analysis. Our results allowed us to identify areas of excessive heating. Proposals for cooling channel relocation were made. Additional suggestions were made to improve hydraulic performance of the cooling system.
Numerical Analysis of Coolant Flow and Heat Transfer in ITER Diagnostic First Wall
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Khodak, A.; Loesser, G.; Zhai, Y.; Udintsev, V.; Klabacha, J.; Wang, W.; Johnson, D.; Feder, R.
2015-07-24
We performed numerical simulations of the ITER Diagnostic First Wall (DFW) using ANSYS workbench. During operation DFW will include solid main body as well as liquid coolant. Thus thermal and hydraulic analysis of the DFW was performed using conjugated heat transfer approach, in which heat transfer was resolved in both solid and liquid parts, and simultaneously fluid dynamics analysis was performed only in the liquid part. This approach includes interface between solid and liquid part of the systemAnalysis was performed using ANSYS CFX software. CFX software allows solution of heat transfer equations in solid and liquid part, and solution ofmoreĀ Ā» the flow equations in the liquid part. Coolant flow in the DFW was assumed turbulent and was resolved using Reynolds averaged Navier-Stokes equations with Shear Stress Transport turbulence model. Meshing was performed using CFX method available within ANSYS. The data cloud for thermal loading consisting of volumetric heating and surface heating was imported into CFX Volumetric heating source was generated using Attila software. Surface heating was obtained using radiation heat transfer analysis. Our results allowed us to identify areas of excessive heating. Proposals for cooling channel relocation were made. Additional suggestions were made to improve hydraulic performance of the cooling system.Ā«Ā less
The impact of small-scale turbulence on laminar magnetic reconnection
Watson, P. G.; Oughton, S.; Craig, I. J. D. [School of Physics, University of Sydney, NSW 2006 (Australia); Department of Mathematics, University of Waikato, Private Bag 3105, Hamilton (New Zealand)
2007-03-15
Initial states in incompressible two-dimensional magnetohydrodynamics that are known to lead to strong current sheets and (laminar) magnetic reconnection are modified by the addition of small-scale turbulent perturbations of various energies. The evolution of these states is computed with the aim of ascertaining the influence of the turbulence on the underlying laminar solution. Two main questions are addressed here: (1) What effect does small-scale turbulence have on the energy dissipation rate of the underlying solution? (2) What is the threshold turbulent perturbation level above which the original laminar reconnective dynamics is no longer recognizable. The simulations show that while the laminar dynamics persist the dissipation rates are largely unaffected by the turbulence, other than modest increases attributable to the additional small length scales present in the new initial condition. The solutions themselves are also remarkably insensitive to small-scale turbulent perturbations unless the perturbations are large enough to undermine the integrity of the underlying cellular flow pattern. Indeed, even initial states that lead to the evolution of small-scale microscopic sheets can survive the addition of modest turbulence. The role of a large-scale organizing background magnetic field is also addressed.
Terascale Direct Numerical Simulations of Turbulent Combustion: Capabilities and Limits (PReSS Talk)
Yoo, Chun Sang
2009-03-26
The rapid growth in computational capabilities has provided great opportunities for direct numerical simulations (DNS) of turbulent combustion, a type of simulations without any turbulence model. With the help of terascale high performance supercomputing (HPC) resources, we are now able to provide fundamental insight into turbulence-chemistry interaction in simple laboratory-scale turbulent flames with detailed chemistry using three-dimensional (3D) DNS. However, the actual domain size of 3D-DNS is still limited within {approx} O(10 cm{sup 3}) due to its tremendously high grid resolution required to resolve the smallest turbulent length scale as well as flame structures. Moreover, 3D-DNS will require more computing powers to investigate next-generation engines, of which operating conditions will be characterized by higher pressures, lower temperatures, and higher levels of dilution. In this talk, I will discuss the capabilities and limits of DNS of turbulent combustion and present some results of ignition/extinction characteristics of a highly diluted hydrogen flame counter-flowing against heated air. The results of our recent 3D-DNS of a spatially-developing turbulent lifted hydrogen jet flame in heated coflow will also be presented. The 3D-DNS was performed at a jet Reynolds number of 11,000 with {approx} 1 billion grid points, which required 3.5 million CPU hours on Cray XT3/XT4 at Oak Ridge National Laboratories.
Turbulence elasticityA new mechanism for transport barrier dynamics
Guo, Z. B.; Diamond, P. H.; Kosuga, Y.; Gürcan, Ö. D.
2014-09-15
We present a new, unified model of transport barrier formation in elastic drift wave-zonal flow (DW-ZF) turbulence. A new physical quantitythe delay time (i.e., the mixing time for the DW turbulence)is demonstrated to parameterize each stage of the transport barrier formation. Quantitative predictions for the onset of limit-cycle-oscillation (LCO) among DW and ZF intensities (also denoted as I-mode) and I-mode to high-confinement mode (H-mode) transition are also given. The LCO occurs when the ZF shearing rate (|?v?{sub ZF}{sup ?}|) enters the regime ??{sub k}<|?V?{sub ZF}{sup ?}|turbulence decorrelation rate and ?{sub cr} is the threshold delay time. In the basic predator-prey feedback system, ?{sub cr} is also derived. The I-H transition occurs when |?V?{sub E×B}{sup ?}|>?{sub cr}{sup ?1}, where the mean E?×?B shear flow driven by ion pressure locks the DW-ZF system to the H-mode by reducing the delay time below the threshold value.
Vorticity dynamics after the shockāturbulence interaction
Livescu, Daniel; Ryu, Jaiyoung
2015-07-23
In this article, the interaction of a shock wave with quasi-vortical isotropic turbulence (IT) represents a basic problem for studying some of the phenomena associated with high speed flows, such as hypersonic flight, supersonic combustion and Inertial Confinement Fusion (ICF). In general, in practical applications, the shock width is much smaller than the turbulence scales and the upstream turbulent Mach number is modest. In this case, recent high resolution shock-resolved Direct Numerical Simulations (DNS) (Ryu and Livescu, J Fluid Mech 756, R1, 2014) show that the interaction can be described by the Linear Interaction Approximation (LIA). Using LIA to alleviate the need to resolve the shock, DNS post-shock data can be generated at much higher Reynolds numbers than previously possible. Here, such results with Taylor Reynolds number approximately 180 are used to investigate the changes in the vortical structure as a function of the shock Mach number, M_{s}, up to M_{s} = 10. It is shown that, as M_{s} increases, the shock interaction induces a tendency towards a local axisymmetric state perpendicular to the shock front, which has a profound influence on the vortex-stretching mechanism and divergence of the Lamb vector and, ultimately, on the flow evolution away from the shock.
Vorticity dynamics after the shockāturbulence interaction
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Livescu, Daniel; Ryu, Jaiyoung
2015-07-23
In this article, the interaction of a shock wave with quasi-vortical isotropic turbulence (IT) represents a basic problem for studying some of the phenomena associated with high speed flows, such as hypersonic flight, supersonic combustion and Inertial Confinement Fusion (ICF). In general, in practical applications, the shock width is much smaller than the turbulence scales and the upstream turbulent Mach number is modest. In this case, recent high resolution shock-resolved Direct Numerical Simulations (DNS) (Ryu and Livescu, J Fluid Mech 756, R1, 2014) show that the interaction can be described by the Linear Interaction Approximation (LIA). Using LIA to alleviatemoreĀ Ā» the need to resolve the shock, DNS post-shock data can be generated at much higher Reynolds numbers than previously possible. Here, such results with Taylor Reynolds number approximately 180 are used to investigate the changes in the vortical structure as a function of the shock Mach number, Ms, up to Ms = 10. It is shown that, as Ms increases, the shock interaction induces a tendency towards a local axisymmetric state perpendicular to the shock front, which has a profound influence on the vortex-stretching mechanism and divergence of the Lamb vector and, ultimately, on the flow evolution away from the shock.Ā«Ā less
Stagnation region heat transfer augmentation at very high turbulence levels
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Kingery, Joseph E.; Ames, Forrest E.
2016-08-01
Current land-based gas turbines are growing in size producing higher approach flow Reynolds numbers at the leading edge of turbine nozzles. These vanes are subjected to high intensity large scale turbulence. This present paper reports on the research which significantly expands the parameter range for stagnation region heat transfer augmenta-tion due to high intensity turbulence. Heat transfer measurements were acquired over two constant heat flux test surfaces with large diameter leading edges (10.16 cm and 40.64 cm). The test surfaces were placed downstream from a new high intensity (17.4%) mock combustor and tested over an eight to one range inmoreĀ Ā» approach flow Reynolds number for each test surface. Stagnation region heat transfer augmentation for the smaller (ReD = 15,625ā125,000) and larger (ReD = 62,500ā500,000) leading edge regions ranged from 45% to 81% and 80% to 136%, respectively. Furthermore, these data also include heat transfer distributions over the full test surface compared with the earlier data acquired at six additional inlet turbulence conditions. These surfaces exhibit continued but more moderate acceleration downstream from the stagnation regions and these data are expected to be useful in testing bypass transition predictive approaches. This database will be useful to gas turbine heat transfer design engineers. [DOI: 10.1115/1.4032677]Ā«Ā less
The turbulent cascade and proton heating in the solar wind during solar minimum
Coburn, Jesse T.; Smith, Charles W.; Vasquez, Bernard J.; Stawarz, Joshua E.; Forman, Miriam A.
2013-06-13
Solar wind measurements at 1 AU during the recent solar minimum and previous studies of solar maximum provide an opportunity to study the effects of the changing solar cycle on in situ heating. Our interest is to compare the levels of activity associated with turbulence and proton heating. Large-scale shears in the flow caused by transient activity are a source that drives turbulence that heats the solar wind, but as the solar cycle progresses the dynamics that drive the turbulence and heat the medium are likely to change. The application of third-moment theory to Advanced Composition Explorer (ACE) data gives the turbulent energy cascade rate which is not seen to vary with the solar cycle. Likewise, an empirical heating rate shows no significan changes in proton heating over the cycle.
Bailey, James L.; Vresk, Josip
1989-01-01
A thermal transient anemometer having a thermocouple probe which is utilized to measure the change in temperature over a period of time to provide a measure of fluid flow velocity. The thermocouple probe is located in the fluid flow path and pulsed to heat or cool the probe. The cooling of the heated probe or the heating of the cooled probe from the fluid flow over a period of time is measured to determine the fluid flow velocity. The probe is desired to be locally heated near the tip to increase the efficiency of devices incorporating the probe.
Bailey, J.L.; Vresk, J.
1989-07-18
A thermal transient anemometer is disclosed having a thermocouple probe which is utilized to measure the change in temperature over a period of time to provide a measure of fluid flow velocity. The thermocouple probe is located in the fluid flow path and pulsed to heat or cool the probe. The cooling of the heated probe or the heating of the cooled probe from the fluid flow over a period of time is measured to determine the fluid flow velocity. The probe is desired to be locally heated near the tip to increase the efficiency of devices incorporating the probe. 12 figs.
Primordial magnetic field amplification from turbulent reheating
Calzetta, Esteban; Kandus, Alejandra E-mail: kandus@uesc.br
2010-08-01
We analyze the possibility of primordial magnetic field amplification by a stochastic large scale kinematic dynamo during reheating. We consider a charged scalar field minimally coupled to gravity. During inflation this field is assumed to be in its vacuum state. At the transition to reheating the state of the field changes to a many particle/anti-particle state. We characterize that state as a fluid flow of zero mean velocity but with a stochastic velocity field. We compute the scale-dependent Reynolds number Re(k), and the characteristic times for decay of turbulence, t{sub d} and pair annihilation t{sub a}, finding t{sub a} << t{sub d}. We calculate the rms value of the kinetic helicity of the flow over a scale L and show that it does not vanish. We use this result to estimate the amplification factor of a seed field from the stochastic kinematic dynamo equations. Although this effect is weak, it shows that the evolution of the cosmic magnetic field from reheating to galaxy formation may well be more complex than as dictated by simple flux freezing.
Computational analysis of a three-dimensional High-Velocity Oxygen-Fuel (HVOF) Thermal Spray torch
Hassan, B.; Lopez, A.R.; Oberkampf, W.L.
1995-07-01
An analysis of a High-Velocity Oxygen-Fuel Thermal Spray torch is presented using computational fluid dynamics (CFD). Three-dimensional CFD results are presented for a curved aircap used for coating interior surfaces such as engine cylinder bores. The device analyzed is similar to the Metco Diamond Jet Rotating Wire torch, but wire feed is not simulated. To the authors` knowledge, these are the first published 3-D results of a thermal spray device. The feed gases are injected through an axisymmetric nozzle into the curved aircap. Argon is injected through the center of the nozzle. Pre-mixed propylene and oxygen are introduced from an annulus in the nozzle, while cooling air is injected between the nozzle and the interior wall of the aircap. The combustion process is modeled assuming instantaneous chemistry. A standard, two-equation, K-{var_epsilon} turbulence model is employed for the turbulent flow field. An implicit, iterative, finite volume numerical technique is used to solve the coupled conservation of mass, momentum, and energy equations for the gas in a sequential manner. Flow fields inside and outside the aircap are presented and discussed.
PROTON KINETIC EFFECTS IN VLASOV AND SOLAR WIND TURBULENCE
Servidio, S.; Valentini, F.; Perrone, D.; Veltri, P.; Osman, K. T.; Chapman, S.; Califano, F.; Matthaeus, W. H.
2014-02-01
Kinetic plasma processes are investigated in the framework of solar wind turbulence, employing hybrid Vlasov-Maxwell (HVM) simulations. Statistical analysis of spacecraft observation data relates proton temperature anisotropy T /T {sub ā„} and parallel plasma beta Ī²{sub ā„}, where subscripts refer to the ambient magnetic field direction. Here, this relationship is recovered using an ensemble of HVM simulations. By varying plasma parameters, such as plasma beta and fluctuation level, the simulations explore distinct regions of the parameter space given by T /T {sub ā„} and Ī²{sub ā„}, similar to solar wind sub-datasets. Moreover, both simulation and solar wind data suggest that temperature anisotropy is not only associated with magnetic intermittent events, but also with gradient-type structures in the flow and in the density. This connection between non-Maxwellian kinetic effects and various types of intermittency may be a key point for understanding the complex nature of plasma turbulence.
Glenn E McCreery; Keith G Condie
2006-09-01
The Very High Temperature Reactor (VHTR) is the leading candidate for the Next Generation Nuclear Power (NGNP) Project in the U.S. which has the goal of demonstrating the production of emissions free electricity and hydrogen by 2015. The present document addresses experimental modeling of flow and thermal mixing phenomena of importance during normal or reduced power operation and during a loss of forced reactor cooling (pressurized conduction cooldown) scenario. The objectives of the experiments are, 1), provide benchmark data for assessment and improvement of codes proposed for NGNP designs and safety studies, and, 2), obtain a better understanding of related phenomena, behavior and needs. Physical models of VHTR vessel upper and lower plenums which use various working fluids to scale phenomena of interest are described. The models may be used to both simulate natural convection conditions during pressurized conduction cooldown and turbulent lower plenum flow during normal or reduced power operation.
Thermal Storage and Advanced Heat Transfer Fluids (Fact Sheet)
Not Available
2010-08-01
Fact sheet describing NREL CSP Program capabilities in the area of thermal storage and advanced heat transfer fluids: measuring thermophysical properties, measuring fluid flow and heat transfer, and simulating flow of thermal energy and fluid.
Evaluation of the Effects of Turbulence on the Behavior of Migratory Fish, 2002 Final Report.
Odeh, Mufeed.
2002-03-01
The fundamental influence of fluid dynamics on aquatic organisms is receiving increasing attention among aquatic ecologists. For example, the importance of turbulence to ocean plankton has long been a subject of investigation (Peters and Redondo 1997). More recently, studies have begun to emerge that explicitly consider the effects of shear and turbulence on freshwater invertebrates (Statzner et al. 1988; Hart et al. 1996) and fishes (Pavlov et al. 1994, 1995). Hydraulic shear stress and turbulence are interdependent natural hydraulic phenomena that are important to fish, and consequently it is important to develop an understanding of how fish sense, react to, and perhaps utilize these phenomena under normal river flows. The appropriate reaction to turbulence may promote movement of migratory fish (Coutant 1998) or prevent displacement of resident fish. It has been suggested that one of the adverse effects of flow regulation by hydroelectric projects is the reduction of normal turbulence, particularly in the headwaters of reservoirs, which can lead to disorientation and slowing of migration (Williams et al. 1996; Coutant et al. 1997; Coutant 1998). On the other hand, greatly elevated levels of shear and turbulence may be injurious to fish; injuries can range from removal of the mucous layer on the body surface to descaling to torn opercula, popped eyes, and decapitation (Neitzel et al. 2000a,b). Damaging levels of fluid stress, such turbulence, can occur in a variety of circumstances in both natural and man-made environments. This report discusses the effects of shear stress and turbulence on fish, with an emphasis on potentially damaging levels in man-made environments. It defines these phenomena, describes studies that have been conducted to understand their effects, and identifies gaps in our knowledge. In particular, this report reviews the available information on the levels of turbulence that can occur within hydroelectric power plants, and the associated
Turbulent equipartitions in two dimensional drift convection
Isichenko, M.B.; Yankov, V.V.
1995-07-25
Unlike the thermodynamic equipartition of energy in conservative systems, turbulent equipartitions (TEP) describe strongly non-equilibrium systems such as turbulent plasmas. In turbulent systems, energy is no longer a good invariant, but one can utilize the conservation of other quantities, such as adiabatic invariants, frozen-in magnetic flux, entropy, or combination thereof, in order to derive new, turbulent quasi-equilibria. These TEP equilibria assume various forms, but in general they sustain spatially inhomogeneous distributions of the usual thermodynamic quantities such as density or temperature. This mechanism explains the effects of particle and energy pinch in tokamaks. The analysis of the relaxed states caused by turbulent mixing is based on the existence of Lagrangian invariants (quantities constant along fluid-particle or other orbits). A turbulent equipartition corresponds to the spatially uniform distribution of relevant Lagrangian invariants. The existence of such turbulent equilibria is demonstrated in the simple model of two dimensional electrostatically turbulent plasma in an inhomogeneous magnetic field. The turbulence is prescribed, and the turbulent transport is assumed to be much stronger than the classical collisional transport. The simplicity of the model makes it possible to derive the equations describing the relaxation to the TEP state in several limits.
Catalytic reaction in confined flow channel
Van Hassel, Bart A.
2016-03-29
A chemical reactor comprises a flow channel, a source, and a destination. The flow channel is configured to house at least one catalytic reaction converting at least a portion of a first nanofluid entering the channel into a second nanofluid exiting the channel. The flow channel includes at least one turbulating flow channel element disposed axially along at least a portion of the flow channel. A plurality of catalytic nanoparticles is dispersed in the first nanofluid and configured to catalytically react the at least one first chemical reactant into the at least one second chemical reaction product in the flow channel.
Recurrent flow analysis in spatiotemporally chaotic 2-dimensional Kolmogorov flow
Lucas, Dan Kerswell, Rich R.
2015-04-15
Motivated by recent success in the dynamical systems approach to transitional flow, we study the efficiency and effectiveness of extracting simple invariant sets (recurrent flows) directly from chaotic/turbulent flows and the potential of these sets for providing predictions of certain statistics of the flow. Two-dimensional Kolmogorov flow (the 2D Navier-Stokes equations with a sinusoidal body force) is studied both over a square [0, 2?]{sup 2} torus and a rectangular torus extended in the forcing direction. In the former case, an order of magnitude more recurrent flows are found than previously [G. J. Chandler and R. R. Kerswell, Invariant recurrent solutions embedded in a turbulent two-dimensional Kolmogorov flow, J. Fluid Mech. 722, 554595 (2013)] and shown to give improved predictions for the dissipation and energy pdfs of the chaos via periodic orbit theory. Analysis of the recurrent flows shows that the energy is largely trapped in the smallest wavenumbers through a combination of the inverse cascade process and a feature of the advective nonlinearity in 2D. Over the extended torus at low forcing amplitudes, some extracted states mimic the statistics of the spatially localised chaos present surprisingly well recalling the findings of Kawahara and Kida [Periodic motion embedded in plane Couette turbulence: Regeneration cycle and burst, J. Fluid Mech. 449, 291 (2001)] in low-Reynolds-number plane Couette flow. At higher forcing amplitudes, however, success is limited highlighting the increased dimensionality of the chaos and the need for larger data sets. Algorithmic developments to improve the extraction procedure are discussed.
Chang, Chih-Wei; Majumdar, Arunava; Zettl, Alexander K.
2014-07-15
Disclosed is a device whereby the thermal conductance of a multiwalled nanostructure such as a multiwalled carbon nanotube (MWCNT) can be controllably and reversibly tuned by sliding one or more outer shells with respect to the inner core. As one example, the thermal conductance of an MWCNT dropped to 15% of the original value after extending the length of the MWCNT by 190 nm. The thermal conductivity returned when the tube was contracted. The device may comprise numbers of multiwalled nanotubes or other graphitic layers connected to a heat source and a heat drain and various means for tuning the overall thermal conductance for applications in structure heat management, heat flow in nanoscale or microscale devices and thermal logic devices.
Wilson, W.R.; Chapman, D.S.
1980-10-01
Separate abstracts were prepared for the three reports included in this volume on the interpretation of heat flow data in a geothermal area. (MHR)
Peterson, J. L.; Hammet, G. W.; Mikkelsen, D. R.; Yuh, H. Y.; Candy, J.; Guttenfelder, W.; Kaye, S. M.; LeBlanc, B.
2011-05-11
The first nonlinear gyrokinetic simulations of electron internal transport barriers (e-ITBs) in the National Spherical Torus Experiment show that reversed magnetic shear can suppress thermal transport by increasing the nonlinear critical gradient for electron-temperature-gradient-driven turbulence to three times its linear critical value. An interesting feature of this turbulence is non- linearly driven off-midplane radial streamers. This work reinforces the experimental observation that magnetic shear is likely an effective way of triggering and sustaining e-ITBs in magnetic fusion devices.
Turbulence and waves in the solar wind
Roberts, D.A.; Goldstein, M.L. )
1991-01-01
Studies of turbulence and waves in the solar wind is discussed. Consideration is given to the observations and theory concerning the origin and evolution of interplanetary MHD fluctuations and to the observations, theory, and simulations of compressive fluctuations. Particular attention is given to extrapolations to near-sun and polar fields regions. Results obtained on turbulence at comets and magnetic turbulence of low-frequency waves excited by unstable distributions of ions are discussed. 230 refs.