Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 199 records · Page 11

Clear air turbulence forecasting techniques

A method to improve clear air turbulence (CAT) forecasting by more effectively using the currently operational rawinsonde (RW) system is discussed. The method is called the Diagnostic Richardson Number Tendency (DRT) technique. The technique does not attempt to use the RW as a direct detector of the turbulent motion or even of the CAT mechanism structure but rather senses the synoptic scale centers of action which provide the energy to the CAT mechanism at the mesoscale level. The DRT algorithm is deterministic rather than statistical in nature, using the hydrodynamic equations (equations of motion) relevant to the synoptic scale. However, interpretation, by necessity, is probabilistic. What is most important with respect to its operational implementation is that this method uses the same input data as currently used by the operational National Meteorological Center prognostic models.

Keller, J. L.↗

Methods of separation of variables in turbulence theory

Two schemes of closing turbulent moment equations are proposed both of which make double correlation equations separated into single-point equations. The first is based on neglected triple correlation, leading to an equation differing from small perturbed gasdynamic equations where the separation constant appears as the frequency. Grid-produced turbulence is described in this light as time-independent, cylindrically-isotropic turbulence. Application to wall turbulence guided by a new asymptotic method for the Orr-Sommerfeld equation reveals a neutrally stable mode of essentially three dimensional nature. The second closure scheme is based on an assumption of identity of the separated variables through which triple and quadruple correlations are formed. The resulting equation adds, to its equivalent of the first scheme, an integral of nonlinear convolution in the frequency describing a role due to triple correlation of direct energy-cascading.

Tsuge, S.↗

Numerical investigation of an internal layer in turbulent flow over a curved hill

The development of an internal layer in a turbulent boundary layer flow over a curved hill is investigated numerically. The turbulence field of the boundary layer flow over the curved hill is compared with that of a turbulent flow over a symmetric airfoil (which has the same geometry as the curved hill except that the leading and trailing edge plates were removed) to study the influence of the strongly curved surface on the turbulence field. The turbulent flow equations are solved by a control-volume based finite difference method. The turbulence is described by a multiple-time-scale turbulence model supplemented with a near-wall turbulence model. Computational results for the mean flow field (pressure distributions on the walls, wall shearing stresses and mean velocity profiles), the turbulence structure (Reynolds stress and turbulent kinetic energy profiles), and the integral parameters (displacement and momentum thicknesses) compared favorably with the measured data. Computational results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. Development of the internal layer was more obviously observed in the Reynolds stress profiles and in the turbulent kinetic energy profiles than in the mean velocity profiles. In this regard, the internal layers is significantly different from wall-bounded simple shear layers in which the mean velocity profile characterizes the boundary layer most distinguishably. Development of such an internal layer, characterized by an intense turbulence field, is attributed to the enormous mean flow strain rate caused by the streamline curvature and the strong pressure gradient. In the turbulent flow over the curved hill, the internal layer begin to form near the forward corner of the hill, merges with the external boundary layer, and develops into a new fully turbulent boundary layer as the fluid flows in the downstream direction. For the flow over the symmetric airfoil, the boundary layer began to form from almost the same location as that of the curved hill, grew in its strength, and formed a fully turbulent boundary layer from mid-part of the airfoil and in the downstream region. Computational results also show that the detailed turbulence structure in the region very close to the wall of the curved hill is almost the same as that of the airfoil in most of the curved regions except near the leading edge. Thus the internal layer of the curved hill and the boundary layer of the airfoil were also almost the same. Development of the wall shearing stress and separation of the boundary layer at the rear end of the curved hill mostly depends on the internal layer and is only slightly influenced by the external boundary layer flow.

Kim, S-W.↗

Scaling and modeling of turbulent suspension flows

Scaling factors determining various aspects of particle-fluid interactions and the development of physical models to predict gas-solid turbulent suspension flow fields are discussed based on two-fluid, continua formulation. The modes of particle-fluid interactions are discussed based on the length and time scale ratio, which depends on the properties of the particles and the characteristics of the flow turbulence. For particle size smaller than or comparable with the Kolmogorov length scale and concentration low enough for neglecting direct particle-particle interaction, scaling rules can be established in various parameter ranges. The various particle-fluid interactions give rise to additional mechanisms which affect the fluid mechanics of the conveying gas phase. These extra mechanisms are incorporated into a turbulence modeling method based on the scaling rules. A multiple-scale two-phase turbulence model is developed, which gives reasonable predictions for dilute suspension flow. Much work still needs to be done to account for the poly-dispersed effects and the extension to dense suspension flows.

