Numerical study of turbulence-influence mechanism on arc characteristics in an air direct current circuit breaker

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wu, Mingliang; Yang, Fei, E-mail: yfei2007@mail.xjtu.edu.cn; Rong, Mingzhe

This paper focuses on the numerical investigation of arc characteristics in an air direct current circuit breaker (air DCCB). Using magneto-hydrodynamics (MHD) theory, 3D laminar model and turbulence model are constructed and calculated. The standard k-epsilon model is utilized to consider the turbulence effect in the arc chamber of the DCCB. Several important phenomena are found: the arc column in the turbulence-model case is more extensive, moves much more slowly than the counterpart in the laminar-model case, and shows stagnation at the entrance of the chamber, unlike in the laminar-model case. Moreover, the arc voltage in the turbulence-model case ismore » much lower than in the laminar-model case. However, the results in the turbulence-model case show a much better agreement with the results of the breaking experiments under DC condition than in the laminar-model case, which is contradictory to the previous conclusions from the arc researches of both the low-voltage circuit breaker and the sulfur hexafluoride (SF6) nozzle. First, in the previous air-arc research of the low-voltage circuit breaker, it is assumed that the air plasma inside the chamber is in the state of laminar, and the laminar-model application gives quite satisfactory results compared with the experiments, while in this paper, the laminar-model application works badly. Second, the turbulence-model application in the arc research of the SF6-nozzle performs much better and gives higher arc voltage than the laminar-model application does, whereas in this paper, the turbulence-model application predicts lower arc voltage than the laminar-model application does. Based on the analysis of simulation results in detail, the mechanism of the above phenomena is revealed. The transport coefficients are strongly changed by turbulence, which will enhance the arc diffusion and make the arc volume much larger. Consequently, the arc appearance and the distribution of Lorentz force in the turbulence-model case substantially differ from the arc appearance and the distribution of Lorentz force in the laminar-model case. Thus, the moving process of the arc in the turbulence-model case is slowed down and slower than in the laminar-model case. Moreover, the more extensive arc column in the turbulence-model case reduces the total arc resistance, which results in a lower arc voltage, more consistent with the experimental results than the arc voltage in the laminar-model case. Therefore, the air plasma inside this air DCCB is believed to be in the turbulence state, and the turbulence model is more suitable than the laminar model for the arc simulation of this kind of air DCCB.« less

Nonaxisymmetric modelling in BOUT++; toward global edge fluid turbulence in stellarators

NASA Astrophysics Data System (ADS)

Shanahan, Brendan; Hill, Peter; Dudson, Ben

2016-10-01

As Wendelstein 7-X has been optimized for neoclassical transport, turbulent transport could potentially become comparable to neoclassical losses. Furthermore, the imminent installation of an island divertor merits global edge modelling to determine heat flux profiles and the efficacy of the system. Currently, however, nonaxisymmetric edge plasma modelling is limited to either steady state (non-turbulent) transport modelling, or computationally expensive gyrokinetics. The implementation of the Flux Coordinate Independent (FCI) approach to parallel derivatives has allowed the extension of the BOUT++ edge fluid turbulence framework to nonaxisymmetric geometries. Here we first investigate the implementation of the FCI method in BOUT++ by modelling diffusion equations in nonaxisymmetric geometries with and without boundary interaction, and quantify the inherent error. We then present the results of non-turbulent transport modelling and compare with analytical theory. The ongoing extension of BOUT++ to nonaxisymmetric configurations, and the prospects of stellarator edge fluid turbulence simulations will be discussed.

An assessment and application of turbulence models for hypersonic flows

NASA Technical Reports Server (NTRS)

Coakley, T. J.; Viegas, J. R.; Huang, P. G.; Rubesin, M. W.

1990-01-01

The current approach to the Accurate Computation of Complex high-speed flows is to solve the Reynolds averaged Navier-Stokes equations using finite difference methods. An integral part of this approach consists of development and applications of mathematical turbulence models which are necessary in predicting the aerothermodynamic loads on the vehicle and the performance of the propulsion plant. Computations of several high speed turbulent flows using various turbulence models are described and the models are evaluated by comparing computations with the results of experimental measurements. The cases investigated include flows over insulated and cooled flat plates with Mach numbers ranging from 2 to 8 and wall temperature ratios ranging from 0.2 to 1.0. The turbulence models investigated include zero-equation, two-equation, and Reynolds-stress transport models.

Algorithm for Simulating Atmospheric Turbulence and Aeroelastic Effects on Simulator Motion Systems

NASA Technical Reports Server (NTRS)

Ercole, Anthony V.; Cardullo, Frank M.; Kelly, Lon C.; Houck, Jacob A.

2012-01-01

Atmospheric turbulence produces high frequency accelerations in aircraft, typically greater than the response to pilot input. Motion system equipped flight simulators must present cues representative of the aircraft response to turbulence in order to maintain the integrity of the simulation. Currently, turbulence motion cueing produced by flight simulator motion systems has been less than satisfactory because the turbulence profiles have been attenuated by the motion cueing algorithms. This report presents a new turbulence motion cueing algorithm, referred to as the augmented turbulence channel. Like the previous turbulence algorithms, the output of the channel only augments the vertical degree of freedom of motion. This algorithm employs a parallel aircraft model and an optional high bandwidth cueing filter. Simulation of aeroelastic effects is also an area where frequency content must be preserved by the cueing algorithm. The current aeroelastic implementation uses a similar secondary channel that supplements the primary motion cue. Two studies were conducted using the NASA Langley Visual Motion Simulator and Cockpit Motion Facility to evaluate the effect of the turbulence channel and aeroelastic model on pilot control input. Results indicate that the pilot is better correlated with the aircraft response, when the augmented channel is in place.

Using Reconstructed POD Modes as Turbulent Inflow for LES Wind Turbine Simulations

NASA Astrophysics Data System (ADS)

Nielson, Jordan; Bhaganagar, Kiran; Juttijudata, Vejapong; Sirisup, Sirod

2016-11-01

Currently, in order to get realistic atmospheric effects of turbulence, wind turbine LES simulations require computationally expensive precursor simulations. At times, the precursor simulation is more computationally expensive than the wind turbine simulation. The precursor simulations are important because they capture turbulence in the atmosphere and as stated above, turbulence impacts the power production estimation. On the other hand, POD analysis has been shown to be capable of capturing turbulent structures. The current study was performed to determine the plausibility of using lower dimension models from POD analysis of LES simulations as turbulent inflow to wind turbine LES simulations. The study will aid the wind energy community by lowering the computational cost of full scale wind turbine LES simulations, while maintaining a high level of turbulent information and being able to quickly apply the turbulent inflow to multi turbine wind farms. This will be done by comparing a pure LES precursor wind turbine simulation with simulations that use reduced POD mod inflow conditions. The study shows the feasibility of using lower dimension models as turbulent inflow of LES wind turbine simulations. Overall the power production estimation and velocity field of the wind turbine wake are well captured with small errors.

The evolution equation for the flame surface density in turbulent premixed combustion

NASA Technical Reports Server (NTRS)

Trouve, Arnaud

1993-01-01

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Modeling of Turbulence Effects on Liquid Jet Atomization and Breakup

NASA Technical Reports Server (NTRS)

Trinh, Huu P.; Chen, C. P.

2005-01-01

Recent experimental investigations and physical modeling studies have indicated that turbulence behaviors within a liquid jet have considerable effects on the atomization process. This study aims to model the turbulence effect in the atomization process of a cylindrical liquid jet. Two widely used models, the Kelvin-Helmholtz (KH) instability of Reitz (blob model) and the Taylor-Analogy-Breakup (TAB) secondary droplet breakup by O Rourke et al, are further extended to include turbulence effects. In the primary breakup model, the level of the turbulence effect on the liquid breakup depends on the characteristic scales and the initial flow conditions. For the secondary breakup, an additional turbulence force acted on parent drops is modeled and integrated into the TAB governing equation. The drop size formed from this breakup regime is estimated based on the energy balance before and after the breakup occurrence. This paper describes theoretical development of the current models, called "T-blob" and "T-TAB", for primary and secondary breakup respectivety. Several assessment studies are also presented in this paper.

The effect of capturing the correct turbulence dissipation rate in BHR

DOE Office of Scientific and Technical Information (OSTI.GOV)

Schwarzkopf, John Dennis; Ristorcelli, Raymond

In this manuscript, we discuss the shortcoming of a quasi-equilibrium assumption made in the BHR closure model. Turbulence closure models generally assume fully developed turbulence, which is not applicable to 1) non-equilibrium turbulence (e.g. change in mean pressure gradient) or 2) laminar-turbulence transition flows. Based on DNS data, we show that the current BHR dissipation equation [modeled based on the fully developed turbulence phenomenology] does not capture important features of nonequilibrium flows. To demonstrate our thesis, we use the BHR equations to predict a non-equilibrium flow both with the BHR dissipation and the dissipation from DNS. We find that themore » prediction can be substantially improved, both qualitatively and quantitatively, with the correct dissipation rate. We conclude that a new set of nonequilibrium phenomenological assumptions must be used to develop a new model equation for the dissipation to accurately predict the turbulence time scale used by other models.« less

Multi-Scale Modeling and the Eddy-Diffusivity/Mass-Flux (EDMF) Parameterization

NASA Astrophysics Data System (ADS)

Teixeira, J.

2015-12-01

Turbulence and convection play a fundamental role in many key weather and climate science topics. Unfortunately, current atmospheric models cannot explicitly resolve most turbulent and convective flow. Because of this fact, turbulence and convection in the atmosphere has to be parameterized - i.e. equations describing the dynamical evolution of the statistical properties of turbulence and convection motions have to be devised. Recently a variety of different models have been developed that attempt at simulating the atmosphere using variable resolution. A key problem however is that parameterizations are in general not explicitly aware of the resolution - the scale awareness problem. In this context, we will present and discuss a specific approach, the Eddy-Diffusivity/Mass-Flux (EDMF) parameterization, that not only is in itself a multi-scale parameterization but it is also particularly well suited to deal with the scale-awareness problems that plague current variable-resolution models. It does so by representing small-scale turbulence using a classic Eddy-Diffusivity (ED) method, and the larger-scale (boundary layer and tropospheric-scale) eddies as a variety of plumes using the Mass-Flux (MF) concept.

Center for modeling of turbulence and transition: Research briefs, 1993

NASA Technical Reports Server (NTRS)

Liou, William W. (Editor)

1994-01-01

This research brief contains the progress reports of the research staff of the Center for Modeling of Turbulence and Transition (CMOTT) from June 1992 to July 1993. It is also an annual report to the Institute for Computational Mechanics in Propulsion located at Ohio Aerospace Institute and NASA Lewis Research Center. The main objectives of the research activities at CMOTT are to develop, validate, and implement turbulence and transition models for flows of interest in propulsion systems. Currently, our research covers eddy viscosity one- and two-equation models, Reynolds-stress algebraic equation models, Reynolds-stress transport equation models, nonequilibrium multiple-scale models, bypass transition models, joint scalar probability density function models, and Renormalization Group Theory and Direct Interaction Approximation methods. Some numerical simulations (LES and DNS) have also been carried out to support the development of turbulence modeling. Last year was CMOTT's third year in operation. During this period, in addition to the above mentioned research, CMOTT has also hosted the following programs: an eighteen-hour short course on 'Turbulence--Fundamentals and Computational Modeling (Part I)' given by CMOTT at the NASA Lewis Research Center; a productive summer visitor research program that has generated many encouraging results; collaborative programs with industry customers to help improve their turbulent flow calculations for propulsion system designs; a biweekly CMOTT seminar series with speakers from within and without the NASA Lewis Research Center including foreign speakers. In addition, CMOTT members have been actively involved in the national and international turbulence research activities. The current CMOTT roster and organization are listed in Appendix A. Listed in Appendix B are the abstracts of the biweekly CMOTT seminar. Appendix C lists the papers contributed by CMOTT members.

Mass and momentum turbulent transport experiments with confined swirling coaxial jets

NASA Technical Reports Server (NTRS)

Roback, R.; Johnson, B. V.

1983-01-01

Swirling coaxial jets mixing downstream, discharging into an expanded duct was conducted to obtain data for the evaluation and improvement of turbulent transport models currently used in a variety of computational procedures throughout the combustion community. A combination of laser velocimeter (LV) and laser induced fluorescence (LIF) techniques was employed to obtain mean and fluctuating velocity and concentration distributions which were used to derive mass and momentum turbulent transport parameters currently incorporated into various combustor flow models. Flow visualization techniques were also employed to determine qualitatively the time dependent characteristics of the flow and the scale of turbulence. The results of these measurements indicated that the largest momentum turbulent transport was in the r-z plane. Peak momentum turbulent transport rates were approximately the same as those for the nonswirling flow condition. The mass turbulent transport process for swirling flow was complicated. Mixing occurred in several steps of axial and radial mass transport and was coupled with a large radial mean convective flux. Mixing for swirling flow was completed in one-third the length required for nonswirling flow.

Automated turbulence forecasts for aviation hazards

NASA Astrophysics Data System (ADS)

Sharman, R.; Frehlich, R.; Vandenberghe, F.

2010-09-01

An operational turbulence forecast system for commercial and aviation use is described that is based on an ensemble of turbulence diagnostics derived from standard NWP model outputs. In the U. S. this forecast product is named GTG (Graphical Turbulence Guidance) and has been described in detail in Sharman et al., WAF 2006. Since turbulence has many sources in the atmosphere, the ensemble approach of combining diagnostics has been shown to provide greater statistical accuracy than the use of a single diagnostic, or of a subgrid tke parameterization. GTG is sponsored by the FAA, and has undergone rigorous accuracy, safety, and usability evaluations. The GTG product is now hosted on NOAA's Aviation Data Service (ADDS), web site (http://aviationweather.gov/), for access by pilots, air traffic controllers, and dispatchers. During this talk the various turbulence diagnostics, their statistical properties, and their relative performance (based on comparisons to observations) will be presented. Importantly, the model output is ɛ1/3 (where ɛ is the eddy dissipation rate), so is aircraft independent. The diagnostics are individually and collectively calibrated so that their PDFs satisfy the expected log normal distribution of ɛ^1/3. Some of the diagnostics try to take into account the role of gravity waves and inertia-gravity waves in the turbulence generation process. Although the current GTG product is based on the RUC forecast model running over the CONUS, it is transitioning to a WRF based product, and in fact WRF-based versions are currently running operationally over Taiwan and has also been implemented for use by the French Navy in climatological studies. Yet another version has been developed which uses GFS model output to provide global turbulence forecasts. Thus the forecast product is available as a postprocessing program for WRF or other model output and provides 3D maps of turbulence likelihood of any region where NWP model data is available. Although the current GTG has been used mainly for large commercial aircraft, since the output is aircraft independent it could readily be scaled to smaller aircraft such as UAVs. Further, the ensemble technique allows the diagnostics to be used to form probabilistic forecasts, in a manner similar to ensemble NWP forecasts.

Status of Turbulence Modeling for Hypersonic Propulsion Flowpaths

NASA Technical Reports Server (NTRS)

Georgiadis, Nicholas J.; Yoder, Dennis A.; Vyas, Manan A.; Engblom, William A.

2012-01-01

This report provides an assessment of current turbulent flow calculation methods for hypersonic propulsion flowpaths, particularly the scramjet engine. Emphasis is placed on Reynolds-averaged Navier-Stokes (RANS) methods, but some discussion of newer meth- ods such as Large Eddy Simulation (LES) is also provided. The report is organized by considering technical issues throughout the scramjet-powered vehicle flowpath including laminar-to-turbulent boundary layer transition, shock wave / turbulent boundary layer interactions, scalar transport modeling (specifically the significance of turbulent Prandtl and Schmidt numbers) and compressible mixing. Unit problems are primarily used to conduct the assessment. In the combustor, results from calculations of a direct connect supersonic combustion experiment are also used to address the effects of turbulence model selection and in particular settings for the turbulent Prandtl and Schmidt numbers. It is concluded that RANS turbulence modeling shortfalls are still a major limitation to the accuracy of hypersonic propulsion simulations, whether considering individual components or an overall system. Newer methods such as LES-based techniques may be promising, but are not yet at a maturity to be used routinely by the hypersonic propulsion community. The need for fundamental experiments to provide data for turbulence model development and validation is discussed.

Status of turbulence modeling for hypersonic propulsion flowpaths

NASA Astrophysics Data System (ADS)

Georgiadis, Nicholas J.; Yoder, Dennis A.; Vyas, Manan A.; Engblom, William A.

2014-06-01

This report provides an assessment of current turbulent flow calculation methods for hypersonic propulsion flowpaths, particularly the scramjet engine. Emphasis is placed on Reynolds-averaged Navier-Stokes (RANS) methods, but some discussion of newer methods such as large eddy simulation (LES) is also provided. The report is organized by considering technical issues throughout the scramjet-powered vehicle flowpath, including laminar-to-turbulent boundary layer transition, shock wave/turbulent boundary layer interactions, scalar transport modeling (specifically the significance of turbulent Prandtl and Schmidt numbers), and compressible mixing. Unit problems are primarily used to conduct the assessment. In the combustor, results from calculations of a direct connect supersonic combustion experiment are also used to address the effects of turbulence model selection and in particular settings for the turbulent Prandtl and Schmidt numbers. It is concluded that RANS turbulence modeling shortfalls are still a major limitation to the accuracy of hypersonic propulsion simulations, whether considering individual components or an overall system. Newer methods such as LES-based techniques may be promising, but are not yet at a maturity to be used routinely by the hypersonic propulsion community. The need for fundamental experiments to provide data for turbulence model development and validation is discussed.

Dissipation Rate of Turbulent Kinetic Energy in Diel Vertical Migrations: Comparison of ANSYS Fluent Model to Measurements

NASA Astrophysics Data System (ADS)

Dean, Cayla; Soloviev, Alexander; Hirons, Amy; Frank, Tamara; Wood, Jon

2015-04-01

Recent studies suggest that diel vertical migrations of zooplankton may have an impact on ocean mixing, though details are not completely clear. A strong sound scattering layer of zooplankton undergoing diel vertical migrations was observed in Saanich Inlet, British Colombia, Canada by Kunze et al. (2006). In this study, a shipboard 200-kHz echosounder was used to track vertical motion of the sound scattering layer, and microstructure profiles were collected to observe turbulence. An increase of dissipation rate of turbulent kinetic energy by four to five orders of magnitude was measured during diel vertical migrations of zooplankton in one case (but not observed during other cases). A strong sound scattering layer undergoing diel vertical migration was also observed in the Straits of Florida via a bottom mounted acoustic Doppler current profiler at 244 m isobath. A 3-D non-hydrostatic computational fluid dynamics model with Lagrangian particle injections (a proxy for migrating zooplankton) via a discrete phase model was used to simulate the effect of diel vertical migrations on the turbulence for both Saanich Inlet and the Straits of Florida. The model was initialized with idealized (but based on observation) density and velocity profiles. Particles, with buoyancy adjusted to serve as a proxy for vertically swimming zooplankton, were injected to simulate diel vertical migration cycles. Results of models run with extreme concentrations of particles showed an increase in dissipation rate of turbulent kinetic energy of approximately five orders of magnitude over background turbulence during migration of particles in both Saanich Inlet and the Straits of Florida cases (though direct relation of the turbulence produced by buoyant particles and swimming organisms isn't straightforward). This increase was quantitatively consistent, with turbulence measurements by Kunze et al. (2006). When 10 times fewer particles were injected into the model, the effect on dissipation rate of turbulent kinetic energy was an order of magnitude smaller than that from the extreme concentration. At a concentration of particles 100 times smaller than the extreme concentration, there was no longer an observable effect. In the Straits of Florida, direct turbulence measurements were not available to make a quantitative comparison. However, a small, but statistically significant decrease in northward current velocity profiles during migration times were observed after averaging these profiles over 11 months. A small decrease of current velocity connected to the vertical migrations of particles was reproduced in the Straits of Florida model case. The deviations in the velocity profiles can be explained by the increase in turbulent mixing during vertical migration periods.

Turbulent mixing in stably stratified flows of the environment: The current state of the problem (Review)

NASA Astrophysics Data System (ADS)

Vasil'Ev, O. F.; Voropaeva, O. F.; Kurbatskii, A. F.

2011-06-01

Specific features of the turbulent transfer of the momentum and heat in stably stratified geophysical flows, as well as possibilities for including them into RANS turbulence models, are analyzed. The momentum (but not heat) transfer by internal gravity waves under conditions of strong stability is, for example, one such feature. Laboratory data and measurements in the atmosphere fix a clear dropping trend of the inverse turbulent Prandtl number with an increasing gradient Richardson number, which must be reproduced by turbulence models. Ignoring this feature can cause a false diffusion of heat under conditions of strong stability and lead, in particular, to noticeable errors in calculations of the temperature in the atmospheric boundary layer. Therefore, models of turbulent transfer must include the effect of the action of buoyancy and internal gravity waves on turbulent flows of the momentum. Such a strategy of modeling the stratified turbulence is presented in the review by a concrete RANS model and original results obtained during the modeling of stratified flows in the environment. Semiempirical turbulence models used for calculations of complex turbulent flows in deep stratified bodies of water are also analyzed. This part of the review is based on the data of investigations within the framework of the large international scientific Comparative Analysis and Rationalization of Second-Moment Turbulence Models (CARTUM) project and other publications of leading specialists. The most economical and effective approach associated with modified two-parameter turbulence models is a real alternative to classical variants of these models. A class of test problems and laboratory and full-scale experiments used by the participants of the CARTUM project for the approbation of numerical models are considered.

Development of a two-equation turbulence model for hypersonic flows. Volume 1; Evaluation of a low Reynolds number correction to the Kappa - epsilon two equation compressible turbulence model

NASA Technical Reports Server (NTRS)

Knight, Doyle D.; Becht, Robert J.

1995-01-01

The objective of the current research is the development of an improved k-epsilon two-equation compressible turbulence model for turbulent boundary layer flows experiencing strong viscous-inviscid interactions. The development of an improved model is important in the design of hypersonic vehicles such as the National Aerospace Plane (NASP) and the High Speed Civil Transport (HSCT). Improvements have been made to the low Reynolds number functions in the eddy viscosity and dissipation of solenoidal dissipation of the k-epsilon turbulence mode. These corrections offer easily applicable modifications that may be utilized for more complex geometries. The low Reynolds number corrections are functions of the turbulent Reynolds number and are therefore independent of the coordinate system. The proposed model offers advantages over some current models which are based upon the physical distance from the wall, that modify the constants of the standard model, or that make more corrections than are necessary to the governing equations. The code has been developed to solve the Favre averaged, boundary layer equations for mass, momentum, energy, turbulence kinetic energy, and dissipation of solenoidal dissipation using Keller's box scheme and the Newton spatial marching method. The code has been validated by removing the turbulent terms and comparing the solution with the Blasius solution, and by comparing the turbulent solution with an existing k-epsilon model code using wall function boundary conditions. Excellent agreement is seen between the computed solution and the Blasius solution, and between the two codes. The model has been tested for both subsonic and supersonic flat-plate turbulent boundary layer flow by comparing the computed skin friction with the Van Driest II theory and the experimental data of Weighardt; by comparing the transformed velocity profile with the data of Weighardt, and the Law of the Wall and the Law of the Wake; and by comparing the computed results of an adverse pressure gradient with the experimental data of Fernando and Smits. Good agreement is obtained with the experimental correlations for all flow conditions.

The evolution equation for the flame surface density in turbulent premixed combustion

NASA Technical Reports Server (NTRS)

Trouve, A.; Poinsot, T.

1992-01-01

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Electron beam generation in the turbulent plasma of Z-pinch discharges

NASA Astrophysics Data System (ADS)

Vikhrev, Victor V.; Baronova, Elena O.

1997-05-01

Numerical modeling of the process of electron beam generation in z-pinch discharges are presented. The proposed model represents the electron beam generation under turbulent plasma conditions. Strong current distribution inhomogeneity in the plasma column has been accounted for the adequate generation process investigation. Electron beam is generated near the maximum of compression due to run away mechanism and it is not related with the current break effect.

Electron magnetic reconnection without ion coupling in Earth's turbulent magnetosheath

NASA Astrophysics Data System (ADS)

Phan, T. D.; Eastwood, J. P.; Shay, M. A.; Drake, J. F.; Sonnerup, B. U. Ö.; Fujimoto, M.; Cassak, P. A.; Øieroset, M.; Burch, J. L.; Torbert, R. B.; Rager, A. C.; Dorelli, J. C.; Gershman, D. J.; Pollock, C.; Pyakurel, P. S.; Haggerty, C. C.; Khotyaintsev, Y.; Lavraud, B.; Saito, Y.; Oka, M.; Ergun, R. E.; Retino, A.; Le Contel, O.; Argall, M. R.; Giles, B. L.; Moore, T. E.; Wilder, F. D.; Strangeway, R. J.; Russell, C. T.; Lindqvist, P. A.; Magnes, W.

2018-05-01

Magnetic reconnection in current sheets is a magnetic-to-particle energy conversion process that is fundamental to many space and laboratory plasma systems. In the standard model of reconnection, this process occurs in a minuscule electron-scale diffusion region1,2. On larger scales, ions couple to the newly reconnected magnetic-field lines and are ejected away from the diffusion region in the form of bi-directional ion jets at the ion Alfvén speed3-5. Much of the energy conversion occurs in spatially extended ion exhausts downstream of the diffusion region6. In turbulent plasmas, which contain a large number of small-scale current sheets, reconnection has long been suggested to have a major role in the dissipation of turbulent energy at kinetic scales7-11. However, evidence for reconnection plasma jetting in small-scale turbulent plasmas has so far been lacking. Here we report observations made in Earth's turbulent magnetosheath region (downstream of the bow shock) of an electron-scale current sheet in which diverging bi-directional super-ion-Alfvénic electron jets, parallel electric fields and enhanced magnetic-to-particle energy conversion were detected. Contrary to the standard model of reconnection, the thin reconnecting current sheet was not embedded in a wider ion-scale current layer and no ion jets were detected. Observations of this and other similar, but unidirectional, electron jet events without signatures of ion reconnection reveal a form of reconnection that can drive turbulent energy transfer and dissipation in electron-scale current sheets without ion coupling.

Electron magnetic reconnection without ion coupling in Earth's turbulent magnetosheath.

PubMed

Phan, T D; Eastwood, J P; Shay, M A; Drake, J F; Sonnerup, B U Ö; Fujimoto, M; Cassak, P A; Øieroset, M; Burch, J L; Torbert, R B; Rager, A C; Dorelli, J C; Gershman, D J; Pollock, C; Pyakurel, P S; Haggerty, C C; Khotyaintsev, Y; Lavraud, B; Saito, Y; Oka, M; Ergun, R E; Retino, A; Le Contel, O; Argall, M R; Giles, B L; Moore, T E; Wilder, F D; Strangeway, R J; Russell, C T; Lindqvist, P A; Magnes, W

2018-05-01

Magnetic reconnection in current sheets is a magnetic-to-particle energy conversion process that is fundamental to many space and laboratory plasma systems. In the standard model of reconnection, this process occurs in a minuscule electron-scale diffusion region 1,2 . On larger scales, ions couple to the newly reconnected magnetic-field lines and are ejected away from the diffusion region in the form of bi-directional ion jets at the ion Alfvén speed 3-5 . Much of the energy conversion occurs in spatially extended ion exhausts downstream of the diffusion region 6 . In turbulent plasmas, which contain a large number of small-scale current sheets, reconnection has long been suggested to have a major role in the dissipation of turbulent energy at kinetic scales 7-11 . However, evidence for reconnection plasma jetting in small-scale turbulent plasmas has so far been lacking. Here we report observations made in Earth's turbulent magnetosheath region (downstream of the bow shock) of an electron-scale current sheet in which diverging bi-directional super-ion-Alfvénic electron jets, parallel electric fields and enhanced magnetic-to-particle energy conversion were detected. Contrary to the standard model of reconnection, the thin reconnecting current sheet was not embedded in a wider ion-scale current layer and no ion jets were detected. Observations of this and other similar, but unidirectional, electron jet events without signatures of ion reconnection reveal a form of reconnection that can drive turbulent energy transfer and dissipation in electron-scale current sheets without ion coupling.

On the Conditioning of Machine-Learning-Assisted Turbulence Modeling

NASA Astrophysics Data System (ADS)

Wu, Jinlong; Sun, Rui; Wang, Qiqi; Xiao, Heng

2017-11-01

Recently, several researchers have demonstrated that machine learning techniques can be used to improve the RANS modeled Reynolds stress by training on available database of high fidelity simulations. However, obtaining improved mean velocity field remains an unsolved challenge, restricting the predictive capability of current machine-learning-assisted turbulence modeling approaches. In this work we define a condition number to evaluate the model conditioning of data-driven turbulence modeling approaches, and propose a stability-oriented machine learning framework to model Reynolds stress. Two canonical flows, the flow in a square duct and the flow over periodic hills, are investigated to demonstrate the predictive capability of the proposed framework. The satisfactory prediction performance of mean velocity field for both flows demonstrates the predictive capability of the proposed framework for machine-learning-assisted turbulence modeling. With showing the capability of improving the prediction of mean flow field, the proposed stability-oriented machine learning framework bridges the gap between the existing machine-learning-assisted turbulence modeling approaches and the demand of predictive capability of turbulence models in real applications.

Effect of double layers on magnetosphere-ionosphere coupling

NASA Technical Reports Server (NTRS)

Lysak, Robert L.; Hudson, Mary K.

1987-01-01

The dynamic aspects of auroral current structures are reviewed with emphasis on consequences for models of microscopic turbulence (MT). A number of models of MT are introduced into a large-scale model of Alfven wave propagation to determine the effect of various models on the overall structure of auroral currents. The effect of a double layer (DL) electric field which scales with the plasma temperature and the Debye length is compared with the effect of anomalous resistivity due to electrostatic ion cyclotron turbulence in which the electric field scales with the magnetic field strength. It is shown that the DL model is less diffusive than the resistive model, indicating the possibility of narrow intense current structures.

Wake-Driven Dynamics of Finite-Sized Buoyant Spheres in Turbulence

NASA Astrophysics Data System (ADS)

Mathai, Varghese; Prakash, Vivek N.; Brons, Jon; Sun, Chao; Lohse, Detlef

2015-09-01

Particles suspended in turbulent flows are affected by the turbulence and at the same time act back on the flow. The resulting coupling can give rise to rich variability in their dynamics. Here we report experimental results from an investigation of finite-sized buoyant spheres in turbulence. We find that even a marginal reduction in the particle's density from that of the fluid can result in strong modification of its dynamics. In contrast to classical spatial filtering arguments and predictions of particle models, we find that the particle acceleration variance increases with size. We trace this reversed trend back to the growing contribution from wake-induced forces, unaccounted for in current particle models in turbulence. Our findings highlight the need for improved multiphysics based models that account for particle wake effects for a faithful representation of buoyant-sphere dynamics in turbulence.

Turbulence: The chief outstanding difficulty of our subject

NASA Technical Reports Server (NTRS)

Bradshaw, Peter

1992-01-01

A review of interesting current topics in turbulence research is decorated with examples of popular fallacies about the behavior of turbulence. Topics include the status of the Law of the Wall, especially in compressible flow; analogies between the effects of Reynolds numbers, pressure gradient, unsteadiness and roughness change; the status of Kolmogorov's universal equilibrium theory and local isotropy of the small eddies; turbulence modelling, with reference to universality, pressure-strain modelling and the dissipation equation; and chaos. Fallacies include the mixing-length concept; the effect of pressure gradient on Reynolds shear stress; the separability of time and space derivatives; models of the dissipation equation; and chaos.

Model for compressible turbulence in hypersonic wall boundary and high-speed mixing layers

NASA Astrophysics Data System (ADS)

Bowersox, Rodney D. W.; Schetz, Joseph A.

1994-07-01

The most common approach to Navier-Stokes predictions of turbulent flows is based on either the classical Reynolds-or Favre-averaged Navier-Stokes equations or some combination. The main goal of the current work was to numerically assess the effects of the compressible turbulence terms that were experimentaly found to be important. The compressible apparent mass mixing length extension (CAMMLE) model, which was based on measured experimental data, was found to produce accurate predictions of the measured compressible turbulence data for both the wall bounded and free mixing layer. Hence, that model was incorporated into a finite volume Navier-Stokes code.

Atmospheric turbulence simulation for Shuttle orbiter

NASA Technical Reports Server (NTRS)

Tatom, F. B.; Smith, S. R.

1979-01-01

An improved non-recursive model for atmospheric turbulence along the flight path of the Shuttle Orbiter is developed which provides for simulation of instantaneous vertical and horizontal gusts at the vehicle center-of-gravity, and also for simulation of instantaneous gust gradients. Based on this model the time series for both gusts and gust gradients are generated and stored on a series of magnetic tapes. Section 2 provides a description of the various technical considerations associated with the turbulence simulation model. Included in this section are descriptions of the digital filter simulation model, the von Karman spectra with finite upper limits, and the final non recursive turbulence simulation model which was used to generate the time series. Section 2 provides a description of the various technical considerations associated with the turbulence simulation model. Included in this section are descriptions of the digial filter simulation model, the von Karman spectra with finite upper limits, and the final non recursive turbulence simulation model which was used to generate the time series. Section 3 provides a description of the time series as currently recorded on magnetic tape. Conclusions and recommendations are presented in Section 4.

Consistency, Verification, and Validation of Turbulence Models for Reynolds-Averaged Navier-Stokes Applications

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.

2009-01-01

In current practice, it is often difficult to draw firm conclusions about turbulence model accuracy when performing multi-code CFD studies ostensibly using the same model because of inconsistencies in model formulation or implementation in different codes. This paper describes an effort to improve the consistency, verification, and validation of turbulence models within the aerospace community through a website database of verification and validation cases. Some of the variants of two widely-used turbulence models are described, and two independent computer codes (one structured and one unstructured) are used in conjunction with two specific versions of these models to demonstrate consistency with grid refinement for several representative problems. Naming conventions, implementation consistency, and thorough grid resolution studies are key factors necessary for success.

Turbulence modelling of flow fields in thrust chambers

NASA Technical Reports Server (NTRS)

Chen, C. P.; Kim, Y. M.; Shang, H. M.

1993-01-01

Following the consensus of a workshop in Turbulence Modelling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows, and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data to account for the non-isotropic turbulence effects.

Helicopter-V/STOL dynamic wind and turbulence design methodology

NASA Technical Reports Server (NTRS)

Bailey, J. Earl

1987-01-01

Aircraft and helicopter accidents due to severe dynamic wind and turbulence continue to present challenging design problems. The development of the current set of design analysis tools for a aircraft wind and turbulence design began in the 1940's and 1950's. The areas of helicopter dynamic wind and turbulence modeling and vehicle response to severe dynamic wind inputs (microburst type phenomena) during takeoff and landing remain as major unsolved design problems from a lack of both environmental data and computational methodology. The development of helicopter and V/STOL dynamic wind and turbulence response computation methology is reviewed, the current state of the design art in industry is outlined, and comments on design methodology are made which may serve to improve future flight vehicle design.

An Investigation of a Hybrid Mixing Model for PDF Simulations of Turbulent Premixed Flames

NASA Astrophysics Data System (ADS)

Zhou, Hua; Li, Shan; Wang, Hu; Ren, Zhuyin

2015-11-01

Predictive simulations of turbulent premixed flames over a wide range of Damköhler numbers in the framework of Probability Density Function (PDF) method still remain challenging due to the deficiency in current micro-mixing models. In this work, a hybrid micro-mixing model, valid in both the flamelet regime and broken reaction zone regime, is proposed. A priori testing of this model is first performed by examining the conditional scalar dissipation rate and conditional scalar diffusion in a 3-D direct numerical simulation dataset of a temporally evolving turbulent slot jet flame of lean premixed H2-air in the thin reaction zone regime. Then, this new model is applied to PDF simulations of the Piloted Premixed Jet Burner (PPJB) flames, which are a set of highly shear turbulent premixed flames and feature strong turbulence-chemistry interaction at high Reynolds and Karlovitz numbers. Supported by NSFC 51476087 and NSFC 91441202.

Sediment Transport Capacity of Turbidity Currents: from Microscale to Geological Scale.

NASA Astrophysics Data System (ADS)

Eggenhuisen, J. T.; Tilston, M.; Cartigny, M.; Pohl, F.; de Leeuw, J.; van der Grind, G. J.

2016-12-01

A big question in sedimentology concerns the magnitude of fluxes of sediment particles, solute matter and dissolved gasses from shallow marine waters to deep basins by turbidity current flow. Here we establish sediment transport capacity of turbidity current flow on three levels. The most elementary level is set by the maximum amount of sediment that can be contained at the base of turbidity currents without causing complete extinction of boundary layer turbulence. The second level concerns the capacity in a vertical column within turbidity currents. The third level involves the amount of sediment that can be transported in turbidite systems on geological timescales. The capacity parameter Γ compares turbulent forces near the boundary of a turbulent suspension to gravity and buoyancy forces acting on suspended particles. The condition of Γ>1 coincides with complete suppression of coherent boundary layer turbulence in Direct Numerical Simulations of sediment-laden turbulent flow. Γ=1 coincides with the upper limit of observed suspended particle concentrations in flume and field measurements. Γ is grainsize independent, yet capacity of the full vertical structure of turbidity currents becomes grainsize dependent. This is due to the appearance of grainsize dependent vertical motions within turbulence as a primary control on the shape of the vertical concentration profile. We illustrate this dependence with experiments and theory and conclude that capacity depends on the competence of prevailing turbulence to suspend particle sizes. The concepts of capacity and competence are thus tangled. Finally, the capacity of turbidity current flow structure is coupled to geological constraints on recurrence times, channel and lobe life cycles, and allogenic forcing on system activity to arrive at system scale sediment transport capacity. We demonstrate a simple model that uses the fundamental process insight described above to estimate geological sediment budgets from architectural information. These predictions are tied to existing S2S analyses to constrain submarine channel and fan dimensions in ancient and subsurface systems. Predictions of sediment budgets in deep marine systems rely on integration of fundamental issues in turbulent particle suspension into geological models of turbidite systems.

Implementation of the Baldwin-Barth turbulence model into the ZETA code and its diagnosis. M.S. Thesis

NASA Technical Reports Server (NTRS)

Low, Scott L.

1993-01-01

The Baldwin-Barth turbulence model was implemented into Zeta, a time-accurate, zonal, integro-differential code for incompressible laminar and turbulent flows. The implementation procedure is patterned after the model subroutine in ARC2D. The results of ZETA with the Baldwin-Barth turbulence model were compared with experimental data, with ZETA using Baldwin-Lomax model, and with ARC2D using the Baldwin-Barth model. The Baldwin-Barth model subroutine was tested by inputting an ARC2D velocity solution of an NACA-0012 airfoil at R(sub e) = 3.9 x 10(exp 6) and alpha = 5 deg. The resultant turbulent viscosity and Reynolds stresses compared favorably with the original data. For the same grid having grid points inside the laminar sublayer, which is necessary due to the one-equation nature of the model, ZETA however predicts early separation. It was found that the current ZETA has problem with such a fine grid. Further work is in progress to solve this problem.

Inference of turbulence parameters from a ROMS simulation using the k-ε closure scheme

NASA Astrophysics Data System (ADS)

Thyng, Kristen M.; Riley, James J.; Thomson, Jim

2013-12-01

Comparisons between high resolution turbulence data from Admiralty Inlet, WA (USA), and a 65-meter horizontal grid resolution simulation using the hydrostatic ocean modelling code, Regional Ocean Modeling System (ROMS), show that the model's k-ε turbulence closure scheme performs reasonably well. Turbulent dissipation rates and Reynolds stresses agree within a factor of two, on average. Turbulent kinetic energy (TKE) also agrees within a factor of two, but only for motions within the observed inertial sub-range of frequencies (i.e., classic approximately isotropic turbulence). TKE spectra from the observations indicate that there is significant energy at lower frequencies than the inertial sub-range; these scales are not captured by the model closure scheme nor the model grid resolution. To account for scales not present in the model, the inertial sub-range is extrapolated to lower frequencies and then integrated to obtain an inferred, diagnostic total TKE, with improved agreement with the observed total TKE. The realistic behavior of the dissipation rate and Reynolds stress, combined with the adjusted total TKE, imply that ROMS simulations can be used to understand and predict spatial and temporal variations in turbulence. The results are suggested for application to siting tidal current turbines.

Simple Analytical Forms of the Perpendicular Diffusion Coefficient for Two-component Turbulence. III. Damping Model of Dynamical Turbulence

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gammon, M.; Shalchi, A., E-mail: andreasm4@yahoo.com

2017-10-01

In several astrophysical applications one needs analytical forms of cosmic-ray diffusion parameters. Some examples are studies of diffusive shock acceleration and solar modulation. In the current article we explore perpendicular diffusion based on the unified nonlinear transport theory. While we focused on magnetostatic turbulence in Paper I, we included the effect of dynamical turbulence in Paper II of the series. In the latter paper we assumed that the temporal correlation time does not depend on the wavenumber. More realistic models have been proposed in the past, such as the so-called damping model of dynamical turbulence. In the present paper wemore » derive analytical forms for the perpendicular diffusion coefficient of energetic particles in two-component turbulence for this type of time-dependent turbulence. We present new formulas for the perpendicular diffusion coefficient and we derive a condition for which the magnetostatic result is recovered.« less

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reacting flows

NASA Technical Reports Server (NTRS)

Chinitz, W.; Foy, E.; Rowan, G.; Goldstein, D.

1982-01-01

The use of probability theory to determine the effects of turbulent fluctuations on reaction rates in turbulent combustion systems is briefly reviewed. Results are presented for the effect of species fluctuations in particular. It is found that turbulent fluctuations of species act to reduce the reaction rates, in contrast with the temperature fluctuations previously determined to increase Arrhenius reaction rate constants. For the temperature fluctuations, a criterion is set forth for determining if, in a given region of a turbulent flow field, the temperature can be expected to exhibit ramp like fluctuations. Using the above results, along with results previously obtained, a model is described for testing the effects of turbulent fluctuations of temperature and species on reaction rates in computer programs dealing with turbulent reacting flows. An alternative model which employs three variable probability density functions (temperature and two species) and is currently being formulated is discussed as well.

Proceedings of the U.S. Air Force and The Federal Republic of Germany Data Exchange Agreement Meeting (9th), Viscous and Interacting Flow Field Effects Held at Silver Spring, Maryland on 9-10 May 1984,

DTIC Science & Technology

1984-08-01

found in References 1-3. 2. Modeling of Roughness Effects on Turbulent Flow In turbulent flow analysis , use of time-averaged equations leads to the...eddy viscosity and the mixing length which are important parameters used in current algebraic modeling of the turbulence shear term. Two different ...surfaces with three-dimensional (distributed) roughness elements. Calculations using the present model have been compared with experimental data from

Analysis of the characteristics of DC nozzle arcs in air and guidance for the search of SF6 replacement gas

NASA Astrophysics Data System (ADS)

Liu, J.; Zhang, Q.; Yan, J. D.; Zhong, J.; Fang, M. T. C.

2016-11-01

It is shown that the arc model based on laminar flow cannot predict satisfactorily the voltage of an air arc burning in a supersonic nozzle. The Prandtl mixing length model (PML) and a modified k-epsilon turbulence model (MKE) are used to introduce turbulence enhanced momentum and energy transport. Arc voltages predicted by these two turbulence models are in good agreement with experiments at the stagnation pressure (P 0) of 10 bar. The predicted arc voltages by MKE for P 0  =  13 bar and 7 bar are in better agreement with experiments than those predicted by PML. MKE is therefore a preferred turbulence model for an air nozzle arc. There are two peaks in ρC P of air at 4000 K and 7000 K due, respectively, to the dissociation of oxygen and that of nitrogen. These peaks produce corresponding peaks in turbulent thermal conductivity, which results in very broad radial temperature profile and a large arc radius. Thus, turbulence indirectly enhances axial enthalpy transport, which becomes the dominant energy transport process for the overall energy balance of the arc column at high currents. When the current reduces, turbulent thermal conduction gradually becomes dominant. The temperature dependence of ρC P has a decisive influence on the radial temperature profile of a turbulent arc, thus the thermal interruption capability of a gas. Comparison between ρC P for air and SF6 shows that ρC P for SF6 has peaks below 4000 K. This renders a distinctive arc core and a small arc radius for turbulent SF6, thus superior arc quenching capability. It is suggested, for the first time, that ρC P provides guidance for the search of a replacement switching gas for SF6.

Vertical structure of mean cross-shore currents across a barred surf zone

USGS Publications Warehouse

Haines, John W.; Sallenger, Asbury H.

1994-01-01

Mean cross-shore currents observed across a barred surf zone are compared to model predictions. The model is based on a simplified momentum balance with a turbulent boundary layer at the bed. Turbulent exchange is parameterized by an eddy viscosity formulation, with the eddy viscosity Aυ independent of time and the vertical coordinate. Mean currents result from gradients due to wave breaking and shoaling, and the presence of a mean setup of the free surface. Descriptions of the wave field are provided by the wave transformation model of Thornton and Guza [1983]. The wave transformation model adequately reproduces the observed wave heights across the surf zone. The mean current model successfully reproduces the observed cross-shore flows. Both observations and predictions show predominantly offshore flow with onshore flow restricted to a relatively thin surface layer. Successful application of the mean flow model requires an eddy viscosity which varies horizontally across the surf zone. Attempts are made to parameterize this variation with some success. The data does not discriminate between alternative parameterizations proposed. The overall variability in eddy viscosity suggested by the model fitting should be resolvable by field measurements of the turbulent stresses. Consistent shortcomings of the parameterizations, and the overall modeling effort, suggest avenues for further development and data collection.

Effects of including electrojet turbulence in LFM-RCM simulations of geospace storms

NASA Astrophysics Data System (ADS)

Oppenheim, M. M.; Wiltberger, M. J.; Merkin, V. G.; Zhang, B.; Toffoletto, F.; Wang, W.; Lyon, J.; Liu, J.; Dimant, Y. S.

2016-12-01

Global geospace system simulations need to incorporate nonlinear and small-scale physical processes in order to accurately model storms and other intense events. During times of strong magnetospheric disturbances, large-amplitude electric fields penetrate from the Earth's magnetosphere to the E-region ionosphere where they drive Farley-Buneman instabilities (FBI) that create small-scale plasma density turbulence. This induces nonlinear currents and leads to anomalous electron heating. Current global Magnetosphere-Ionosphere-Thermosphere (MIT) models disregard these effects by assuming simple laminar ionospheric currents. This paper discusses the effects of incorporating accurate turbulent conductivities into MIT models. Recently, we showed in Liu et al. (2016) that during storm-time, turbulence increases the electron temperatures and conductivities more than precipitation. In this talk, we present the effect of adding these effects to the combined Lyon-Fedder-Mobarry (LFM) global MHD magnetosphere simulator and the Rice Convection Model (RCM). The LFM combines a magnetohydrodynamic (MHD) simulation of the magnetosphere with a 2D electrostatic solution of the ionosphere. The RCM uses drift physics to accurately model the inner magnetosphere, including a storm enhanced ring current. The LFM and coupled LFM-RCM simulations have previously shown unrealistically high cross-polar-cap potentials during strong solar wind driving conditions. We have recently implemented an LFM module that modifies the ionospheric conductivity to account for FBI driven anomalous electron heating and non-linear cross-field current enhancements as a function of the predicted ionospheric electric field. We have also improved the LFM-RCM code by making it capable of handling dipole tilts and asymmetric ionospheric solutions. We have tested this new LFM version by simulating the March 17, 2013 geomagnetic storm. These simulations showed a significant reduction in the cross-polar-cap potential during the strongest driving conditions, significant increases in the ionospheric conductivity in the auroral oval, and better agreement with DMSP observations of sub-auroral polarization streams. We conclude that accurate MIT simulations of geospace storms require the inclusion of turbulent conductivities.

Exploring a Method for Improving Turbulent Separated-Flow Predictions with kappa-omega Models

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.

2009-01-01

A particular failing of Reynolds-averaged Navier-Stokes separated turbulent flow computations is addressed within the context of a kappa-omega two-equation turbulence model. The failing is the tendency for turbulence models to under-predict turbulent shear stress in the shear layers of some separation bubbles, yielding late boundary layer reattachment and recovery. Inspired by unpublished work of Volker, Langtry, and Menter, the author undertook an independent investigation in an attempt to improve the ability of the Menter shear stress transport (SST) model to predict flowfield characteristics in and downstream of separation bubbles. The fix is an ad hoc term that is a function of the local ratio of turbulent production to dissipation; it is used to multiply the omega-destruction term, increasing eddy viscosity in separated regions. With this fix, several flowfields are investigated. Results show that, although the "separation fix" can provide dramatic improvement in some cases, it is not consistently good for all flows. Thus, although it may prove helpful in many situations in its current form, this model may benefit from further refinements, including better sensitization to the energetics of turbulence in the separated region.

Numerical Simulation of Liquid Jet Atomization Including Turbulence Effects

NASA Technical Reports Server (NTRS)

Trinh, Huu P.; Chen, C. P.; Balasubramanyam, M. S.

2005-01-01

This paper describes numerical implementation of a newly developed hybrid model, T-blob/T-TAB, into an existing computational fluid dynamics (CFD) program for primary and secondary breakup simulation of liquid jet atomization. This model extend two widely used models, the Kelvin-Helmholtz (KH) instability of Reitz (blob model) and the Taylor-Analogy-Breakup (TAB) secondary droplet breakup by O'Rourke and Amsden to include turbulence effects. In the primary breakup model, the level of the turbulence effect on the liquid breakup depends on the characteristic scales and the initial flow conditions. For the secondary breakup, an additional turbulence force acted on parent drops is modeled and integrated into the TAB governing equation. Several assessment studies are presented and the results indicate that the existing KH and TAB models tend to under-predict the product drop size and spray angle, while the current model provides superior results when compared with the measured data.

Integrated modeling of plasma ramp-up in DIII-D ITER-like and high bootstrap current scenario discharges

NASA Astrophysics Data System (ADS)

Wu, M. Q.; Pan, C. K.; Chan, V. S.; Li, G. Q.; Garofalo, A. M.; Jian, X.; Liu, L.; Ren, Q. L.; Chen, J. L.; Gao, X.; Gong, X. Z.; Ding, S. Y.; Qian, J. P.; Cfetr Physics Team

2018-04-01

Time-dependent integrated modeling of DIII-D ITER-like and high bootstrap current plasma ramp-up discharges has been performed with the equilibrium code EFIT, and the transport codes TGYRO and ONETWO. Electron and ion temperature profiles are simulated by TGYRO with the TGLF (SAT0 or VX model) turbulent and NEO neoclassical transport models. The VX model is a new empirical extension of the TGLF turbulent model [Jian et al., Nucl. Fusion 58, 016011 (2018)], which captures the physics of multi-scale interaction between low-k and high-k turbulence from nonlinear gyro-kinetic simulation. This model is demonstrated to accurately model low Ip discharges from the EAST tokamak. Time evolution of the plasma current density profile is simulated by ONETWO with the experimental current ramp-up rate. The general trend of the predicted evolution of the current density profile is consistent with that obtained from the equilibrium reconstruction with Motional Stark effect constraints. The predicted evolution of βN , li , and βP also agrees well with the experiments. For the ITER-like cases, the predicted electron and ion temperature profiles using TGLF_Sat0 agree closely with the experimental measured profiles, and are demonstrably better than other proposed transport models. For the high bootstrap current case, the predicted electron and ion temperature profiles perform better in the VX model. It is found that the SAT0 model works well at high IP (>0.76 MA) while the VX model covers a wider range of plasma current ( IP > 0.6 MA). The results reported in this paper suggest that the developed integrated modeling could be a candidate for ITER and CFETR ramp-up engineering design modeling.

Prediction of the structure of fuel sprays in gas turbine combustors

NASA Technical Reports Server (NTRS)

Shuen, J. S.

1985-01-01

The structure of fuel sprays in a combustion chamber is theoretically investigated using computer models of current interest. Three representative spray models are considered: (1) a locally homogeneous flow (LHF) model, which assumes infinitely fast interphase transport rates; (2) a deterministic separated flow (DSF) model, which considers finite rates of interphase transport but ignores effects of droplet/turbulence interactions; and (3) a stochastic separated flow (SSF) model, which considers droplet/turbulence interactions using random sampling for turbulence properties in conjunction with random-walk computations for droplet motion and transport. Two flow conditions are studied to investigate the influence of swirl on droplet life histories and the effects of droplet/turbulence interactions on flow properties. Comparison of computed results with the experimental data show that general features of the flow structure can be predicted with reasonable accuracy using the two separated flow models. In contrast, the LHF model overpredicts the rate of development of the flow. While the SSF model provides better agreement with measurements than the DSF model, definitive evaluation of the significance of droplet/turbulence interaction is not achieved due to uncertainties in the spray initial conditions.

Characterisation of turbulence downstream of a linear compressor cascade

NASA Astrophysics Data System (ADS)

di Mare, Luca; Jelly, Thomas; Day, Ivor

2014-11-01

Characterisation of turbulence in turbomachinery remains one of the most complex tasks in fluid mechanics. In addition, current closure models required for Reynolds-averaged Navier-Stokes computations do not accurately represent the action of turbulent forces against the mean flow. Therefore, the statistical properties of turbulence in turbomachinery are of significant interest. In the current work, single- and two-point hot-wire measurements have been acquired downstream of a linear compressor cascade in order to examine the properties of large-scale turbulent structures and to assess how they affect turbulent momentum and energy transfer in compressor passages. The cascade has seven controlled diffusion which are representative of high-pressure stator blades found in turbofan engines. Blade chord, thickness and camber are 0.1515 m, 9.3% and 42 degrees, respectively. Measurements were acquired at a chord Reynolds number of 6 . 92 ×105 . Single-point statistics highlight differences in turbulence structure when comparing mid-span and end-wall regions. Evaluation of two-point correlations and their corresponding spectra reveal the length-scales of the energy-bearing eddies in the cascade. Ultimately, these measurements can be used to calibrate future computational models. The authors gratefully acknowledge Rolls-Royce plc for funding this work and granting permission for its publication.

The structure of evaporating and combusting sprays: Measurements and predictions

NASA Technical Reports Server (NTRS)

Shuen, J. S.; Solomon, A. S. P.; Faeth, G. M.

1982-01-01

An apparatus was constructed to provide measurements in open sprays with no zones of recirculation, in order to provide well-defined conditions for use in evaluating spray models. Measurements were completed in a gas jet, in order to test experimental methods, and are currently in progress for nonevaporating sprays. A locally homogeneous flow (LHF) model where interphase transport rates are assumed to be infinitely fast; a separated flow (SF) model which allows for finite interphase transport rates but neglects effects of turbulent fluctuations on drop motion; and a stochastic SF model which considers effects of turbulent fluctuations on drop motion were evaluated using existing data on particle-laden jets. The LHF model generally overestimates rates of particle dispersion while the SF model underestimates dispersion rates. The stochastic SF flow yield satisfactory predictions except at high particle mass loadings where effects of turbulence modulation may have caused the model to overestimate turbulence levels.

Military specifications

NASA Technical Reports Server (NTRS)

Reynolds, Philip

1987-01-01

The current situation relative to the military specification is that there is not one specific model of turbulence which people are using. Particular disagreement exists on how turbulence levels will vary with qualitative analysis. It does not tie one down to specifics. When it comes to flying quality specifications, many feel that one should stay with the definitions of the Cooper-Harper rating scale but allow the levels to shift depending on the level of turbulence. There is a ride quality specification in the MIL-SPEC having to do with flight control systems design that is related to a turbulence model. This spec (MIL-F8785C) and others are discussed.

Exploring reconnection, current sheets, and dissipation in a laboratory MHD turbulence experiment

NASA Astrophysics Data System (ADS)

Schaffner, D. A.

2015-12-01

The Swarthmore Spheromak Experiment (SSX) can serve as a testbed for studying MHD turbulence in a controllable laboratory setting, and in particular, explore the phenomena of reconnection, current sheets and dissipation in MHD turbulence. Plasma with turbulently fluctuating magnetic and velocity fields can be generated using a plasma gun source and launched into a flux-conserving cylindrical tunnel. No background magnetic field is applied so internal fields are allowed to evolve dynamically. Point measurements of magnetic and velocity fluctuations yield broadband power-law spectra with a steepening breakpoint indicative of the onset of a dissipation scale. The frequency range at which this steepening occurs can be correlated to the ion inertial scale of the plasma, a length which is characteristic of the size of current sheets in MHD plasmas and suggests a connection to dissipation. Observation of non-Gaussian intermittent jumps in magnetic field magnitude and angle along with measurements of ion temperature bursts suggests the presence of current sheets embedded within the turbulent plasma, and possibly even active reconnection sites. Additionally, structure function analysis coupled with appeals to fractal scaling models support the hypothesis that current sheets are associated with dissipation in this system.

Turbulent particulate transportation during electrostatic precipitation

NASA Astrophysics Data System (ADS)

Choi, Bum Seog

The generation of secondary flows and turbulence by a corona discharge influences particle transport in an electrostatic precipitator (ESP), and is known to play an important role in the particle collection process. However, it is difficult to characterise theoretically and experimentally the ``turbulent'' fluctuations of the gas flow produced by negative tuft corona. Because of this difficulty, only limited studies have been undertaken previously to understand the structure of corona-induced turbulence and its influence on particle transport in ESPs. The present study is aimed at modelling electrohydrodynamic turbulent flows and particle transport, and at establishing an unproved understanding of them. For a multiply interactive coupling of electrostatics, fluid dynamics and particle dynamics, a strongly coupled system of the governing equations has been solved. The present computer model has considered the most important interaction mechanisms including an ionic wind, corona- induced turbulence and the particle space charge effect. Numerical simulations have been performed for the extensive validation of the numerical and physical models. To account for electrically excited turbulence associated with the inhomogeneous and unsteady characteristics of negative corona discharges, a new turbulence model (called the electrostatic turbulence model) has been developed. In this, an additional production or destruction term is included into the turbulent kinetic energy and dissipation rate equations. It employs a gradient type model of the current density and an electrostatic diffusivity concept. The results of the computation show that the electrostatic turbulence model gives much better agreement with the experimental data than the conventional RNG k-ɛ turbulence model when predicting turbulent gas flows and particle distributions in an ESP. Computations of turbulent particulate two-phase flows for both mono-dispersed and poly-dispersed particles have been performed. The effects of coriona-induced turbulence and the particle space charge on particle transport and the collection process have been investigated. The calculated results for the poly-dispersed particulate flow were compared with those of the mono-dispersed particulate flow, and significant differences were demonstrated. It is established that effective particle- particle interaction occurs, due to the influence of the particle space charge, even for dilute gas-particle flows that occur in ESPs.

Observations of turbulent mixing in a shallow coral reef

NASA Astrophysics Data System (ADS)

Huang, Z. C.

2016-02-01

In situ measurements of waves, currents, and turbulence are presented to study turbulence properties within a depression that is surrounded by multiple coral-reef colonies in a fringing reef in Hobihu, Nan-Wan Bay, southern Taiwan. Turbulence was measured using a dual velocimetry technique, and wave bias contamination in the turbulence is controlled using ogive curve testing of the turbulent shear stress. The observed turbulent dissipation rate is approximately five times greater than simultaneous observations over the nearby sandy bottom site, which indicates stronger mixing within the coral reef than on sandy bottoms. Energetic downward momentum flux exists due to sweeping process; the turbulent kinetic energy is transported downward into the depression through the mechanisms of vertical turbulent transport and advection. The observed turbulent dissipation rate exceeds the shear production rate, which suggests that transport terms or other source terms might be important. The wake flow caused by the resistance force of coral colonies is examined. The form drag coefficient was estimated from the time-averaged alongshore linear momentum between two sites upstream and within the coral reef. The work done due to the form drag, which is termed the wake production, is found to strongly correlate and approximate well to the observed turbulent dissipation rate. The effects of waves and currents on the wake production are discussed. The observed TSS can be described well by classic turbulence closure model when the empirical stability function is adjusted. This study suggests that the complex canopy structure of multiple colonies and the coexistence of the wave-induced and current flows are significant factors for energetic turbulence in the coral reef, which could have positive effects to the health of the coral reefs.

Market Assessment of Forward-Looking Turbulence Sensing Systems

NASA Technical Reports Server (NTRS)

Kauffmann, Paul; Sousa-Poza, Andres

2001-01-01

In recognition of the importance of turbulence mitigation as a tool to improve aviation safety, NASA's Aviation Safety Program developed a Turbulence Detection and Mitigation Sub-element. The objective of this effort is to develop highly reliable turbulence detection technologies for commercial transport aircraft to sense dangerous turbulence with sufficient time warning so that defensive measures can be implemented and prevent passenger and crew injuries. Current research involves three forward sensing products to improve the cockpit awareness of possible turbulence hazards. X-band radar enhancements will improve the capabilities of current weather radar to detect turbulence associated with convective activity. LIDAR (Light Detection and Ranging) is a laser-based technology that is capable of detecting turbulence in clear air. Finally, a possible Radar-LIDAR hybrid sensor is envisioned to detect the full range of convective and clear air turbulence. To support decisions relating to the development of these three forward-looking turbulence sensor technologies, the objective of this study was defined as examination of cost and implementation metrics. Tasks performed included the identification of cost factors and certification issues, the development and application of an implementation model, and the development of cost budget/targets for installing the turbulence sensor and associated software devices into the commercial transport fleet.

Computation of turbulent flows-state-of-the-art, 1970

NASA Technical Reports Server (NTRS)

Reynolds, W. C.

1972-01-01

The state-of-the-art of turbulent flow computation is surveyed. The formulations were generalized to increase the range of their applicability, and the excitement of current debate on equation models was brought into the review. Some new ideas on the modeling of the pressure-strain term in the Reynolds stress equations are also suggested.

Inside pyroclastic density currents - uncovering the enigmatic flow structure and transport behaviour in large-scale experiments

NASA Astrophysics Data System (ADS)

Breard, Eric C. P.; Lube, Gert

2017-01-01

Pyroclastic density currents (PDCs) are the most lethal threat from volcanoes. While there are two main types of PDCs (fully turbulent, fully dilute pyroclastic surges and more concentrated pyroclastic flows encompassing non-turbulent to turbulent transport) pyroclastic flows, which are the subject of the present study, are far more complex than dilute pyroclastic surges and remain the least understood type despite their far greater hazard, greater runout length and ability to transport vast quantities of material across the Earth's surface. Here we present large-scale experiments of natural volcanic material and gas in order to provide the missing quantitative view of the internal structure and gas-particle transport mechanisms in pyroclastic flows. We show that the outer flow structure with head, body and wake regions broadly resembles current PDC analogues of dilute gravity currents. However, the internal structure, in which lower levels consist of a concentrated granular fluid and upper levels are more dilute, contrasts significantly with the internal structure of fully dilute gravity currents. This bipartite vertical structure shows strong analogy to current conceptual models of high-density turbidity currents, which are responsible for the distribution of coarse sediment in marine basins and of great interest to the hydrocarbon industry. The lower concentrated and non-turbulent levels of the PDC (granular-fluid basal flow) act as a fast-flowing carrier for the more dilute and turbulent upper levels of the current (ash-cloud surge). Strong kinematic coupling between these flow parts reduces viscous dissipation and entrainment of ambient air into the lower part of the ash-cloud surge. This leads to a state of forced super-criticality whereby fast and destructive PDCs can endure even at large distances from volcanoes. Importantly, the basal flow/ash-cloud surge coupling yields a characteristically smooth rheological boundary across the non-turbulent/turbulent interface, as well as vertical velocity and density profiles in the ash-cloud surge, which strongly differ from current theoretical predictions. Observed generation of successive pulses of high dynamic pressure within the upper dilute levels of the PDC may be important to understand the destructive potential of PDCs. The experiments further show that a wide range in the degree of coupling between particle and gas phases is critical to the vertical and longitudinal segregation of the currents into reaches that have starkly contrasting sediment transport capacities. In particular, the formation of mesoscale turbulence clusters under strong particle-gas feedback controls vertical stratification inside the turbulent upper levels of the current (ash-cloud surge) and triggers significant transfers of mass and momentum from the ash-cloud surge onto the granular-fluid basal flow. These results open up new pathways to advance current computational PDC hazard models and to describe and interpret PDCs as well as other types of high-density gravity currents transported across the surfaces of Earth and other planets and across marine basins.

The distribution of density in supersonic turbulence

NASA Astrophysics Data System (ADS)

Squire, Jonathan; Hopkins, Philip F.

2017-11-01

We propose a model for the statistics of the mass density in supersonic turbulence, which plays a crucial role in star formation and the physics of the interstellar medium (ISM). The model is derived by considering the density to be arranged as a collection of strong shocks of width ˜ M^{-2}, where M is the turbulent Mach number. With two physically motivated parameters, the model predicts all density statistics for M>1 turbulence: the density probability distribution and its intermittency (deviation from lognormality), the density variance-Mach number relation, power spectra and structure functions. For the proposed model parameters, reasonable agreement is seen between model predictions and numerical simulations, albeit within the large uncertainties associated with current simulation results. More generally, the model could provide a useful framework for more detailed analysis of future simulations and observational data. Due to the simple physical motivations for the model in terms of shocks, it is straightforward to generalize to more complex physical processes, which will be helpful in future more detailed applications to the ISM. We see good qualitative agreement between such extensions and recent simulations of non-isothermal turbulence.

Modeling the transition region

NASA Technical Reports Server (NTRS)

Singer, Bart A.

1993-01-01

The current status of transition-region models is reviewed in this report. To understand modeling problems, various flow features that influence the transition process are discussed first. Then an overview of the different approaches to transition-region modeling is given. This is followed by a detailed discussion of turbulence models and the specific modifications that are needed to predict flows undergoing laminar-turbulent transition. Methods for determining the usefulness of the models are presented, and an outlook for the future of transition-region modeling is suggested.

Solving Navier-Stokes Equations with Advanced Turbulence Models on Three-Dimensional Unstructured Grids

NASA Technical Reports Server (NTRS)

Wang, Qun-Zhen; Massey, Steven J.; Abdol-Hamid, Khaled S.; Frink, Neal T.

1999-01-01

USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flows. We have implemented two modified versions of the original Jones and Launder k-epsilon two-equation turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for two flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those of empirical formulae, theoretical results and the existing Spalart-Allmaras one-equation model.

Validation and optimization of SST k-ω turbulence model for pollutant dispersion within a building array

NASA Astrophysics Data System (ADS)

Yu, Hesheng; Thé, Jesse

2016-11-01

The prediction of the dispersion of air pollutants in urban areas is of great importance to public health, homeland security, and environmental protection. Computational Fluid Dynamics (CFD) emerges as an effective tool for pollutant dispersion modelling. This paper reports and quantitatively validates the shear stress transport (SST) k-ω turbulence closure model and its transitional variant for pollutant dispersion under complex urban environment for the first time. Sensitivity analysis is performed to establish recommendation for the proper use of turbulence models in urban settings. The current SST k-ω simulation is validated rigorously by extensive experimental data using hit rate for velocity components, and the "factor of two" of observations (FAC2) and fractional bias (FB) for concentration field. The simulation results show that current SST k-ω model can predict flow field nicely with an overall hit rate of 0.870, and concentration dispersion with FAC2 = 0.721 and FB = 0.045. The flow simulation of the current SST k-ω model is slightly inferior to that of a detached eddy simulation (DES), but better than that of standard k-ε model. However, the current study is the best among these three model approaches, when validated against measurements of pollutant dispersion in the atmosphere. This work aims to provide recommendation for proper use of CFD to predict pollutant dispersion in urban environment.

A Variable Turbulent Schmidt Number Formulation for Scramjet Application

NASA Technical Reports Server (NTRS)

Xiao, X.; Edwards, J. R.; Hassan, H. A.; Cutler, A. D.

2004-01-01

In high speed engines, thorough turbulent mixing of fuel and air is required to obtain high performance and high efficiency. Thus, the ability to predict turbulent mixing is crucial in obtaining accurate numerical simulation of an engine and its performance. Current state of the art in CFD simulation is to assume both turbulent Prandtl number and Schmidt numbers to be constants. However, since the mixing of fuel and air is inversely proportional to the Schmidt number, a value of 0.45 for the Schmidt number will produce twice as much diffusion as that with a value of 0.9. Because of this, current CFD tools and models have not been able to provide the needed guidance required for the efficient design of a scramjet engine. The goal of this investigation is to develop the framework needed to calculate turbulent Prandtl and Schmidt numbers as part of the solution. This requires four additional equations: two for the temperature variance and its dissipation rate and two for the concentration variance and its dissipation rate. In the current investigation emphasis will be placed on studying mixing without reactions. For such flows, variable Prandtl number does not play a major role in determining the flow. This, however, will have to be addressed when combustion is present. The approach to be used is similar to that used to develop the k-zeta model. In this approach, relevant equations are derived from the exact Navier-Stokes equations and each individual correlation is modeled. This ensures that relevant physics is incorporated into the model equations. This task has been accomplished. The final set of equations have no wall or damping functions. Moreover, they are tensorially consistent and Galilean invariant. The derivation of the model equations is rather lengthy and thus will not be incorporated into this abstract, but will be included in the final paper. As a preliminary to formulating the proposed model, the original k-zeta model with constant turbulent Prandtl and Schmidt numbers is used to model the supersonic coaxial jet mixing experiments involving He, O2 and air.

Numerical simulation of turbulent combustion: Scientific challenges

NASA Astrophysics Data System (ADS)

Ren, ZhuYin; Lu, Zhen; Hou, LingYun; Lu, LiuYan

2014-08-01

Predictive simulation of engine combustion is key to understanding the underlying complicated physicochemical processes, improving engine performance, and reducing pollutant emissions. Critical issues as turbulence modeling, turbulence-chemistry interaction, and accommodation of detailed chemical kinetics in complex flows remain challenging and essential for high-fidelity combustion simulation. This paper reviews the current status of the state-of-the-art large eddy simulation (LES)/prob-ability density function (PDF)/detailed chemistry approach that can address the three challenging modelling issues. PDF as a subgrid model for LES is formulated and the hybrid mesh-particle method for LES/PDF simulations is described. Then the development need in micro-mixing models for the PDF simulations of turbulent premixed combustion is identified. Finally the different acceleration methods for detailed chemistry are reviewed and a combined strategy is proposed for further development.

A domain analysis approach to clear-air turbulence forecasting using high-density in-situ measurements

NASA Astrophysics Data System (ADS)

Abernethy, Jennifer A.

Pilots' ability to avoid clear-air turbulence (CAT) during flight affects the safety of the millions of people who fly commercial airlines and other aircraft, and turbulence costs millions in injuries and aircraft maintenance every year. Forecasting CAT is not straightforward, however; microscale features like the turbulence eddies that affect aircraft (100m) are below the current resolution of operational numerical weather prediction (NWP) models, and the only evidence of CAT episodes, until recently, has been sparse, subjective reports from pilots known as PIREPs. To forecast CAT, researchers use a simple weighted sum of top-performing turbulence indicators derived from NWP model outputs---termed diagnostics---based on their agreement with current PIREPs. However, a new, quantitative source of observation data---high-density measurements made by sensor equipment and software on aircraft, called in-situ measurements---is now available. The main goal of this thesis is to develop new data analysis and processing techniques to apply to the model and new observation data, in order to improve CAT forecasting accuracy. This thesis shows that using in-situ data improves forecasting accuracy and that automated machine learning algorithms such as support vector machines (SVM), logistic regression, and random forests, can match current performance while eliminating almost all hand-tuning. Feature subset selection is paired with the new algorithms to choose diagnostics that predict well as a group rather than individually. Specializing forecasts and choice of diagnostics by geographic region further improves accuracy because of the geographic variation in turbulence sources. This work uses random forests to find climatologically-relevant regions based on these variations and implements a forecasting system testbed which brings these techniques together to rapidly prototype new, regionalized versions of operational CAT forecasting systems.

Aggregates: The Fundamental Building Blocks of Planetesimals?

NASA Technical Reports Server (NTRS)

Cuzzi, J. N.; Hartlep, T.; Simon, J. I.; Cato, M. J.

2017-01-01

The initial accretion of primitive asteroids (meteorite parent bodies) from freely-floating nebula particles remains problematic. Traditional growth-by-sticking models in turbulent nebulae encounter a formidable "meter-size barrier" due to both drift and destruction, or even a mmtocmsize "bouncing" barrier. Even if growth by sticking could somehow breach these barriers (perhaps if the actual sticking or strength is larger than current estimates, which are based on pure ice or pure silicate), turbulent nebulae present further obstacles through the 1-10km size range. On the other hand, nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids. Thus, the intensity of nebula turbulence (or "alpha ") is critical to the entire process. Theoretical understanding of nebula turbulence continues to evolve; while recent models of MRI (magneticallydriven) turbulence favor lowornoturbulence environments, purely hydrodynamic turbulence is making a comeback with three recently discovered mechanisms generating turbulence of moderate which do not rely on magnetic fields at all.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Xiao, Heng; Gustafson, William I.; Wang, Hailong

Subgrid-scale interactions between turbulence and radiation are potentially important for accurately reproducing marine low clouds in climate models. To better understand the impact of these interactions, the Weather Research and Forecasting (WRF) model is configured for large eddy simulation (LES) to study the stratocumulus-to-trade cumulus (Sc-to-Cu) transition. Using the GEWEX Atmospheric System Studies (GASS) composite Lagrangian transition case and the Atlantic Trade Wind Experiment (ATEX) case, it is shown that the lack of subgrid-scale turbulence-radiation interaction, as is the case in current generation climate models, accelerates the Sc-to-Cu transition. Our analysis suggests that in cloud-topped boundary layers subgrid-scale turbulence-radiation interactionsmore » contribute to stronger production of temperature variance, which in turn leads to stronger buoyancy production of turbulent kinetic energy and helps to maintain the Sc cover.« less

Second order closure modeling of turbulent buoyant wall plumes

NASA Technical Reports Server (NTRS)

Zhu, Gang; Lai, Ming-Chia; Shih, Tsan-Hsing

1992-01-01

Non-intrusive measurements of scalar and momentum transport in turbulent wall plumes, using a combined technique of laser Doppler anemometry and laser-induced fluorescence, has shown some interesting features not present in the free jet or plumes. First, buoyancy-generation of turbulence is shown to be important throughout the flow field. Combined with low-Reynolds-number turbulence and near-wall effect, this may raise the anisotropic turbulence structure beyond the prediction of eddy-viscosity models. Second, the transverse scalar fluxes do not correspond only to the mean scalar gradients, as would be expected from gradient-diffusion modeling. Third, higher-order velocity-scalar correlations which describe turbulent transport phenomena could not be predicted using simple turbulence models. A second-order closure simulation of turbulent adiabatic wall plumes, taking into account the recent progress in scalar transport, near-wall effect and buoyancy, is reported in the current study to compare with the non-intrusive measurements. In spite of the small velocity scale of the wall plumes, the results showed that low-Reynolds-number correction is not critically important to predict the adiabatic cases tested and cannot be applied beyond the maximum velocity location. The mean and turbulent velocity profiles are very closely predicted by the second-order closure models. but the scalar field is less satisfactory, with the scalar fluctuation level underpredicted. Strong intermittency of the low-Reynolds-number flow field is suspected of these discrepancies. The trends in second- and third-order velocity-scalar correlations, which describe turbulent transport phenomena, are also predicted in general, with the cross-streamwise correlations better than the streamwise one. Buoyancy terms modeling the pressure-correlation are shown to improve the prediction slightly. The effects of equilibrium time-scale ratio and boundary condition are also discussed.

Effect of Turbulence Models on Two Massively-Separated Benchmark Flow Cases

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.

2003-01-01

Two massively-separated flow cases (the 2-D hill and the 3-D Ahmed body) were computed with several different turbulence models in the Reynolds-averaged Navier-Stokes code CFL3D as part of participation in a turbulence modeling workshop held in Poitiers, France in October, 2002. Overall, results were disappointing, but were consistent with results from other RANS codes and other turbulence models at the workshop. For the 2-D hill case, those turbulence models that predicted separation location accurately ended up yielding a too-long separation extent downstream. The one model that predicted a shorter separation extent in better agreement with LES data did so only by coincidence: its prediction of earlier reattachment was due to a too-late prediction of the separation location. For the Ahmed body, two slant angles were computed, and CFD performed fairly well for one of the cases (the larger slant angle). Both turbulence models tested in this case were very similar to each other. For the smaller slant angle, CFD predicted massive separation, whereas the experiment showed reattachment about half-way down the center of the face. These test cases serve as reminders that state- of-the-art CFD is currently not a reliable predictor of massively-separated flow physics, and that further validation studies in this area would be beneficial.

A new approach for turbulent simulations in complex geometries

NASA Astrophysics Data System (ADS)

Israel, Daniel M.

Historically turbulence modeling has been sharply divided into Reynolds averaged Navier-Stokes (RANS), in which all the turbulent scales of motion are modeled, and large-eddy simulation (LES), in which only a portion of the turbulent spectrum is modeled. In recent years there have been numerous attempts to couple these two approaches either by patching RANS and LES calculations together (zonal methods) or by blending the two sets of equations. In order to create a proper bridging model, that is, a single set of equations which captures both RANS and LES like behavior, it is necessary to place both RANS and LES in a more general framework. The goal of the current work is threefold: to provide such a framework, to demonstrate how the Flow Simulation Methodology (FSM) fits into this framework, and to evaluate the strengths and weaknesses of the current version of the FSM. To do this, first a set of filtered Navier-Stokes (FNS) equations are introduced in terms of an arbitrary generalized filter. Additional exact equations are given for the second order moments and the generalized subfilter dissipation rate tensor. This is followed by a discussion of the role of implicit and explicit filters in turbulence modeling. The FSM is then described with particular attention to its role as a bridging model. In order to evaluate the method a specific implementation of the FSM approach is proposed. Simulations are presented using this model for the case of a separating flow over a "hump" with and without flow control. Careful attention is paid to error estimation, and, in particular, how using flow statistics and time series affects the error analysis. Both mean flow and Reynolds stress profiles are presented, as well as the phase averaged turbulent structures and wall pressure spectra. Using the phase averaged data it is possible to examine how the FSM partitions the energy between the coherent resolved scale motions, the random resolved scale fluctuations, and the subfilter quantities. The method proves to be qualitatively successful at reproducing large turbulent structures. However, like other hybrid methods, it has difficulty in the region where the model behavior transitions from RANS to LES. Consequently the phase averaged structures reproduce the experiments quite well, and the forcing does significantly reduce the length of the separated region. Nevertheless, the recirculation length is significantly too large for all the cases. Overall the current results demonstrate the promise of bridging models in general and the FSM in particular. However, current bridging techniques are still in their infancy. There is still important progress to be made and it is hoped that this work points out the more important avenues for exploration.

Nearshore Current Model Workshop Summary.

DTIC Science & Technology

1983-09-01

dissipation , and wave-current interaction. b. Incorporation into models of wave-breaking. c. Parameterization of turbulence in models. d. Incorporation...into models of surf zone energy dissipation . e. Methods to specify waves and currents on the boundaries of the grid. f. Incorporation into models of...also recommended. Improvements should include nonlinear and irregular wave effects and improved models of wave-breaking and wave energy dissipation in

A 3D unstructured grid nearshore hydrodynamic model based on the vortex force formalism

NASA Astrophysics Data System (ADS)

Zheng, Peng; Li, Ming; van der A, Dominic A.; van der Zanden, Joep; Wolf, Judith; Chen, Xueen; Wang, Caixia

2017-08-01

A new three-dimensional nearshore hydrodynamic model system is developed based on the unstructured-grid version of the third generation spectral wave model SWAN (Un-SWAN) coupled with the three-dimensional ocean circulation model FVCOM to enable the full representation of the wave-current interaction in the nearshore region. A new wave-current coupling scheme is developed by adopting the vortex-force (VF) scheme to represent the wave-current interaction. The GLS turbulence model is also modified to better reproduce wave-breaking enhanced turbulence, together with a roller transport model to account for the effect of surface wave roller. This new model system is validated first against a theoretical case of obliquely incident waves on a planar beach, and then applied to three test cases: a laboratory scale experiment of normal waves on a beach with a fixed breaker bar, a field experiment of oblique incident waves on a natural, sandy barred beach (Duck'94 experiment), and a laboratory study of normal-incident waves propagating around a shore-parallel breakwater. Overall, the model predictions agree well with the available measurements in these tests, illustrating the robustness and efficiency of the present model for very different spatial scales and hydrodynamic conditions. Sensitivity tests indicate the importance of roller effects and wave energy dissipation on the mean flow (undertow) profile over the depth. These tests further suggest to adopt a spatially varying value for roller effects across the beach. In addition, the parameter values in the GLS turbulence model should be spatially inhomogeneous, which leads to better prediction of the turbulent kinetic energy and an improved prediction of the undertow velocity profile.

Learning to soar in turbulent environments

NASA Astrophysics Data System (ADS)

Reddy, Gautam; Celani, Antonio; Sejnowski, Terrence; Vergassola, Massimo

Birds and gliders exploit warm, rising atmospheric currents (thermals) to reach heights comparable to low-lying clouds with a reduced expenditure of energy. Soaring provides a remarkable instance of complex decision-making in biology and requires a long-term strategy to effectively use the ascending thermals. Furthermore, the problem is technologically relevant to extend the flying range of autonomous gliders. The formation of thermals unavoidably generates strong turbulent fluctuations, which make deriving an efficient policy harder and thus constitute an essential element of soaring. Here, we approach soaring flight as a problem of learning to navigate highly fluctuating turbulent environments. We simulate the atmospheric boundary layer by numerical models of turbulent convective flow and combine them with model-free, experience-based, reinforcement learning algorithms to train the virtual gliders. For the learned policies in the regimes of moderate and strong turbulence levels, the virtual glider adopts an increasingly conservative policy as turbulence levels increase, quantifying the degree of risk affordable in turbulent environments. Reinforcement learning uncovers those sensorimotor cues that permit effective control over soaring in turbulent environments.

Current Status on the use of Parallel Computing in Turbulent Reacting Flow Computations Involving Sprays, Monte Carlo PDF and Unstructured Grids. Chapter 4

NASA Technical Reports Server (NTRS)

Raju, M. S.

1998-01-01

The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Liu, nan-Suey

2010-01-01

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

Implementation of Advanced Two Equation Turbulence Models in the USM3D Unstructured Flow Solver

NASA Technical Reports Server (NTRS)

Wang, Qun-Zhen; Massey, Steven J.; Abdol-Hamid, Khaled S.

2000-01-01

USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flow. We have implemented two modified versions of the original Jones and Launder k-epsilon "two-equation" turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for three flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those from direct numerical simulation, empirical formulae, theoretical results, and the existing Spalart-Allmaras one-equation model.

Exploiting similarity in turbulent shear flows for turbulence modeling

NASA Technical Reports Server (NTRS)

Robinson, David F.; Harris, Julius E.; Hassan, H. A.

1992-01-01

It is well known that current k-epsilon models cannot predict the flow over a flat plate and its wake. In an effort to address this issue and other issues associated with turbulence closure, a new approach for turbulence modeling is proposed which exploits similarities in the flow field. Thus, if we consider the flow over a flat plate and its wake, then in addition to taking advantage of the log-law region, we can exploit the fact that the flow becomes self-similar in the far wake. This latter behavior makes it possible to cast the governing equations as a set of total differential equations. Solutions of this set and comparison with measured shear stress and velocity profiles yields the desired set of model constants. Such a set is, in general, different from other sets of model constants. The rational for such an approach is that if we can correctly model the flow over a flat plate and its far wake, then we can have a better chance of predicting the behavior in between. It is to be noted that the approach does not appeal, in any way, to the decay of homogeneous turbulence. This is because the asymptotic behavior of the flow under consideration is not representative of the decay of homogeneous turbulence.

Exploiting similarity in turbulent shear flows for turbulence modeling

NASA Astrophysics Data System (ADS)

Robinson, David F.; Harris, Julius E.; Hassan, H. A.

1992-12-01

It is well known that current k-epsilon models cannot predict the flow over a flat plate and its wake. In an effort to address this issue and other issues associated with turbulence closure, a new approach for turbulence modeling is proposed which exploits similarities in the flow field. Thus, if we consider the flow over a flat plate and its wake, then in addition to taking advantage of the log-law region, we can exploit the fact that the flow becomes self-similar in the far wake. This latter behavior makes it possible to cast the governing equations as a set of total differential equations. Solutions of this set and comparison with measured shear stress and velocity profiles yields the desired set of model constants. Such a set is, in general, different from other sets of model constants. The rational for such an approach is that if we can correctly model the flow over a flat plate and its far wake, then we can have a better chance of predicting the behavior in between. It is to be noted that the approach does not appeal, in any way, to the decay of homogeneous turbulence. This is because the asymptotic behavior of the flow under consideration is not representative of the decay of homogeneous turbulence.

Contributions of numerical simulation data bases to the physics, modeling and measurement of turbulence

NASA Technical Reports Server (NTRS)

Moin, Parviz; Spalart, Philippe R.

1987-01-01

The use of simulation data bases for the examination of turbulent flows is an effective research tool. Studies of the structure of turbulence have been hampered by the limited number of probes and the impossibility of measuring all desired quantities. Also, flow visualization is confined to the observation of passive markers with limited field of view and contamination caused by time-history effects. Computer flow fields are a new resource for turbulence research, providing all the instantaneous flow variables in three-dimensional space. Simulation data bases also provide much-needed information for phenomenological turbulence modeling. Three dimensional velocity and pressure fields from direct simulations can be used to compute all the terms in the transport equations for the Reynolds stresses and the dissipation rate. However, only a few, geometrically simple flows have been computed by direct numerical simulation, and the inventory of simulation does not fully address the current modeling needs in complex turbulent flows. The availability of three-dimensional flow fields also poses challenges in developing new techniques for their analysis, techniques based on experimental methods, some of which are used here for the analysis of direct-simulation data bases in studies of the mechanics of turbulent flows.

Model simulations of dense bottom currents in the Western Baltic Sea

NASA Astrophysics Data System (ADS)

Burchard, Hans; Janssen, Frank; Bolding, Karsten; Umlauf, Lars; Rennau, Hannes

2009-01-01

Only recently, medium intensity inflow events into the Baltic Sea have gained more awareness because of their potential to ventilate intermediate layers in the Southern Baltic Sea basins. With the present high-resolution model study of the Western Baltic Sea a first attempt is made to obtain model based realistic estimates of turbulent mixing in this area where dense bottom currents resulting from medium intensity inflow events are weakened by turbulent entrainment. The numerical model simulation which is carried out using the General Estuarine Transport Model (GETM) during nine months in 2003 and 2004 is first validated by means of three automatic stations at the Drogden and Darss Sills and in the Arkona Sea. In order to obtain good agreement between observations and model results, the 0.5×0.5 nautical mile bathymetry had to be adjusted in order to account for the fact that even at that scale many relevant topographic features are not resolved. Current velocity, salinity and turbulence observations during a medium intensity inflow event through the Øresund are then compared to the model results. Given the general problems of point to point comparisons between observations and model simulations, the agreement is fairly good with the characteristic features of the inflow event well represented by the model simulations. Two different bulk measures for mixing activity are then introduced, the vertically integrated decay of salinity variance, which is equal to the production of micro-scale salinity variance, and the vertically integrated turbulent salt flux, which is related to an increase of potential energy due to vertical mixing of stably stratified flow. Both measures give qualitatively similar results and identify the Drogden and Darss Sills as well as the Bornholm Channel as mixing hot spots. Further regions of strong mixing are the dense bottom current pathways from these sills into the Arkona Sea, areas around Kriegers Flak (a shoal in the western Arkona Sea) and north-west of the island of Rügen.

An error reduction algorithm to improve lidar turbulence estimates for wind energy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Newman, Jennifer F.; Clifton, Andrew

Remote-sensing devices such as lidars are currently being investigated as alternatives to cup anemometers on meteorological towers for the measurement of wind speed and direction. Although lidars can measure mean wind speeds at heights spanning an entire turbine rotor disk and can be easily moved from one location to another, they measure different values of turbulence than an instrument on a tower. Current methods for improving lidar turbulence estimates include the use of analytical turbulence models and expensive scanning lidars. While these methods provide accurate results in a research setting, they cannot be easily applied to smaller, vertically profiling lidarsmore » in locations where high-resolution sonic anemometer data are not available. Thus, there is clearly a need for a turbulence error reduction model that is simpler and more easily applicable to lidars that are used in the wind energy industry. In this work, a new turbulence error reduction algorithm for lidars is described. The Lidar Turbulence Error Reduction Algorithm, L-TERRA, can be applied using only data from a stand-alone vertically profiling lidar and requires minimal training with meteorological tower data. The basis of L-TERRA is a series of physics-based corrections that are applied to the lidar data to mitigate errors from instrument noise, volume averaging, and variance contamination. These corrections are applied in conjunction with a trained machine-learning model to improve turbulence estimates from a vertically profiling WINDCUBE v2 lidar. The lessons learned from creating the L-TERRA model for a WINDCUBE v2 lidar can also be applied to other lidar devices. L-TERRA was tested on data from two sites in the Southern Plains region of the United States. The physics-based corrections in L-TERRA brought regression line slopes much closer to 1 at both sites and significantly reduced the sensitivity of lidar turbulence errors to atmospheric stability. The accuracy of machine-learning methods in L-TERRA was highly dependent on the input variables and training dataset used, suggesting that machine learning may not be the best technique for reducing lidar turbulence intensity (TI) error. Future work will include the use of a lidar simulator to better understand how different factors affect lidar turbulence error and to determine how these errors can be reduced using information from a stand-alone lidar.« less

An error reduction algorithm to improve lidar turbulence estimates for wind energy

DOE PAGES

Newman, Jennifer F.; Clifton, Andrew

2017-02-10

Remote-sensing devices such as lidars are currently being investigated as alternatives to cup anemometers on meteorological towers for the measurement of wind speed and direction. Although lidars can measure mean wind speeds at heights spanning an entire turbine rotor disk and can be easily moved from one location to another, they measure different values of turbulence than an instrument on a tower. Current methods for improving lidar turbulence estimates include the use of analytical turbulence models and expensive scanning lidars. While these methods provide accurate results in a research setting, they cannot be easily applied to smaller, vertically profiling lidarsmore » in locations where high-resolution sonic anemometer data are not available. Thus, there is clearly a need for a turbulence error reduction model that is simpler and more easily applicable to lidars that are used in the wind energy industry. In this work, a new turbulence error reduction algorithm for lidars is described. The Lidar Turbulence Error Reduction Algorithm, L-TERRA, can be applied using only data from a stand-alone vertically profiling lidar and requires minimal training with meteorological tower data. The basis of L-TERRA is a series of physics-based corrections that are applied to the lidar data to mitigate errors from instrument noise, volume averaging, and variance contamination. These corrections are applied in conjunction with a trained machine-learning model to improve turbulence estimates from a vertically profiling WINDCUBE v2 lidar. The lessons learned from creating the L-TERRA model for a WINDCUBE v2 lidar can also be applied to other lidar devices. L-TERRA was tested on data from two sites in the Southern Plains region of the United States. The physics-based corrections in L-TERRA brought regression line slopes much closer to 1 at both sites and significantly reduced the sensitivity of lidar turbulence errors to atmospheric stability. The accuracy of machine-learning methods in L-TERRA was highly dependent on the input variables and training dataset used, suggesting that machine learning may not be the best technique for reducing lidar turbulence intensity (TI) error. Future work will include the use of a lidar simulator to better understand how different factors affect lidar turbulence error and to determine how these errors can be reduced using information from a stand-alone lidar.« less

Effects of turbulence on a kinetic auroral arc model

NASA Technical Reports Server (NTRS)

Cornwall, J. M.; Chiu, Y. T.

1981-01-01

A plasma kinetic model of an inverted-V auroral arc structure which includes the effects of electrostatic turbulence is proposed. In the absence of turbulence, a parallel potential drop is supported by magnetic mirror forces and charge quasi neutrality, with energetic auroral ions penetrating to low altitudes; relative to the electrons, the ions' pitch angle distribution is skewed toward smaller pitch angles. The electrons energized by the potential drop form a current which excites electrostatic turbulence. In equilibrium the plasma is marginally stable. The conventional anomalous resistivity contribution to the potential drop is very small. Anomalous resistivity processes are far too dissipative to be powered by auroral particles. It is concluded that under certain circumstances equilibrium may be impossible and relaxation oscillations set in.

Current Trends in Modeling Research for Turbulent Aerodynamic Flows

NASA Technical Reports Server (NTRS)

Gatski, Thomas B.; Rumsey, Christopher L.; Manceau, Remi

2007-01-01

The engineering tools of choice for the computation of practical engineering flows have begun to migrate from those based on the traditional Reynolds-averaged Navier-Stokes approach to methodologies capable, in theory if not in practice, of accurately predicting some instantaneous scales of motion in the flow. The migration has largely been driven by both the success of Reynolds-averaged methods over a wide variety of flows as well as the inherent limitations of the method itself. Practitioners, emboldened by their ability to predict a wide-variety of statistically steady, equilibrium turbulent flows, have now turned their attention to flow control and non-equilibrium flows, that is, separation control. This review gives some current priorities in traditional Reynolds-averaged modeling research as well as some methodologies being applied to a new class of turbulent flow control problems.

Counterintuitive effects of substrate roughness on PDCs

NASA Astrophysics Data System (ADS)

Andrews, B. J.; Manga, M.

2012-12-01

We model dilute pyroclastic density currents (PDCs) using scaled, warm, particle-laden density currents in a 6 m long, 0.6 m wide, 1.8 m tall air-filled tank. In this set of experiments, we run currents over substrates with characteristic roughness scales, hr, ranging over ~3 orders of magnitude from smooth, through 250 μm sandpaper, 0.1-, 1-, 2-, 5-, and 10 cm hemispheres. As substrate roughness increases, runout distance increases until a critical roughness height, hrc, is reached; further increases in roughness height decrease runout. The critical roughness height appears to be 0.25-0.5 htb, the thickness of the turbulent lower layer of the density currents. The dependence of runout on hr is most likely the result of increases in substrate roughness decreasing the average current velocity and converting that energy into increased turbulence intensity. Small values of hr thus result in increased runout as sedimentation is inhibited by the increased turbulence intensity. At larger values of hr current behavior is controlled by much larger decreases in average current velocity, even though sedimentation decreases. Scaling our experiments up to the size of real volcanic eruptions suggests that landscapes must have characteristic roughness hr>10 m to reduce the runout of natural PDCs, smaller roughness scales can increase runout. Comparison of relevant bulk (Reynolds number, densimetric and thermal Richardson numbers, excess buoyant thermal energy density) and turbulent (Stokes and settling numbers) between our experiments and natural dilute PDCs indicates that we are accurately modeling at least the large scale behaviors and dynamics of dilute PDCs.

Mixing parametrizations for ocean climate modelling

NASA Astrophysics Data System (ADS)

Gusev, Anatoly; Moshonkin, Sergey; Diansky, Nikolay; Zalesny, Vladimir

2016-04-01

The algorithm is presented of splitting the total evolutionary equations for the turbulence kinetic energy (TKE) and turbulence dissipation frequency (TDF), which is used to parameterize the viscosity and diffusion coefficients in ocean circulation models. The turbulence model equations are split into the stages of transport-diffusion and generation-dissipation. For the generation-dissipation stage, the following schemes are implemented: the explicit-implicit numerical scheme, analytical solution and the asymptotic behavior of the analytical solutions. The experiments were performed with different mixing parameterizations for the modelling of Arctic and the Atlantic climate decadal variability with the eddy-permitting circulation model INMOM (Institute of Numerical Mathematics Ocean Model) using vertical grid refinement in the zone of fully developed turbulence. The proposed model with the split equations for turbulence characteristics is similar to the contemporary differential turbulence models, concerning the physical formulations. At the same time, its algorithm has high enough computational efficiency. Parameterizations with using the split turbulence model make it possible to obtain more adequate structure of temperature and salinity at decadal timescales, compared to the simpler Pacanowski-Philander (PP) turbulence parameterization. Parameterizations with using analytical solution or numerical scheme at the generation-dissipation step of the turbulence model leads to better representation of ocean climate than the faster parameterization using the asymptotic behavior of the analytical solution. At the same time, the computational efficiency left almost unchanged relative to the simple PP parameterization. Usage of PP parametrization in the circulation model leads to realistic simulation of density and circulation with violation of T,S-relationships. This error is majorly avoided with using the proposed parameterizations containing the split turbulence model. The high sensitivity of the eddy-permitting circulation model to the definition of mixing is revealed, which is associated with significant changes of density fields in the upper baroclinic ocean layer over the total considered area. For instance, usage of the turbulence parameterization instead of PP algorithm leads to increasing circulation velocity in the Gulf Stream and North Atlantic Current, as well as the subpolar cyclonic gyre in the North Atlantic and Beaufort Gyre in the Arctic basin are reproduced more realistically. Consideration of the Prandtl number as a function of the Richardson number significantly increases the modelling quality. The research was supported by the Russian Foundation for Basic Research (grant № 16-05-00534) and the Council on the Russian Federation President Grants (grant № MK-3241.2015.5)

Turbulence simulation mechanization for Space Shuttle Orbiter dynamics and control studies

NASA Technical Reports Server (NTRS)

Tatom, F. B.; King, R. L.

1977-01-01

The current version of the NASA turbulent simulation model in the form of a digital computer program, TBMOD, is described. The logic of the program is discussed and all inputs and outputs are defined. An alternate method of shear simulation suitable for incorporation into the model is presented. The simulation is based on a von Karman spectrum and the assumption of isotropy. The resulting spectral density functions for the shear model are included.

Numerical modeling of the transitional boundary layer over a flat plate

NASA Astrophysics Data System (ADS)

Ivanov, Dimitry; Chorny, Andrei

2015-11-01

Our example is connected with fundamental research on understanding how an initially laminar boundary layer becomes turbulent. We have chosen the flow over a flat plate as a prototype for boundary-layer flows around bodies. Special attention was paid to the near-wall region in order to capture all levels of the boundary layer. In this study, the numerical software package OpenFOAM has been used in order to solve the flow field. The results were used in a comparative study with data obtained from Large Eddy Simulation (LES). The composite SGS-wall model is presently incorporated into a computer code suitable for the LES of developing flat-plate boundary layers. Presently this model is extended to the LES of the zero-pressure gradient, flat-plate turbulent boundary layer. In current study the time discretization is based on a second order Crank-Nicolson/Adams-Bashforth method. LES solver using Smagorinsky and the one-equation LES turbulence models. The transition models significantly improve the prediction of the onset location compared to the fully turbulent models.LES methods appear to be the most promising new tool for the design and analysis of flow devices including transition regions of the turbulent flow.

A Simplified Ab Initio Cosmic-ray Modulation Model with Simulated Time Dependence and Predictive Capability

NASA Astrophysics Data System (ADS)

Moloto, K. D.; Engelbrecht, N. E.; Burger, R. A.

2018-06-01

A simplified ab initio approach is followed to model cosmic-ray proton modulation, using a steady-state three-dimensional stochastic solver of the Parker transport equation that simulates some effects of time dependence. Standard diffusion coefficients based on Quasilinear Theory and Nonlinear Guiding Center Theory are employed. The spatial and temporal dependences of the various turbulence quantities required as inputs for the diffusion, as well as the turbulence-reduced drift coefficients, follow from parametric fits to results from a turbulence transport model as well as from spacecraft observations of these turbulence quantities. Effective values are used for the solar wind speed, magnetic field magnitude, and tilt angle in the modulation model to simulate temporal effects due to changes in the large-scale heliospheric plasma. The unusually high cosmic-ray intensities observed during the 2009 solar minimum follow naturally from the current model for most of the energies considered. This demonstrates that changes in turbulence contribute significantly to the high intensities during that solar minimum. We also discuss and illustrate how this model can be used to predict future cosmic-ray intensities, and comment on the reliability of such predictions.

Turbulence vs Self-organized criticality: A hybrid approach, with implications for substorm dynamics of the magnetosphere

NASA Astrophysics Data System (ADS)

Milovanov, Alexander

Plasmas in astrophysics, cosmical geophysics, and laboratory plasmas are often found in far-from-equilibrium dynamical state usually described as ``turbulence". It has been argued and discussed in the literature that the typical signatures of turbulent systems including power-law power spectral density and the scale-free statistics of fluctuating observable quantities can more or less successfully be reproduced by complex systems in the state of self-organized criticality (SOC). An obvious distinction between the theoretical concepts of turbulence and SOC has not been obtained, though (for the challenges that lie ahead, and current scientific debate, see the recently announced book on ``Self-Organized Criticality Systems" - available in open access from Open Academic Press, http://www.openacademicpress.de/). Here we discuss these issues further and show that the behavior crucially depends on the type of boundary conditions, feedback mechanisms, and the role of nonlinearity. We then apply this approach to the dynamics of Earth's geomagnetic tail and propose a hybrid model of ``turbulent" current sheet, which explicitly takes into account the self-organization processes taking place. The model yields the slope of magnetic fluctuation spectra in the near-Earth stretched magnetotail prior to the substorm below the characteristic frequency turnover scales posed by convection. A comparison between the model theoretical predictions and the available data of in situ satellite observations is given.

Characterizing Ocean Turbulence from Argo, Acoustic Doppler, and Simulation Data

NASA Astrophysics Data System (ADS)

McCaffrey, Katherine

Turbulence is inherently chaotic and unsteady, so observing it and modeling it are no easy tasks. The ocean's sheer size makes it even more difficult to observe, and its unpredictable and ever-changing forcings introduce additional complexities. Turbulence in the oceans ranges from basin scale to the scale of the molecular viscosity. The method of energy transfer between scales is, however, an area of active research, so observations of the ocean at all scales are crucial to understanding the basic dynamics of its motions. In this collection of work, I use a variety of datasets to characterize a wide range of scales of turbulence, including observations from multiple instruments and from models with different governing equations. I analyzed the largest scales of the turbulent range using the global salinity data of the Argo profiling float network. Taking advantage of the scattered and discontinuous nature of this dataset, the second-order structure function was calculated down to 2000m depth, and shown to be useful for predicting spectral slopes. Results showed structure function slopes of 2/3 at small scales, and 0 at large scales, which corresponds with spectral slopes of -5/3 at small scales, and -1 at large scales. Using acoustic Doppler velocity measurements, I characterized the meter- to kilometer-scale turbulence at a potential tidal energy site in the Puget Sound, WA. Acoustic Doppler current profiler (ADCP) and acoustic Doppler velocimeter (ADV) observations provided the data for an analysis that includes coherence, anisotropy, and intermittency. In order to more simply describe these features, a parameterization was done with four turbulence metrics, and the anisotropy magnitude, introduced here, was shown to most closely capture the coherent events. Then, using both the NREL TurbSim stochastic turbulence generator and the NCAR large-eddy simulation (LES) model, I calculated turbulence statistics to validate the accuracy of these methods in reproducing the tidal channel. TurbSim models statistics at the height of a turbine hub (5m) well, but do not model coherent events, while the LES does create these events, but not realistically in this configuration, based on comparisons with observations. Each of the datasets have disadvantages when it comes to observing turbulence. The Argo network is sparse in space, and few measurements are taken simultaneously in time. Therefore spatial and temporal averaging is needed, which requires the turbulence to be homogeneous and stationary if it is to be generalized. Though the acoustic Doppler current profiler provides a vertical profile of velocities, the fluctuations are dominated by instrument noise and beam spread, preventing it from being used for most turbulence metrics. ADV measurements have much less noise, and no beam spread, but the observations are made at one point in space, limiting us to temporal statistics or an assumption of "frozen turbulence" to infer spatial scales. As for the models, TurbSim does not have any real-world forcing, and uses parameterized spectra, and coherence functions and randomizes phase information, while LES models must make assumptions about sub-grid scales, which may be inaccurate. Additionally, all models are set up with idealizations of the forcing and domain, which may make the results unlike observations in a particular location and time. Despite these difficulties in observing and characterizing turbulence, I present several quantities that use the imperfect, yet still valuable observations, to attain a better description of the turbulence in the oceans.

A NEW SIMPLE DYNAMO MODEL FOR STELLAR ACTIVITY CYCLE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yokoi, N.; Hamba, F.; Schmitt, D.

2016-06-20

A new simple dynamo model for stellar activity cycle is proposed. By considering an inhomogeneous flow effect on turbulence, it is shown that turbulent cross helicity (velocity–magnetic-field correlation) enters the expression of turbulent electromotive force as the coupling coefficient for the mean absolute vorticity. This makes the present model different from the current α –Ω-type models in two main ways. First, in addition to the usual helicity ( α ) and turbulent magnetic diffusivity ( β ) effects, we consider the cross-helicity effect as a key ingredient of the dynamo process. Second, the spatiotemporal evolution of cross helicity is solvedmore » simultaneously with the mean magnetic fields. The basic scenario is as follows. In the presence of turbulent cross helicity, the toroidal field is induced by the toroidal rotation. Then, as in usual models, the α effect generates the poloidal field from the toroidal one. This induced poloidal field produces a turbulent cross helicity whose sign is opposite to the original one (negative production). With this cross helicity of the reversed sign, a reversal in field configuration starts. Eigenvalue analyses of the simplest possible model give a butterfly diagram, which confirms the above scenario and the equatorward migrations, the phase relationship between the cross helicity and magnetic fields. These results suggest that the oscillation of the turbulent cross helicity is a key for the activity cycle. The reversal of the cross helicity is not the result of the magnetic-field reversal, but the cause of the latter. This new model is expected to open up the possibility of the mean-field or turbulence closure dynamo approaches.« less

Wind Shear/Turbulence Inputs to Flight Simulation and Systems Certification

NASA Technical Reports Server (NTRS)

Bowles, Roland L. (Editor); Frost, Walter (Editor)

1987-01-01

The purpose of the workshop was to provide a forum for industry, universities, and government to assess current status and likely future requirements for application of flight simulators to aviation safety concerns and system certification issues associated with wind shear and atmospheric turbulence. Research findings presented included characterization of wind shear and turbulence hazards based on modeling efforts and quantitative results obtained from field measurement programs. Future research thrusts needed to maximally exploit flight simulators for aviation safety application involving wind shear and turbulence were identified. The conference contained sessions on: Existing wind shear data and simulator implementation initiatives; Invited papers regarding wind shear and turbulence simulation requirements; and Committee working session reports.

Learning to soar in turbulent environments

PubMed Central

Reddy, Gautam; Celani, Antonio; Sejnowski, Terrence J.; Vergassola, Massimo

2016-01-01

Birds and gliders exploit warm, rising atmospheric currents (thermals) to reach heights comparable to low-lying clouds with a reduced expenditure of energy. This strategy of flight (thermal soaring) is frequently used by migratory birds. Soaring provides a remarkable instance of complex decision making in biology and requires a long-term strategy to effectively use the ascending thermals. Furthermore, the problem is technologically relevant to extend the flying range of autonomous gliders. Thermal soaring is commonly observed in the atmospheric convective boundary layer on warm, sunny days. The formation of thermals unavoidably generates strong turbulent fluctuations, which constitute an essential element of soaring. Here, we approach soaring flight as a problem of learning to navigate complex, highly fluctuating turbulent environments. We simulate the atmospheric boundary layer by numerical models of turbulent convective flow and combine them with model-free, experience-based, reinforcement learning algorithms to train the gliders. For the learned policies in the regimes of moderate and strong turbulence levels, the glider adopts an increasingly conservative policy as turbulence levels increase, quantifying the degree of risk affordable in turbulent environments. Reinforcement learning uncovers those sensorimotor cues that permit effective control over soaring in turbulent environments. PMID:27482099

Learning to soar in turbulent environments.

PubMed

Reddy, Gautam; Celani, Antonio; Sejnowski, Terrence J; Vergassola, Massimo

2016-08-16

Birds and gliders exploit warm, rising atmospheric currents (thermals) to reach heights comparable to low-lying clouds with a reduced expenditure of energy. This strategy of flight (thermal soaring) is frequently used by migratory birds. Soaring provides a remarkable instance of complex decision making in biology and requires a long-term strategy to effectively use the ascending thermals. Furthermore, the problem is technologically relevant to extend the flying range of autonomous gliders. Thermal soaring is commonly observed in the atmospheric convective boundary layer on warm, sunny days. The formation of thermals unavoidably generates strong turbulent fluctuations, which constitute an essential element of soaring. Here, we approach soaring flight as a problem of learning to navigate complex, highly fluctuating turbulent environments. We simulate the atmospheric boundary layer by numerical models of turbulent convective flow and combine them with model-free, experience-based, reinforcement learning algorithms to train the gliders. For the learned policies in the regimes of moderate and strong turbulence levels, the glider adopts an increasingly conservative policy as turbulence levels increase, quantifying the degree of risk affordable in turbulent environments. Reinforcement learning uncovers those sensorimotor cues that permit effective control over soaring in turbulent environments.

Twelve mortal sins of the turbulence propagation science

NASA Astrophysics Data System (ADS)

Charnotskii, Mikhail

2011-09-01

In this review paper we discuss a series of typical mistakes and omissions that are made by engineers and scientists involved in the theoretical research and modeling of the optical propagation through atmospheric turbulence. We show how the use of the oversimplified Gaussian spectral model of turbulence delivers the completely erroneous results for the beam wander. We address a series of common omissions related to calculations of the average beam intensity: unnecessary use of the approximations when rigorous result is available, invalid application of the RMS beam size to the turbulence-distorted beams, overlooking the simple theoretical result - average beam intensity is a convolution with the turbulent Point Spread Function (PSF). We discuss the meaning and potential dangers of the use of the quadratic structure function for modeling of the turbulent perturbations. We will also address the issues related to the energy conservation principle and reciprocity that have very important consequences for the turbulence propagation, but are frequently overlooked in the current literature. We discuss a series of misconceptions that very common in of the Scintillation Index (SI) calculations. We will clarify the infamous misunderstanding of the Rytov's approximation: vanishing scintillation at the beam focus, and show the correct weak and strong scintillation solutions for the SI at the beam focus. We discuss the flaws of the Fried model of the short-term PSF, and direct to the more accurate PSF model. We will briefly review the propagation of the polarized optical waves through turbulence and discuss the inadequacy of the recently published calculations of the electromagnetic beams calculations. We discuss some common errors in representation of the calculation results for the non-Kolmogorov turbulence.

Decay of Far-Flowfield in Trailing Vortices

NASA Technical Reports Server (NTRS)

Baldwin, B. S.; Chigier, N. A.; Sheaffer, Y. S.

1973-01-01

Methods for reduction of velocities in trailing vortices of large aircraft are of current interest for the purpose of shortening the waiting time between landings at central airports. We have made finite-difference calculations of the flow in turbulent wake vortices as an aid to interpretation of wind-tunnel and flight experiments directed toward that end. Finite-difference solutions are capable of adding flexibility to such investigations if they are based on an adequate model of turbulence. Interesting developments have been taking place in the knowledge of turbulence that may lead to a complete theory in the future. In the meantime, approximate methods that yield reasonable agreement with experiment are appropriate. The simplified turbulence model we have selected contains features that account for the major effects disclosed by more sophisticated models in which the parameters are not yet established. Several puzzles are thereby resolved that arose in previous theoretical investigations of wake vortices.

Atmospheric Turbulence Relative to Aviation, Missile, and Space Programs

NASA Technical Reports Server (NTRS)

Camp, Dennis W. (Editor); Frost, Walter (Editor)

1987-01-01

The purpose of the workshop was to bring together various disciplines of the aviation, missile, and space programs involved in predicting, measuring, modeling, and understanding the processes of atmospheric turbulence. Working committees re-examined the current state of knowledge, identified present and future needs, and documented and prioritized integrated and cooperative research programs.

Turbulent Mixing Chemistry in Disks

NASA Astrophysics Data System (ADS)

Semenov, D.; Wiebe, D.

2006-11-01

A gas-grain chemical model with surface reaction and 1D/2D turbulent mixing is available for protoplanetary disks and molecular clouds. Current version is based on the updated UMIST'95 database with gas-grain interactions (accretion, desorption, photoevaporation, etc.) and modified rate equation approach to surface chemistry (see also abstract for the static chemistry code).

Phenomenological Nusselt-Rayleigh Scaling of Turbulent Thermal Convection

NASA Astrophysics Data System (ADS)

Liu, Chien-Chia

2017-12-01

Natural convection between the hot floor and the cool ceiling, so called Rayleigh-Bénard convection, is pervasive and of both fundamental and industrial interests. One key issue is how heat transfer varies with increasing thermal potential, or equivalently how the Nusselt number (Nu) scales with the Rayleigh number (Ra). The overview of experimental findings remains to show the need of extra explanation complemental to the current theories. Here we present a model based on the phenomenological theory of turbulence, where the power-law spectral exponent of the energy spectrum is the only input parameter required. The goal aims to elucidate the unexplained aspect in the Nu-Ra scaling. We find that Kolmogorov turbulence in the current model leads to Nu ˜ Ra0.3, in good agreement with the modern experimental results. We hope that this model could stimulate the discussion as to the effects of the spectral phenomena on the Nu-Ra scaling, and thus augment our understanding of buoyancy-driven thermal convection.

Application of RANS Simulations for Contact Time Predictions in Turbulent Reactor Tanks for Water Purification Process

NASA Astrophysics Data System (ADS)

Nickles, Cassandra; Goodman, Matthew; Saez, Jose; Issakhanian, Emin

2016-11-01

California's current drought has renewed public interest in recycled water from Water Reclamation Plants (WRPs). It is critical that the recycled water meets public health standards. This project consists of simulating the transport of an instantaneous conservative tracer through the WRP chlorine contact tanks. Local recycled water regulations stipulate a minimum 90-minute modal contact time during disinfection at peak dry weather design flow. In-situ testing is extremely difficult given flowrate dependence on real world sewage line supply and recycled water demand. Given as-built drawings and operation parameters, the chlorine contact tanks are modeled to simulate extreme situations, which may not meet regulatory standards. The turbulent flow solutions are used as the basis to model the transport of a turbulently diffusing conservative tracer added instantaneously to the inlet of the reactors. This tracer simulates the transport through advection and dispersion of chlorine in the WRPs. Previous work validated the models against experimental data. The current work shows the predictive value of the simulations.

SOME EXPERIMENTS ON THE TURBULENCE AND MUTUAL FRICTION IN LIQUID HELIUM (in French)

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kramers, H.C.

1962-05-01

Recent experiments on the hydrodynamic properties of liquid helium II, performed in the Kamerlingh Onnes Laboratorium (Leiden), are discussed on the basis of the twofluid model. Staas and Taconis measured temperature and pressure gradients occurring in a narrow circular capillary in the presence of a heat current. They discovered a turbulent phenomenon, very analogous to ordinary classical turbulence. The attenuation of second sound was studied by Wiarda and the present author in its dependence on a simultaneous continuous heat current. The method employed proves to be very suitable for a study of the so-called mutual friction and of phenomena occurringmore » in the critical velocity region. (auth)« less

Reduced model simulations of the scrape-off-layer heat-flux width and comparison with experiment

DOE PAGES

Myra, J. R.; Russell, D. A.; D’Ippolito, D. A.; ...

2011-01-01

Reduced model simulations of turbulence in the edge and scrape-off-layer (SOL) region of a spherical torus or tokamak plasma are employed to address the physics of the scrape-off-layer heat flux width. The simulation model is an electrostatic two-dimensional fluid turbulence model, applied in the plane perpendicular to the magnetic field at the outboard midplane of the torus. The model contains curvature-driven-interchange modes, sheath losses, and both perpendicular turbulent diffusive and convective (blob) transport. These transport processes compete with classical parallel transport to set the SOL width. Midplane SOL profiles of density, temperature and parallel heat flux are obtained from themore » simulation and compared with experimental results from the National Spherical Torus Experiment (NSTX) to study the scaling of the heat flux width with power and plasma current. It is concluded that midplane turbulence is the main contributor to the SOL heat flux width for the low power H-mode discharges studied, while additional physics is required to explain the plasma current scaling of the SOL heat flux width observed experimentally in higher power discharges. Intermittent separatrix spanning convective cells are found to be the main mechanism that sets the near-SOL width in the simulations. The roles of sheared flows and blob trapping vs. emission are discussed.« less

Reynolds Stress Distributions and the Measurement and Calculation of Eddy Viscosity in Gravity Currents

NASA Astrophysics Data System (ADS)

Kelly, R. W.; Chalk, C.; Dorrell, R. M.; Peakall, J.; Burns, A. D.; Keevil, G. M.; Thomas, R. E.; Williams, G.

2016-12-01

In the natural environment, gravity currents transport large volumes of sediment great distances and are often considered one of the most important mechanisms for sediment transport in ocean basins. Deposits from many individual submarine gravity currents, turbidites, ultimately form submarine fan systems. These are the largest sedimentary systems on the planet and contain valuable hydrocarbon reserves. Moreover, the impact of these currents on submarine technologies and seafloor infrastructure can be devastating and therefore they are of significant interest to a wide range of industries. Here we present experimental, numerical and theoretical models of time-averaged turbulent shear stresses, i.e. Reynolds stresses. Reynolds stresses can be conceptually parameterised by an eddy viscosity parameter that relates chaotic fluid motion to diffusive type processes. As such, it is a useful parameter for indicating the extent of internal mixing and is used extensively in both numerical and analytical modelling of both open-channel and gravity driven flows. However, a lack of knowledge of the turbulent structure of gravity currents limits many hydro- and morphodynamic models. High resolution 3-dimensional experimental velocity data, gathered using acoustic Doppler profiling velocimetry, enabled direct calculation of stresses and eddy viscosity. Comparison of experimental data to CFD and analytical models allowed the testing of eddy viscosity-based turbulent mixing models. The calculated eddy viscosity profile is parabolic in nature in both the upper and lower shear layers. However, an apparent breakdown in the Boussinesq hypothesis (used to calculate the eddy viscosity and upon which many numerical models are based) is observed in the region of the current around the velocity maximum. With the use of accompanying density data it is suggested that the effect of stratification on eddy viscosity is significant and alternative formulations may be required.

Turbulent structures in cylindrical density currents in a rotating frame of reference

NASA Astrophysics Data System (ADS)

Salinas, Jorge S.; Cantero, Mariano I.; Dari, Enzo A.; Bonometti, Thomas

2018-06-01

Gravity currents are flows generated by the action of gravity on fluids with different densities. In some geophysical applications, modeling such flows makes it necessary to account for rotating effects, modifying the dynamics of the flow. While previous works on rotating stratified flows focused on currents of large Coriolis number, the present work focuses on flows with small Coriolis numbers (i.e. moderate-to-large Rossby numbers). In this work, cylindrical rotating gravity currents are investigated by means of highly resolved simulations. A brief analysis of the mean flow evolution to the final state is presented to provide a complete picture of the flow dynamics. The numerical results, showing the well-known oscillatory behavior of the flow (inertial waves) and a final state lens shape (geostrophic adjustment), are in good agreement with experimental observations and theoretical models. The turbulent structures in the flow are visualized and described using, among others, a stereoscopic visualization and videos as supplementary material. In particular, the structure of the lobes and clefts at the front of the current is presented in association to local turbulent structures. In rotating gravity currents, the vortices observed at the lobes front are not of hairpin type but are rather of Kelvin-Helmholtz type.

Optical turbulence on underwater image degradation in natural environments.

PubMed

Hou, Weilin; Woods, Sarah; Jarosz, Ewa; Goode, Wesley; Weidemann, Alan

2012-05-10

It is a well-known fact that the major degradation source on electro-optical imaging underwater is from scattering by particles of various origins and sizes. Recent research indicates that, under certain conditions, the apparent degradation could also be caused by the variations of index of refraction associated with temperature and salinity microstructures in the ocean and lakes. The combined impact has been modeled previously through the simple underwater imaging model. The current study presents the first attempts in quantifying the level of image degradation due to optical turbulence in natural waters in terms of modulation transfer functions using measured turbulence dissipation rates. Image data collected from natural environments during the Skaneateles Optical Turbulence Exercise are presented. Accurate assessments of the turbulence conditions are critical to the model validation and were measured by two instruments to ensure consistency and accuracy. Optical properties of the water column in the field were also measured in coordination with temperature, conductivity, and depth. The results show that optical turbulence degrades the image quality as predicted and on a level comparable to that caused by the particle scattering just above the thermocline. Other contributing elements involving model closure, including temporal and spatial measurement scale differences among sensors and mitigation efforts, are discussed.

Evolution of three-dimensional relativistic current sheets and development of self-generated turbulence

NASA Astrophysics Data System (ADS)

Takamoto, M.

2018-05-01

In this paper, the temporal evolution of three-dimensional relativistic current sheets in Poynting-dominated plasma is studied for the first time. Over the past few decades, a lot of efforts have been conducted on studying the evolution of current sheets in two-dimensional space, and concluded that sufficiently long current sheets always evolve into the so-called plasmoid chain, which provides a fast reconnection rate independent of its resistivity. However, it is suspected that plasmoid chain can exist only in the case of two-dimensional approximation, and would show transition to turbulence in three-dimensional space. We performed three-dimensional numerical simulation of relativistic current sheet using resistive relativistic magnetohydrodynamic approximation. The results showed that the three-dimensional current sheets evolve not into plasmoid chain but turbulence. The resulting reconnection rate is 0.004, which is much smaller than that of plasmoid chain. The energy conversion from magnetic field to kinetic energy of turbulence is just 0.01 per cent, which is much smaller than typical non-relativistic cases. Using the energy principle, we also showed that the plasmoid is always unstable for a displacement in the opposite direction to its acceleration, probably interchange-type instability, and this always results in seeds of turbulence behind the plasmoids. Finally, the temperature distribution along the sheet is discussed, and it is found that the sheet is less active than plasmoid chain. Our finding can be applied for many high-energy astrophysical phenomena, and can provide a basic model of the general current sheet in Poynting-dominated plasma.

Development of a Hybrid RANS/LES Method for Turbulent Mixing Layers

NASA Technical Reports Server (NTRS)

Georgiadis, Nicholas J.; Alexander, J. Iwan D.; Reshotko, Eli

2001-01-01

Significant research has been underway for several years in NASA Glenn Research Center's nozzle branch to develop advanced computational methods for simulating turbulent flows in exhaust nozzles. The primary efforts of this research have concentrated on improving our ability to calculate the turbulent mixing layers that dominate flows both in the exhaust systems of modern-day aircraft and in those of hypersonic vehicles under development. As part of these efforts, a hybrid numerical method was recently developed to simulate such turbulent mixing layers. The method developed here is intended for configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. Interest in Large Eddy Simulation (LES) methods have increased in recent years, but applying an LES method to calculate the wide range of turbulent scales from small eddies in the wall-bounded regions to large eddies in the mixing region is not yet possible with current computers. As a result, the hybrid method developed here uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall-bounded regions entering a mixing section and uses a LES procedure to calculate the mixing-dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. With this technique, closure for the RANS equations is obtained by using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The LES equations are closed using the Smagorinsky subgrid scale model. Although the function of the Cebeci-Smith model to replace all of the turbulent stresses is quite different from that of the Smagorinsky subgrid model, which only replaces the small subgrid turbulent stresses, both are eddy viscosity models and both are derived at least in part from mixing-length theory. The similar formulation of these two models enables the RANS and LES equations to be solved with a single solution scheme and computational grid. The hybrid RANS-LES method has been applied to a benchmark compressible mixing layer experiment in which two isolated supersonic streams, separated by a splitter plate, provide the flows to a constant-area mixing section. Although the configuration is largely two dimensional in nature, three-dimensional calculations were found to be necessary to enable disturbances to develop in three spatial directions and to transition to turbulence. The flow in the initial part of the mixing section consists of a periodic vortex shedding downstream of the splitter plate trailing edge. This organized vortex shedding then rapidly transitions to a turbulent structure, which is very similar to the flow development observed in the experiments. Although the qualitative nature of the large-scale turbulent development in the entire mixing section is captured well by the LES part of the current hybrid method, further efforts are planned to directly calculate a greater portion of the turbulence spectrum and to limit the subgrid scale modeling to only the very small scales. This will be accomplished by the use of higher accuracy solution schemes and more powerful computers, measured both in speed and memory capabilities.

Transport modeling of the DIII-D high $${{\\beta}_{p}}$$ scenario and extrapolations to ITER steady-state operation

DOE PAGES

McClenaghan, Joseph; Garofalo, Andrea M.; Meneghini, Orso; ...

2017-08-03

In this study, transport modeling of a proposed ITER steady-state scenario based on DIII-D high poloidal-beta (more » $${{\\beta}_{p}}$$ ) discharges finds that ITB formation can occur with either sufficient rotation or a negative central shear q-profile. The high $${{\\beta}_{p}}$$ scenario is characterized by a large bootstrap current fraction (80%) which reduces the demands on the external current drive, and a large radius internal transport barrier which is associated with excellent normalized confinement. Modeling predictions of the electron transport in the high $${{\\beta}_{p}}$$ scenario improve as $${{q}_{95}}$$ approaches levels similar to typical existing models of ITER steady-state and the ion transport is turbulence dominated. Typical temperature and density profiles from the non-inductive high $${{\\beta}_{p}}$$ scenario on DIII-D are scaled according to 0D modeling predictions of the requirements for achieving a $Q=5$ steady-state fusion gain in ITER with 'day one' heating and current drive capabilities. Then, TGLF turbulence modeling is carried out under systematic variations of the toroidal rotation and the core q-profile. A high bootstrap fraction, high $${{\\beta}_{p}}$$ scenario is found to be near an ITB formation threshold, and either strong negative central magnetic shear or rotation in a high bootstrap fraction are found to successfully provide the turbulence suppression required to achieve $Q=5$.« less

Turbulence Modeling Workshop

NASA Technical Reports Server (NTRS)

Rubinstein, R. (Editor); Rumsey, C. L. (Editor); Salas, M. D. (Editor); Thomas, J. L. (Editor); Bushnell, Dennis M. (Technical Monitor)

2001-01-01

Advances in turbulence modeling are needed in order to calculate high Reynolds number flows near the onset of separation and beyond. To this end, the participants in this workshop made the following recommendations. (1) A national/international database and standards for turbulence modeling assessment should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be taken of other efforts already under-way, such as that of the European Research Community on Flow, Turbulence, and Combustion (ERCOFTAC) consortium. Carefully selected "unit" experiments will be needed, as well as advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least understood and is a key component of two-equation and higher models. Second priority should be given to the development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes (RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in turbulence modeling development, validation, and implementation should include these approaches as well.

Flare Energy Release: Internal Conflict, Contradiction with High Resolution Observations, Possible Solutions

NASA Astrophysics Data System (ADS)

Pustilnik, L.

2017-06-01

All accepted paradigm of solar and stellar flares energy release based on 2 whales: 1. Source of energy is free energy of non-potential force free magnetic field in atmosphere above active region; 2. Process of ultrafast dissipation of magnetic fields is Reconnection in Thin Turbulent Current Sheet (RTTCS). Progress in observational techniques in last years provided ultra-high spatial resolution and in physics of turbulent plasma showed that real situation is much more complicated and standard approach is in contradiction both with observations and with problem of RTTCS stability. We present critical analysis of classic models of pre-flare energy accumulation and its dissipation during flare energy release from pioneer works Giovanelli (1939, 1947) up to topological reconnection. We show that all accepted description of global force-free fields as source of future flare cannot be agreed with discovered in last years fine and ultra-fine current-magnetic structure included numerouse arcs-threads with diameters up to 100 km with constant sequence from photosphere to corona. This magnetic skeleton of thin current magnetic threads with strong interaction between them is main source of reserved magnetic energy insolar atmosphere. Its dynamics will be controlled by percolation of magnetic stresses through network of current-magnetic threads with transition to flare state caused by critical value of global current. We show that thin turbulent current sheet is absolutely unstable configuration both caused by splitting to numerous linear currents by dissipative modes like to tearing, and as sequence of suppress of plasma turbulence caused by anomalous heating of turbulent plasma. In result of these factors primary RTTCS will be disrupted in numerous turbulent and normal plasma domains like to resistors network. Current propagation through this network will have percolation character with all accompanied properties of percolated systems: self-organization with formation power spectrum of distribution of flares and micro-flares, and possibility of phase transition to flare energy release with huge increasing of energy release.

Boundary layer integral matrix procedure: Verification of models

NASA Technical Reports Server (NTRS)

Bonnett, W. S.; Evans, R. M.

1977-01-01

The three turbulent models currently available in the JANNAF version of the Aerotherm Boundary Layer Integral Matrix Procedure (BLIMP-J) code were studied. The BLIMP-J program is the standard prediction method for boundary layer effects in liquid rocket engine thrust chambers. Experimental data from flow fields with large edge-to-wall temperature ratios are compared to the predictions of the three turbulence models contained in BLIMP-J. In addition, test conditions necessary to generate additional data on a flat plate or in a nozzle are given. It is concluded that the Cebeci-Smith turbulence model be the recommended model for the prediction of boundary layer effects in liquid rocket engines. In addition, the effects of homogeneous chemical reaction kinetics were examined for a hydrogen/oxygen system. Results show that for most flows, kinetics are probably only significant for stoichiometric mixture ratios.

Physical consistency of subgrid-scale models for large-eddy simulation of incompressible turbulent flows

NASA Astrophysics Data System (ADS)

Silvis, Maurits H.; Remmerswaal, Ronald A.; Verstappen, Roel

2017-01-01

We study the construction of subgrid-scale models for large-eddy simulation of incompressible turbulent flows. In particular, we aim to consolidate a systematic approach of constructing subgrid-scale models, based on the idea that it is desirable that subgrid-scale models are consistent with the mathematical and physical properties of the Navier-Stokes equations and the turbulent stresses. To that end, we first discuss in detail the symmetries of the Navier-Stokes equations, and the near-wall scaling behavior, realizability and dissipation properties of the turbulent stresses. We furthermore summarize the requirements that subgrid-scale models have to satisfy in order to preserve these important mathematical and physical properties. In this fashion, a framework of model constraints arises that we apply to analyze the behavior of a number of existing subgrid-scale models that are based on the local velocity gradient. We show that these subgrid-scale models do not satisfy all the desired properties, after which we explain that this is partly due to incompatibilities between model constraints and limitations of velocity-gradient-based subgrid-scale models. However, we also reason that the current framework shows that there is room for improvement in the properties and, hence, the behavior of existing subgrid-scale models. We furthermore show how compatible model constraints can be combined to construct new subgrid-scale models that have desirable properties built into them. We provide a few examples of such new models, of which a new model of eddy viscosity type, that is based on the vortex stretching magnitude, is successfully tested in large-eddy simulations of decaying homogeneous isotropic turbulence and turbulent plane-channel flow.

A detached eddy simulation model for the study of lateral separation zones along a large canyon-bound river

NASA Astrophysics Data System (ADS)

Alvarez, Laura V.; Schmeeckle, Mark W.; Grams, Paul E.

2017-01-01

Lateral flow separation occurs in rivers where banks exhibit strong curvature. In canyon-bound rivers, lateral recirculation zones are the principal storage of fine-sediment deposits. A parallelized, three-dimensional, turbulence-resolving model was developed to study the flow structures along lateral separation zones located in two pools along the Colorado River in Marble Canyon. The model employs the detached eddy simulation (DES) technique, which resolves turbulence structures larger than the grid spacing in the interior of the flow. The DES-3D model is validated using Acoustic Doppler Current Profiler flow measurements taken during the 2008 controlled flood release from Glen Canyon Dam. A point-to-point validation using a number of skill metrics, often employed in hydrological research, is proposed here for fluvial modeling. The validation results show predictive capabilities of the DES model. The model reproduces the pattern and magnitude of the velocity in the lateral recirculation zone, including the size and position of the primary and secondary eddy cells, and return current. The lateral recirculation zone is open, having continuous import of fluid upstream of the point of reattachment and export by the recirculation return current downstream of the point of separation. Differences in magnitude and direction of near-bed and near-surface velocity vectors are found, resulting in an inward vertical spiral. Interaction between the recirculation return current and the main flow is dynamic, with large temporal changes in flow direction and magnitude. Turbulence structures with a predominately vertical axis of vorticity are observed in the shear layer becoming three-dimensional without preferred orientation downstream.

Development of a coupled wave-flow-vegetation interaction model

USGS Publications Warehouse

Beudin, Alexis; Kalra, Tarandeep S.; Ganju, Neil K.; Warner, John C.

2017-01-01

Emergent and submerged vegetation can significantly affect coastal hydrodynamics. However, most deterministic numerical models do not take into account their influence on currents, waves, and turbulence. In this paper, we describe the implementation of a wave-flow-vegetation module into a Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system that includes a flow model (ROMS) and a wave model (SWAN), and illustrate various interacting processes using an idealized shallow basin application. The flow model has been modified to include plant posture-dependent three-dimensional drag, in-canopy wave-induced streaming, and production of turbulent kinetic energy and enstrophy to parameterize vertical mixing. The coupling framework has been updated to exchange vegetation-related variables between the flow model and the wave model to account for wave energy dissipation due to vegetation. This study i) demonstrates the validity of the plant posture-dependent drag parameterization against field measurements, ii) shows that the model is capable of reproducing the mean and turbulent flow field in the presence of vegetation as compared to various laboratory experiments, iii) provides insight into the flow-vegetation interaction through an analysis of the terms in the momentum balance, iv) describes the influence of a submerged vegetation patch on tidal currents and waves separately and combined, and v) proposes future directions for research and development.