Chen, C. P.↗

Turbulence Modeling: A NASA Perspective

This paper presents turbulence modeling from NASA's perspective. The topics include: 1) Hierarchy of Solution Methods; 2) Turbulence Modeling Focus; 3) Linear Eddy Viscosity Models; and 4) Nonlinear Eddy Viscosity Algebraic Stress Models.

Gatski, T. B.↗

Propagation of sound in turbulent media

Perturbation methods commonly used to study the propagation of acoustic waves in turbulent media are reviewed. Emphasis is on those techniques which are applicable to problems involving long-range propagation in the atmosphere and ocean. Characteristic features of the various methods are illustrated by applying them to particular problems. It is shown that conventional perturbation techniques, such as the Born approximation, yield solutions which contain secular terms, and which therefore have a relatively limited range of validity. In contrast, it is found that solutions obtained with the aid of the Rytov method or the smoothing method do not contain secular terms, and consequently have a much greater range of validity.

Wenzel, A. R.↗

Modeling and simulation of a turbulent far wake

Work continued on two projects which had been started during previous years. Both projects involve calculations of the subsonic, turbulent far wake of a two-dimensional object at a Reynolds number of 1000 (based on wake momentum thickness). This flow was used as a test case for direct comparison of various turbulence models and a direct numerical simulations (DNS) of this flow were undertaken. In the turbulence model comparison studies, for any particular model tested, a unique self-similar solution was obtained far enough downstream, regardless of inlet conditions. Furthermore, different turbulence models led to different far-wake equilibrium solutions. No turbulence model could correctly predict all features of the turbulent far wake. For example, the spreading rate and turbulent shear stresses were underpredicted by all the standard models (both two-equation and full Reynolds stress models). In cases where a more correct spreading rate was achieved, it was at the expense of the turbulent kinetic energy, which was overpredicted. In general, the Algebraic Dissipation Rate Model of Gatski and Speziale, 1992, when added to any of the standard models, improved the results dramatically. Also, full Reynolds stress closure models did a much better job at predicting the shapes of both the mean and turbulence profiles, but the spreading rate was not significantly improved over that predicted by the simpler two-equation models. There are two main conclusions from these studies: First, in a comparison such as this, it is not enough to compare just one parameter, like the spreading rate. A good prediction for one parameter does not necessarily imply good predictions for all parameters in a flow. Second, since no turbulence model could correctly predict the turbulent far wake, much of the important physics of turbulent free shear flows is apparently lost by the assumptions inherent in today's methods of turbulence modeling; turbulence models must be improved. Direct simulations of this flow were begun last year in order to provide a data base through which some of the deficiencies of the existing turbulence models could be identified. Quantities such as the pressure-strain correlation, turbulent diffusion, and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. During this summer, diagnostics codes were written which will calculate the parameters mentioned above, along with other single-point and multi-point statistics. The DNS calculations are still in progress at the time of this writing. When these calculations are complete, the diagnostics codes will be applied so that the results can aid turbulence modelers. In addition, the results will show whether or not there exists a universal equilibrium turbulent far wake, independent of initial conditions.

Cimbala, John M.↗

Helicopter response to atmospheric turbulence

A new time-domain method for simulating cyclostationary turbulence as seen by a translating and rotating blade element has recently been developed for the case of one-dimensional spectral distribution. This paper extends the simulation method to the cases of two- and three-dimensional spectral distributions and presents validation results for the two-dimensional case. The statistics of an isolated rigid blade flapping response to turbulence are computed using a two-dimensional spectral representation of the von Karman turbulence model, and the results are compared with those obtained using the conventional space-fixed turbulence analysis. The new turbulence simulation method is used for predicting the Black Hawk helicopter response to atmospheric turbulence.

Riaz, J.↗

Optimization of active control systems to suppress flutter and minimize turbulence response