Turbulence and sediment transport over sand dunes and ripples

NASA Astrophysics Data System (ADS)

Bennis, A.; Le Bot, S.; lafite, R.; Bonneton, P.; Ardhuin, F.

2013-12-01

Several bedforms are present near to the surfzone of natural beaches. Dunes and ripples are frequently observed. Understanding the turbulence over these forms is essential for the sediment transport. The turbulent flow and the suspended sand particles interact with each other. At the moment, the modelling strategy for turbulence is still a challenge. According to the spatial scales, some different methods to model the turbulence are employed, in particular the RANS (Reynolds Averaged Navier-Stokes) and the LES (Large Eddy Simulation). A hybrid method combining both RANS and LES is set up here. We have adapted this method, initially developed for atmospheric flow, to the oceanic flow. This new method is implemented inside the 3D hydrodynamic model, MARS 3D, which is forced by waves. LES is currently the best way to simulate turbulent flow but its higher cost prevents it from being used for large scale applications. So, here we use RANS near the bottom while LES is set elsewhere. It allows us minimize the computational cost and ensure a better accuracy of the results than with a fully RANS model. In the case of megaripples, the validation step was performed with two sets of field data (Sandy Duck'97 and Forsoms'13) but also with the data from Dune2D model which uses only RANS for turbulence. The main findings are: a) the vertical profiles of the velocity are similar throughout the data b) the turbulent kinetic energy, which was underestimated by Dune2D, is in line with the observations c) the concentration of the suspended sediment is simulated with a better accuracy than with Dune2D but this remains lower than the observations.