A method for optimizing an active control system which will suppress flutter and minimize the response of a lifting surface to atmospheric turbulence is presented. A mathematical search method is developed which will find a control law which will cause an active control system to flutter at a specified freestream velocity, air density, and Mach number. With the flutter velocity of the system held constant, the control law is then modified in such a way that the peak output power spectral density function of the angular response of a lifting surface (as a result of atmospheric turbulence) is minimized for a specified flight velocity which is less than the flutter velocity. The von Karman generalized power spectrum for the transverse components of turbulence is used in an example problem to increase the flutter velocity and minimize the turbulence response of a simplified delta-wing model which has leading and trailing edge control surfaces.

Rudisill, C. S.↗

Electrochemical Control for Corrosion in Molten Chlorides During CSP Plant Operation

The Liquid Pathway of the Concentrating Solar Power Generation 3 (CSP Gen3) program proposed low-cost molten chloride salt for energy storage. However, online corrosion control was identified as a remining major risk of the Liquid Pathway approach. This project addressed that risk. Electrochemical solutions for corrosion mitigation during CSP plant operation were investigated and their feasibility and scalability were evaluated. The leading cause of corrosion in molten chloride salt systems was identified as corrosive impurities that form within the salt upon exposure to trace amounts of air and moisture. Leveraging electrochemistry, reduction/oxidation reactions can be employed to remove these corrosive impurities. In Phase 1 of this project, a bench-scale batch electrochemical reactor was designed, fabricated, and used to assess the kinetics and thermodynamics of electrochemical salt purification. In Phase 2, a laboratory-scale flow reactor was designed, fabricated, and used to assess the efficacy of the electrochemical method under flowing conditions. Results show that under proposed operating conditions for the Liquid Pathway Gen3 Pilot Plant, the electrochemical method is significantly more effective at removing impurities than alternative chemical and thermal methods, and that the electrochemical method produces less harmful byproducts. A key advance made in the course of this project was the development of a 2-electrode method for electrochemical purification that is more scalable than previously developed 3 electrode methods. This novel method is based on Magnesium (Mg) electrowinning. A provisional patent based on this invention has been submitted (USPTO Application No. 63/480,355). Additional key advances made during this project include assessment of the effect of dissimilar alloys on corrosion, kinetic and thermodynamic evaluation of thermolysis reactions of impurities within the molten salt, characterization of byproducts of purification reactions, and generation of IP focused on isolating value-added products using molten salt-based electrochemistry that could be deployed to valorize the process (USPTO Application No. 63/478,806). Ultimately, this project represented a step toward feasibility of Liquid Pathway Gen3 CSP. The method developed under this project could significantly reduce capital expenses and operating costs and increase plant profitability by enabling use of less expensive alloys, decreasing maintenance, and increasing plant longevity. Key focus areas for follow-on work have been identified as 1) evaluation of the efficacy of the electrochemical method under turbulent conditions in a larger flow system, such as the FASTR loop, 2) development of methods for removal of purification byproducts, 3) modeling pilot and industrial scale performance of electrochemical salt purification during plant operation and 4) further assessment of the effect of impurities on salt vapor phase.

14 SOLAR ENERGY↗

Characteristics of a trailing flap flow with small separation

Detailed measurements of pressure and velocity characteristics are reported for the flow on and downstream of a trailing flap, whose 16 degree angle of incidence resulted in boundary layer separation, a small region of recirculating flow, and a curved downstream wake. Emphasis is placed on the region of recirculating flow and on the downstream wake. The characteristics were quantified by the use of pressure, hot-wire, and flying hot-wire probes. The flow characteristics are compared to those of a similar flow but with a larger region of recirculation reported by Thompson and Whitelaw (1985). The relative importance of the terms in the transport equations for mean momentum and turbulence energy are quantified and the implications of the results for viscous-inviscid interaction methods and turbulence modeling are discussed.

Adair, D.↗

Analysis and testing for rotordynamic coefficients of turbulent annular seals with different, directionally-homogeneous surface-roughness treatment for rotor and stator elements

A combined analytical-computational method has been developed to calculate the transient pressure field and rotordynamic coefficients for high-pressure annular seals. The solution procedure applies to constant-clearance or the convergent tapered geometries of multistage centrifugal pumps which may have different surface roughness conditions on the stator and rotor seal elements. In experimental calculations with the method, the turbulent equations of Hirs (1973) are modified slightly to account for different surface roughness conditions, and a perturbation analysis is employed to develop zeroth and first order perturbation equations. Zeroth equations are also used to define both the leakage and circumferential flow due to shear stresses around the stator and rotor surfaces. The solution to analytical equations for four different surface roughnesses confirm the predicted net damping for the seals. A round-holed stator pattern yielded the highest net damping and lowest leakage of all the seals tested.