Characterizing Turbulent Events at a Tidal Energy Site from Acoustic Doppler Velocity Observations

NASA Astrophysics Data System (ADS)

McCaffrey, Katherine; Fox-Kemper, Baylor; Hamlington, Peter

2013-11-01

As interest in marine renewable energy increases, observations are crucial to understanding the environments encountered by energy conversion devices. Data obtained from an acoustic Doppler current profiler and an acoustic Doppler velocimeter at two locations in the Puget Sound, WA are used to perform a detailed analysis of the turbulent environment that is expected to be present at a turbine placed in a tidal strait. Metrics such as turbulence intensity, structure functions, probability density functions, intermittency, coherent turbulence kinetic energy, anisotropy invariants, and linear combinations of eigenvalues are used to characterize the turbulence. The results indicate that coherent turbulence kinetic energy and turbulence intensity can be used to identify and parameterize different turbulent events in the flow. An analysis of the anisotropy characteristics leads to a physical description of turbulent events (defined using both turbulence intensity and coherent turbulent kinetic energy) as being dominated by one component of the Reynolds stresses. During non-turbulent events, the flow is dominated by two Reynolds stress components. The importance of these results for the development of realistic models of energy conversion devices is outlined. Cooperative Institute for Research in Environmental Sciences, Department of Atmospheric and Oceanic Sciences.

ANALYTIC FORMS OF THE PERPENDICULAR DIFFUSION COEFFICIENT IN NRMHD TURBULENCE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Shalchi, A., E-mail: andreasm4@yahoo.com

2015-02-01

In the past different analytic limits for the perpendicular diffusion coefficient of energetic particles interacting with magnetic turbulence were discussed. These different limits or cases correspond to different transport modes describing how the particles are diffusing across the large-scale magnetic field. In the current paper we describe a new transport regime by considering the model of noisy reduced magnetohydrodynamic turbulence. We derive different analytic forms of the perpendicular diffusion coefficient, and while we do this, we focus on the aforementioned new transport mode. We show that for this turbulence model a small perpendicular diffusion coefficient can be obtained so thatmore » the latter diffusion coefficient is more than hundred times smaller than the parallel diffusion coefficient. This result is relevant to explain observations in the solar system where such small perpendicular diffusion coefficients have been reported.« less

The Numerical Analysis of a Turbulent Compressible Jet. Degree awarded by Ohio State Univ., 2000

NASA Technical Reports Server (NTRS)

DeBonis, James R.

2001-01-01

A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).

Development of an algebraic stress/two-layer model for calculating thrust chamber flow fields

NASA Technical Reports Server (NTRS)

Chen, C. P.; Shang, H. M.; Huang, J.

1993-01-01

Following the consensus of a workshop in Turbulence Modeling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data, to account for the non-isotropic turbulence effects.

Turbulence Modeling Verification and Validation

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.

2014-01-01

Computational fluid dynamics (CFD) software that solves the Reynolds-averaged Navier-Stokes (RANS) equations has been in routine use for more than a quarter of a century. It is currently employed not only for basic research in fluid dynamics, but also for the analysis and design processes in many industries worldwide, including aerospace, automotive, power generation, chemical manufacturing, polymer processing, and petroleum exploration. A key feature of RANS CFD is the turbulence model. Because the RANS equations are unclosed, a model is necessary to describe the effects of the turbulence on the mean flow, through the Reynolds stress terms. The turbulence model is one of the largest sources of uncertainty in RANS CFD, and most models are known to be flawed in one way or another. Alternative methods such as direct numerical simulations (DNS) and large eddy simulations (LES) rely less on modeling and hence include more physics than RANS. In DNS all turbulent scales are resolved, and in LES the large scales are resolved and the effects of the smallest turbulence scales are modeled. However, both DNS and LES are too expensive for most routine industrial usage on today's computers. Hybrid RANS-LES, which blends RANS near walls with LES away from walls, helps to moderate the cost while still retaining some of the scale-resolving capability of LES, but for some applications it can still be too expensive. Even considering its associated uncertainties, RANS turbulence modeling has proved to be very useful for a wide variety of applications. For example, in the aerospace field, many RANS models are considered to be reliable for computing attached flows. However, existing turbulence models are known to be inaccurate for many flows involving separation. Research has been ongoing for decades in an attempt to improve turbulence models for separated and other nonequilibrium flows. When developing or improving turbulence models, both verification and validation are important steps in the process. Verification insures that the CFD code is solving the equations as intended (no errors in the implementation). This is typically done either through the method of manufactured solutions (MMS) or through careful step-by-step comparisons with other verified codes. After the verification step is concluded, validation is performed to document the ability of the turbulence model to represent different types of flow physics. Validation can involve a large number of test case comparisons with experiments, theory, or DNS. Organized workshops have proved to be valuable resources for the turbulence modeling community in its pursuit of turbulence modeling verification and validation. Workshop contributors using different CFD codes run the same cases, often according to strict guidelines, and compare results. Through these comparisons, it is often possible to (1) identify codes that have likely implementation errors, and (2) gain insight into the capabilities and shortcomings of different turbulence models to predict the flow physics associated with particular types of flows. These are valuable lessons because they help bring consistency to CFD codes by encouraging the correction of faulty programming and facilitating the adoption of better models. They also sometimes point to specific areas needed for improvement in the models. In this paper, several recent workshops are summarized primarily from the point of view of turbulence modeling verification and validation. Furthermore, the NASA Langley Turbulence Modeling Resource website is described. The purpose of this site is to provide a central location where RANS turbulence models are documented, and test cases, grids, and data are provided. The goal of this paper is to provide an abbreviated survey of turbulence modeling verification and validation efforts, summarize some of the outcomes, and give some ideas for future endeavors in this area.

Effects of anisotropic turbulent thermal diffusion on spherical magnetoconvection in the Earth's core

NASA Astrophysics Data System (ADS)

Ivers, D. J.; Phillips, C. G.

2018-03-01

We re-consider the plate-like model of turbulence in the Earth's core, proposed by Braginsky and Meytlis (1990), and show that it is plausible for core parameters not only in polar regions but extends to mid- and low-latitudes where rotation and gravity are not parallel, except in a very thin equatorial layer. In this model the turbulence is highly anisotropic with preferred directions imposed by the Earth's rotation and the magnetic field. Current geodynamo computations effectively model sub-grid scale turbulence by using isotropic viscous and thermal diffusion values significantly greater than the molecular values of the Earth's core. We consider a local turbulent dynamo model for the Earth's core in which the mean magnetic field, velocity and temperature satisfy the Boussinesq induction, momentum and heat equations with an isotropic turbulent Ekman number and Roberts number. The anisotropy is modelled only in the thermal diffusion tensor with the Earth's rotation and magnetic field as preferred directions. Nonlocal organising effects of gravity and rotation (but not aspect ratio in the Earth's core) such as an inverse cascade and nonlocal transport are assumed to occur at longer length scales, which computations may accurately capture with sufficient resolution. To investigate the implications of this anisotropy for the proposed turbulent dynamo model we investigate the linear instability of turbulent magnetoconvection on length scales longer than the background turbulence in a rotating sphere with electrically insulating exterior for no-slip and isothermal boundary conditions. The equations are linearised about an axisymmetric basic state with a conductive temperature, azimuthal magnetic field and differential rotation. The basic state temperature is a function of the anisotropy and the spherical radius. Elsasser numbers in the range 1-20 and turbulent Roberts numbers 0.01-1 are considered for both equatorial symmetries of the magnetic basic state. It is found that anisotropic turbulent thermal diffusivity has a strong destabilising effect on magneto-convective instabilities, which may relax the tight energy budget constraining geodynamo models. The enhanced instability is not due to a reduction of the total diffusivity. The anisotropy also strengthens instabilities which break the symmetry of the underlying state, which may facilitate magnetic field reversal. Geostrophic flow appears to suppress the symmetry breaking modes and magnetic instabilities. Through symmetry breaking and the geostrophic flow the anisotropy may provide a mechanism of magnetic field reversal and its suppression in computational dynamo models.

Effect of roughness formulation on the performance of a coupled wave, hydrodynamic, and sediment transport model

USGS Publications Warehouse

Ganju, Neil K.; Sherwood, Christopher R.

2010-01-01

A variety of algorithms are available for parameterizing the hydrodynamic bottom roughness associated with grain size, saltation, bedforms, and wave–current interaction in coastal ocean models. These parameterizations give rise to spatially and temporally variable bottom-drag coefficients that ostensibly provide better representations of physical processes than uniform and constant coefficients. However, few studies have been performed to determine whether improved representation of these variable bottom roughness components translates into measurable improvements in model skill. We test the hypothesis that improved representation of variable bottom roughness improves performance with respect to near-bed circulation, bottom stresses, or turbulence dissipation. The inner shelf south of Martha’s Vineyard, Massachusetts, is the site of sorted grain-size features which exhibit sharp alongshore variations in grain size and ripple geometry over gentle bathymetric relief; this area provides a suitable testing ground for roughness parameterizations. We first establish the skill of a nested regional model for currents, waves, stresses, and turbulent quantities using a uniform and constant roughness; we then gauge model skill with various parameterization of roughness, which account for the influence of the wave-boundary layer, grain size, saltation, and rippled bedforms. We find that commonly used representations of ripple-induced roughness, when combined with a wave–current interaction routine, do not significantly improve skill for circulation, and significantly decrease skill with respect to stresses and turbulence dissipation. Ripple orientation with respect to dominant currents and ripple shape may be responsible for complicating a straightforward estimate of the roughness contribution from ripples. In addition, sediment-induced stratification may be responsible for lower stresses than predicted by the wave–current interaction model.

In situ measurements of wind and current speed and relationship between output power and turbulence

NASA Astrophysics Data System (ADS)

Duran Medina, Olmo; Schmitt, François G.; Sentchev, Alexei; Calif, Rudy

2015-04-01

In a context of energy transition, wind and tidal energy are sources of clean energy with the potential of partially satisfying the growing demand. The main problem of this type of energy, and other types of renewable energy remains the discontinuity of the electric power produced in different scales, inducing large fluctuations also called intermittency. This intermittency of wind and tidal energy is inherent to the turbulent nature of wind and marine currents. We consider this intermittent power production in strong relation with the turbulent intermittency of the resource. The turbulence theory is multifractal energy cascades models, a classic in physics of turbulence. From earlier studies in atmospheric sciences, we learn that wind speed and the aggregate power output are intermittent and multifractal over a wide range of scales [Calif and Schmitt 2014]. We want to extend this study to a marine current turbine and compare the scaling properties for those renewable energy sources. We consider here coupling between simultaneous velocity time series and output power from a wind turbine and a marine current turbine. Wind turbine data were obtained from Denmark and marine current data from Western Scheldt, Belgium where a prototype of a vertical and horizontal marine current turbines are tested. After an estimation of their Fourier density power spectra, we study their scaling properties in Kolmogorov's theory and the framework of fully developed turbulence. Hence, we employ a Hilbert-based methodology, namely arbitrary-order Hilbert spectral analysis [Calif et al. 2013a, 2013b] to characterize the intermittent property of the wind and marine current velocity in order to characterize the intermittent nature of the fluid. This method is used in order to obtain the spectrum and the corresponding power law for non-linear and non-stationary time series. The goal is to study the non-linear transfer characteristics in a multi-scale and multi-intensity framework.

Preliminary testing of turbulence and radionuclide transport modeling in deep ocean environment

DOE Office of Scientific and Technical Information (OSTI.GOV)

Onishi, Y.; Dummuller, D.C.; Trent, D.S.

Pacific Northwest Laboratory (PNL) performed a study for the US Environmental Protection Agency's Office of Radiation Programs to (1) identify candidate models for regional modeling of low-level waste ocean disposal sites in the mid-Atlantic ocean; (2) evaluate mathematical representation of the model's eddy viscosity/dispersion coefficients; and (3) evaluate the adequacy of the k-{epsilon} turbulence model and the feasibility of one of the candidate models, TEMPEST{copyright}/FLESCOT{copyright}, to deep-ocean applications on a preliminary basis. PNL identified the TEMPEST{copyright}/FLESCOT{copyright}, FLOWER, Blumberg's, and RMA 10 models as appropriate candidates for the regional radionuclide modeling. Among these models, TEMPEST/FLESCOT is currently the only model thatmore » solves distributions of flow, turbulence (with the k-{epsilon} model), salinity, water temperature, sediment, dissolved contaminants, and sediment-sorbed contaminants. Solving the Navier-Stokes equations using higher order correlations is not practical for regional modeling because of the prohibitive computational requirements; therefore, the turbulence modeling is a more practical approach. PNL applied the three-dimensional code, TEMPEST{copyright}/FLESCOT{copyright} with the k-{epsilon} model, to a very simple, hypothetical, two-dimensional, deep-ocean case, producing at least qualitatively appropriate results. However, more detailed testing should be performed for the further testing of the code. 46 refs., 39 figs., 6 tabs.« less

Turbulent mixing and removal of ozone within an Amazon rainforest canopy

NASA Astrophysics Data System (ADS)

Freire, L. S.; Gerken, T.; Ruiz-Plancarte, J.; Wei, D.; Fuentes, J. D.; Katul, G. G.; Dias, N. L.; Acevedo, O. C.; Chamecki, M.

2017-03-01

Simultaneous profiles of turbulence statistics and mean ozone mixing ratio are used to establish a relation between eddy diffusivity and ozone mixing within the Amazon forest. A one-dimensional diffusion model is proposed and used to infer mixing time scales from the eddy diffusivity profiles. Data and model results indicate that during daytime conditions, the upper (lower) half of the canopy is well (partially) mixed most of the time and that most of the vertical extent of the forest can be mixed in less than an hour. During nighttime, most of the canopy is predominantly poorly mixed, except for periods with bursts of intermittent turbulence. Even though turbulence is faster than chemistry during daytime, both processes have comparable time scales in the lower canopy layers during nighttime conditions. Nonchemical loss time scales (associated with stomatal uptake and dry deposition) for the entire forest are comparable to turbulent mixing time scale in the lower canopy during the day and in the entire canopy during the night, indicating a tight coupling between turbulent transport and dry deposition and stomatal uptake processes. Because of the significant time of day and height variability of the turbulent mixing time scale inside the canopy, it is important to take it into account when studying chemical and biophysical processes happening in the forest environment. The method proposed here to estimate turbulent mixing time scales is a reliable alternative to currently used models, especially for situations in which the vertical distribution of the time scale is relevant.

Applying an economical scale-aware PDF-based turbulence closure model in NOAA NCEP GCMs.

NASA Astrophysics Data System (ADS)

Krueger, S. K.; Belochitski, A.; Moorthi, S.; Bogenschutz, P.; Pincus, R.

2015-12-01

A novel unified representation of sub-grid scale (SGS) turbulence, cloudiness, and shallow convection is being implemented into the NOAA NCEP Global Forecasting System (GFS) general circulation model. The approach, known as Simplified High Order Closure (SHOC), is based on predicting a joint PDF of SGS thermodynamic variables and vertical velocity and using it to diagnose turbulent diffusion coefficients, SGS fluxes, condensation and cloudiness. Unlike other similar methods, only one new prognostic variable, turbulent kinetic energy (TKE), needs to be intoduced, making the technique computationally efficient.SHOC code was adopted for a global model environment from its origins in a cloud resolving model, and incorporated into NCEP GFS. SHOC was first tested in a non-interactive mode, a configuration where SHOC receives inputs from the host model, but its outputs are not returned to the GFS. In this configuration: a) SGS TKE values produced by GFS SHOC are consistent with those produced by SHOC in a CRM, b) SGS TKE in GFS SHOC exhibits a well defined diurnal cycle, c) there's enhanced boundary layer turbulence in the subtropical stratocumulus and tropical transition-to-cumulus areas d) buoyancy flux diagnosed from the assumed PDF is consistent with independently calculated Brunt-Vaisala frequency in identifying stable and unstable regions.Next, SHOC was coupled to GFS, namely turbulent diffusion coefficients computed by SHOC are now used in place of those currently produced by the GFS boundary layer and shallow convection schemes (Han and Pan, 2011), as well as condensation and cloud fraction diagnosed from the SGS PDF replace those calculated in the current large-scale cloudines scheme (Zhao and Carr, 1997). Ongoing activities consist of debugging the fully coupled GFS/SHOC.Future work will consist of evaluating model performance and tuning the physics if necessary, by performing medium-range NWP forecasts with prescribed initial conditions, and AMIP-type climate tests with prescribed SSTs. Depending on the results, the model will be tuned or parameterizations modified. Next, SHOC will be implemented in the NCEP CFS, and tuned and evaluated for climate applications - seasonal prediction and long coupled climate runs. Impact of new physics on ENSO, MJO, ISO, monsoon variability, etc will be examined.

A Hybrid RANS/LES Approach for Predicting Jet Noise

NASA Technical Reports Server (NTRS)

Goldstein, Marvin E.

2006-01-01

Hybrid acoustic prediction methods have an important advantage over the current Reynolds averaged Navier-Stokes (RANS) based methods in that they only involve modeling of the relatively universal subscale motion and not the configuration dependent larger scale turbulence. Unfortunately, they are unable to account for the high frequency sound generated by the turbulence in the initial mixing layers. This paper introduces an alternative approach that directly calculates the sound from a hybrid RANS/LES flow model (which can resolve the steep gradients in the initial mixing layers near the nozzle lip) and adopts modeling techniques similar to those used in current RANS based noise prediction methods to determine the unknown sources in the equations for the remaining unresolved components of the sound field. The resulting prediction method would then be intermediate between the current noise prediction codes and previously proposed hybrid noise prediction methods.

Influence of turbulence on the wake of a marine current turbine simulator.

PubMed

Blackmore, T; Batten, W M J; Bahaj, A S

2014-10-08

Marine current turbine commercial prototypes have now been deployed and arrays of multiple turbines under design. The tidal flows in which they operate are highly turbulent, but the characteristics of the inflow turbulence have not being considered in present design methods. This work considers the effects of inflow turbulence on the wake behind an actuator disc representation of a marine current turbine. Different turbulence intensities and integral length scales were generated in a large eddy simulation using a gridInlet, which produces turbulence from a grid pattern on the inlet boundary. The results highlight the significance of turbulence on the wake profile, with a different flow regime occurring for the zero turbulence case. Increasing the turbulence intensity reduced the velocity deficit and shifted the maximum deficit closer to the turbine. Increasing the integral length scale increased the velocity deficit close to the turbine due to an increased production of turbulent energy. However, the wake recovery was increased due to the higher rate of turbulent mixing causing the wake to expand. The implication of this work is that marine current turbine arrays could be further optimized, increasing the energy yield of the array when the site-specific turbulence characteristics are considered.

Influence of turbulence on the wake of a marine current turbine simulator

PubMed Central

Blackmore, T.; Batten, W. M. J.; Bahaj, A. S.

2014-01-01

Marine current turbine commercial prototypes have now been deployed and arrays of multiple turbines under design. The tidal flows in which they operate are highly turbulent, but the characteristics of the inflow turbulence have not being considered in present design methods. This work considers the effects of inflow turbulence on the wake behind an actuator disc representation of a marine current turbine. Different turbulence intensities and integral length scales were generated in a large eddy simulation using a gridInlet, which produces turbulence from a grid pattern on the inlet boundary. The results highlight the significance of turbulence on the wake profile, with a different flow regime occurring for the zero turbulence case. Increasing the turbulence intensity reduced the velocity deficit and shifted the maximum deficit closer to the turbine. Increasing the integral length scale increased the velocity deficit close to the turbine due to an increased production of turbulent energy. However, the wake recovery was increased due to the higher rate of turbulent mixing causing the wake to expand. The implication of this work is that marine current turbine arrays could be further optimized, increasing the energy yield of the array when the site-specific turbulence characteristics are considered. PMID:25294966

MODELING AND EXPERIMENTAL EVALUATION OF AN AEROSOL GENERATOR FOR VERY HIGH NUMBER CURRENTS BASED ON A FREE TURBULENT JET. (R827354C008)

EPA Science Inventory

In this paper we report on theoretical and experimental work on aerosol formation in a free turbulent jet. A hot DEHS vapor issues through a circular nozzle into slowly moving cold air. Vapor concentration and temperatures are such that particles are formed via homogeneous nuc...

The structure of the solution obtained with Reynolds-stress-transport models at the free-stream edges of turbulent flows

NASA Astrophysics Data System (ADS)

Cazalbou, J.-B.; Chassaing, P.

2002-02-01

The behavior of Reynolds-stress-transport models at the free-stream edges of turbulent flows is investigated. Current turbulent-diffusion models are found to produce propagative (possibly weak) solutions of the same type as those reported earlier by Cazalbou, Spalart, and Bradshaw [Phys. Fluids 6, 1797 (1994)] for two-equation models. As in the latter study, an analysis is presented that provides qualitative information on the flow structure predicted near the edge if a condition on the values of the diffusion constants is satisfied. In this case, the solution appears to be fairly insensitive to the residual free-stream turbulence levels needed with conventional numerical methods. The main specific result is that, depending on the diffusion model, the propagative solution can force turbulence toward definite and rather extreme anisotropy states at the edge (one- or two-component limit). This is not the case with the model of Daly and Harlow [Phys. Fluids 13, 2634 (1970)]; it may be one of the reasons why this "old" scheme is still the most widely used, even in recent Reynolds-stress-transport models. In addition, the analysis helps us to interpret some difficulties encountered in computing even very simple flows with Lumley's pressure-diffusion model [Adv. Appl. Mech. 18, 123 (1978)]. A new realizability condition, according to which the diffusion model should not globally become "anti-diffusive," is introduced, and a recalibration of Lumley's model satisfying this condition is performed using information drawn from the analysis.

Study of the origin and structure of a nocturnal atmospheric density current from observations and numerical simulations

NASA Astrophysics Data System (ADS)

Ander Arrillaga, Jon; Yagüe, Carlos; Román-Cascón, Carlos; Sastre, Mariano

2016-04-01

Density currents are flows generated when a dense fluid passes through a less dense surrounding, under the influence of gravity. They usually appear as a consequence of sea-breeze circulations, thunderstorm outflows or katabatic flows. Density currents acquire a particular relevance during nocturnal stable situations, as their onset causes a significant turbulence increase (both from buoyancy and shear) and they occasionally produce turbulence intermittency through the formation of gravity waves. In this work, the arrival of a density current on 23 September 2015 is analysed in the CIBA site (Spain), which is located in the Spanish Northern Plateau, approximately 200 km away from the sea and 100 km away from the closest mountain ranges. Previous studies at this location associated similar nocturnal events with daytime sea breeze in the eastern Cantabrian coast [1]. Micrometeorological measurements from sonic anemometers and different sensors at multiple levels up to 100 m agl provide a solid database. In this specific case, the outbreak of the density current occurs 2 hours after sunset, causing an abrupt increase of the wind speed and a significant weakening of the surface-based thermal inversion. Besides, turbulent parameters and fluxes such as the friction velocity, the sensible heat flux and the Turbulent Kinetic Energy (TKE) are sharply altered with its arrival. The latter, indeed, increases by two orders of magnitude and the Multi Resolution Flux Decomposition (MRFD) of this and other turbulent variables gives the approximate size of the contributing eddies. Furthermore, simulations with the WRF model, which is tested for different Planetary Boundary Layer (PBL) schemes and the topo_wind option for complex topography [2], give meaningful information about the vertical structure and origin of this density current. [1] Udina, M., Soler, M.R., Viana, S. & Yagüe, C. (2013). Model simulation of gravity waves triggered by a density current. Q J R Meteorol Soc, 139, 701-714. [2] Jiménez, P.A., Dudhia, J., González-Rouco, J.F., Navarro, J., Montávez, J.P., García-Bustamante, E. (2012). A revised scheme for the wrf surface layer formulation. Mon Weather Rev, 140, 898-918.

An Error-Reduction Algorithm to Improve Lidar Turbulence Estimates for Wind Energy

DOE Office of Scientific and Technical Information (OSTI.GOV)

Newman, Jennifer F.; Clifton, Andrew

2016-08-01

Currently, cup anemometers on meteorological (met) towers are used to measure wind speeds and turbulence intensity to make decisions about wind turbine class and site suitability. However, as modern turbine hub heights increase and wind energy expands to complex and remote sites, it becomes more difficult and costly to install met towers at potential sites. As a result, remote sensing devices (e.g., lidars) are now commonly used by wind farm managers and researchers to estimate the flow field at heights spanned by a turbine. While lidars can accurately estimate mean wind speeds and wind directions, there is still a largemore » amount of uncertainty surrounding the measurement of turbulence with lidars. This uncertainty in lidar turbulence measurements is one of the key roadblocks that must be overcome in order to replace met towers with lidars for wind energy applications. In this talk, a model for reducing errors in lidar turbulence estimates is presented. Techniques for reducing errors from instrument noise, volume averaging, and variance contamination are combined in the model to produce a corrected value of the turbulence intensity (TI), a commonly used parameter in wind energy. In the next step of the model, machine learning techniques are used to further decrease the error in lidar TI estimates.« less

Sedimentological regimes for turbidity currents: Depth-averaged theory

NASA Astrophysics Data System (ADS)

Halsey, Thomas C.; Kumar, Amit; Perillo, Mauricio M.

2017-07-01

Turbidity currents are one of the most significant means by which sediment is moved from the continents into the deep ocean; their properties are interesting both as elements of the global sediment cycle and due to their role in contributing to the formation of deep water oil and gas reservoirs. One of the simplest models of the dynamics of turbidity current flow was introduced three decades ago, and is based on depth-averaging of the fluid mechanical equations governing the turbulent gravity-driven flow of relatively dilute turbidity currents. We examine the sedimentological regimes of a simplified version of this model, focusing on the role of the Richardson number Ri [dimensionless inertia] and Rouse number Ro [dimensionless sedimentation velocity] in determining whether a current is net depositional or net erosional. We find that for large Rouse numbers, the currents are strongly net depositional due to the disappearance of local equilibria between erosion and deposition. At lower Rouse numbers, the Richardson number also plays a role in determining the degree of erosion versus deposition. The currents become more erosive at lower values of the product Ro × Ri, due to the effect of clear water entrainment. At higher values of this product, the turbulence becomes insufficient to maintain the sediment in suspension, as first pointed out by Knapp and Bagnold. We speculate on the potential for two-layer solutions in this insufficiently turbulent regime, which would comprise substantial bedload flow with an overlying turbidity current.

Laboratory evidence of dynamo amplification of magnetic fields in a turbulent plasma

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tzeferacos, P.; Rigby, A.; Bott, A. F. A.

Magnetic fields are ubiquitous in the Universe. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. However, experimental demonstration of the turbulent dynamo mechanism has remained elusive, since it requires plasma conditions that are extremely hard to re-create in terrestrial laboratories. Heremore » in this paper, we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization.« less

Laboratory evidence of dynamo amplification of magnetic fields in a turbulent plasma

DOE PAGES

Tzeferacos, P.; Rigby, A.; Bott, A. F. A.; ...

2018-02-09

Magnetic fields are ubiquitous in the Universe. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. However, experimental demonstration of the turbulent dynamo mechanism has remained elusive, since it requires plasma conditions that are extremely hard to re-create in terrestrial laboratories. Heremore » in this paper, we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization.« less

Laboratory evidence of dynamo amplification of magnetic fields in a turbulent plasma.

PubMed

Tzeferacos, P; Rigby, A; Bott, A F A; Bell, A R; Bingham, R; Casner, A; Cattaneo, F; Churazov, E M; Emig, J; Fiuza, F; Forest, C B; Foster, J; Graziani, C; Katz, J; Koenig, M; Li, C-K; Meinecke, J; Petrasso, R; Park, H-S; Remington, B A; Ross, J S; Ryu, D; Ryutov, D; White, T G; Reville, B; Miniati, F; Schekochihin, A A; Lamb, D Q; Froula, D H; Gregori, G

2018-02-09

Magnetic fields are ubiquitous in the Universe. The energy density of these fields is typically comparable to the energy density of the fluid motions of the plasma in which they are embedded, making magnetic fields essential players in the dynamics of the luminous matter. The standard theoretical model for the origin of these strong magnetic fields is through the amplification of tiny seed fields via turbulent dynamo to the level consistent with current observations. However, experimental demonstration of the turbulent dynamo mechanism has remained elusive, since it requires plasma conditions that are extremely hard to re-create in terrestrial laboratories. Here we demonstrate, using laser-produced colliding plasma flows, that turbulence is indeed capable of rapidly amplifying seed fields to near equipartition with the turbulent fluid motions. These results support the notion that turbulent dynamo is a viable mechanism responsible for the observed present-day magnetization.

Numerical Modeling of Turbulence Effects within an Evaporating Droplet in Atomizing Sprays

NASA Technical Reports Server (NTRS)

Balasubramanyam, M. S.; Chen, C. P.; Trinh, H. P.

2006-01-01

A new approach to account for finite thermal conductivity and turbulence effects within atomizing liquid sprays is presented in this paper. The model is an extension of the T-blob and T-TAB atomization/spray model of Trinh and Chen (2005). This finite conductivity model is based on the two-temperature film theory, where the turbulence characteristics of the droplet are used to estimate the effective thermal diffhsivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. The current evaporation model is incorporated into the T-blob atomization model of Trinh and Chen (2005) and implemented in an existing CFD Eulerian-Lagrangian two-way coupling numerical scheme. Validation studies were carried out by comparing with available evaporating atomization spray experimental data in terms of jet penetration, temperature field, and droplet SMD distribution within the spray. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating spray.

Large eddy simulation of dust-uplift by haboob density currents

NASA Astrophysics Data System (ADS)

Huang, Q.

2017-12-01

Cold pool outflows have been shown from both observations and convection-permitting models to be a dominant source of dust uplift ("haboobs") in the summertime Sahel and Sahara, and to cause dust uplift over deserts across the world. In this paper large eddy model (LEM) simulations, which resolve the turbulence within the cold-pools much better than previous studies of haboobs which have used convection-permitting models, are used to investigate the winds that cause dust uplift in cold pools, and the resultant dust uplift and transport. Dust uplift largely occurs in the head of the density current, consistent with the few existing observations. In the modeled density current dust is largely restricted to the lowest coldest and well mixed layer of the cold pool outflow (below around 400 m), except above the head of the cold pool where some dust reaches 2.5 km. This rapid transport to high altitude will contribute to long atmospheric lifetimes of large dust particles from haboobs. Decreasing the model horizontal grid-spacing from 1.0 km to 100 m resolves more turbulence, locally increasing winds, increasing mixing and reducing the propagation speed of the density current. Total accumulated dust uplift is approximately twice as large in 1.0 km runs compared with 100 m runs, suggesting that for studying haboobs in convection-permitting runs the representation of turbulence and mixing is significant. Simulations with surface sensible heat fluxes representative of those from a desert region in daytime show that increasing surface fluxes slow the density current due to increased mixing, but increase dust uplift rates, due to increased downward transport of momentum to the surface.

Laser Doppler systems in atmospheric turbulence

NASA Technical Reports Server (NTRS)

Murty, S. S. R.

1976-01-01

The loss of heterodyne signal power for the Marshall Space Flight Center laser Doppler system due to the random changes in the atmospheric index of refraction is investigated. The current status in the physics of low energy laser propagation through turbulent atmosphere is presented. The analysis and approximate evaluation of the loss of the heterodyne signal power due to the atmospheric absorption, scattering, and turbulence are estimated for the conditions of the January 1973 flight tests. Theoretical and experimental signal to noise values are compared. Maximum and minimum values of the atmospheric attenuation over a two way path of 20 km range are calculated as a function of altitude using models of atmosphere, aerosol concentration, and turbulence.

Comprehensive Validation of an Intermittency Transport Model for Transitional Low-Pressure Turbine Flows

NASA Technical Reports Server (NTRS)

Suzen, Y. B.; Huang, P. G.

2005-01-01

A transport equation for the intermittency factor is employed to predict transitional flows under the effects of pressure gradients, freestream turbulence intensities, Reynolds number variations, flow separation and reattachment. and unsteady wake-blade interactions representing diverse operating conditions encountered in low-pressure turbines. The intermittent behaviour of the transitional flows is taken into account and incorporated into computations by modifying the eddy viscosity, Mu(sub t), with the intermittency factor, gamma. Turbulent quantities are predicted by using Menter's two-equation turbulence model (SST). The onset location of transition is obtained from correlations based on boundary-layer momentum thickness, acceleration parameter, and turbulence intensity. The intermittency factor is obtained from a transport model which can produce both the experimentally observed streamwise variation of intermittency and a realistic profile in the cross stream direction. The intermittency transport model is tested and validated against several well documented low pressure turbine experiments ranging from flat plate cases to unsteady wake-blade interaction experiments. Overall, good agreement between the experimental data and computational results is obtained illustrating the predicting capabilities of the model and the current intermittency transport modelling approach for transitional flow simulations.

Numerical model of turbulence, sediment transport, and morphodynamics tested in the Colorado River at Grand Canyon

NASA Astrophysics Data System (ADS)

Alvarez, L. V.; Grams, P.

2017-12-01

We present a parallelized, three-dimensional, turbulence-resolving model using the Detached-Eddy Simulation (DES) technique, tested at the scale of the river-reach in the Colorado River. DES is a hybrid large eddy simulation (LES) and Reynolds-averaged Navier Stokes (RANS). RANS is applied to the near-bed grid cells, where grid resolution is not sufficient to fully resolve wall turbulence. LES is applied in the flow interior. We utilize the Spalart-Allmaras one equation turbulence closure with a rough wall extension. The model resolves large-scale turbulence using DES and simultaneously integrates the suspended sediment advection-diffusion equation. The Smith and McLean suspended sediment boundary condition is used to calculate the upward and downward settling of sediment fluxes in the grid cells attached to the bed. Model results compare favorably with ADCP measurements of flow taken on the Colorado River in Grand Canyon during the High Flow Experiment (HFE) of 2008. The model accurately reproduces the size and position of the major recirculation currents, and the error in velocity magnitude was found to be less than 17% or 0.22 m/s absolute error. The mean deviation of the direction of velocity with respect to the measured velocity was found to be 20 degrees. Large-scale turbulence structures with vorticity predominantly in the vertical direction are produced at the shear layer between the main channel and the separation zone. However, these structures rapidly become three-dimensional with no preferred orientation of vorticity. Cross-stream velocities, into the main recirculation zone just upstream of the point of reattachment and out of the main recirculation region just downstream of the point of separation, are highest near the bed. Lateral separation eddies are more efficient at storing and exporting sediment than previously modeled. The input of sediment to the eddy recirculation zone occurs in the interface of the eddy and main channel. Pulsation of the strength of the return current becomes a key factor to determine the rates of erosion and deposition in the main recirculation zone.

Turbulence evolution and transport behavior during current ramp-up in ITER-like plasmas on DIII-D

DOE PAGES

McKee, George R.; Austin, Max E.; Boedo, Jose A.; ...

2017-07-12

Low-wavenumber density fluctuations exhibit unique characteristics during the current ramp-up phase of ITER-like discharges that can partially explain the challenges of correctly modeling transport behavior and predicting global plasma parameters during this period. A strong interaction takes place between the evolving transport, safety factor (q) and kinetic profiles as well as the appearance and evolution of low-order rational surfaces. Density fluctuations from 0.75 < ρ < 0.9 are transiently reduced to exceptionally low levels during early times and from 0.8 < ρ < 0.9 at late times in the ramp-up in a manner that is different from behavior observed duringmore » steady-state plasma conditions with similar values of q 95. Turbulence is suppressed as low-order-rational q-surfaces enter the plasma; the local electron temperature likewise exhibits transient increases during these periods of reduced fluctuations indicating changes in transport that impact temperature and consequently the evolution of current density and plasma inductance. These observations can explain discrepancies between CORSICA modelling and the higher electron temperature found previously over the outer half radius. Comparison of turbulence properties with time-varying linear growth rates with GYRO and GENE demonstrate qualitative consistency with measured fluctuation levels, but calculations don’t exhibit reduced growth rates near low-order rational surfaces, which is inconsistent with experimental observations. Here, this indicates a mechanism that can contribute to reconciling observed turbulence behavior with transport models, allowing for the development of more accurate predictive tools.« less

Turbulence evolution and transport behavior during current ramp-up in ITER-like plasmas on DIII-D

NASA Astrophysics Data System (ADS)

McKee, G. R.; Austin, M.; Boedo, J.; Bravenec, R.; Holland, C.; Jackson, G.; Luce, T. C.; Rhodes, T. L.; Rudakov, D.; Wang, G.; Yan, Z.; Zeng, L.; Zhao, Y.

2017-08-01

Low-wavenumber density fluctuations exhibit unique characteristics during the current ramp-up phase of ITER-like discharges that can partially explain the challenges of correctly modeling transport behavior and predicting global plasma parameters during this period. A strong interaction takes place between the evolving transport, safety factor (q) and kinetic profiles as well as the appearance and evolution of low-order rational surfaces. Density fluctuations from 0.75  <  ρ  <  0.9 are transiently reduced to exceptionally low levels during early times and from 0.8  <  ρ  <  0.9 at late times in the ramp-up in a manner that is different from behavior observed during steady-state plasma conditions with similar values of q 95. Turbulence is suppressed as low-order-rational q-surfaces enter the plasma; the local electron temperature likewise exhibits transient increases during these periods of reduced fluctuations indicating changes in transport that impact temperature and consequently the evolution of current density and plasma inductance. These observations can explain discrepancies between CORSICA modelling and the higher electron temperature found previously over the outer half radius. Comparison of turbulence properties with time-varying linear growth rates with GYRO and GENE demonstrate qualitative consistency with measured fluctuation levels, but calculations don’t exhibit reduced growth rates near low-order rational surfaces, which is inconsistent with experimental observations. This indicates a mechanism that can contribute to reconciling observed turbulence behavior with transport models, allowing for the development of more accurate predictive tools.

On the interaction between fluid turbulence and particle loading: numerical simulation of turbidity currents and prediction of deep-sea arenites

NASA Astrophysics Data System (ADS)

Doronzo, D. M.; Dufek, J.

2012-04-01

Turbidity currents are water-particle flows able to move large distance over the seafloor, and the deep-sea arenitic facies of their deposits often represents an important class of hydrocarbon reservoirs. Coupling flow behavior and the resulting deposits may thus help finding new reservoirs, as well as reconstructing the sediment transport mechanisms from the continental shelf to the abyssal plain. There is a broad literature of turbidity currents, which includes field, theoretical, experimental, and numerical studies on flow dynamics and associated deposits. Generally, the field and theoretical approaches focus on the scale of actual deposits and currents, respectively, whereas experimental and numerical approaches are often restricted to the laboratory scale and relatively low-Reynolds number, respectively. Fully resolved simulations that incorporate complex bathymetry, large-scale flow, multiphase and 3D effects, are computationally expensive and require closure schemes. Here, a 2D numerical model of turbidity current is proposed, which is based on the Euler-Lagrange formulation of multiphase physics, and on the Reynolds-averaged Navier-Stokes closure of turbulence. This strategy has been recently used in volcanology to simulate the gas-particle flow of pyroclastic density currents, in order to predict their deposits. The incompressible conservation equations of mass and momentum are solved for the water, and the equation of particle motion is solved for the sediment, which for this example, has an initial concentration of 1 % of 0.5 mm sand particles. The equations are solved numerically with the finite-volume method of Ansys Fluent software, and particle and fluid motion are two-way coupled during calculation, which means that the particles are tracked on the basis of water solution, then are allowed to affect the liquid turbulence through a momentum exchange. The Reynolds (turbulent) stresses, which dominate over the viscous ones in the turbidity current, are calculated with a two-equation model (RNG k-ɛ) solving for the turbulent kinetic energy and the turbulent dissipation rate. The simulated seafloor is represented by a ramp 8 km long and 3° steep, over which the particles rebound inelastically, in order to capture the bed-load of the current. Although the sediment is mainly transported as suspended-load (this makes the flow "turbid"), the ground-hugging processes play a fundamental role in the emplacement of deposits, as well as in the flow behavior. A highly refined grid of 0.2 m at the base is thus used to solve for these processes. After 6750 s of flow time, sedimentation rates of 4 and 0.5 kg/m2 s are calculated over the seafloor in proximal (1 km) and medial (4 km) regions, respectively. These values are converted to deposit thickness, resulting in arenitic turbidite sequences of 14.5 and 1.8 m, respectively. Turbulence intensities of 54 and 66 %, respectively in the same areas, indicate the water is further made turbulent by the sediment (water-sand interaction), so the sand moves in suspension toward the deep-sea, where is able to deposit. Richardson numbers of 0.79 and 0.58, respectively, show how the water entrainment in the current increases with distance.

Final Scientific Report: Experimental Investigation of Reconnection in a Line-tied Plasma

DOE Office of Scientific and Technical Information (OSTI.GOV)

Forest, Cary

This grant used funding from the NSF/DoE Partnership on Plasma Science to investigate magnetic reconnection phenomena in a line-tied pinch experiment. The experiment was upgraded from a previous device intended to study fusion plasma-related instabilities to a new configuration capable of studying a number of new, previously unstudied configurations. A high spatial and time resolution array of magnetic probes was constructed to measure time evolving structures present as instability and turbulence developed. The most important new equilibrium made possible by this grant was a Zero-Net-Current equilibrium that models the footpoint twisting of solar flux tubes that occurs prior to solarmore » eruptions (flares and coronal mass ejections). This new equilibrium was successfully created in the lab, and it exhibited a host of instabilities. In particular, at low current when the equilibrium was not overly stressed, a saturated internal kink mode oscillation was observed. At high current, 2 D magnetic turbulence developed which we attribute to the lack of a equilibrium brought about by a subcritical transition to turbulence. A second set of experiments involved the turbulent interactions of a collection of flux tubes all being twisted independently, a problem known as the Parker Problem. Current profiles consisting of 2, 3 and 4 guns were used to impose a fine scale drive, and resulted in a new experimental platform in which the injection scale of the magnetic turbulence could be controlled. First experiments in this configuration support the conclusion that an inverse cascade of magnetic energy occurred which self-organized the plasma into a nearly axisymmetric current distribution.« less

Development of a Localized Low-Dimensional Approach to Turbulence Simulation

NASA Astrophysics Data System (ADS)

Juttijudata, Vejapong; Rempfer, Dietmar; Lumley, John

2000-11-01

Our previous study has shown that the localized low-dimensional model derived from a projection of Navier-Stokes equations onto a set of one-dimensional scalar POD modes, with boundary conditions at y^+=40, can predict wall turbulence accurately for short times while failing to give a stable long-term solution. The structures obtained from the model and later studies suggest our boundary conditions from DNS are not consistent with the solution from the localized model resulting in an injection of energy at the top boundary. In the current study, we develop low-dimensional models using one-dimensional scalar POD modes derived from an explicitly filtered DNS. This model problem has exact no-slip boundary conditions at both walls while the locality of the wall layer is still retained. Furthermore, the interaction between wall and core region is attenuated via an explicit filter which allows us to investigate the quality of the model without requiring complicated modeling of the top boundary conditions. The full-channel model gives reasonable wall turbulence structures as well as long-term turbulent statistics while still having difficulty with the prediction of the mean velocity profile farther from the wall. We also consider a localized model with modified boundary conditions in the last part of our study.

The structure of particle-laden jets and nonevaporating sprays

NASA Technical Reports Server (NTRS)

Shuen, J. S.; Solomon, A. S. P.; Zhang, Q. F.; Faeth, G. M.

1983-01-01

Mean and fluctuating gas velocities, liquid mass fluxes and drop sizes were in nonevaporating sprays. These results, as well as existing measurements in solid particle-laden jets, were used to evaluate models of these processes. The following models were considered: (1) a locally homogeneous flow (LHF) model, where slip between the phases was neglected; (2) a deterministic separated flow (DSF) model, where slip was considered but effects of particle dispersion by turbulence were ignored; and (3) a stochastic separated flow (SSF) model, where effects of interphase slip and turbulent dispersion were considered using random-walk computations for particle motion. The LHF and DSF models did not provide very satisfactory predictions over the present data base. In contrast, the SSF model performed reasonably well - including conditions in nonevaporating sprays where enhanced dispersion of particles by turbulence caused the spray to spread more rapidly than single-phase jets for comparable conditions. While these results are encouraging, uncertainties in initial conditions limit the reliability of the evaluation. Current work is seeking to eliminate this deficiency.

Validation of the kinetic-turbulent-neoclassical theory for edge intrinsic rotation in DIII-D

NASA Astrophysics Data System (ADS)

Ashourvan, Arash; Grierson, B. A.; Battaglia, D. J.; Haskey, S. R.; Stoltzfus-Dueck, T.

2018-05-01

In a recent kinetic model of edge main-ion (deuterium) toroidal velocity, intrinsic rotation results from neoclassical orbits in an inhomogeneous turbulent field [T. Stoltzfus-Dueck, Phys. Rev. Lett. 108, 065002 (2012)]. This model predicts a value for the toroidal velocity that is co-current for a typical inboard X-point plasma at the core-edge boundary (ρ ˜ 0.9). Using this model, the velocity prediction is tested on the DIII-D tokamak for a database of L-mode and H-mode plasmas with nominally low neutral beam torque, including both signs of plasma current. Values for the flux-surface-averaged main-ion rotation velocity in the database are obtained from the impurity carbon rotation by analytically calculating the main-ion—impurity neoclassical offset. The deuterium rotation obtained in this manner has been validated by direct main-ion measurements for a limited number of cases. Key theoretical parameters of ion temperature and turbulent scale length are varied across a wide range in an experimental database of discharges. Using a characteristic electron temperature scale length as a proxy for a turbulent scale length, the predicted main-ion rotation velocity has a general agreement with the experimental measurements for neutral beam injection (NBI) powers in the range PNBI < 4 MW. At higher NBI power, the experimental rotation is observed to saturate and even degrade compared to theory. TRANSP-NUBEAM simulations performed for the database show that for discharges with nominally balanced—but high powered—NBI, the net injected torque through the edge can exceed 1 Nm in the counter-current direction. The theory model has been extended to compute the rotation degradation from this counter-current NBI torque by solving a reduced momentum evolution equation for the edge and found the revised velocity prediction to be in agreement with experiment. Using the theory modeled—and now tested—velocity to predict the bulk plasma rotation opens up a path to more confidently projecting the confinement and stability in ITER.

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

NASA Technical Reports Server (NTRS)

Lillard, Randolph P.

2011-01-01

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Sub-grid-scale description of turbulent magnetic reconnection in magnetohydrodynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Widmer, F., E-mail: widmer@mps.mpg.de; Institut für Astrophysik, Georg-August-Universität, Friedrich-Hund-Platz 1, 37077 Göttingen; Büchner, J.

Magnetic reconnection requires, at least locally, a non-ideal plasma response. In collisionless space and astrophysical plasmas, turbulence could transport energy from large to small scales where binary particle collisions are rare. We have investigated the influence of small scale magnetohydrodynamics (MHD) turbulence on the reconnection rate in the framework of a compressible MHD approach including sub-grid-scale (SGS) turbulence. For this sake, we considered Harris-type and force-free current sheets with finite guide magnetic fields directed out of the reconnection plane. The goal is to find out whether unresolved by conventional simulations MHD turbulence can enhance the reconnection process in high-Reynolds-number astrophysicalmore » plasmas. Together with the MHD equations, we solve evolution equations for the SGS energy and cross-helicity due to turbulence according to a Reynolds-averaged turbulence model. The SGS turbulence is self-generated and -sustained through the inhomogeneities of the mean fields. By this way, the feedback of the unresolved turbulence into the MHD reconnection process is taken into account. It is shown that the turbulence controls the regimes of reconnection by its characteristic timescale τ{sub t}. The dependence on resistivity was investigated for large-Reynolds-number plasmas for Harris-type as well as force-free current sheets with guide field. We found that magnetic reconnection depends on the relation between the molecular and apparent effective turbulent resistivity. We found that the turbulence timescale τ{sub t} decides whether fast reconnection takes place or whether the stored energy is just diffused away to small scale turbulence. If the amount of energy transferred from large to small scales is enhanced, fast reconnection can take place. Energy spectra allowed us to characterize the different regimes of reconnection. It was found that reconnection is even faster for larger Reynolds numbers controlled by the molecular resistivity η, as long as the initial level of turbulence is not too large. This implies that turbulence plays an important role to reach the limit of fast reconnection in large Reynolds number plasmas even for smaller amounts of turbulence.« less