Childs, D. W.↗

Evaluation of four subcritical response methods for on-line prediction flutter onset in wind-tunnel tests

The methods were evaluated for use in tests where the flutter model is excited solely by airstream turbulence. The methods were: randomdec, power-spectral-density, peak-hold, and cross-spectrum. The test procedure was to maintain a constant Mach number (M) and increase the dynamic pressure (g) in incremental steps. The test Mach numbers were 0.65, 0.75, 0.82, 0.90, and 1.15. The four methods provided damping trends by which the flutter mode could be tracked and extrapolated to a flutter-onset q. A hard flutter point was obtained at M = 0.82. The peak-hold and cross-spectrum methods gave reliable results and could be most readily used for on-line testing. At M = 0.82, a p-k analysis predicted the same flutter mode as the experiment but a 6-percent lower flutter q. At the subcritical dynamic pressures, calculated damping values were appreciably lower than measured data.

Ruhlin, C. L.↗

The Two-Energies Turbulence Scheme Coupled to the Assumed PDF Method

An update of the two-energy turbulence scheme is presented, the 2TE + APDF scheme. The original version of the two-energy scheme is able to successfully model shallow convection without the need of an additional parameterization for non-local fluxes. However, the performance of the two-energy scheme is worse in stratocumulus cases, where it tends to overestimate the erosion of the stable layers. We have identified the causes: the non-local stability parameter does not consider local stratification, the scheme lacks an internal parameter that could distinguish between a shallow convection regime and a stratocumulus regime, and it uses an inflexible turbulence length scale formulation. To alleviate this problem, we propose several modifications: an update of the stability parameter, a modified computation of the turbulence length scale, and the introduction of the entropy potential temperature to distinguish between a shallow convection and a stratocumulus regime. In addition, the two-energy scheme is coupled to a simplified assumed probability density function method in order to achieve a more universal representation of the cloudy regimes. The updated turbulence scheme is evaluated for several idealized cases and one selected real case in the ICOsahedral Nonhydrostatic (ICON) modeling framework. The results show that the updated scheme corrects the overmixing problem in the stratocumulus cases. The performance of the updated scheme is comparable to the operational setup, and can be thus used instead of the operational turbulence and shallow convection scheme in ICON. Additionally, the updated scheme improves the coupling with dynamics, which is beneficial for the modeling of coherent flow structures in the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Updated users' guide for TAWFIVE with multigrid

A program for the Transonic Analysis of a Wing and Fuselage with Interacted Viscous Effects (TAWFIVE) was improved by the incorporation of multigrid and a method to specify lift coefficient rather than angle-of-attack. A finite volume full potential multigrid method is used to model the outer inviscid flow field. First order viscous effects are modeled by a 3-D integral boundary layer method. Both turbulent and laminar boundary layers are treated. Wake thickness effects are modeled using a 2-D strip method. A brief discussion of the engineering aspects of the program is given. The input, output, and use of the program are covered in detail. Sample results are given showing the effects of boundary layer corrections and the capability of the lift specification method.

Melson, N. Duane↗

Jet Noise Modeling Using an Acoustic Analogy: From Theoretical Formulation to Practical Noise Predictions

In this presentation an overview of the development of an acoustic analogy formulation, and its implementation into a reduced-order prediction method for turbulent jet mixing noise, will be given. The definition and rationale of an acoustic analogy approach to the prediction of flow-generated noise will be provided. A high-level description of the mathematical processes involved will be outlined in the context of a particular formulation developed at NASA Glenn Research Center. The modelling approximations needed to implement the general theoretical formulation into a practical noise prediction tool will be discussed. Sample results obtained using a noise prediction code based on the NASA Glenn formulation, and comparisons with experimental data, will be presented for a series of turbulent jet test cases. A view on the prospects for continued use of acoustic analogy-based methods for the prediction of flow-generated noise and the assessment of noise-reduction technologies will conclude the presentation.

Jet noise↗