Development of wavelet analysis tools for turbulence

NASA Technical Reports Server (NTRS)

Bertelrud, A.; Erlebacher, G.; Dussouillez, PH.; Liandrat, M. P.; Liandrat, J.; Bailly, F. Moret; Tchamitchian, PH.

1992-01-01

Presented here is the general framework and the initial results of a joint effort to derive novel research tools and easy to use software to analyze and model turbulence and transition. Given here is a brief review of the issues, a summary of some basic properties of wavelets, and preliminary results. Technical aspects of the implementation, the physical conclusions reached at this time, and current developments are discussed.

Modeling of turbulence effects on the heat and mass transfer of evaporating sprays

NASA Astrophysics Data System (ADS)

Madhanabharatam, Balasubramanyam

A large diversity of two-phase gas-liquid flows of both scientific and practical interest involves the evaporation of near spherical liquid droplets in high temperature turbulent environments. Current numerical modeling approaches are predominantly focused towards the effects of continuous phase (gas phase) turbulence on the evaporation rates of liquid fuel sprays during the evaporation process, failing to account for the inherent turbulence present in the dispersed phase (liquid phase), due to the injection of sprays at high velocities. Existing models accounting for internal turbulence effects use Direct Numerical Simulations and Large Eddy Simulations that are computationally intensive. This research provides an alternative phenomenological approach of modeling droplet internal turbulence effects through the mass and heat transfer between the droplet surface and the external gas phase within a thin film inside the droplet. This finite conductivity (F-C) model was based on the two-temperature film theory, where the turbulence characteristics of the droplet are used to estimate the effective thermal diffusivity (alphaeff) within the droplet phase. The alphaeff is estimated from the physical properties of the flow within the droplet rather than from a 'curve-fit' as done conventionally. The results of the one-way coupled study indicated that the equilibrium drop temperature predictions were higher than calculations by the infinite conductivity (I-C) model. The liquid internal turbulence has a considerable effect on the diffusivity in the primary atomization regime. The thermal boundary layer was found to be substantially thick initially, decreasing quickly to a small value, exhibiting a reasonable physical trend. The two-way coupled studies (CFD) indicated that the F-C model, slowed down the evaporation process, produced larger droplets and longer tip penetration lengths during the initial stages of injection. For a jet in a supersonic cross-flow, results indicated that jet penetration increased rapidly in the vicinity of the injector exit and then gradually increased due to increase in the drag of the air stream. A modified drag coefficient was incorporated to improve model accuracy in predictions. Overall the results obtained from the numerical calculations during this study were reasonably comparable to measured data and showed more accurate comparisons to that of the I-C model.

Direct numerical simulations and spectral modeling of premixed turbulent combustion in the flamelet regime

NASA Astrophysics Data System (ADS)

Ulitsky, Mark

1997-11-01

A model for premixed turbulent combustion in the so called 'flamelet regime' has been developed. This regime, often referred to as the fast chemistry or high Damkohler number regime, is characterized by turbulent length and time scales that are much larger and slower than the flame thickness and reaction time scales respectively. There is currently great interest in trying to better understand flamelet combustion, as many practical devices (i.e., spark ignition engines, gas turbines, etc.) have been found to operate in this regime. Before a model could be developed however, it was first necessary to ascertain which part of the turbulence (either the nearly Gaussian background turbulence or the tube-like coherent vortical structures) was responsible for the multi-scale wrinkling of the flame surface. This question motivated a DNS study of flames passing through both structure containing the structure free isotropic turbulence. After it was determined that the presence of the coherent structures was merely ancillary in terms of increasing the surface area of the flame, a spectral model based on the EDQNM (Eddy Damped Quasi Normal Markovian) theory of turbulence was developed. This theory implicitly assumes that joint distributions of the fluctuating velocity components are nearly Gaussian, and as only spectra are transported in this model, there is no direct information about any of the coherent structures which might be embedded within the flow field. One of the advantages of this model is that both the Reynolds number and the ratio of the rms fluctuating velocity to the laminar flame speed can be varied independently. To test the model's ability to capture the nonlinear dynamics of the governing field equation a DNS study was performed and both steady-state and transient single- and two-point statistics were compared. Finally, the model was compared to two-point experimental measurements taken from a lean premixed methane-air flame.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Newman, Jennifer F.; Clifton, Andrew

Currently, cup anemometers on meteorological towers are used to measure wind speeds and turbulence intensity to make decisions about wind turbine class and site suitability; however, as modern turbine hub heights increase and wind energy expands to complex and remote sites, it becomes more difficult and costly to install meteorological towers at potential sites. As a result, remote-sensing devices (e.g., lidars) are now commonly used by wind farm managers and researchers to estimate the flow field at heights spanned by a turbine. Although lidars can accurately estimate mean wind speeds and wind directions, there is still a large amount ofmore » uncertainty surrounding the measurement of turbulence using these devices. Errors in lidar turbulence estimates are caused by a variety of factors, including instrument noise, volume averaging, and variance contamination, in which the magnitude of these factors is highly dependent on measurement height and atmospheric stability. As turbulence has a large impact on wind power production, errors in turbulence measurements will translate into errors in wind power prediction. The impact of using lidars rather than cup anemometers for wind power prediction must be understood if lidars are to be considered a viable alternative to cup anemometers.In this poster, the sensitivity of power prediction error to typical lidar turbulence measurement errors is assessed. Turbulence estimates from a vertically profiling WINDCUBE v2 lidar are compared to high-resolution sonic anemometer measurements at field sites in Oklahoma and Colorado to determine the degree of lidar turbulence error that can be expected under different atmospheric conditions. These errors are then incorporated into a power prediction model to estimate the sensitivity of power prediction error to turbulence measurement error. Power prediction models, including the standard binning method and a random forest method, were developed using data from the aeroelastic simulator FAST for a 1.5 MW turbine. The impact of lidar turbulence error on the predicted power from these different models is examined to determine the degree of turbulence measurement accuracy needed for accurate power prediction.« less

Current and turbulence measurements at the FINO1 offshore wind energy site: analysis using 5-beam ADCPs

NASA Astrophysics Data System (ADS)

Bakhoday-Paskyabi, Mostafa; Fer, Ilker; Reuder, Joachim

2018-01-01

We report concurrent measurements of ocean currents and turbulence at two sites in the North Sea, one site at upwind of the FINO1 platform and the other 200-m downwind of the Alpha Ventus wind farm. At each site, mean currents, Reynolds stresses, turbulence intensity and production of turbulent kinetic energy are obtained from two bottom-mounted 5-beam Nortek Signature1000s, high-frequency Doppler current profiler, at a water depth of approximately 30 m. Measurements from the two sites are compared to statistically identify the effects of wind farm and waves on ocean current variability and the turbulent structure in the water column. Profiles of Reynolds stresses are found to be sensible to both environmental forcing and the wind farm wake-induced distortions in both boundary layers near the surface and the seabed. Production of turbulent kinetic energy and turbulence intensity exhibit approximately similar, but less pronounced, patterns in the presence of farm wake effects.

Turbulence measurements in axisymmetric jets of air and helium. I - Air jet. II - Helium jet

NASA Technical Reports Server (NTRS)

Panchapakesan, N. R.; Lumley, J. L.

1993-01-01

Results are presented of measurements on turbulent round jets of air and of helium of the same nozzle momentum efflux, using, for the air jets, x-wire hot-wire probes mounted on a moving shuttle and, for He jets, a composite probe consisting of an interference probe of the Way-Libby type and an x-probe. Current models for scalar triple moments were evaluated. It was found that the performance of the model termed the Full model, which includes all terms except advection, was very good for both the air and the He jets.

The self-consistent parallel electric field due to electrostatic ion-cyclotron turbulence in downward auroral-current regions of the Earth's magnetosphere. IV

NASA Astrophysics Data System (ADS)

Jasperse, John R.; Basu, Bamandas; Lund, Eric J.; Grossbard, Neil

2010-06-01

The physical processes that determine the self-consistent electric field (E∥) parallel to the magnetic field have been an unresolved problem in magnetospheric physics for over 40 years. Recently, a new multimoment fluid theory was developed for inhomogeneous, nonuniformly magnetized plasma in the guiding-center and gyrotropic approximation that includes the effect of electrostatic, turbulent, wave-particle interactions (see Jasperse et al. [Phys. Plasmas 13, 072903 (2006); Jasperse et al., Phys. Plasmas13, 112902 (2006)]). In the present paper and its companion paper [Jasperse et al., Phys. Plasmas 17, 062903 (2010)], which are intended as sequels to the earlier work, a fundamental model for downward, magnetic field-aligned (Birkeland) currents for quasisteady conditions is presented. The model includes the production of electrostatic ion-cyclotron turbulence in the long-range potential region by an electron, bump-on-tail-driven ion-cyclotron instability. Anomalous momentum transfer (anomalous resistivity) by itself is found to produce a very small contribution to E∥; however, the presence of electrostatic, ion-cyclotron turbulence has a very large effect on the altitude dependence of the entire quasisteady solution. Anomalous energy transfer (anomalous heating and cooling) modifies the density, drift, and temperature altitude profiles and hence the generalized parallel-pressure gradients and mirror forces in the electron and ion momentum-balance equations. As a result, |E∥| is enhanced by nearly a factor of 40 compared to its value when turbulence is absent. The space-averaged potential increase associated with the strong double layer at the bottom of the downward-current sheet is estimated using the FAST satellite data and the multimoment fluid theory.

Dynamic properties of ionospheric plasma turbulence driven by high-power high-frequency radiowaves

NASA Astrophysics Data System (ADS)

Grach, S. M.; Sergeev, E. N.; Mishin, E. V.; Shindin, A. V.

2016-11-01

A review is given of the current state-of-the-art of experimental studies and the theoretical understanding of nonlinear phenomena that occur in the ionospheric F-layer irradiated by high-power high-frequency ground-based transmitters. The main focus is on the dynamic features of high-frequency turbulence (plasma waves) and low-frequency turbulence (density irregularities of various scales) that have been studied in experiments at the Sura and HAARP heating facilities operated in temporal and frequency regimes specially designed with consideration of the characteristic properties of nonlinear processes in the perturbed ionosphere using modern radio receivers and optical instruments. Experimental results are compared with theoretical turbulence models for a magnetized collisional plasma in a high-frequency electromagnetic field, allowing the identification of the processes responsible for the observed features of artificial ionospheric turbulence.

On the Origin of the Saturnian Satellite System: Did Iapetus Form In-Situ?

NASA Technical Reports Server (NTRS)

Mosqueira, I.; Estrada, P. R.

2005-01-01

Current models of planet and satellite formation are marred by our lack of understanding regarding the turbulent state of accretion disks. According to the Rayleigh criterion, Keplerian disks are hydrodynamically stable. Indeed, it has been argued that a carefully designed Taylor Couette experiment shows stability in the case of positive radial gradients in specific angular momentum even for high Reynolds numbers [1], in agreement with numerical simulations which consistently show turbulence decay [2]. Other possible sources of turbulence may fail due to low ionization, may decay as the optical depth decreases due to dust coagulation, may involve unrealistic boundary conditions, or result in limited transport. The difficulty stems not only from the degree of turbulence, but also from the kind of turbulence, and whether it may be characterized by an alpha parameter.

Stochastic modelling of turbulent combustion for design optimization of gas turbine combustors

NASA Astrophysics Data System (ADS)

Mehanna Ismail, Mohammed Ali

The present work covers the development and the implementation of an efficient algorithm for the design optimization of gas turbine combustors. The purpose is to explore the possibilities and indicate constructive suggestions for optimization techniques as alternative methods for designing gas turbine combustors. The algorithm is general to the extent that no constraints are imposed on the combustion phenomena or on the combustor configuration. The optimization problem is broken down into two elementary problems: the first is the optimum search algorithm, and the second is the turbulent combustion model used to determine the combustor performance parameters. These performance parameters constitute the objective and physical constraints in the optimization problem formulation. The examination of both turbulent combustion phenomena and the gas turbine design process suggests that the turbulent combustion model represents a crucial part of the optimization algorithm. The basic requirements needed for a turbulent combustion model to be successfully used in a practical optimization algorithm are discussed. In principle, the combustion model should comply with the conflicting requirements of high fidelity, robustness and computational efficiency. To that end, the problem of turbulent combustion is discussed and the current state of the art of turbulent combustion modelling is reviewed. According to this review, turbulent combustion models based on the composition PDF transport equation are found to be good candidates for application in the present context. However, these models are computationally expensive. To overcome this difficulty, two different models based on the composition PDF transport equation were developed: an improved Lagrangian Monte Carlo composition PDF algorithm and the generalized stochastic reactor model. Improvements in the Lagrangian Monte Carlo composition PDF model performance and its computational efficiency were achieved through the implementation of time splitting, variable stochastic fluid particle mass control, and a second order time accurate (predictor-corrector) scheme used for solving the stochastic differential equations governing the particles evolution. The model compared well against experimental data found in the literature for two different configurations: bluff body and swirl stabilized combustors. The generalized stochastic reactor is a newly developed model. This model relies on the generalization of the concept of the classical stochastic reactor theory in the sense that it accounts for both finite micro- and macro-mixing processes. (Abstract shortened by UMI.)

A detached eddy simulation model for the study of lateral separation zones along a large canyon-bound river

USGS Publications Warehouse

Alvarez, Laura V.; Schmeeckle, Mark W.; Grams, Paul E.

2017-01-01

Lateral flow separation occurs in rivers where banks exhibit strong curvature. In canyon-boundrivers, lateral recirculation zones are the principal storage of fine-sediment deposits. A parallelized,three-dimensional, turbulence-resolving model was developed to study the flow structures along lateralseparation zones located in two pools along the Colorado River in Marble Canyon. The model employs thedetached eddy simulation (DES) technique, which resolves turbulence structures larger than the grid spacingin the interior of the flow. The DES-3D model is validated using Acoustic Doppler Current Profiler flowmeasurements taken during the 2008 controlled flood release from Glen Canyon Dam. A point-to-pointvalidation using a number of skill metrics, often employed in hydrological research, is proposed here forfluvial modeling. The validation results show predictive capabilities of the DES model. The model reproducesthe pattern and magnitude of the velocity in the lateral recirculation zone, including the size and position ofthe primary and secondary eddy cells, and return current. The lateral recirculation zone is open, havingcontinuous import of fluid upstream of the point of reattachment and export by the recirculation returncurrent downstream of the point of separation. Differences in magnitude and direction of near-bed andnear-surface velocity vectors are found, resulting in an inward vertical spiral. Interaction between therecirculation return current and the main flow is dynamic, with large temporal changes in flow direction andmagnitude. Turbulence structures with a predominately vertical axis of vorticity are observed in the shearlayer becoming three-dimensional without preferred orientation downstream.

Assessment of RANS and LES Turbulence Modeling for Buoyancy-Aided/Opposed Forced and Mixed Convection

NASA Astrophysics Data System (ADS)

Clifford, Corey; Kimber, Mark

2017-11-01

Over the last 30 years, an industry-wide shift within the nuclear community has led to increased utilization of computational fluid dynamics (CFD) to supplement nuclear reactor safety analyses. One such area that is of particular interest to the nuclear community, specifically to those performing loss-of-flow accident (LOFA) analyses for next-generation very-high temperature reactors (VHTR), is the capacity of current computational models to predict heat transfer across a wide range of buoyancy conditions. In the present investigation, a critical evaluation of Reynolds-averaged Navier-Stokes (RANS) and large-eddy simulation (LES) turbulence modeling techniques is conducted based on CFD validation data collected from the Rotatable Buoyancy Tunnel (RoBuT) at Utah State University. Four different experimental flow conditions are investigated: (1) buoyancy-aided forced convection; (2) buoyancy-opposed forced convection; (3) buoyancy-aided mixed convection; (4) buoyancy-opposed mixed convection. Overall, good agreement is found for both forced convection-dominated scenarios, but an overly-diffusive prediction of the normal Reynolds stress is observed for the RANS-based turbulence models. Low-Reynolds number RANS models perform adequately for mixed convection, while higher-order RANS approaches underestimate the influence of buoyancy on the production of turbulence.

Applying Turbulence Models to Hydroturbine Flows: A Sensitivity Analysis Using the GAMM Francis Turbine

NASA Astrophysics Data System (ADS)

Lewis, Bryan; Cimbala, John; Wouden, Alex

2011-11-01

Turbulence models are generally developed to study common academic geometries, such as flat plates and channels. Creating quality computational grids for such geometries is trivial, and allows stringent requirements to be met for boundary layer grid refinement. However, engineering applications, such as flow through hydroturbines, require the analysis of complex, highly curved geometries. To produce body-fitted grids for such geometries, the mesh quality requirements must be relaxed. Relaxing these requirements, along with the complexity of rotating flows, forces turbulence models to be employed beyond their developed scope. This study explores the solution sensitivity to boundary layer grid quality for various turbulence models and boundary conditions currently implemented in OpenFOAM. The following models are resented: k-omega, k-omega SST, k-epsilon, realizable k-epsilon, and RNG k-epsilon. Standard wall functions, adaptive wall functions, and sub-grid integration are compared using various grid refinements. The chosen geometry is the GAMM Francis Turbine because experimental data and comparison computational results are available for this turbine. This research was supported by a grant from the DoE and a National Defense Science and Engineering Graduate Fellowship.

Experience with turbulence interaction and turbulence-chemistry models at Fluent Inc.

NASA Technical Reports Server (NTRS)

Choudhury, D.; Kim, S. E.; Tselepidakis, D. P.; Missaghi, M.

1995-01-01

This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.

Water velocity, turbulence, and migration rate of subyearling fall Chinook salmon in the free-flowing and impounded Snake River

USGS Publications Warehouse

Tiffan, Kenneth F.; Kock, Tobias J.; Haskell, Craig A.; Connor, William P.; Steinhorst, R. Kirk

2009-01-01

We studied the migratory behavior of subyearling fall Chinook salmon Oncorhynchus tshawytscha in free-flowing and impounded reaches of the Snake River to evaluate the hypothesis that velocity and turbulence are the primary causal mechanisms of downstream migration. The hypothesis states that impoundment reduces velocity and turbulence and alters the migratory behavior of juvenile Chinook salmon as a result of their reduced perception of these cues. At a constant flow (m3 /s), both velocity (km/d) and turbulence (the SD of velocity) decreased from riverine to impounded habitat as cross-sectional areas increased. We found evidence for the hypothesis that subyearling Chinook salmon perceive velocity and turbulence cues and respond to these cues by varying their behavior. The percentage of the subyearlings that moved faster than the average current speed decreased as fish made the transition from riverine reaches with high velocities and turbulence to upper reservoir reaches with low velocities and turbulence but increased to riverine levels again as the fish moved further down in the reservoir, where velocity and turbulence remained low. The migration rate (km/d) decreased in accordance with longitudinal reductions in velocity and turbulence, as predicted by the hypothesis. The variation in migration rate was better explained by a repeatedmeasures regression model containing velocity (Akaike’s information criterion ¼ 1,769.0) than a model containing flow (2,232.6). We conclude that subyearling fall Chinook salmon respond to changes in water velocity and turbulence, which work together to affect the migration rate.

Experimental results on current-driven turbulence in plasmas - a survey

NASA Astrophysics Data System (ADS)

de Kluiver, H.; Perepelkin, N. F.; Hirose, A.

1991-01-01

The experimental consequences of plasma turbulence driven by a current parallel to a magnetic field and concurrent anomalous plasma heating are reviewed, with an attempt to deduce universalities in key parameters such as the anomalous electrical conductivities observed in diverse devices. It has been found that the nature of plasma turbulence and turbulent heating depends on several parameters including the electric field, current and magnetic fields. A classification of turbulence regimes based on these parameters has been made. Experimental observations of the anomalous electrical conductivity, plasma heating, skin effect, runaway electron braking and turbulent fluctuations are surveyed, and current theoretical understanding is briefly reviewed. Experimental results recently obtained in stellarators (SIRIUS, URAGAN at Kharkov), and in tokamaks (TORTUR at Nieuwegein, STOR-1M at Saskatoon) are presented in some detail in the light of investigating the feasibility of using turbulent heating as a means of injecting a large power into toroidal devices.

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

NASA Technical Reports Server (NTRS)

Abdol-Hamid, Khaled S.; Elmiligui, Alaa; Giriamaji, Sharath S.

2008-01-01

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

Pulsating Instability of Turbulent Thermonuclear Flames in Type Ia Supernovae

NASA Astrophysics Data System (ADS)

Poludnenko, Alexei Y.

2014-01-01

Presently, one of the main explosion scenarios of type Ia supernovae (SNIa), aimed at explaining both "normal" and subluminous events, is the thermonuclear incineration of a white-dwarf in a single-degenerate system. The underlying engine of such explosions is the turbulent thermonuclear flame. Modern, large-scale, multidimensional simulations of SNIa cannot resolve the internal flame structure, and instead must include a subgrid-scale prescription for the turbulent-flame properties. As a result, development of robust, parameter-free, large-scale models of SNIa crucially relies on the detailed understanding of the turbulent flame properties during each stage of the flame evolution. Due to the complexity of the flame dynamics, such understanding must be validated by the first-principles direct numerical simulations (DNS). In our previous work, we showed that sufficiently fast turbulent flames are inherently susceptible to the development of detonations, which may provide the mechanism for the deflagration-to-detonation transition (DDT) in the delayed-detonation model of SNIa. Here we extend this study by performing detailed analysis of the turbulent flame properties at turbulent intensities below the critical threshold for DDT. We carried out a suite of 3D DNS of turbulent flames for a broad range of turbulent intensities and system sizes using a simplified, single-step, Arrhenius-type reaction kinetics. Our results show that at the later stages of the explosion, as the turbulence intensity increases prior to the possible onset of DDT, the flame front will become violently unstable. We find that the burning rate exhibits periodic pulsations with the energy release rate varying by almost an order of magnitude. Furthermore, such flame pulsations can produce pressure waves and shocks as the flame speed approaches the critical Chapman-Jouguet deflagration speed. Finally, in contrast with the current theoretical understanding, such fast turbulent flames can propagate at speeds, which are much higher than the characteristic speeds of turbulent fluctuations. These effects can qualitatively change the dynamics of the explosion and, therefore, must be properly accounted for in the turbulent-flame subgrid-scale models.

SIERRA/Aero Theory Manual Version 4.46.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sierra Thermal/Fluid Team

2017-09-01

SIERRA/Aero is a two and three dimensional, node-centered, edge-based finite volume code that approximates the compressible Navier-Stokes equations on unstructured meshes. It is applicable to inviscid and high Reynolds number laminar and turbulent flows. Currently, two classes of turbulence models are provided: Reynolds Averaged Navier-Stokes (RANS) and hybrid methods such as Detached Eddy Simulation (DES). Large Eddy Simulation (LES) models are currently under development. The gas may be modeled either as ideal, or as a non-equilibrium, chemically reacting mixture of ideal gases. This document describes the mathematical models contained in the code, as well as certain implementation details. First, themore » governing equations are presented, followed by a description of the spatial discretization. Next, the time discretization is described, and finally the boundary conditions. Throughout the document, SIERRA/ Aero is referred to simply as Aero for brevity.« less

SIERRA/Aero Theory Manual Version 4.44

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sierra Thermal /Fluid Team

2017-04-01

SIERRA/Aero is a two and three dimensional, node-centered, edge-based finite volume code that approximates the compressible Navier-Stokes equations on unstructured meshes. It is applicable to inviscid and high Reynolds number laminar and turbulent flows. Currently, two classes of turbulence models are provided: Reynolds Averaged Navier-Stokes (RANS) and hybrid methods such as Detached Eddy Simulation (DES). Large Eddy Simulation (LES) models are currently under development. The gas may be modeled either as ideal, or as a non-equilibrium, chemically reacting mixture of ideal gases. This document describes the mathematical models contained in the code, as well as certain implementation details. First, themore » governing equations are presented, followed by a description of the spatial discretization. Next, the time discretization is described, and finally the boundary conditions. Throughout the document, SIERRA/ Aero is referred to simply as Aero for brevity.« less

Anomalous pinch of turbulent plasmas driven by the magnetic-drift-induced Lorentz force through the Stokes-Einstein relation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Shaojie, E-mail: wangsj@ustc.edu.cn

It is found that the Lorentz force generated by the magnetic drift drives a generic plasma pinch flux of particle, energy and momentum through the Stokes-Einstein relation. The proposed theoretical model applies for both electrons and ions, trapped particles, and passing particles. An anomalous parallel current pinch due to the electrostatic turbulence with long parallel wave-length is predicted.

Laminar Flow in the Ocean Ekman Layer

NASA Astrophysics Data System (ADS)

Woods, J. T. H.

INTRODUCTION THE EFFECT OF A STABLE DENSITY GRADIENT THE FATAL FLAW FLOW VISUALIZATION THE DISCOVERY OF LAMINAR FLOW FINE STRUCTURE WAVE-INDUCED SHEAR INSTABILITY BILLOW TURBULENCE REVERSE TRANSITION REVISED PARADIGM ONE-DIMENSIONAL MODELLING OF THE UPPER OCEAN DIURNAL VARIATION BUOYANT CONVECTION BILLOW TURBULENCE IN THE DIURNAL THERMOCLINE CONSEQUENCES FOR THE EKMAN CURRENT PROFILE SOLAR RADIATION APPLICATIONS Slippery Seas of Acapulco Pollution Afternoon Effect in Sonar Patchiness Fisheries Climate DISCUSSION CONCLUSION REFERENCES

Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary-Layer Interaction

NASA Technical Reports Server (NTRS)

Vyas, Manan A.; Yoder, Dennis A.; Gaitonde, Datta V.

2018-01-01

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Comparisons with experimental data showed a sensitivity of the current prediction to the modeling of the sidewalls. This was found to be common among various computational studies in the literature where periodic boundary conditions were used in the spanwise direction, as was the case in the present work. Thus, although the experiment was quasi-two-dimensional, the present simulation was determined to be two-dimensional. Quantities present in the exact equation of the Reynolds-stress transport, i.e., production, molecular diffusion, turbulent transport, pressure diffusion, pressure strain, dissipation, and turbulent mass flux were calculated. Reynolds-stress budgets were compared with past large-eddy simulation and direct numerical simulation datasets in the undisturbed portion of the turbulent boundary layer to validate the current approach. The budgets in SBLI showed the growth in the production term for the primary normal stress and energy transfer mechanism was led by the pressure strain term in the secondary normal stresses. The pressure diffusion term, commonly assumed as negligible by turbulence model developers, was shown to be small but non-zero in the normal stress budgets, however it played a key role in the primary shear stress budget.

Reconnection-Driven Magnetohydrodynamic Turbulence in a Simulated Coronal-Hole Jet

NASA Technical Reports Server (NTRS)

Uritskiy, Vadim M.; Roberts, Merrill A.; DeVore, C. Richard; Karpen, Judith T.

2017-01-01

Extreme-ultraviolet and X-ray jets occur frequently in magnetically open coronal holes on the Sun, especially at high solar latitudes. Some of these jets are observed by white-light coronagraphs as they propagate through the outer corona toward the inner heliosphere, and it has been proposed that they give rise to microstreams and torsional Alfven waves detected in situ in the solar wind. To predict and understand the signatures of coronal-hole jets, we have performed a detailed statistical analysis of such a jet simulated with an adaptively refined magnetohydrodynamics model. The results confirm the generation and persistence of three-dimensional, reconnection-driven magnetic turbulence in the simulation. We calculate the spatial correlations of magnetic fluctuations within the jet and find that they agree best with the Meuller - Biskamp scaling model including intermittent current sheets of various sizes coupled via hydrodynamic turbulent cascade. The anisotropy of the magnetic fluctuations and the spatial orientation of the current sheets are consistent with an ensemble of nonlinear Alfven waves. These properties also reflect the overall collimated jet structure imposed by the geometry of the reconnecting magnetic field. A comparison with Ulysses observations shows that turbulence in the jet wake is in quantitative agreement with that in the fast solar wind.

Extracting atmospheric turbulence and aerosol characteristics from passive imagery

NASA Astrophysics Data System (ADS)

Reinhardt, Colin N.; Wayne, D.; McBryde, K.; Cauble, G.

2013-09-01

Obtaining accurate, precise and timely information about the local atmospheric turbulence and extinction conditions and aerosol/particulate content remains a difficult problem with incomplete solutions. It has important applications in areas such as optical and IR free-space communications, imaging systems performance, and the propagation of directed energy. The capability to utilize passive imaging data to extract parameters characterizing atmospheric turbulence and aerosol/particulate conditions would represent a valuable addition to the current piecemeal toolset for atmospheric sensing. Our research investigates an application of fundamental results from optical turbulence theory and aerosol extinction theory combined with recent advances in image-quality-metrics (IQM) and image-quality-assessment (IQA) methods. We have developed an algorithm which extracts important parameters used for characterizing atmospheric turbulence and extinction along the propagation channel, such as the refractive-index structure parameter C2n , the Fried atmospheric coherence width r0 , and the atmospheric extinction coefficient βext , from passive image data. We will analyze the algorithm performance using simulations based on modeling with turbulence modulation transfer functions. An experimental field campaign was organized and data were collected from passive imaging through turbulence of Siemens star resolution targets over several short littoral paths in Point Loma, San Diego, under conditions various turbulence intensities. We present initial results of the algorithm's effectiveness using this field data and compare against measurements taken concurrently with other standard atmospheric characterization equipment. We also discuss some of the challenges encountered with the algorithm, tasks currently in progress, and approaches planned for improving the performance in the near future.

Impact of E × B flow shear on turbulence and resulting power fall-off width in H-mode plasmas in experimental advanced superconducting tokamak

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Q. Q., E-mail: yangqq@ipp.ac.cn; Zhong, F. C., E-mail: gsxu@ipp.ac.cn, E-mail: fczhong@dhu.edu.cn; Jia, M. N.

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 lengthmore » 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.« less

THE DYNAMICAL GENERATION OF CURRENT SHEETS IN ASTROPHYSICAL PLASMA TURBULENCE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Howes, Gregory G.

2016-08-20

Turbulence profoundly affects particle transport and plasma heating in many astrophysical plasma environments, from galaxy clusters to the solar corona and solar wind to Earth's magnetosphere. Both fluid and kinetic simulations of plasma turbulence ubiquitously generate coherent structures, in the form of current sheets, at small scales, and the locations of these current sheets appear to be associated with enhanced rates of dissipation of the turbulent energy. Therefore, illuminating the origin and nature of these current sheets is critical to identifying the dominant physical mechanisms of dissipation, a primary aim at the forefront of plasma turbulence research. Here, we presentmore » evidence from nonlinear gyrokinetic simulations that strong nonlinear interactions between counterpropagating Alfvén waves, or strong Alfvén wave collisions, are a natural mechanism for the generation of current sheets in plasma turbulence. Furthermore, we conceptually explain this current sheet development in terms of the nonlinear dynamics of Alfvén wave collisions, showing that these current sheets arise through constructive interference among the initial Alfvén waves and nonlinearly generated modes. The properties of current sheets generated by strong Alfvén wave collisions are compared to published observations of current sheets in the Earth's magnetosheath and the solar wind, and the nature of these current sheets leads to the expectation that Landau damping of the constituent Alfvén waves plays a dominant role in the damping of turbulently generated current sheets.« less

A Laboratory Study of River Discharges into Shallow Seas

NASA Astrophysics Data System (ADS)

Crawford, T. J.; Linden, P. F.

2016-02-01

We present an experimental study that aims to simulate the buoyancy driven coastal currents produced by estuarine freshwater discharges into the ocean. The currents are generated inside a rotating tank filled with saltwater by the continuous release of buoyant freshwater from a source structure located at the fluid surface. The freshwater is discharged horizontally from a finite-depth source, giving rise to significant momentum-flux effects and a non-zero potential vorticity. We perform a parametric study in which we vary the rotation rate, freshwater discharge magnitude, the density difference and the source cross-sectional area. The parameter values are chosen to match the regimes appropriate to the River Rhine and River Elbe when entering the North Sea. Persistent features of an anticyclonic outflow vortex and a propagating boundary current were identified and their properties quantified. We also present a finite potential vorticity, geostrophic model that provides theoretical predictions for the current height, width and velocity as functions of the experimental parameters. The experiments and model are compared with each other in terms of a set of non-dimensional parameters identified in the theoretical analysis of the problem. Good agreement between the model and the experimental data is found. The effect of mixing in the turbulent ocean is also addressed with the addition of an oscillating grid to the experimental setup. The grid generates turbulence in the saltwater ambient that is designed to represent the mixing effects of the wind, tides and bathymetry in a shallow shelf sea. The impact of the addition of turbulence is discussed in terms of the experimental data and through modifications to the theoretical model to include mixing. Once again, good agreement is seen between the experiments and the model.

Propagation of a Gaussian-beam wave in general anisotropic turbulence

NASA Astrophysics Data System (ADS)

Andrews, L. C.; Phillips, R. L.; Crabbs, R.

2014-10-01

Mathematical models for a Gaussian-beam wave propagating through anisotropic non-Kolmogorov turbulence have been developed in the past by several researchers. In previous publications, the anisotropic spatial power spectrum model was based on the assumption that propagation was in the z direction with circular symmetry maintained in the orthogonal xy-plane throughout the path. In the present analysis, however, the anisotropic spectrum model is no longer based on a single anisotropy parameter—instead, two such parameters are introduced in the orthogonal xyplane so that circular symmetry in this plane is no longer required. In addition, deviations from the 11/3 power-law behavior in the spectrum model are allowed by assuming power-law index variations 3 < α < 4 . In the current study we develop theoretical models for beam spot size, spatial coherence, and scintillation index that are valid in weak irradiance fluctuation regimes as well as in deep turbulence, or strong irradiance fluctuation regimes. These new results are compared with those derived from the more specialized anisotropic spectrum used in previous analyses.

Spatiotemporal structure of wind farm-atmospheric boundary layer interactions

NASA Astrophysics Data System (ADS)

Cervarich, Matthew; Baidya Roy, Somnath; Zhou, Liming

2013-04-01

Wind power is currently one of the fastest growing energy sources in the world. Most of the growth is in the utility sector consisting of large wind farms with numerous industrial-scale wind turbines. Wind turbines act as a sink of mean kinetic energy and a source of turbulent kinetic energy in the atmospheric boundary layer (ABL). In doing so, they modify the ABL profiles and land-atmosphere exchanges of energy, momentum, mass and moisture. This project explores theses interactions using remote sensing data and numerical model simulations. The domain is central Texas where 4 of the world's largest wind farms are located. A companion study of seasonally-averaged Land Surface Temperature data from the Moderate Resolution Imaging Spectroradiometer (MODIS) on TERRA and AQUA satellites shows a warming signal at night and a mixed cooling/warming signal during the daytime within the wind farms. In the present study, wind farm-ABL interactions are simulated with the Weather Research and Forecasting (WRF) model. The simulations show that the model is capable of replicating the observed signal in land surface temperature. Moreover, similar warming/cooling effect, up to 1C, was observed in seasonal mean 2m air temperature as well. Further analysis show that enhanced turbulent mixing in the rotor wakes is responsible for the impacts on 2m and surface air temperatures. The mixing is due to 2 reasons: (i) turbulent momentum transport to compensate the momentum deficit in the wakes of the turbines and (ii) turbulence generated due to motion of turbine rotors. Turbulent mixing also alters vertical profiles of moisture. Changes in land-atmosphere temperature and moisture gradient and increase in turbulent mixing leads to more than 10% change in seasonal mean surface sensible and latent heat flux. Given the current installed capacity and the projected installation across the world, wind farms are likely becoming a major driver of anthropogenic land use change on Earth. Hence, understanding WF-ABL interactions and its effects is of significant scientific and societal importance.

Large-eddy simulation of dust-uplift by a haboob density current

NASA Astrophysics Data System (ADS)

Huang, Qian; Marsham, John H.; Tian, Wenshou; Parker, Douglas J.; Garcia-Carreras, Luis

2018-04-01

Cold pool outflows have been shown from both observations and convection-permitting models to be a dominant source of dust emissions ("haboobs") in the summertime Sahel and Sahara, and to cause dust uplift over deserts across the world. In this paper Met Office Large Eddy Model (LEM) simulations, which resolve the turbulence within the cold-pools much better than previous studies of haboobs with convection-permitting models, are used to investigate the winds that uplift dust in cold pools, and the resultant dust transport. In order to simulate the cold pool outflow, an idealized cooling is added in the model during the first 2 h of 5.7 h run time. Given the short duration of the runs, dust is treated as a passive tracer. Dust uplift largely occurs in the "head" of the density current, consistent with the few existing observations. In the modeled density current dust is largely restricted to the lowest, coldest and well mixed layers of the cold pool outflow (below around 400 m), except above the "head" of the cold pool where some dust reaches 2.5 km. This rapid transport to above 2 km will contribute to long atmospheric lifetimes of large dust particles from haboobs. Decreasing the model horizontal grid-spacing from 1.0 km to 100 m resolves more turbulence, locally increasing winds, increasing mixing and reducing the propagation speed of the density current. Total accumulated dust uplift is approximately twice as large in 1.0 km runs compared with 100 m runs, suggesting that for studying haboobs in convection-permitting runs the representation of turbulence and mixing is significant. Simulations with surface sensible heat fluxes representative of those from a desert region during daytime show that increasing surface fluxes slows the density current due to increased mixing, but increase dust uplift rates, due to increased downward transport of momentum to the surface.

Pressure oscillations and instability of working processes in the combustion chambers of solid rocket motors

NASA Astrophysics Data System (ADS)

Emelyanov, V. N.; Teterina, I. V.; Volkov, K. N.; Garkushev, A. U.

2017-06-01

Metal particles are widely used in space engineering to increase specific impulse and to supress acoustic instability of intra-champber processes. A numerical analysis of the internal injection-driven turbulent gas-particle flows is performed to improve the current understanding and modeling capabilities of the complex flow characteristics in the combustion chambers of solid rocket motors (SRMs) in presence of forced pressure oscillations. The two-phase flow is simulated with a combined Eulerian-Lagrangian approach. The Reynolds-averaged Navier-Stokes equations and transport equations of k - ε model are solved numerically for the gas. The particulate phase is simulated through a Lagrangian deterministic and stochastic tracking models to provide particle trajectories and particle concentration. The results obtained highlight the crucial significance of the particle dispersion in turbulent flowfield and high potential of statistical methods. Strong coupling between acoustic oscillations, vortical motion, turbulent fluctuations and particle dynamics is observed.

Simulation of two-dimensional turbulent flows in a rotating annulus

NASA Astrophysics Data System (ADS)

Storey, Brian D.

2004-05-01

Rotating water tank experiments have been used to study fundamental processes of atmospheric and geophysical turbulence in a controlled laboratory setting. When these tanks are undergoing strong rotation the forced turbulent flow becomes highly two dimensional along the axis of rotation. An efficient numerical method has been developed for simulating the forced quasi-geostrophic equations in an annular geometry to model current laboratory experiments. The algorithm employs a spectral method with Fourier series and Chebyshev polynomials as basis functions. The algorithm has been implemented on a parallel architecture to allow modelling of a wide range of spatial scales over long integration times. This paper describes the derivation of the model equations, numerical method, testing and performance of the algorithm. Results provide reasonable agreement with the experimental data, indicating that such computations can be used as a predictive tool to design future experiments.

Two-Equation Low-Reynolds-Number Turbulence Modeling of Transitional Boundary Layer Flows Characteristic of Gas Turbine Blades. Ph.D. Thesis. Final Contractor Report

NASA Technical Reports Server (NTRS)

Schmidt, Rodney C.; Patankar, Suhas V.

1988-01-01

The use of low Reynolds number (LRN) forms of the k-epsilon turbulence model in predicting transitional boundary layer flow characteristic of gas turbine blades is developed. The research presented consists of: (1) an evaluation of two existing models; (2) the development of a modification to current LRN models; and (3) the extensive testing of the proposed model against experimental data. The prediction characteristics and capabilities of the Jones-Launder (1972) and Lam-Bremhorst (1981) LRN k-epsilon models are evaluated with respect to the prediction of transition on flat plates. Next, the mechanism by which the models simulate transition is considered and the need for additional constraints is discussed. Finally, the transition predictions of a new model are compared with a wide range of different experiments, including transitional flows with free-stream turbulence under conditions of flat plate constant velocity, flat plate constant acceleration, flat plate but strongly variable acceleration, and flow around turbine blade test cascades. In general, calculational procedure yields good agreement with most of the experiments.

Turbulent Flame Propagation Characteristics of High Hydrogen Content Fuels

DOE Office of Scientific and Technical Information (OSTI.GOV)

Seitzman, Jerry; Lieuwen, Timothy

2014-09-30

This final report describes the results of an effort to better understand turbulent flame propagation, especially at conditions relevant to gas turbines employing fuels with syngas or hydrogen mixtures. Turbulent flame speeds were measured for a variety of hydrogen/carbon monoxide (H2/CO) and hydrogen/methane (H2/CH4) fuel mixtures with air as the oxidizer. The measurements include global consumption speeds (ST,GC) acquired in a turbulent jet flame at pressures of 1-10 atm and local displacement speeds (ST,LD) acquired in a low-swirl burner at atmospheric pressure. The results verify the importance of fuel composition in determining turbulent flame speeds. For example, different fuel-air mixturesmore » having the same unstretched laminar flame speed (SL,0) but different fuel compositions resulted in significantly different ST,GC for the same turbulence levels (u'). This demonstrates the weakness of turbulent flame speed correlations based simply on u'/SL,0. The results were analyzed using a steady-steady leading points concept to explain the sensitivity of turbulent burning rates to fuel (and oxidizer) composition. Leading point theories suggest that the premixed turbulent flame speed is controlled by the flame front characteristics at the flame brush leading edge, or, in other words, by the flamelets that advance farthest into the unburned mixture (the so-called leading points). For negative Markstein length mixtures, this is assumed to be close to the maximum stretched laminar flame speed (SL,max) for the given fuel-oxidizer mixture. For the ST,GC measurements, the data at a given pressure were well-correlated with an SL,max scaling. However the variation with pressure was not captured, which may be due to non-quasi-steady effects that are not included in the current model. For the ST,LD data, the leading points model again faithfully captured the variation of turbulent flame speed over a wide range of fuel-compositions and turbulence intensities. These results provide evidence that the leading points model can provide useful predictions of turbulent flame speed over a wide range of operating conditions and flow geometries.« less

Meso-scale turbulence in living fluids

PubMed Central

Wensink, Henricus H.; Dunkel, Jörn; Heidenreich, Sebastian; Drescher, Knut; Goldstein, Raymond E.; Löwen, Hartmut; Yeomans, Julia M.

2012-01-01

Turbulence is ubiquitous, from oceanic currents to small-scale biological and quantum systems. Self-sustained turbulent motion in microbial suspensions presents an intriguing example of collective dynamical behavior among the simplest forms of life and is important for fluid mixing and molecular transport on the microscale. The mathematical characterization of turbulence phenomena in active nonequilibrium fluids proves even more difficult than for conventional liquids or gases. It is not known which features of turbulent phases in living matter are universal or system-specific or which generalizations of the Navier–Stokes equations are able to describe them adequately. Here, we combine experiments, particle simulations, and continuum theory to identify the statistical properties of self-sustained meso-scale turbulence in active systems. To study how dimensionality and boundary conditions affect collective bacterial dynamics, we measured energy spectra and structure functions in dense Bacillus subtilis suspensions in quasi-2D and 3D geometries. Our experimental results for the bacterial flow statistics agree well with predictions from a minimal model for self-propelled rods, suggesting that at high concentrations the collective motion of the bacteria is dominated by short-range interactions. To provide a basis for future theoretical studies, we propose a minimal continuum model for incompressible bacterial flow. A detailed numerical analysis of the 2D case shows that this theory can reproduce many of the experimentally observed features of self-sustained active turbulence. PMID:22908244

Meso-scale turbulence in living fluids.

PubMed

Wensink, Henricus H; Dunkel, Jörn; Heidenreich, Sebastian; Drescher, Knut; Goldstein, Raymond E; Löwen, Hartmut; Yeomans, Julia M

2012-09-04

Turbulence is ubiquitous, from oceanic currents to small-scale biological and quantum systems. Self-sustained turbulent motion in microbial suspensions presents an intriguing example of collective dynamical behavior among the simplest forms of life and is important for fluid mixing and molecular transport on the microscale. The mathematical characterization of turbulence phenomena in active nonequilibrium fluids proves even more difficult than for conventional liquids or gases. It is not known which features of turbulent phases in living matter are universal or system-specific or which generalizations of the Navier-Stokes equations are able to describe them adequately. Here, we combine experiments, particle simulations, and continuum theory to identify the statistical properties of self-sustained meso-scale turbulence in active systems. To study how dimensionality and boundary conditions affect collective bacterial dynamics, we measured energy spectra and structure functions in dense Bacillus subtilis suspensions in quasi-2D and 3D geometries. Our experimental results for the bacterial flow statistics agree well with predictions from a minimal model for self-propelled rods, suggesting that at high concentrations the collective motion of the bacteria is dominated by short-range interactions. To provide a basis for future theoretical studies, we propose a minimal continuum model for incompressible bacterial flow. A detailed numerical analysis of the 2D case shows that this theory can reproduce many of the experimentally observed features of self-sustained active turbulence.

Gyrokinetic simulation of residual turbulence in transport barriers

NASA Astrophysics Data System (ADS)

Jenko, Frank; Told, Daniel; Goerler, Tobias; Brunner, Stephan; Sautter, Olivier

2011-10-01

One of the ultimate aims for gyrokinetic simulation is to describe the formation and evolution of transport barriers. An important step in that direction is the study of the residual turbulence in established barriers - a challenging task in itself, given that a wide range of spatio-temporal scales can be involved. In the present work, we employ the physically comprehensive, nonlocal gyrokinetic turbulence code GENE to study turbulence in both core and edge transport barriers. First, we apply GENE to a set of discharges in the TCV tokamak which exhibit electron ITBs. Nonlinear gyrokinetic simulations are used to examine the influence of a varying current profile on the strength of the barrier. For each case, the transport spectra reveal how much transport (for each channel) is done in the low-k, medium-k, and high-k regimes, respectively. The role of ETG turbulence is discussed. Second, we explore the role of ETG turbulence in a typical ASDEX Upgrade H-mode discharge. Numerical convergence is carefully examined, and new insights on the characteristics of ETG turbulence in the edge will be discussed, focusing particularly on the role of streamers, which had been found to be a necessary ingredient for experimentally relevant ETG transport in core plasmas. The radial dependence of the resulting electron heat diffusivity is also examined and a simple ETG model is presented which can be used in future edge modeling efforts.

Air-flow distortion and turbulence statistics near an animal facility

NASA Astrophysics Data System (ADS)

Prueger, J. H.; Eichinger, W. E.; Hipps, L. E.; Hatfield, J. L.; Cooper, D. I.

The emission and dispersion of particulates and gases from concentrated animal feeding operations (CAFO) at local to regional scales is a current issue in science and society. The transport of particulates, odors and toxic chemical species from the source into the local and eventually regional atmosphere is largely determined by turbulence. Any models that attempt to simulate the dispersion of particles must either specify or assume various statistical properties of the turbulence field. Statistical properties of turbulence are well documented for idealized boundary layers above uniform surfaces. However, an animal production facility is a complex surface with structures that act as bluff bodies that distort the turbulence intensity near the buildings. As a result, the initial release and subsequent dispersion of effluents in the region near a facility will be affected by the complex nature of the surface. Previous Lidar studies of plume dispersion over the facility used in this study indicated that plumes move in complex yet organized patterns that would not be explained by the properties of turbulence generally assumed in models. The objective of this study was to characterize the near-surface turbulence statistics in the flow field around an array of animal confinement buildings. Eddy covariance towers were erected in the upwind, within the building array and downwind regions of the flow field. Substantial changes in turbulence intensity statistics and turbulence-kinetic energy (TKE) were observed as the mean wind flow encountered the building structures. Spectra analysis demonstrated unique distribution of the spectral energy in the vertical profile above the buildings.

Evaluation of subgrid-scale models in large-eddy simulations of turbulent flow in a centrifugal pump impeller

NASA Astrophysics Data System (ADS)

Yang, Zhengjun; Wang, Fujun; Zhou, Peijian

2012-09-01

The current research of large eddy simulation (LES) of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale (SGS) model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry (PIV) and laser Doppler velocimetry (LDV). The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.

Some practical turbulence modeling options for Reynolds-averaged full Navier-Stokes calculations of three-dimensional flows

NASA Technical Reports Server (NTRS)

Bui, Trong T.

1993-01-01

New turbulence modeling options recently implemented for the 3-D version of Proteus, a Reynolds-averaged compressible Navier-Stokes code, are described. The implemented turbulence models include: the Baldwin-Lomax algebraic model, the Baldwin-Barth one-equation model, the Chien k-epsilon model, and the Launder-Sharma k-epsilon model. Features of this turbulence modeling package include: well documented and easy to use turbulence modeling options, uniform integration of turbulence models from different classes, automatic initialization of turbulence variables for calculations using one- or two-equation turbulence models, multiple solid boundaries treatment, and fully vectorized L-U solver for one- and two-equation models. Validation test cases include the incompressible and compressible flat plate turbulent boundary layers, turbulent developing S-duct flow, and glancing shock wave/turbulent boundary layer interaction. Good agreement is obtained between the computational results and experimental data. Sensitivity of the compressible turbulent solutions with the method of y(sup +) computation, the turbulent length scale correction, and some compressibility corrections are examined in detail. The test cases show that the highly optimized one-and two-equation turbulence models can be used in routine 3-D Navier-Stokes computations with no significant increase in CPU time as compared with the Baldwin-Lomax algebraic model.

Flow-combustion interactions in ducted flameholder-stabilized premixed flames

NASA Astrophysics Data System (ADS)

Soteriou, Marios; Arienti, Marco; Erickson, Robert

2006-11-01

Turbulent premixed combustion is present in many power generation and propulsion systems due to its large energy conversion rate (as compared to non-premixed combustion) and its potential for reduced emissions (at the lean limit). As a result, the study of turbulent premixed flames has received substantial attention in the past through experiment, analysis and simulation. In the recent past, unsteady Computational Fluid Dynamics (CFD) based models have been increasingly leveraged towards the in depth study of the physics of turbulent premixed flames. The bulk of this effort focuses on the response of the flame to turbulence. In contrast, we focus on the opposite problem, i.e. the modification of the turbulent flowfield by the flame. This topic has also received some attention but with a strong emphasis on planar (in the mean), flames propagating normal to the flow. Instead, we focus on flameholder-stabilized ducted flames, i.e. ones in which the flame is confined and substantially inclined to the incoming flow. The fundamental mechanisms by which the flame impacts the flow, i.e. dilatation, baroclinic vorticity generation and molecular diffusion enhancement are discussed in detail and their relative impact quantified. Limitations of modeling these mechanisms in current state of the art CFD models are also addressed.

Aqueous turbulence structure immediately adjacent to the air - water interface and interfacial gas exchange

NASA Astrophysics Data System (ADS)

Wang, Binbin

Air-sea interaction and the interfacial exchange of gas across the air-water interface are of great importance in coupled atmospheric-oceanic environmental systems. Aqueous turbulence structure immediately adjacent to the air-water interface is the combined result of wind, surface waves, currents and other environmental forces and plays a key role in energy budgets, gas fluxes and hence the global climate system. However, the quantification of turbulence structure sufficiently close to the air-water interface is extremely difficult. The physical relationship between interfacial gas exchange and near surface turbulence remains insufficiently investigated. This dissertation aims to measure turbulence in situ in a complex environmental forcing system on Lake Michigan and to reveal the relationship between turbulent statistics and the CO2 flux across the air-water interface. The major objective of this dissertation is to investigate the physical control of the interfacial gas exchange and to provide a universal parameterization of gas transfer velocity from environmental factors, as well as to propose a mechanistic model for the global CO2 flux that can be applied in three dimensional climate-ocean models. Firstly, this dissertation presents an advanced measurement instrument, an in situ free floating Particle Image Velocimetry (FPIV) system, designed and developed to investigate the small scale turbulence structure immediately below the air-water interface. Description of hardware components, design of the system, measurement theory, data analysis procedure and estimation of measurement error were provided. Secondly, with the FPIV system, statistics of small scale turbulence immediately below the air-water interface were investigated under a variety of environmental conditions. One dimensional wave-number spectrum and structure function sufficiently close to the water surface were examined. The vertical profiles of turbulent dissipation rate were intensively studied. Comparison between the turbulence structures measured during the wind wave initiation period and those obtained during the growing period was presented. Significant wave effects on near surface turbulence were found. A universal scaling law was proposed to parameterize turbulent dissipation rate immediately below the air-water interface with friction velocity, significant wave height and wave age. Finally, the gas transfer velocity was measured with a floating chamber (FC) system, along with simultaneously FPIV measurements. Turbulent dissipation rate both at the interface and at a short distance away from the interface (~ 10 cm) were analyzed and used to examine the small scale eddy model. The model coefficient was found to be dependent on the level of turbulence, instead of being a constant. An empirical relationship between the model coefficient and turbulent dissipation rate was provided, which improved the accuracy of the gas transfer velocity estimation by more than 100% for data acquired. Other data from the literature also supported this empirical relation. Furthermore, the relationship between model coefficient and turbulent Reynolds number was also investigated. In addition to physical control of gas exchange, the disturbance on near surface hydrodynamics by the FC was also discussed. Turbulent dissipation rates are enhanced at the short distance away from the interface, while the surface dissipation rates do not change significantly.

Design of Bi-Directional Hydrofoils for Tidal Current Turbines

NASA Astrophysics Data System (ADS)

Nedyalkov, Ivaylo; Wosnik, Martin

2015-11-01

Tidal Current Turbines operate in flows which reverse direction. Bi-directional hydrofoils have rotational symmetry and allow such turbines to operate without the need for pitch or yaw control, decreasing the initial and maintenance costs. A numerical test-bed was developed to automate the simulations of hydrofoils in OpenFOAM and was utilized to simulate the flow over eleven classes of hydrofoils comprising a total of 700 foil shapes at different angles of attack. For promising candidate foil shapes physical models of 75 mm chord and 150 mm span were fabricated and tested in the University of New Hampshire High-Speed Cavitation Tunnel (HiCaT). The experimental results were compared to the simulations for model validation. The numerical test-bed successfully generated simulations for a wide range of foil shapes, although, as expected, the k - ω - SST turbulence model employed here was not adequate for some of the foils and for large angles of attack at which separation occurred. An optimization algorithm is currently being coupled with the numerical test-bed and additional turbulence models will be implemented in the future.

Comprehensive Approaches to Multiphase Flows in Geophysics - Application to nonisothermal, nonhomogenous, unsteady, large-scale, turbulent dusty clouds I. Hydrodynamic and Thermodynamic RANS and LES Models

DOE Office of Scientific and Technical Information (OSTI.GOV)

S. Dartevelle

2005-09-05

The objective of this manuscript is to fully derive a geophysical multiphase model able to ''accommodate'' different multiphase turbulence approaches; viz., the Reynolds Averaged Navier-Stokes (RANS), the Large Eddy Simulation (LES), or hybrid RANSLES. This manuscript is the first part of a larger geophysical multiphase project--lead by LANL--that aims to develop comprehensive modeling tools for large-scale, atmospheric, transient-buoyancy dusty jets and plume (e.g., plinian clouds, nuclear ''mushrooms'', ''supercell'' forest fire plumes) and for boundary-dominated geophysical multiphase gravity currents (e.g., dusty surges, diluted pyroclastic flows, dusty gravity currents in street canyons). LES is a partially deterministic approach constructed on either amore » spatial- or a temporal-separation between the large and small scales of the flow, whereas RANS is an entirely probabilistic approach constructed on a statistical separation between an ensemble-averaged mean and higher-order statistical moments (the so-called ''fluctuating parts''). Within this specific multiphase context, both turbulence approaches are built up upon the same phasic binary-valued ''function of presence''. This function of presence formally describes the occurrence--or not--of any phase at a given position and time and, therefore, allows to derive the same basic multiphase Navier-Stokes model for either the RANS or the LES frameworks. The only differences between these turbulence frameworks are the closures for the various ''turbulence'' terms involving the unknown variables from the fluctuating (RANS) or from the subgrid (LES) parts. Even though the hydrodynamic and thermodynamic models for RANS and LES have the same set of Partial Differential Equations, the physical interpretations of these PDEs cannot be the same, i.e., RANS models an averaged field, while LES simulates a filtered field. In this manuscript, we also demonstrate that this multiphase model fully fulfills the second law of thermodynamics and fulfills the necessary requirements for a well-posed initial-value problem. In the next manuscripts, we will further develop specific closures for multiphase RANS, LES, and hybrid-LES.« less

Flamelet Model Application for Non-Premixed Turbulent Combustion

NASA Technical Reports Server (NTRS)

Secundov, A.; Bezgin, L.; Buriko, Yu.; Guskov, O.; Kopchenov, V.; Laskin, I.; Lomkov, K.; Tshepin, S.; Volkov, D.; Zaitsev, S.

1996-01-01

The current Final Report contains results of the study which was performed in Scientific Research Center 'ECOLEN' (Moscow, Russia). The study concerns the development and verification of non-expensive approach for modeling of supersonic turbulent diffusion flames based on flamelet consideration of the chemistry/turbulence interaction (FL approach). Research work included: development of the approach and CFD tests of the flamelet model for supersonic jet flames; development of the simplified procedure for solution of the flamelet equations based on partial equilibrium chemistry assumption; study of the flame ignition/extinction predictions provided by flamelet model. The performed investigation demonstrated that FL approach allowed to describe satisfactory main features of supersonic H 2/air jet flames. Model demonstrated also high capabilities for reduction of the computational expenses in CFD modeling of the supersonic flames taking into account detailed oxidation chemistry. However, some disadvantages and restrictions of the existing version of approach were found in this study. They were: (1) inaccuracy in predictions of the passive scalar statistics by our turbulence model for one of the considered test cases; and (2) applicability of the available version of the flamelet model to flames without large ignition delay distance only. Based on the results of the performed investigation, we formulated and submitted to the National Aeronautics and Space Administration our Project Proposal for the next step research directed toward further improvement of the FL approach.

Cancellation exponent and multifractal structure in two-dimensional magnetohydrodynamics: direct numerical simulations and Lagrangian averaged modeling.

PubMed

Graham, Jonathan Pietarila; Mininni, Pablo D; Pouquet, Annick

2005-10-01

We present direct numerical simulations and Lagrangian averaged (also known as alpha model) simulations of forced and free decaying magnetohydrodynamic turbulence in two dimensions. The statistics of sign cancellations of the current at small scales is studied using both the cancellation exponent and the fractal dimension of the structures. The alpha model is found to have the same scaling behavior between positive and negative contributions as the direct numerical simulations. The alpha model is also able to reproduce the time evolution of these quantities in free decaying turbulence. At large Reynolds numbers, an independence of the cancellation exponent with the Reynolds numbers is observed.

Aerothermal modeling program, phase 2

NASA Technical Reports Server (NTRS)

Mongia, H. C.; Patankar, S. V.; Murthy, S. N. B.; Sullivan, J. P.; Samuelsen, G. S.

1985-01-01

The main objectives of the Aerothermal Modeling Program, Phase 2 are: to develop an improved numerical scheme for incorporation in a 3-D combustor flow model; to conduct a benchmark quality experiment to study the interaction of a primary jet with a confined swirling crossflow and to assess current and advanced turbulence and scalar transport models; and to conduct experimental evaluation of the air swirler interaction with fuel injectors, assessments of current two-phase models, and verification the improved spray evaporation/dispersion models.

Performance of Smagorinsky and dynamic models in near surface turbulence

NASA Astrophysics Data System (ADS)

Brasseur, James G.; Juneja, Anurag

1997-11-01

In LES of high-Reynolds-number wall bounded turbulence such as the atmospheric boundary layer (ABL), a viscous sublayer either does not exist or is within the first grid cell, and some integral scale motions are necessarily under-resolved at the first few grid locations. Here the subgrid terms dominate the evolution of resolved velocity and the SGS model performance becomes crucial. To develop improved closures for surface layer turbulence (under-resolved and anisotropic), we explore (a) why current SGS closures fail and (b) what needs to be fixed. We evaluate the performance of the Smagorinsky and dynamic models using DNS data from shear- and buoyancy-driven turbulence as a function of filter cutoff location. We find that the underlying assumption of good alignment between the subgrid stress and resolved strain-rate tensors is not correct in general. More importantly, the Smagorinsky model incorrectly predicts a strong preference in the direction of the SGS stress divergence vector, a spurious prediction that is directly related to the anisotropic structure of the resolved turbulence field. This, and its under-estimation of the SGS pressure gradient, are likely sources of the errors observed in LES of the ABL. Whereas the dynamic formulations do a better job predicting some SGS dynamics, the model fails when the filter cutoff is near an integral scale, and predicts unreasonable fluctuation levels-- although performance is sensitive to type of averaging. *supported by ARO grant DAAL03-92-0117.

Inhomogeneous kinetic effects related to intermittent magnetic discontinuities

NASA Astrophysics Data System (ADS)

Greco, A.; Valentini, F.; Servidio, S.; Matthaeus, W. H.

2012-12-01

A connection between kinetic processes and two-dimensional intermittent plasma turbulence is observed using direct numerical simulations of a hybrid Vlasov-Maxwell model, in which the Vlasov equation is solved for protons, while the electrons are described as a massless fluid. During the development of turbulence, the proton distribution functions depart from the typical configuration of local thermodynamic equilibrium, displaying statistically significant non-Maxwellian features. In particular, temperature anisotropy and distortions are concentrated near coherent structures, generated as the result of the turbulent cascade, such as current sheets, which are nonuniformly distributed in space. Here, the partial variance of increments (PVI) method has been employed to identify high magnetic stress regions within a two-dimensional turbulent pattern. A quantitative association between non-Maxwellian features and coherent structures is established.

Local-area simulations of rotating compressible convection and associated mean flows

NASA Technical Reports Server (NTRS)

Hurlburt, Neal E.; Brummell, N. H.; Toomre, Juri

1995-01-01

The dynamics of compressible convection within a curved local segment of a rotating spherical shell are considered in relation to the turbulent redistribution of angular momentum within the solar convection zone. Current supercomputers permit fully turbulent flows to be considered within the restricted geometry of local area models. By considering motions in a curvilinear geometry in which the Coriolos parameters vary with latitude, Rossby waves which couple with the turbulent convection are thought of as being possible. Simulations of rotating convection are presented in such a curved local segment of a spherical shell using a newly developed, sixth-order accurate code based on compact finite differences.

Aerodynamic Design on Unstructured Grids for Turbulent Flows

NASA Technical Reports Server (NTRS)

Anderson, W. Kyle; Bonhaus, Daryl L.

1997-01-01

An aerodynamic design algorithm for turbulent flows using unstructured grids is described. The current approach uses adjoint (costate) variables for obtaining derivatives of the cost function. The solution of the adjoint equations is obtained using an implicit formulation in which the turbulence model is fully coupled with the flow equations when solving for the costate variables. The accuracy of the derivatives is demonstrated by comparison with finite-difference gradients and a few example computations are shown. In addition, a user interface is described which significantly reduces the time required for setting up the design problems. Recommendations on directions of further research into the Navier Stokes design process are made.

Influence of pinches on magnetic reconnection in turbulent space plasmas

NASA Astrophysics Data System (ADS)

Olshevsky, Vyacheslav; Lapenta, Giovanni; Markidis, Stefano; Divin, Andrey

A generally accepted scenario of magnetic reconnection in space plasmas is the breakage of magnetic field lines in X-points. In laboratory, reconnection is widely studied in pinches, current channels embedded into twisted magnetic fields. No model of magnetic reconnection in space plasmas considers both null-points and pinches as peers. We have performed a particle-in-cell simulation of magnetic reconnection in a three-dimensional configuration where null-points are present nitially, and Z-pinches are formed during the simulation. The X-points are relatively stable, and no substantial energy dissipation is associated with them. On contrary, turbulent magnetic reconnection in the pinches causes the magnetic energy to decay at a rate of approximately 1.5 percent per ion gyro period. Current channels and twisted magnetic fields are ubiquitous in turbulent space plasmas, so pinches can be responsible for the observed high magnetic reconnection rates.

Torque fluctuations caused by upstream mean flow and turbulence

NASA Astrophysics Data System (ADS)

Farr, T. D.; Hancock, P. E.

2014-12-01

A series of studies are in progress investigating the effects of turbine-array-wake interactions for a range of atmospheric boundary layer states by means of the EnFlo meteorological wind tunnel. The small, three-blade model wind turbines drive 4-quadrant motor-generators. Only a single turbine in neutral flow is considered here. The motor-generator current can be measured with adequate sensitivity by means of a current sensor allowing the mean and fluctuating torque to be inferred. Spectra of torque fluctuations and streamwise velocity fluctuations ahead of the rotor, between 0.1 and 2 diameters, show that only the large-scale turbulent motions contribute significantly to the torque fluctuations. Time-lagged cross-correlation between upstream velocity and torque fluctuations are largest over the inner part of the blade. They also show the turbulence to be frozen in behaviour over the 2 diameters upstream of the turbine.

Radar Reflectivity in Wingtip-Generated Wake Vortices

NASA Technical Reports Server (NTRS)

Marshall, Robert E.; Mudukutore, Ashok; Wissel, Vicki

1997-01-01

This report documents new predictive models of radar reflectivity, with meter-scale resolution, for aircraft wakes in clear air and fog. The models result from a radar design program to locate and quantify wake vortices from commercial aircraft in support of the NASA Aircraft Vortex Spacing System (AVOSS). The radar reflectivity model for clear air assumes: 1) turbulent eddies in the wake produce small discontinuities in radar refractive index; and 2) these turbulent eddies are in the 'inertial subrange' of turbulence. From these assumptions, the maximum radar frequency for detecting a particular aircraft wake, as well as the refractive index structure constant and radar volume reflectivity in the wake can be obtained from the NASA Terminal Area Simulation System (TASS) output. For fog conditions, an empirical relationship is used to calculate radar reflectivity factor from TASS output of bulk liquid water. Currently, two models exist: 1) Atlas-based on observations of liquid water and radar reflectivity factor in clouds; and 2) de Wolf- specifically tailored to a specific measured dataset (1992 Vandenberg Air Force Base).

Wind-Wave Effects on Vertical Mixing in Chesapeake Bay, USA: comparing observations to second-moment closure predictions.

NASA Astrophysics Data System (ADS)

Fisher, A. W.; Sanford, L. P.; Scully, M. E.

2016-12-01

Coherent wave-driven turbulence generated through wave breaking or nonlinear wave-current interactions, e.g. Langmuir turbulence (LT), can significantly enhance the downward transfer of momentum, kinetic energy, and dissolved gases in the oceanic surface layer. There are few observations of these processes in the estuarine or coastal environments, where wind-driven mixing may co-occur with energetic tidal mixing and strong density stratification. This presents a major challenge for evaluating vertical mixing parameterizations used in modeling estuarine and coastal dynamics. We carried out a large, multi-investigator study of wind-driven estuarine dynamics in the middle reaches of Chesapeake Bay, USA, during 2012-2013. The center of the observational array was an instrumented turbulence tower with both atmospheric and marine turbulence sensors as well as rapidly sampled temperature and conductivity sensors. For this paper, we examined the impacts of surface gravity waves on vertical profiles of turbulent mixing and compared our results to second-moment turbulence closure predictions. Wave and turbulence measurements collected from the vertical array of Acoustic Doppler Velocimeters (ADVs) provided direct estimates of the dominant terms in the TKE budget and the surface wave field. Observed dissipation rates, TKE levels, and turbulent length scales are compared to published scaling relations and used in the calculation of second-moment nonequilibrium stability functions. Results indicate that in the surface layer of the estuary, where elevated dissipation is balanced by vertical divergence in TKE flux, existing nonequilibrium stability functions underpredict observed eddy viscosities. The influences of wave breaking and coherent wave-driven turbulence on modeled and observed stability functions will be discussed further in the context of turbulent length scales, TKE and dissipation profiles, and the depth at which the wave-dominated turbulent transport layer transitions to a turbulent log layer. The influences of fetch-limited wind waves, density stratification, and surface buoyancy fluxes will also be discussed.

Turbulent current drive mechanisms

NASA Astrophysics Data System (ADS)

McDevitt, Christopher J.; Tang, Xian-Zhu; Guo, Zehua

2017-08-01

Mechanisms through which plasma microturbulence can drive a mean electron plasma current are derived. The efficiency through which these turbulent contributions can drive deviations from neoclassical predictions of the electron current profile is computed by employing a linearized Coulomb collision operator. It is found that a non-diffusive contribution to the electron momentum flux as well as an anomalous electron-ion momentum exchange term provide the most efficient means through which turbulence can modify the mean electron current for the cases considered. Such turbulent contributions appear as an effective EMF within Ohm's law and hence provide an ideal means for driving deviations from neoclassical predictions.

Modeling of turbulent chemical reaction

NASA Technical Reports Server (NTRS)

Chen, J.-Y.

1995-01-01

Viewgraphs are presented on modeling turbulent reacting flows, regimes of turbulent combustion, regimes of premixed and regimes of non-premixed turbulent combustion, chemical closure models, flamelet model, conditional moment closure (CMC), NO(x) emissions from turbulent H2 jet flames, probability density function (PDF), departures from chemical equilibrium, mixing models for PDF methods, comparison of predicted and measured H2O mass fractions in turbulent nonpremixed jet flames, experimental evidence of preferential diffusion in turbulent jet flames, and computation of turbulent reacting flows.

Measurements of Turbulence at Two Tidal Energy Sites in Puget Sound, WA

DOE Office of Scientific and Technical Information (OSTI.GOV)

Thomson, Jim; Polagye, Brian; Durgesh, Vibhav

2012-06-05

Field measurements of turbulence are pre- sented from two sites in Puget Sound, WA (USA) that are proposed for electrical power generation using tidal current turbines. Rapidly sampled data from multiple acoustic Doppler instruments are analyzed to obtain statistical mea- sures of fluctuations in both the magnitude and direction of the tidal currents. The resulting turbulence intensities (i.e., the turbulent velocity fluctuations normalized by the harmonic tidal currents) are typically 10% at the hub- heights (i.e., the relevant depth bin) of the proposed turbines. Length and time scales of the turbulence are also analyzed. Large-scale, anisotropic eddies dominate the energymore » spectra, which may be the result of proximity to headlands at each site. At small scales, an isotropic turbulent cascade is observed and used to estimate the dissipation rate of turbulent kinetic energy. Data quality and sampling parameters are discussed, with an emphasis on the removal of Doppler noise from turbulence statistics.« less

Improved Modeling of Finite-Rate Turbulent Combustion Processes in Research Combustors

NASA Technical Reports Server (NTRS)

VanOverbeke, Thomas J.

1998-01-01

The objective of this thesis is to further develop and test a stochastic model of turbulent combustion in recirculating flows. There is a requirement to increase the accuracy of multi-dimensional combustion predictions. As turbulence affects reaction rates, this interaction must be more accurately evaluated. In this work a more physically correct way of handling the interaction of turbulence on combustion is further developed and tested. As turbulence involves randomness, stochastic modeling is used. Averaged values such as temperature and species concentration are found by integrating the probability density function (pdf) over the range of the scalar. The model in this work does not assume the pdf type, but solves for the evolution of the pdf using the Monte Carlo solution technique. The model is further developed by including a more robust reaction solver, by using accurate thermodynamics and by more accurate transport elements. The stochastic method is used with Semi-Implicit Method for Pressure-Linked Equations. The SIMPLE method is used to solve for velocity, pressure, turbulent kinetic energy and dissipation. The pdf solver solves for temperature and species concentration. Thus, the method is partially familiar to combustor engineers. The method is compared to benchmark experimental data and baseline calculations. The baseline method was tested on isothermal flows, evaporating sprays and combusting sprays. Pdf and baseline predictions were performed for three diffusion flames and one premixed flame. The pdf method predicted lower combustion rates than the baseline method in agreement with the data, except for the premixed flame. The baseline and stochastic predictions bounded the experimental data for the premixed flame. The use of a continuous mixing model or relax to mean mixing model had little effect on the prediction of average temperature. Two grids were used in a hydrogen diffusion flame simulation. Grid density did not effect the predictions except for peak temperature and tangential velocity. The hybrid pdf method did take longer and required more memory, but has a theoretical basis to extend to many reaction steps which cannot be said of current turbulent combustion models.

Turbulent Mixing in Gravity Currents with Transverse Shear

NASA Astrophysics Data System (ADS)

White, Brian; Helfrich, Karl; Scotti, Alberto

2010-11-01

A parallel flow with horizontal shear and horizontal density gradient undergoes an intensification of the shear by gravitational tilting and stretching, rapidly breaking down into turbulence. Such flows have the potential for substantial mixing in estuaries and the coastal ocean. We present high-resolution numerical results for the mixing efficiency of these flows, which can be viewed as gravity currents with transverse shear, and contrast them with the well-studied case of stably stratified, homogeneous turbulence (uniform vertical density and velocity gradients). For a sheared gravity current, the buoyancy flux, turbulent Reynolds stress, and dissipation are well out of equilibrium. The total kinetic energy first increases as potential energy is transferred to the gravity current, but rapidly decays once turbulence sets in. Despite the non-equilibrium character, mixing efficiencies are slightly higher but qualitatively similar to homogeneous stratified turbulence. Efficiency decreases in the highly energetic regime where the dissipation rate is large compared with viscosity and stratification, ɛ/(νN^2)>100, further declining as turbulence decays and kinetic energy dissipation dominates the buoyancy flux. In general, the mixing rate, parameterized by a turbulent eddy diffusivity, increases with the strength of the transverse shear.

Scaling laws in magnetized plasma turbulence

DOE Office of Scientific and Technical Information (OSTI.GOV)

Boldyrev, Stanislav

2015-06-28

Interactions of plasma motion with magnetic fields occur in nature and in the laboratory in an impressively broad range of scales, from megaparsecs in astrophysical systems to centimeters in fusion devices. The fact that such an enormous array of phenomena can be effectively studied lies in the existence of fundamental scaling laws in plasma turbulence, which allow one to scale the results of analytic and numerical modeling to the sized of galaxies, velocities of supernovae explosions, or magnetic fields in fusion devices. Magnetohydrodynamics (MHD) provides the simplest framework for describing magnetic plasma turbulence. Recently, a number of new features ofmore » MHD turbulence have been discovered and an impressive array of thought-provoking phenomenological theories have been put forward. However, these theories have conflicting predictions, and the currently available numerical simulations are not able to resolve the contradictions. MHD turbulence exhibits a variety of regimes unusual in regular hydrodynamic turbulence. Depending on the strength of the guide magnetic field it can be dominated by weakly interacting Alfv\\'en waves or strongly interacting wave packets. At small scales such turbulence is locally anisotropic and imbalanced (cross-helical). In a stark contrast with hydrodynamic turbulence, which tends to ``forget'' global constrains and become uniform and isotropic at small scales, MHD turbulence becomes progressively more anisotropic and unbalanced at small scales. Magnetic field plays a fundamental role in turbulent dynamics. Even when such a field is not imposed by external sources, it is self-consistently generated by the magnetic dynamo action. This project aims at a comprehensive study of universal regimes of magnetic plasma turbulence, combining the modern analytic approaches with the state of the art numerical simulations. The proposed study focuses on the three topics: weak MHD turbulence, which is relevant for laboratory devices, the solar wind, solar corona heating, and planetary magnetospheres; strong MHD turbulence, which is relevant for fusion devices, star formation, cosmic rays acceleration, scattering and trapping in galaxies, as well as many aspects of dynamics, distribution and composition of space plasmas, and the process of magnetic dynamo action, which explains the generation and the structure of magnetic fields in turbulent plasmas. The planned work will aim at developing new analytic approaches, conducting new numerical simulations with currently unmatched resolution, and training students in the methods of the modern theory of plasma turbulence. The work will be performed at the University of Wisconsin--Madison.« less

Validation metrics for turbulent plasma transport

DOE PAGES

Holland, C.

2016-06-22

Developing accurate models of plasma dynamics is essential for confident predictive modeling of current and future fusion devices. In modern computer science and engineering, formal verification and validation processes are used to assess model accuracy and establish confidence in the predictive capabilities of a given model. This paper provides an overview of the key guiding principles and best practices for the development of validation metrics, illustrated using examples from investigations of turbulent transport in magnetically confined plasmas. Particular emphasis is given to the importance of uncertainty quantification and its inclusion within the metrics, and the need for utilizing synthetic diagnosticsmore » to enable quantitatively meaningful comparisons between simulation and experiment. As a starting point, the structure of commonly used global transport model metrics and their limitations is reviewed. An alternate approach is then presented, which focuses upon comparisons of predicted local fluxes, fluctuations, and equilibrium gradients against observation. Furthermore, the utility of metrics based upon these comparisons is demonstrated by applying them to gyrokinetic predictions of turbulent transport in a variety of discharges performed on the DIII-D tokamak, as part of a multi-year transport model validation activity.« less

Validation metrics for turbulent plasma transport

DOE Office of Scientific and Technical Information (OSTI.GOV)

Holland, C.

Developing accurate models of plasma dynamics is essential for confident predictive modeling of current and future fusion devices. In modern computer science and engineering, formal verification and validation processes are used to assess model accuracy and establish confidence in the predictive capabilities of a given model. This paper provides an overview of the key guiding principles and best practices for the development of validation metrics, illustrated using examples from investigations of turbulent transport in magnetically confined plasmas. Particular emphasis is given to the importance of uncertainty quantification and its inclusion within the metrics, and the need for utilizing synthetic diagnosticsmore » to enable quantitatively meaningful comparisons between simulation and experiment. As a starting point, the structure of commonly used global transport model metrics and their limitations is reviewed. An alternate approach is then presented, which focuses upon comparisons of predicted local fluxes, fluctuations, and equilibrium gradients against observation. Furthermore, the utility of metrics based upon these comparisons is demonstrated by applying them to gyrokinetic predictions of turbulent transport in a variety of discharges performed on the DIII-D tokamak, as part of a multi-year transport model validation activity.« less

Numerical simulation support to the ESA/THOR mission

NASA Astrophysics Data System (ADS)

Valentini, F.; Servidio, S.; Perri, S.; Perrone, D.; De Marco, R.; Marcucci, M. F.; Daniele, B.; Bruno, R.; Camporeale, E.

2016-12-01

THOR is a spacecraft concept currently undergoing study phase as acandidate for the next ESA medium size mission M4. THOR has been designedto solve the longstanding physical problems of particle heating andenergization in turbulent plasmas. It will provide high resolutionmeasurements of electromagnetic fields and particle distribution functionswith unprecedented resolution, with the aim of exploring the so-calledkinetic scales. We present the numerical simulation framework which is supporting the THOR mission during the study phase. The THOR teamincludes many scientists developing and running different simulation codes(Eulerian-Vlasov, Particle-In-Cell, Gyrokinetics, Two-fluid, MHD, etc.),addressing the physics of plasma turbulence, shocks, magnetic reconnectionand so on.These numerical codes are being used during the study phase, mainly withthe aim of addressing the following points:(i) to simulate the response of real particle instruments on board THOR, byemploying an electrostatic analyser simulator which mimics the response ofthe CSW, IMS and TEA instruments to the particle velocity distributions ofprotons, alpha particle and electrons, as obtained from kinetic numericalsimulations of plasma turbulence.(ii) to compare multi-spacecraft with single-spacecraft configurations inmeasuring current density, by making use of both numerical models ofsynthetic turbulence and real data from MMS spacecraft.(iii) to investigate the validity of the Taylor hypothesis indifferent configurations of plasma turbulence

A multiple-time-scale turbulence model based on variable partitioning of turbulent kinetic energy spectrum

NASA Technical Reports Server (NTRS)

Kim, S.-W.; Chen, C.-P.

1988-01-01

The paper presents a multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method. Consideration is given to a class of turbulent boundary layer flows and of separated and/or swirling elliptic turbulent flows. For the separated and/or swirling turbulent flows, the present turbulence model yielded significantly improved computational results over those obtained with the standard k-epsilon turbulence model.

Using random forests to diagnose aviation turbulence.

PubMed

Williams, John K

Atmospheric turbulence poses a significant hazard to aviation, with severe encounters costing airlines millions of dollars per year in compensation, aircraft damage, and delays due to required post-event inspections and repairs. Moreover, attempts to avoid turbulent airspace cause flight delays and en route deviations that increase air traffic controller workload, disrupt schedules of air crews and passengers and use extra fuel. For these reasons, the Federal Aviation Administration and the National Aeronautics and Space Administration have funded the development of automated turbulence detection, diagnosis and forecasting products. This paper describes a methodology for fusing data from diverse sources and producing a real-time diagnosis of turbulence associated with thunderstorms, a significant cause of weather delays and turbulence encounters that is not well-addressed by current turbulence forecasts. The data fusion algorithm is trained using a retrospective dataset that includes objective turbulence reports from commercial aircraft and collocated predictor data. It is evaluated on an independent test set using several performance metrics including receiver operating characteristic curves, which are used for FAA turbulence product evaluations prior to their deployment. A prototype implementation fuses data from Doppler radar, geostationary satellites, a lightning detection network and a numerical weather prediction model to produce deterministic and probabilistic turbulence assessments suitable for use by air traffic managers, dispatchers and pilots. The algorithm is scheduled to be operationally implemented at the National Weather Service's Aviation Weather Center in 2014.

Anisotropic Stochastic Vortex Structure Method for Simulating Particle Collision in Turbulent Shear Flows

NASA Astrophysics Data System (ADS)

Dizaji, Farzad; Marshall, Jeffrey; Grant, John; Jin, Xing

2017-11-01

Accounting for the effect of subgrid-scale turbulence on interacting particles remains a challenge when using Reynolds-Averaged Navier Stokes (RANS) or Large Eddy Simulation (LES) approaches for simulation of turbulent particulate flows. The standard stochastic Lagrangian method for introducing turbulence into particulate flow computations is not effective when the particles interact via collisions, contact electrification, etc., since this method is not intended to accurately model relative motion between particles. We have recently developed the stochastic vortex structure (SVS) method and demonstrated its use for accurate simulation of particle collision in homogeneous turbulence; the current work presents an extension of the SVS method to turbulent shear flows. The SVS method simulates subgrid-scale turbulence using a set of randomly-positioned, finite-length vortices to generate a synthetic fluctuating velocity field. It has been shown to accurately reproduce the turbulence inertial-range spectrum and the probability density functions for the velocity and acceleration fields. In order to extend SVS to turbulent shear flows, a new inversion method has been developed to orient the vortices in order to generate a specified Reynolds stress field. The extended SVS method is validated in the present study with comparison to direct numerical simulations for a planar turbulent jet flow. This research was supported by the U.S. National Science Foundation under Grant CBET-1332472.

The Kinematics of Turbulent Boundary Layer Structure

NASA Technical Reports Server (NTRS)

Robinson, Stephen Kern

1991-01-01

The long history of research into the internal structure of turbulent boundary layers has not provided a unified picture of the physics responsible for turbulence production and dissipation. The goals of the present research are to: (1) define the current state of boundary layer structure knowledge; and (2) utilize direct numerical simulation results to help close the unresolved issues identified in part A and to unify the fragmented knowledge of various coherent motions into a consistent kinematic model of boundary layer structure. The results of the current study show that all classes of coherent motion in the low Reynolds number turbulent boundary layer may be related to vortical structures, but that no single form of vortex is representative of the wide variety of vortical structures observed. In particular, ejection and sweep motions, as well as entrainment from the free-streem are shown to have strong spatial and temporal relationships with vortical structures. Disturbances of vortex size, location, and intensity show that quasi-streamwise vortices dominate the buffer region, while transverse vortices and vortical arches dominate the wake region. Both types of vortical structure are common in the log region. The interrelationships between the various structures and the population distributions of vortices are combined into a conceptual kinematic model for the boundary layer. Aspects of vortical structure dynamics are also postulated, based on time-sequence animations of the numerically simulated flow.

Description of a Website Resource for Turbulence Modeling Verification and Validation

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Smith, Brian R.; Huang, George P.

2010-01-01

The activities of the Turbulence Model Benchmarking Working Group - which is a subcommittee of the American Institute of Aeronautics and Astronautics (AIAA) Fluid Dynamics Technical Committee - are described. The group s main purpose is to establish a web-based repository for Reynolds-averaged Navier-Stokes turbulence model documentation, including verification and validation cases. This turbulence modeling resource has been established based on feedback from a survey on what is needed to achieve consistency and repeatability in turbulence model implementation and usage, and to document and disseminate information on new turbulence models or improvements to existing models. The various components of the website are described in detail: description of turbulence models, turbulence model readiness rating system, verification cases, validation cases, validation databases, and turbulence manufactured solutions. An outline of future plans of the working group is also provided.

Near-wall turbulence model and its application to fully developed turbulent channel and pipe flows

NASA Technical Reports Server (NTRS)

Kim, S.-W.

1990-01-01

A near-wall turbulence model and its incorporation into a multiple-timescale turbulence model are presented. The near-wall turbulence model is obtained from a k-equation turbulence model and a near-wall analysis. In the method, the equations for the conservation of mass, momentum, and turbulent kinetic energy are integrated up to the wall, and the energy transfer and the dissipation rates inside the near-wall layer are obtained from algebraic equations. Fully developed turbulent channel and pipe flows are solved using a finite element method. The computational results compare favorably with experimental data. It is also shown that the turbulence model can resolve the overshoot phenomena of the turbulent kinetic energy and the dissipation rate in the region very close to the wall.

Criticality and turbulence in a resistive magnetohydrodynamic current sheet

NASA Astrophysics Data System (ADS)

Klimas, Alexander J.; Uritsky, Vadim M.

2017-02-01

Scaling properties of a two-dimensional (2d) plasma physical current-sheet simulation model involving a full set of magnetohydrodynamic (MHD) equations with current-dependent resistivity are investigated. The current sheet supports a spatial magnetic field reversal that is forced through loading of magnetic flux containing plasma at boundaries of the simulation domain. A balance is reached between loading and annihilation of the magnetic flux through reconnection at the current sheet; the transport of magnetic flux from boundaries to current sheet is realized in the form of spatiotemporal avalanches exhibiting power-law statistics of lifetimes and sizes. We identify this dynamics as self-organized criticality (SOC) by verifying an extended set of scaling laws related to both global and local properties of the current sheet (critical susceptibility, finite-size scaling of probability distributions, geometric exponents). The critical exponents obtained from this analysis suggest that the model operates in a slowly driven SOC state similar to the mean-field state of the directed stochastic sandpile model. We also investigate multiscale correlations in the velocity field and find them numerically indistinguishable from certain intermittent turbulence (IT) theories. The results provide clues on physical conditions for SOC behavior in a broad class of plasma systems with propagating instabilities, and suggest that SOC and IT may coexist in driven current sheets which occur ubiquitously in astrophysical and space plasmas.

Criticality and turbulence in a resistive magnetohydrodynamic current sheet.

PubMed

Klimas, Alexander J; Uritsky, Vadim M

2017-02-01

Scaling properties of a two-dimensional (2d) plasma physical current-sheet simulation model involving a full set of magnetohydrodynamic (MHD) equations with current-dependent resistivity are investigated. The current sheet supports a spatial magnetic field reversal that is forced through loading of magnetic flux containing plasma at boundaries of the simulation domain. A balance is reached between loading and annihilation of the magnetic flux through reconnection at the current sheet; the transport of magnetic flux from boundaries to current sheet is realized in the form of spatiotemporal avalanches exhibiting power-law statistics of lifetimes and sizes. We identify this dynamics as self-organized criticality (SOC) by verifying an extended set of scaling laws related to both global and local properties of the current sheet (critical susceptibility, finite-size scaling of probability distributions, geometric exponents). The critical exponents obtained from this analysis suggest that the model operates in a slowly driven SOC state similar to the mean-field state of the directed stochastic sandpile model. We also investigate multiscale correlations in the velocity field and find them numerically indistinguishable from certain intermittent turbulence (IT) theories. The results provide clues on physical conditions for SOC behavior in a broad class of plasma systems with propagating instabilities, and suggest that SOC and IT may coexist in driven current sheets which occur ubiquitously in astrophysical and space plasmas.

Prediction of Transitional Flows in the Low Pressure Turbine

NASA Technical Reports Server (NTRS)

Huang, George; Xiong, Guohua

1998-01-01

Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Direct numerical simulation and reduced-order modeling of the sound-induced flow through a cavity-backed circular under a turbulent boundary layer

NASA Astrophysics Data System (ADS)

Zhang, Qi; Bodony, Daniel

2014-11-01

Commercial jet aircraft generate undesirable noise from several sources, with the engines being the most dominant sources at take-off and major contributors at all other stages of flight. Acoustic liners, which are perforated sheets of metal or composite mounted within the engine, have been an effective means of reducing internal engine noise from the fan, compressor, combustor, and turbine but their performance suffers when subjected to a turbulent grazing flow or to high-amplitude incident sound due to poorly understood interactions between the liner orifices and the exterior flow. Through the use of direct numerical simulations, the flow-orifice interaction is examined numerically, quantified, and modeled over a range of conditions that includes current and envisioned uses of acoustic liners and with detail that exceeds experimental capabilities. A new time-domain model of acoustic liners is developed that extends currently-available reduced-order models to more complex flow conditions but is still efficient for use at the design stage.

Turbulent current drive mechanisms

DOE Office of Scientific and Technical Information (OSTI.GOV)

McDevitt, Christopher J.; Tang, Xian-Zhu; Guo, Zehua

Mechanisms through which plasma microturbulence can drive a mean electron plasma current are derived. The efficiency through which these turbulent contributions can drive deviations from neoclassical predictions of the electron current profile is computed by employing a linearized Coulomb collision operator. It is found that a non-diffusive contribution to the electron momentum flux as well as an anomalous electron-ion momentum exchange term provide the most efficient means through which turbulence can modify the mean electron current for the cases considered. Such turbulent contributions appear as an effective EMF within Ohm’s law, and hence provide an ideal means for driving deviationsmore » from neoclassical predictions.« less

Turbulent current drive mechanisms

DOE PAGES

McDevitt, Christopher J.; Tang, Xian-Zhu; Guo, Zehua

2017-07-01

Mechanisms through which plasma microturbulence can drive a mean electron plasma current are derived. The efficiency through which these turbulent contributions can drive deviations from neoclassical predictions of the electron current profile is computed by employing a linearized Coulomb collision operator. It is found that a non-diffusive contribution to the electron momentum flux as well as an anomalous electron-ion momentum exchange term provide the most efficient means through which turbulence can modify the mean electron current for the cases considered. Such turbulent contributions appear as an effective EMF within Ohm’s law, and hence provide an ideal means for driving deviationsmore » from neoclassical predictions.« less

Modeling the effect of small-scale magnetic turbulence on the X-ray properties of Pulsar Wind Nebulae

NASA Astrophysics Data System (ADS)

Bucciantini, N.; Bandiera, R.; Olmi, B.; Del Zanna, L.

2017-10-01

Pulsar Wind Nebulae (PWNe) constitute an ideal astrophysical environment to test our current understanding of relativistic plasma processes. It is well known that magnetic fields play a crucial role in their dynamics and emission properties. At present, one of the main issues concerns the level of magnetic turbulence present in these systems, which in the absence of space resolved X-ray polarization measures cannot be directly constrained. In this work, we investigate, for the first time using simulated synchrotron maps, the effect of a small-scale fluctuating component of the magnetic field on the emission properties in X-ray. We illustrate how to include the effects of a turbulent component in standard emission models for PWNe and which consequences are expected in terms of net emissivity and depolarization, showing that the X-ray surface brightness maps can provide already some rough constraints. We then apply our analysis to the Crab and Vela nebulae and by comparing our model with Chandra and Vela data, we found that the typical energies in the turbulent component of the magnetic field are about 1.5-3 times the one in the ordered field.

CASL VMA Milestone Report FY16 (L3:VMA.VUQ.P13.08): Westinghouse Mixing with STAR-CCM+

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gilkey, Lindsay Noelle

2016-09-30

STAR-CCM+ (STAR) is a high-resolution computational fluid dynamics (CFD) code developed by CD-adapco. STAR includes validated physics models and a full suite of turbulence models including ones from the k-ε and k-ω families. STAR is currently being developed to be able to do two phase flows, but the current focus of the software is single phase flow. STAR can use imported meshes or use the built in meshing software to create computation domains for CFD. Since the solvers generally require a fine mesh for good computational results, the meshes used with STAR tend to number in the millions of cells,more » with that number growing with simulation and geometry complexity. The time required to model the flow of a full 5x5 Mixing Vane Grid Assembly (5x5MVG) in the current STAR configuration is on the order of hours, and can be very computationally expensive. COBRA-TF (CTF) is a low-resolution subchannel code that can be trained using high fidelity data from STAR. CTF does not have turbulence models and instead uses a turbulent mixing coefficient β. With a properly calibrated β, CTF can be used a low-computational cost alternative to expensive full CFD calculations performed with STAR. During the Hi2Lo work with CTF and STAR, STAR-CCM+ will be used to calibrate β and to provide high-resolution results that can be used in the place of and in addition to experimental results to reduce the uncertainty in the CTF results.« less

Numerical Model of Turbulence, Sediment Transport, and Sediment Cover in a Large Canyon-Bound River

NASA Astrophysics Data System (ADS)

Alvarez, L. V.; Schmeeckle, M. W.

2013-12-01

The Colorado River in Grand Canyon is confined by bedrock and coarse-grained sediments. Finer grain sizes are supply limited, and sandbars primarily occur in lateral separation eddies downstream of coarse-grained tributary debris fans. These sandbars are important resources for native fish, recreational boaters, and as a source of aeolian transport preventing the erosion of archaeological resources by gully extension. Relatively accurate prediction of deposition and, especially, erosion of these sandbar beaches has proven difficult using two- and three-dimensional, time-averaged morphodynamic models. We present a parallelized, three-dimensional, turbulence-resolving model using the Detached-Eddy Simulation (DES) technique. DES is a hybrid large eddy simulation (LES) and Reynolds-averaged Navier Stokes (RANS). RANS is applied to the near-bed grid cells, where grid resolution is not sufficient to fully resolve wall turbulence. LES is applied further from the bed and banks. We utilize the Spalart-Allmaras one equation turbulence closure with a rough wall extension. The model resolves large-scale turbulence using DES and simultaneously integrates the suspended sediment advection-diffusion equation. The Smith and McLean suspended sediment boundary condition is used to calculate the upward and downward settling of sediment fluxes in the grid cells attached to the bed. The model calculates the entrainment of five grain sizes at every time step using a mixing layer model. Where the mixing layer depth becomes zero, the net entrainment is zero or negative. As such, the model is able to predict the exposure and burial of bedrock and coarse-grained surfaces by fine-grained sediments. A separate program was written to automatically construct the computational domain between the water surface and a triangulated surface of a digital elevation model of the given river reach. Model results compare favorably with ADCP measurements of flow taken on the Colorado River in Grand Canyon during the High Flow Experiment (HFE) of 2008. The model accurately reproduces the size and position of the major recirculation currents, and the error in velocity magnitude was found to be less than 17% or 0.22 m/s absolute error. The mean deviation of the direction of velocity with respect to the measured velocity was found to be 20 degrees. Large-scale turbulence structures with vorticity predominantly in the vertical direction are produced at the shear layer between the main channel and the separation zone. However, these structures rapidly become three-dimensional with no preferred orientation of vorticity. Surprisingly, cross-stream velocities, into the main recirculation zone just upstream of the point of reattachment and out of the main recirculation region just downstream of the point of separation, are highest near the bed. Lateral separation eddies are more efficient at storing and exporting sediment than previously modeled. The input of sediment to the eddy recirculation zone occurs near the reattachment zone and is relatively continuous in time. While, the export of sediment to the main channel by the return current occurs in pulses. Pulsation of the strength of the return current becomes a key factor to determine the rates of erosion and deposition in the main recirculation zone.

Predictions of spray combustion interactions

NASA Technical Reports Server (NTRS)

Shuen, J. S.; Solomon, A. S. P.; Faeth, G. M.

1984-01-01

Mean and fluctuating phase velocities; mean particle mass flux; particle size; and mean gas-phase Reynolds stress, composition and temperature were measured in stationary, turbulent, axisymmetric, and flows which conform to the boundary layer approximations while having well-defined initial and boundary conditions in dilute particle-laden jets, nonevaporating sprays, and evaporating sprays injected into a still air environment. Three models of the processes, typical of current practice, were evaluated. The local homogeneous flow and deterministic separated flow models did not provide very satisfactory predictions over the present data base. In contrast, the stochastic separated flow model generally provided good predictions and appears to be an attractive approach for treating nonlinear interphase transport processes in turbulent flows containing particles (drops).

Large Eddy Simulations of a Bottom Boundary Layer Under a Shallow Geostrophic Front

NASA Astrophysics Data System (ADS)

Bateman, S. P.; Simeonov, J.; Calantoni, J.

2017-12-01

The unstratified surf zone and the stratified shelf waters are often separated by dynamic fronts that can strongly impact the character of the Ekman bottom boundary layer. Here, we use large eddy simulations to study the turbulent bottom boundary layer associated with a geostrophic current on a stratified shelf of uniform depth. The simulations are initialized with a spatially uniform vertical shear that is in geostrophic balance with a pressure gradient due to a linear horizontal temperature variation. Superposed on the temperature front is a stable vertical temperature gradient. As turbulence develops near the bottom, the turbulence-induced mixing gradually erodes the initial uniform temperature stratification and a well-mixed layer grows in height until the turbulence becomes fully developed. The simulations provide the spatial distribution of the turbulent dissipation and the Reynolds stresses in the fully developed boundary layer. We vary the initial linear stratification and investigate its effect on the height of the bottom boundary layer and the turbulence statistics. The results are compared to previous models and simulations of stratified bottom Ekman layers.

Large Eddy Simulation of Turbulent Flow in a Ribbed Pipe

NASA Astrophysics Data System (ADS)

Kang, Changwoo; Yang, Kyung-Soo

2011-11-01

Turbulent flow in a pipe with periodically wall-mounted ribs has been investigated by large eddy simulation with a dynamic subgrid-scale model. The value of Re considered is 98,000, based on hydraulic diameter and mean bulk velocity. An immersed boundary method was employed to implement the ribs in the computational domain. The spacing of the ribs is the key parameter to produce the d-type, intermediate and k-type roughness flows. The mean velocity profiles and turbulent intensities obtained from the present LES are in good agreement with the experimental measurements currently available. Turbulence statistics, including budgets of the Reynolds stresses, were computed, and analyzed to elucidate turbulence structures, especially around the ribs. In particular, effects of the ribs are identified by comparing the turbulence structures with those of smooth pipe flow. The present investigation is relevant to the erosion/corrosion that often occurs around a protruding roughness in a pipe system. This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2010-0008457).

In vitro evaluation of flow patterns and turbulent kinetic energy in trans-catheter aortic valve prostheses.

PubMed

Giese, Daniel; Weiss, Kilian; Baeßler, Bettina; Madershahian, Navid; Choi, Yeong-Hoon; Maintz, David; Bunck, Alexander C

2018-02-01

The objective of the current work was to evaluate flow and turbulent kinetic energy in different transcatheter aortic valve implants using highly undersampled time-resolved multi-point 3-directional phase-contrast measurements (4D Flow MRI) in an in vitro setup. A pulsatile flow setup was used with a compliant tubing mimicking a stiff left ventricular outflow tract and ascending aorta. Five different implants were measured using a highly undersampled multi-point 4D Flow MRI sequence. Velocities and turbulent kinetic energy values were analysed and compared. Strong variations of turbulent kinetic energy distributions between the valves were observed. Maximum turbulent kinetic energy values ranged from 100 to over 500 J/m 3 while through-plane velocities were similar between all valves. Highly accelerated 4D Flow MRI for the measurement of velocities and turbulent kinetic energy values allowed for the assessment of hemodynamic parameters in five different implant models. The presented setup, measurement protocol and analysis methods provides an efficient approach to compare different valve implants and could aid future novel valve designs.

Analogy between the Navier-Stokes equations and Maxwell's equations: Application to turbulence

NASA Astrophysics Data System (ADS)

Marmanis, Haralambos

1998-06-01

A new theory of turbulence is initiated, based on the analogy between electromagnetism and turbulent hydrodynamics, for the purpose of describing the dynamical behavior of averaged flow quantities in incompressible fluid flows of high Reynolds numbers. The starting point is the recognition that the vorticity (w=∇×u) and the Lamb vector (l=w×u) should be taken as the kernel of a dynamical theory of turbulence. The governing equations for these fields can be obtained by the Navier-Stokes equations, which underlie the whole evolution. Then whatever parts are not explicitly expressed as a function of w or l only are gathered and treated as source terms. This is done by introducing the concepts of turbulent charge and turbulent current. Thus we are led to a closed set of linear equations for the averaged field quantities. The premise is that the earlier introduced sources will be apt for modeling, in the sense that their distribution will depend only on the geometry and the total energetics of the flow. The dynamics described in the preceding manner is what we call the metafluid dynamics.

Prediction of Turbulence-Generated Noise in Unheated Jets. Part 2; JeNo Users' Manual (Version 1.0)

NASA Technical Reports Server (NTRS)

Khavaran, Abbas; Wolter, John D.; Koch, L. Danielle

2009-01-01

JeNo (Version 1.0) is a Fortran90 computer code that calculates the far-field sound spectral density produced by axisymmetric, unheated jets at a user specified observer location and frequency range. The user must provide a structured computational grid and a mean flow solution from a Reynolds-Averaged Navier Stokes (RANS) code as input. Turbulence kinetic energy and its dissipation rate from a k-epsilon or k-omega turbulence model must also be provided. JeNo is a research code, and as such, its development is ongoing. The goal is to create a code that is able to accurately compute far-field sound pressure levels for jets at all observer angles and all operating conditions. In order to achieve this goal, current theories must be combined with the best practices in numerical modeling, all of which must be validated by experiment. Since the acoustic predictions from JeNo are based on the mean flow solutions from a RANS code, quality predictions depend on accurate aerodynamic input.This is why acoustic source modeling, turbulence modeling, together with the development of advanced measurement systems are the leading areas of research in jet noise research at NASA Glenn Research Center.

Modeling Compressibility Effects in High-Speed Turbulent Flows

NASA Technical Reports Server (NTRS)

Sarkar, S.

2004-01-01

Man has strived to make objects fly faster, first from subsonic to supersonic and then to hypersonic speeds. Spacecraft and high-speed missiles routinely fly at hypersonic Mach numbers, M greater than 5. In defense applications, aircraft reach hypersonic speeds at high altitude and so may civilian aircraft in the future. Hypersonic flight, while presenting opportunities, has formidable challenges that have spurred vigorous research and development, mainly by NASA and the Air Force in the USA. Although NASP, the premier hypersonic concept of the eighties and early nineties, did not lead to flight demonstration, much basic research and technology development was possible. There is renewed interest in supersonic and hypersonic flight with the HyTech program of the Air Force and the Hyper-X program at NASA being examples of current thrusts in the field. At high-subsonic to supersonic speeds, fluid compressibility becomes increasingly important in the turbulent boundary layers and shear layers associated with the flow around aerospace vehicles. Changes in thermodynamic variables: density, temperature and pressure, interact strongly with the underlying vortical, turbulent flow. The ensuing changes to the flow may be qualitative such as shocks which have no incompressible counterpart, or quantitative such as the reduction of skin friction with Mach number, large heat transfer rates due to viscous heating, and the dramatic reduction of fuel/oxidant mixing at high convective Mach number. The peculiarities of compressible turbulence, so-called compressibility effects, have been reviewed by Fernholz and Finley. Predictions of aerodynamic performance in high-speed applications require accurate computational modeling of these "compressibility effects" on turbulence. During the course of the project we have made fundamental advances in modeling the pressure-strain correlation and developed a code to evaluate alternate turbulence models in the compressible shear layer.

Mass and Momentum Turbulent Transport Experiments with Confined Coaxial Jets

NASA Technical Reports Server (NTRS)

Johnson, B. V.; Bennett, J. C.

1981-01-01

Downstream mixing of coaxial jets discharging in an expanded duct was studied to obtain data for the evaluation and improvement of turbulent transport models currently used in a variety of computational procedures throughout the propulsion community for combustor flow modeling. Flow visualization studies showed four major shear regions occurring; a wake region immediately downstream of the inlet jet inlet duct; a shear region further downstream between the inner and annular jets; a recirculation zone; and a reattachment zone. A combination of turbulent momentum transport rate and two velocity component data were obtained from simultaneous measurements with a two color laser velocimeter (LV) system. Axial, radial and azimuthal velocities and turbulent momentum transport rate measurements in the r-z and r-theta planes were used to determine the mean value, second central moment (or rms fluctuation from mean), skewness and kurtosis for each data set probability density function (p.d.f.). A combination of turbulent mass transport rate, concentration and velocity data were obtained system. Velocity and mass transport in all three directions as well as concentration distributions were used to obtain the mean, second central moments, skewness and kurtosis for each p.d.f. These LV/LIF measurements also exposed the existence of a large region of countergradient turbulent axial mass transport in the region where the annular jet fluid was accelerating the inner jet fluid.

Towards predictive simulations of soot formation: from surrogate to turbulence

DOE Office of Scientific and Technical Information (OSTI.GOV)

Blanquart, Guillaume

The combustion of transportation fuels leads to the formation of several kinds of pollutants, among which are soot particles. These particles, also formed during coal combustion and in fires, are the source of several health problems and environmental issues. Unfortunately, our current understanding of the chemical and physical phenomena leading to the formation of soot particles remains incomplete, and as a result, the predictive capability of our numerical tools is lacking. The objective of the work was to reduce the gap in the present understanding and modeling of soot formation both in laminar and turbulent flames. The effort spanned severalmore » length scales from the molecular level to large scale turbulent transport.« less

Turbulent reconnection and its implications

PubMed Central

Lazarian, A.; Eyink, G.; Vishniac, E.; Kowal, G.

2015-01-01

Magnetic reconnection is a process of magnetic field topology change, which is one of the most fundamental processes happening in magnetized plasmas. In most astrophysical environments, the Reynolds numbers corresponding to plasma flows are large and therefore the transition to turbulence is inevitable. This turbulence, which can be pre-existing or driven by magnetic reconnection itself, must be taken into account for any theory of magnetic reconnection that attempts to describe the process in the aforementioned environments. This necessity is obvious as three-dimensional high-resolution numerical simulations show the transition to the turbulence state of initially laminar reconnecting magnetic fields. We discuss ideas of how turbulence can modify reconnection with the focus on the Lazarian & Vishniac (Lazarian & Vishniac 1999 Astrophys. J. 517, 700–718 ()) reconnection model. We present numerical evidence supporting the model and demonstrate that it is closely connected to the experimentally proven concept of Richardson dispersion/diffusion as well as to more recent advances in understanding of the Lagrangian dynamics of magnetized fluids. We point out that the generalized Ohm's law that accounts for turbulent motion predicts the subdominance of the microphysical plasma effects for reconnection for realistically turbulent media. We show that one of the most dramatic consequences of turbulence is the violation of the generally accepted notion of magnetic flux freezing. This notion is a cornerstone of most theories dealing with magnetized plasmas, and therefore its change induces fundamental shifts in accepted paradigms, for instance, turbulent reconnection entails reconnection diffusion process that is essential for understanding star formation. We argue that at sufficiently high Reynolds numbers the process of tearing reconnection should transfer to turbulent reconnection. We discuss flares that are predicted by turbulent reconnection and relate this process to solar flares and γ-ray bursts. With reference to experiments, we analyse solar observations in situ as measurements in the solar wind or heliospheric current sheet and show the correspondence of data with turbulent reconnection predictions. Finally, we discuss first-order Fermi acceleration of particles that is a natural consequence of the turbulent reconnection. PMID:25848076

Turbulent kinetic energy of the ocean winds over the Kuroshio Extension from QuikSCAT winds (1999-2009)

NASA Astrophysics Data System (ADS)

Yu, Kai; Dong, Changming; King, Gregory P.

2017-06-01

We investigate mesoscale turbulence (10-1000 km) in the ocean winds over the Kuroshio Extension (28°N-40°N, 140°E-180°E) using the QuikSCAT data set (November 1999 to October 2009). We calculate the second (Djj) and third-order structure functions (Djjj) and the spatial variance (Vj) as a function of scale r (j=L,T denotes, respectively, the longitudinal (divergent) and transverse (vortical) component). The most interesting results of the analysis follow. Although both Vj>(r>) and Djj>(r>) measure the turbulent kinetic energy (TKE), we find that Vj>(r>) is the more robust measure. The spatial variance density (dVj/dr) has a broad peak near 450 km (close to the midlatitude Rossby radius of deformation). On interannual time scales, TKE correlates well with the El Niño 3.4 index. According to turbulence theory, the kinetic energy cascades downscale (upscale) if DLLL>(r>) (also skewness SL=DLLL/DLL3/2) is negative (positive). Our results for the Kuroshio Extension are consistent with a downscale cascade (indicating convergence dominates). Furthermore, classical turbulence theory predicts that SL=-0.3 and independent of r; however, we find SL varies strongly with r, from -4 at small scales to -0.3 at large scales. This nonclassical behavior implies strong-scale interaction, which we attribute to the rapid, and sometimes explosive, growth of storms in the region through baroclinic instability. Finally, we find that ST (a measure of cyclonic/anticyclonic asymmetry) is positive (cyclonic) and also varies strongly with r, from 4 at small scales to 0.5 at large scales. New turbulence models are needed to explain these results, and that will benefit Weather Prediction and climate modeling.