Turbulence modeling of free shear layers for high-performance aircraft

NASA Technical Reports Server (NTRS)

Sondak, Douglas L.

1993-01-01

The High Performance Aircraft (HPA) Grand Challenge of the High Performance Computing and Communications (HPCC) program involves the computation of the flow over a high performance aircraft. A variety of free shear layers, including mixing layers over cavities, impinging jets, blown flaps, and exhaust plumes, may be encountered in such flowfields. Since these free shear layers are usually turbulent, appropriate turbulence models must be utilized in computations in order to accurately simulate these flow features. The HPCC program is relying heavily on parallel computers. A Navier-Stokes solver (POVERFLOW) utilizing the Baldwin-Lomax algebraic turbulence model was developed and tested on a 128-node Intel iPSC/860. Algebraic turbulence models run very fast, and give good results for many flowfields. For complex flowfields such as those mentioned above, however, they are often inadequate. It was therefore deemed that a two-equation turbulence model will be required for the HPA computations. The k-epsilon two-equation turbulence model was implemented on the Intel iPSC/860. Both the Chien low-Reynolds-number model and a generalized wall-function formulation were included.

Turbulent transport models for scramjet flowfields

NASA Technical Reports Server (NTRS)

Sindir, M. M.; Harsha, P. T.

1984-01-01

Turbulence modeling approaches were examined from the standpoint of their capability to predict the complex flowfield features observed in scramjet combustions. Thus, for example, the accuracy of each turbulence model, with respect to the prediction of recirculating flows, was examined. It was observed that for large diameter ratio axisymmetric sudden expansion flows, a choice of turbulence model was not critical because of the domination of their flowfields by pressure forces. For low diameter ratio axisymmetric sudden expansions and planar backward-facing steps flows, where turbulent shear stresses are of greater significance, the algebraic Reynolds stress approach, modified to increase its sensitivity to streamline curvature, was found to provide the best results. Results of the study also showed that strongly swirling flows provide a stringent test of turbulence model assumptions. Thus, although flows with very high swirl are not of great practical interest, they are useful for turbulence model development. Finally, it was also noted that numerical flowfields solution techniques have a strong interrelation with turbulence models, particularly with the turbulent transport models which involve source-dominated transport equations.

Prediction of Complex Aerodynamic Flows with Explicit Algebraic Stress Models

NASA Technical Reports Server (NTRS)

Abid, Ridha; Morrison, Joseph H.; Gatski, Thomas B.; Speziale, Charles G.

1996-01-01

An explicit algebraic stress equation, developed by Gatski and Speziale, is used in the framework of K-epsilon formulation to predict complex aerodynamic turbulent flows. The nonequilibrium effects are modeled through coefficients that depend nonlinearly on both rotational and irrotational strains. The proposed model was implemented in the ISAAC Navier-Stokes code. Comparisons with the experimental data are presented which clearly demonstrate that explicit algebraic stress models can predict the correct response to nonequilibrium flow.

Predicting NonInertial Effects with Algebraic Stress Models which Account for Dissipation Rate Anisotropies

NASA Technical Reports Server (NTRS)

Jongen, T.; Machiels, L.; Gatski, T. B.

1997-01-01

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coriolis-modified eddy-viscosity model, a realizable algebraic stress model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the turbulent model validation. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.

Algebraic Turbulence-Chemistry Interaction Model

NASA Technical Reports Server (NTRS)

Norris, Andrew T.

2012-01-01

The results of a series of Perfectly Stirred Reactor (PSR) and Partially Stirred Reactor (PaSR) simulations are compared to each other over a wide range of operating conditions. It is found that the PaSR results can be simulated by a PSR solution with just an adjusted chemical reaction rate. A simple expression has been developed that gives the required change in reaction rate for a PSR solution to simulate the PaSR results. This expression is the basis of a simple turbulence-chemistry interaction model. The interaction model that has been developed is intended for use with simple one-step global reaction mechanisms and for steady-state flow simulations. Due to the simplicity of the model there is very little additional computational cost in adding it to existing CFD codes.

Prediction of High-Lift Flows using Turbulent Closure Models

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Gatski, Thomas B.; Ying, Susan X.; Bertelrud, Arild

1997-01-01

The flow over two different multi-element airfoil configurations is computed using linear eddy viscosity turbulence models and a nonlinear explicit algebraic stress model. A subset of recently-measured transition locations using hot film on a McDonnell Douglas configuration is presented, and the effect of transition location on the computed solutions is explored. Deficiencies in wake profile computations are found to be attributable in large part to poor boundary layer prediction on the generating element, and not necessarily inadequate turbulence modeling in the wake. Using measured transition locations for the main element improves the prediction of its boundary layer thickness, skin friction, and wake profile shape. However, using measured transition locations on the slat still yields poor slat wake predictions. The computation of the slat flow field represents a key roadblock to successful predictions of multi-element flows. In general, the nonlinear explicit algebraic stress turbulence model gives very similar results to the linear eddy viscosity models.

An abbreviated Reynolds stress turbulence model for airfoil flows

NASA Technical Reports Server (NTRS)

Gaffney, R. L., Jr.; Hassan, H. A.; Salas, M. D.

1990-01-01

An abbreviated Reynolds stress turbulence model is presented for solving turbulent flow over airfoils. The model consists of two partial differential equations, one for the Reynolds shear stress and the other for the turbulent kinetic energy. The normal stresses and the dissipation rate of turbulent kinetic energy are computed from algebraic relationships having the correct asymptotic near wall behavior. This allows the model to be integrated all the way to the wall without the use of wall functions. Results for a flat plate at zero angle of attack, a NACA 0012 airfoil and a RAE 2822 airfoil are presented.

Prediction of Transonic Vortex Flows Using Linear and Nonlinear Turbulent Eddy Viscosity Models

NASA Technical Reports Server (NTRS)

Bartels, Robert E.; Gatski, Thomas B.

2000-01-01

Three-dimensional transonic flow over a delta wing is investigated with a focus on the effect of transition and influence of turbulence stress anisotropies. The performance of linear eddy viscosity models and an explicit algebraic stress model is assessed at the start of vortex flow, and the results compared with experimental data. To assess the effect of transition location, computations that either fix transition or are fully turbulent are performed. To assess the effect of the turbulent stress anisotropy, comparisons are made between predictions from the algebraic stress model and the linear eddy viscosity models. Both transition location and turbulent stress anisotropy significantly affect the 3D flow field. The most significant effect is found to be the modeling of transition location. At a Mach number of 0.90, the computed solution changes character from steady to unsteady depending on transition onset. Accounting for the anisotropies in the turbulent stresses also considerably impacts the flow, most notably in the outboard region of flow separation.

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.

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.

Assessment of turbulent models for scramjet flowfields

NASA Technical Reports Server (NTRS)

Sindir, M. M.; Harsha, P. T.

1982-01-01

The behavior of several turbulence models applied to the prediction of scramjet combustor flows is described. These models include the basic two equation model, the multiple dissipation length scale variant of the two equation model, and the algebraic stress model (ASM). Predictions were made of planar backward facing step flows and axisymmetric sudden expansion flows using each of these approaches. The formulation of each of these models are discussed, and the application of the different approaches to supersonic flows is described. A modified version of the ASM is found to provide the best prediction of the planar backward facing step flow in the region near the recirculation zone, while the basic ASM provides the best results downstream of the recirculation. Aspects of the interaction of numerica modeling and turbulences modeling as they affect the assessment of turbulence models are discussed.

The Use of DNS in Turbulence Modeling

NASA Technical Reports Server (NTRS)

Mansour, Nagi N.; Merriam, Marshal (Technical Monitor)

1997-01-01

The use of Direct numerical simulations (DNS) data in developing and testing turbulence models is reviewed. The data is used to test turbulence models at all levels: algebraic, one-equation, two-equation and full Reynolds stress models were tested. Particular examples on the development of models for the dissipation rate equation are presented. Homogeneous flows are used to test new scaling arguments for the various terms in the dissipation rate equation. The channel flow data is used to develop modifications to the equation model that take into account near-wall effects. DNS of compressible flows under mean compression are used in testing new compressible modifications to the two-equation models.

A comparative study of turbulence models for overset grids

NASA Technical Reports Server (NTRS)

Renze, Kevin J.; Buning, Pieter G.; Rajagopalan, R. G.

1992-01-01

The implementation of two different types of turbulence models for a flow solver using the Chimera overset grid method is examined. Various turbulence model characteristics, such as length scale determination and transition modeling, are found to have a significant impact on the computed pressure distribution for a multielement airfoil case. No inherent problem is found with using either algebraic or one-equation turbulence models with an overset grid scheme, but simulation of turbulence for multiple-body or complex geometry flows is very difficult regardless of the gridding method. For complex geometry flowfields, modification of the Baldwin-Lomax turbulence model is necessary to select the appropriate length scale in wall-bounded regions. The overset grid approach presents no obstacle to use of a one- or two-equation turbulence model. Both Baldwin-Lomax and Baldwin-Barth models have problems providing accurate eddy viscosity levels for complex multiple-body flowfields such as those involving the Space Shuttle.

Computation of confined coflow jets with three turbulence models

NASA Technical Reports Server (NTRS)

Zhu, J.; Shih, T. H.

1993-01-01

A numerical study of confined jets in a cylindrical duct is carried out to examine the performance of two recently proposed turbulence models: an RNG-based K-epsilon model and a realizable Reynolds stress algebraic equation model. The former is of the same form as the standard K-epsilon model but has different model coefficients. The latter uses an explicit quadratic stress-strain relationship to model the turbulent stresses and is capable of ensuring the positivity of each turbulent normal stress. The flow considered involves recirculation with unfixed separation and reattachment points and severe adverse pressure gradients, thereby providing a valuable test of the predictive capability of the models for complex flows. Calculations are performed with a finite-volume procedure. Numerical credibility of the solutions is ensured by using second-order accurate differencing schemes and sufficiently fine grids. Calculations with the standard K-epsilon model are also made for comparison. Detailed comparisons with experiments show that the realizable Reynolds stress algebraic equation model consistently works better than does the standard K-epsilon model in capturing the essential flow features, while the RNG-based K-epsilon model does not seem to give improvements over the standard K-epsilon model under the flow conditions considered.

A critical evaluation of various turbulence models as applied to internal fluid flows

NASA Technical Reports Server (NTRS)

Nallasamy, M.

1985-01-01

Models employed in the computation of turbulent flows are described and their application to internal flows is evaluated by examining the predictions of various turbulence models in selected flow configurations. The main conclusions are: (1) the k-epsilon model is used in a majority of all the two-dimensional flow calculations reported in the literature; (2) modified forms of the k-epsilon model improve the performance for flows with streamline curvature and heat transfer; (3) for flows with swirl, the k-epsilon model performs rather poorly; the algebraic stress model performs better in this case; and (4) for flows with regions of secondary flow (noncircular duct flows), the algebraic stress model performs fairly well for fully developed flow, for developing flow, the algebraic stress model performance is not good; a Reynolds stress model should be used. False diffusion and inlet boundary conditions are discussed. Countergradient transport and its implications in turbulence modeling is mentioned. Two examples of recirculating flow predictions obtained using PHOENICS code are discussed. The vortex method, large eddy simulation (modeling of subgrid scale Reynolds stresses), and direct simulation, are considered. Some recommendations for improving the model performance are made. The need for detailed experimental data in flows with strong curvature is emphasized.

Turbulence Modeling and Computation of Turbine Aerodynamics and Heat Transfer

NASA Technical Reports Server (NTRS)

Lakshminarayana, B.; Luo, J.

1996-01-01

The objective of the present research is to develop improved turbulence models for the computation of complex flows through turbomachinery passages, including the effects of streamline curvature, heat transfer and secondary flows. Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existance of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flows in strongly curved ducts, including two 180-deg ducts, one 90-deg duct and a strongly concave curved turbulent boundary layer have been carried out with Reynolds stress models (RSM) and algebraic Reynolds stress models (ARSM). An improved near-wall pressure-strain correlation has been developed for capturing the anisotropy of turbulence in the concave region. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several representative low-Reynolds number (LRN) k-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN k-epsilon models employed. The three-dimensional Navier-Stokes procedure developed by the present authors has been used to compute the three-dimensional viscous flow through the turbine nozzle passage of a single stage turbine. A low Reynolds number k-epsilon model and a zonal k-epsilon/ARSM (algebraic Reynolds stress model) are utilized for turbulence closure. An assessment of the performance of the turbulence models has been carried out. The two models are found to provide similar predictions for the mean flow parameters, although slight improvement in the prediction of some secondary flow quantities has been obtained by the

An extended algebraic variational multiscale-multigrid-multifractal method (XAVM4) for large-eddy simulation of turbulent two-phase flow

NASA Astrophysics Data System (ADS)

Rasthofer, U.; Wall, W. A.; Gravemeier, V.

2018-04-01

A novel and comprehensive computational method, referred to as the eXtended Algebraic Variational Multiscale-Multigrid-Multifractal Method (XAVM4), is proposed for large-eddy simulation of the particularly challenging problem of turbulent two-phase flow. The XAVM4 involves multifractal subgrid-scale modeling as well as a Nitsche-type extended finite element method as an approach for two-phase flow. The application of an advanced structural subgrid-scale modeling approach in conjunction with a sharp representation of the discontinuities at the interface between two bulk fluids promise high-fidelity large-eddy simulation of turbulent two-phase flow. The high potential of the XAVM4 is demonstrated for large-eddy simulation of turbulent two-phase bubbly channel flow, that is, turbulent channel flow carrying a single large bubble of the size of the channel half-width in this particular application.

Recent Turbulence Model Advances Applied to Multielement Airfoil Computations

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Gatski, Thomas B.

2000-01-01

A one-equation linear turbulence model and a two-equation nonlinear explicit algebraic stress model (EASM) are applied to the flow over a multielement airfoil. The effect of the K-epsilon and K-omega forms of the two-equation model are explored, and the K-epsilon form is shown to be deficient in the wall-bounded regions of adverse pressure gradient flows. A new K-omega form of EASM is introduced. Nonlinear terms present in EASM are shown to improve predictions of turbulent shear stress behind the trailing edge of the main element and near midflap. Curvature corrections are applied to both the one- and two-equation turbulence models and yield only relatively small local differences in the flap region, where the flow field undergoes the greatest curvature. Predictions of maximum lift are essentially unaffected by the turbulence model variations studied.

Research activities at the Center for Modeling of Turbulence and Transition

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing

1993-01-01

The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.

Center for Modeling of Turbulence and Transition (CMOTT): Research Briefs, 1992

NASA Technical Reports Server (NTRS)

Liou, William W. (Editor)

1992-01-01

The progress is reported of the Center for Modeling of Turbulence and Transition (CMOTT). The main objective of the CMOTT is to develop, validate and implement the turbulence and transition models for practical engineering flows. The flows of interest are three-dimensional, incompressible and compressible flows with chemical reaction. The research covers two-equation (e.g., k-e) and algebraic Reynolds-stress models, second moment closure models, probability density function (pdf) models, Renormalization Group Theory (RNG), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS).

On the validation of a code and a turbulence model appropriate to circulation control airfoils

NASA Technical Reports Server (NTRS)

Viegas, J. R.; Rubesin, M. W.; Maccormack, R. W.

1988-01-01

A computer code for calculating flow about a circulation control airfoil within a wind tunnel test section has been developed. This code is being validated for eventual use as an aid to design such airfoils. The concept of code validation being used is explained. The initial stages of the process have been accomplished. The present code has been applied to a low-subsonic, 2-D flow about a circulation control airfoil for which extensive data exist. Two basic turbulence models and variants thereof have been successfully introduced into the algorithm, the Baldwin-Lomax algebraic and the Jones-Launder two-equation models of turbulence. The variants include adding a history of the jet development for the algebraic model and adding streamwise curvature effects for both models. Numerical difficulties and difficulties in the validation process are discussed. Turbulence model and code improvements to proceed with the validation process are also discussed.

Some Recent Developments in Turbulence Closure Modeling

NASA Astrophysics Data System (ADS)

Durbin, Paul A.

2018-01-01

Turbulence closure models are central to a good deal of applied computational fluid dynamical analysis. Closure modeling endures as a productive area of research. This review covers recent developments in elliptic relaxation and elliptic blending models, unified rotation and curvature corrections, transition prediction, hybrid simulation, and data-driven methods. The focus is on closure models in which transport equations are solved for scalar variables, such as the turbulent kinetic energy, a timescale, or a measure of anisotropy. Algebraic constitutive representations are reviewed for their role in relating scalar closures to the Reynolds stress tensor. Seamless and nonzonal methods, which invoke a single closure model, are reviewed, especially detached eddy simulation (DES) and adaptive DES. Other topics surveyed include data-driven modeling and intermittency and laminar fluctuation models for transition prediction. The review concludes with an outlook.

Turbulence modeling for hypersonic flows

NASA Technical Reports Server (NTRS)

Marvin, J. G.; Coakley, T. J.

1989-01-01

Turbulence modeling for high speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary layer flows, shock wave boundary layer interactions and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.

Comparison of Turbulence Models for Nozzle-Afterbody Flows with Propulsive Jets

NASA Technical Reports Server (NTRS)

Compton, William B., III

1996-01-01

A numerical investigation was conducted to assess the accuracy of two turbulence models when computing non-axisymmetric nozzle-afterbody flows with propulsive jets. Navier-Stokes solutions were obtained for a Convergent-divergent non-axisymmetric nozzle-afterbody and its associated jet exhaust plume at free-stream Mach numbers of 0.600 and 0.938 at an angle of attack of 0 deg. The Reynolds number based on model length was approximately 20 x 10(exp 6). Turbulent dissipation was modeled by the algebraic Baldwin-Lomax turbulence model with the Degani-Schiff modification and by the standard Jones-Launder kappa-epsilon turbulence model. At flow conditions without strong shocks and with little or no separation, both turbulence models predicted the pressures on the surfaces of the nozzle very well. When strong shocks and massive separation existed, both turbulence models were unable to predict the flow accurately. Mixing of the jet exhaust plume and the external flow was underpredicted. The differences in drag coefficients for the two turbulence models illustrate that substantial development is still required for computing very complex flows before nozzle performance can be predicted accurately for all external flow conditions.

A comparison of three algebraic stress closures for combustor flow calculations

NASA Technical Reports Server (NTRS)

Nikjooy, M.; So, R. M. C.; Hwang, B. C.

1985-01-01

A comparison is made of the performance of two locally nonequilibrium and one equilibrium algebraic stress closures in calculating combustor flows. Effects of four different pressure-strain models on these closure models are also analyzed. The results show that the pressure-strain models have a much greater influence on the calculated mean velocity and turbulence field than the algebraic stress closures, and that the best mean strain model for the pressure-strain terms is that proposed by Launder, Reece and Rodi (1975). However, the equilibrium algebraic stress closure with the Rotta return-to-isotropy model (1951) for the pressure-strain terms gives as good a correlation with measurements as when the Launder et al. mean strain model is included in the pressure-strain model. Finally, comparison of the calculations with the standard k-epsilon closure results show that the algebraic stress closures are better suited for simple turbulent flow calculations.

A Study of Grid Resolution, Transition and Turbulence Model Using the Transonic Simple Straked Delta Wing

NASA Technical Reports Server (NTRS)

Bartels, Robert E.

2001-01-01

Three-dimensional transonic flow over a delta wing is investigated using several turbulence models. The performance of linear eddy viscosity models and an explicit algebraic stress model is assessed at the start of vortex flow, and the results compared with experimental data. To assess the effect of transition location, computations that either fix transition aft of the leading edge or are fully turbulent are performed. These computations show that grid resolution, transition location and turbulence model significantly affect the 3D flowfield.

Modifications of the law of the wall and algebraic turbulence modelling for separated boundary layers

NASA Technical Reports Server (NTRS)

Baldwin, B. S.; Maccormack, R. W.

1976-01-01

Various modifications of the conventional algebraic eddy viscosity turbulence model are investigated for application to separated flows. Friction velocity is defined in a way that avoids singular behavior at separation and reattachment but reverts to the conventional definition for flows with small pressure gradients. This leads to a modified law of the wall for separated flows. The effect on the calculated flow field of changes in the model that affect the eddy viscosity at various distances from the wall are determined by (1) switching from Prandtl's form to an inner layer formula due to Clauser at various distances from the wall, (2) varying the constant in the Van Driest damping factor, (3) using Clauser's inner layer formula all the way to the wall, and (4) applying a relaxation procedure in the evaluation of the constant in Clauser's inner layer formula. Numerical solutions of the compressible Navier-Stokes equations are used to determine the effects of the modifications. Experimental results from shock-induced separated flows at Mach numbers 2.93 and 8.45 are used for comparison. For these cases improved predictions of wall pressure distribution and positions of separation and reattachment are obtained from the relaxation version of the Clauser inner layer eddy viscosity formula.

Center for Modeling of Turbulence and Transition (CMOTT). Research briefs: 1990

NASA Technical Reports Server (NTRS)

Povinelli, Louis A. (Compiler); Liou, Meng-Sing (Compiler); Shih, Tsan-Hsing (Compiler)

1991-01-01

Brief progress reports of the Center for Modeling of Turbulence and Transition (CMOTT) research staff from May 1990 to May 1991 are given. The objectives of the CMOTT are to develop, validate, and implement the models for turbulence and boundary layer transition in the practical engineering flows. The flows of interest are three dimensional, incompressible, and compressible flows with chemistry. The schemes being studied include the two-equation and algebraic Reynolds stress models, the full Reynolds stress (or second moment closure) models, the probability density function models, the Renormalization Group Theory (RNG) and Interaction Approximation (DIA), the Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS).

Applications of Analytical Self-Similar Solutions of Reynolds-Averaged Models for Instability-Induced Turbulent Mixing

NASA Astrophysics Data System (ADS)

Hartland, Tucker; Schilling, Oleg

2017-11-01

Analytical self-similar solutions to several families of single- and two-scale, eddy viscosity and Reynolds stress turbulence models are presented for Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instability-induced turbulent mixing. The use of algebraic relationships between model coefficients and physical observables (e.g., experimental growth rates) following from the self-similar solutions to calibrate a member of a given family of turbulence models is shown. It is demonstrated numerically that the algebraic relations accurately predict the value and variation of physical outputs of a Reynolds-averaged simulation in flow regimes that are consistent with the simplifying assumptions used to derive the solutions. The use of experimental and numerical simulation data on Reynolds stress anisotropy ratios to calibrate a Reynolds stress model is briefly illustrated. The implications of the analytical solutions for future Reynolds-averaged modeling of hydrodynamic instability-induced mixing are briefly discussed. This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344.

RAS one-equation turbulence model with non-singular eddy-viscosity coefficient

NASA Astrophysics Data System (ADS)

Rahman, M. M.; Agarwal, R. K.; Siikonen, T.

2016-02-01

A simplified consistency formulation for Pk/ε (production to dissipation ratio) is devised to obtain a non-singular Cμ (coefficient of eddy-viscosity) in the explicit algebraic Reynolds stress model of Gatski and Speziale. The coefficient Cμ depends non-linearly on both rotational/irrotational strains and is used in the framework of an improved RAS (Rahman-Agarwal-Siikonen) one-equation turbulence model to calculate a few well-documented turbulent flows, yielding predictions in good agreement with the direct numerical simulation and experimental data. The strain-dependent Cμ assists the RAS model in constructing the coefficients and functions such as to benefit complex flows with non-equilibrium turbulence. Comparisons with the Spalart-Allmaras one-equation model and the shear stress transport k-ω model demonstrate that Cμ improves the response of RAS model to non-equilibrium effects.

Application of Navier-Stokes code PAB3D with kappa-epsilon turbulence model to attached and separated flows

NASA Technical Reports Server (NTRS)

Abdol-Hamid, Khaled S.; Lakshmanan, B.; Carlson, John R.

1995-01-01

A three-dimensional Navier-Stokes solver was used to determine how accurately computations can predict local and average skin friction coefficients for attached and separated flows for simple experimental geometries. Algebraic and transport equation closures were used to model turbulence. To simulate anisotropic turbulence, the standard two-equation turbulence model was modified by adding nonlinear terms. The effects of both grid density and the turbulence model on the computed flow fields were also investigated and compared with available experimental data for subsonic and supersonic free-stream conditions.

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.

Turbulent mass flux closure modeling for variable density turbulence in the wake of an air-entraining transom stern

NASA Astrophysics Data System (ADS)

Hendrickson, Kelli; Yue, Dick

2016-11-01

This work presents the development and a priori testing of closure models for the incompressible highly-variable density turbulent (IHVDT) flow in the near wake region of a transom stern. This complex, three-dimensional flow includes three regions with distinctly different flow behavior: (i) the convergent corner waves that originate from the body and collide on the ship center plane; (ii) the "rooster tail" that forms from the collision; and (iii) the diverging wave train. The characteristics of these regions involve violent free-surface flows and breaking waves with significant turbulent mass flux (TMF) at Atwood number At = (ρ2 -ρ1) / (ρ2 +ρ1) 1 for which there is little guidance in turbulence closure modeling for the momentum and scalar transport along the wake. Utilizing datasets from high-resolution simulations of the near wake of a canonical three-dimensional transom stern using conservative Volume-of-Fluid (cVOF), implicit Large Eddy Simulation (iLES), and Boundary Data Immersion Method (BDIM), we develop explicit algebraic turbulent mass flux closure models that incorporate the most relevant physical processes. Performance of these models in predicting the turbulent mass flux in all three regions of the wake will be presented. Office of Naval Research.

Separated transonic airfoil flow calculations with a nonequilibrium turbulence model

NASA Technical Reports Server (NTRS)

King, L. S.; Johnson, D. A.

1985-01-01

Navier-Stokes transonic airfoil calculations based on a recently developed nonequilibrium, turbulence closure model are presented for a supercritical airfoil section at transonic cruise conditions and for a conventional airfoil section at shock-induced stall conditions. Comparisons with experimental data are presented which show that this nonequilibrium closure model performs significantly better than the popular Baldwin-Lomax and Cebeci-Smith equilibrium algebraic models when there is boundary-layer separation that results from the inviscid-viscous interactions.

The applicability of turbulence models to aerodynamic and propulsion flowfields at McDonnell-Douglas Aerospace

NASA Technical Reports Server (NTRS)

Kral, Linda D.; Ladd, John A.; Mani, Mori

1995-01-01

The objective of this viewgraph presentation is to evaluate turbulence models for integrated aircraft components such as the forebody, wing, inlet, diffuser, nozzle, and afterbody. The one-equation models have replaced the algebraic models as the baseline turbulence models. The Spalart-Allmaras one-equation model consistently performs better than the Baldwin-Barth model, particularly in the log-layer and free shear layers. Also, the Sparlart-Allmaras model is not grid dependent like the Baldwin-Barth model. No general turbulence model exists for all engineering applications. The Spalart-Allmaras one-equation model and the Chien k-epsilon models are the preferred turbulence models. Although the two-equation models often better predict the flow field, they may take from two to five times the CPU time. Future directions are in further benchmarking the Menter blended k-w/k-epsilon and algorithmic improvements to reduce CPU time of the two-equation model.

A study of the second and third order closure models of turbulence for prediction of separated shear flows

NASA Technical Reports Server (NTRS)

Amano, R. S.

1985-01-01

The hybrid model of the Reynolds-stress turbulence closure is tested for the computation of the flows over a step and disk. Here it is attempted to improve the redistributive action of the turbulence energy among the Reynolds stresses. By evaluating the existing models for the pressure-strain correlation, better coefficients are obtained for the prediction of separating shear flows. Furthermore, the diffusion rate of the Reynolds stresses is reevaluated adopting several algebraic correlations for the triple-velocity products. The models of Cormack et al., Daly-Harlow, Hanjalic-Launder, and Shir were tested for the reattaching shear flows. It was generally observed that all these algebraic models give considerably low values of the triple-velocity products. This is attributed to the fact that none of the algebraic models can take the convective effect of the triple-velocity products into account in the separating shear flows, thus resulting in much lower diffusion rate than Reynolds stresses. In order to improve the evaluation of these quantities correction factors are introduced based on the comparison with some experimental data.

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.

Direct Numerical Simulation of Complex Turbulence

NASA Astrophysics Data System (ADS)

Hsieh, Alan

Direct numerical simulations (DNS) of spanwise-rotating turbulent channel flow were conducted. The data base obtained from these DNS simulations were used to investigate the turbulence generation cycle for simple and complex turbulence. For turbulent channel flow, three theoretical models concerning the formation and evolution of sublayer streaks, three-dimensional hairpin vortices and propagating plane waves were validated using visualizations from the present DNS data. The principal orthogonal decomposition (POD) method was used to verify the existence of the propagating plane waves; a new extension of the POD method was derived to demonstrate these plane waves in a spatial channel model. The analyses of coherent structures was extended to complex turbulence and used to determine the proper computational box size for a minimal flow unit (MFU) at Rob < 0.5. Proper realization of Taylor-Gortler vortices in the highly turbulent pressure region was demonstrated to be necessary for acceptably accurate MFU turbulence statistics, which required a minimum spanwise domain length Lz = pi. A dependence of MFU accuracy on Reynolds number was also discovered and MFU models required a larger domain to accurately approximate higher-Reynolds number flows. In addition, the results obtained from the DNS simulations were utilized to evaluate several turbulence closure models for momentum and thermal transport in rotating turbulent channel flow. Four nonlinear eddy viscosity turbulence models were tested and among these, Explicit Algebraic Reynolds Stress Models (EARSM) obtained the Reynolds stress distributions in best agreement with DNS data for rotational flows. The modeled pressure-strain functions of EARSM were shown to have strong influence on the Reynolds stress distributions near the wall. Turbulent heatflux distributions obtained from two explicit algebraic heat flux models consistently displayed increasing disagreement with DNS data with increasing rotation rate. Results

On the modelling of non-reactive and reactive turbulent combustor flows

NASA Technical Reports Server (NTRS)

Nikjooy, Mohammad; So, Ronald M. C.

1987-01-01

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Numerical prediction of an axisymmetric turbulent mixing layer using two turbulence models

NASA Astrophysics Data System (ADS)

Johnson, Richard W.

1992-01-01

Nuclear power, once considered and then rejected (in the U. S.) for application to space vehicle propulsion, is being reconsidered for powering space rockets, especially for interplanetary travel. The gas core reactor, a high risk, high payoff nuclear engine concept, is one that was considered in the 1960s and 70s. As envisioned then, the gas core reactor would consist of a heavy, slow moving core of fissioning uranium vapor surrounded by a fast moving outer stream of hydrogen propellant. Satisfactory operation of such a configuration would require stable nuclear reaction kinetics to occur simultaneously with a stable, coflowing, probably turbulent fluid system having a dense inner stream and a light outer stream. The present study examines the behavior of two turbulence models in numerically simulating an idealized version of the above coflowing fluid system. The two models are the standard k˜ɛ model and a thin shear algebraic stress model (ASM). The idealized flow system can be described as an axisymmetric mixing layer of constant density. Predictions for the radial distribution of the mean streamwise velocity and shear stress for several axial stations are compared with experiment. Results for the k˜ɛe predictions are broadly satisfactory while those for the ASM are distinctly poorer.

An application of a two-equation model of turbulence to three-dimensional chemically reacting flows

NASA Technical Reports Server (NTRS)

Lee, J.

1994-01-01

A numerical study of three dimensional chemically reacting and non-reacting flowfields is conducted using a two-equation model of turbulence. A generalized flow solver using an implicit Lower-Upper (LU) diagonal decomposition numerical technique and finite-rate chemistry has been coupled with a low-Reynolds number two-equation model of turbulence. This flow solver is then used to study chemically reacting turbulent supersonic flows inside combustors with synergetic fuel injectors. The reacting and non-reacting turbulent combustor solutions obtained are compared with zero-equation turbulence model solutions and with available experimental data. The hydrogen-air chemistry is modeled using a nine-species/eighteen reaction model. A low-Reynolds number k-epsilon model was used to model the effect of turbulence because, in general, the low-Reynolds number k-epsilon models are easier to implement numerically and are far more general than algebraic models. However, low-Reynolds number k-epsilon models require a much finer near-wall grid resolution than high-Reynolds number models to resolve accurately the near-wall physics. This is especially true in complex flowfields, where the stiff nature of the near-wall turbulence must be resolved. Therefore, the limitations imposed by the near-wall characteristics and compressible model corrections need to be evaluated further. The gradient-diffusion hypothesis is used to model the effects of turbulence on the mass diffusion process. The influence of this low-Reynolds number turbulence model on the reacting flowfield predictions was studied parametrically.

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.

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.

A resolvable subfilter-scale model specific to large-eddy simulation of under-resolved turbulence

NASA Astrophysics Data System (ADS)

Zhou, Yong; Brasseur, James G.; Juneja, Anurag

2001-09-01

Large-eddy simulation (LES) of boundary-layer flows has serious deficiencies near the surface when a viscous sublayer either does not exist (rough walls) or is not practical to resolve (high Reynolds numbers). In previous work, we have shown that the near-surface errors arise from the poor performance of algebraic subfilter-scale (SFS) models at the first several grid levels, where integral scales are necessarily under-resolved and the turbulence is highly anisotropic. In under-resolved turbulence, eddy viscosity and similarity SFS models create a spurious feedback loop between predicted resolved-scale (RS) velocity and modeled SFS acceleration, and are unable to simultaneously capture SFS acceleration and RS-SFS energy flux. To break the spurious coupling in a dynamically meaningful manner, we introduce a new modeling strategy in which the grid-resolved subfilter velocity is estimated from a separate dynamical equation containing the essential inertial interactions between SFS and RS velocity. This resolved SFS (RSFS) velocity is then used as a surrogate for the complete SFS velocity in the SFS stress tensor. We test the RSFS model by comparing LES of highly under-resolved anisotropic buoyancy-generated homogeneous turbulence with a corresponding direct numerical simulation (DNS). The new model successfully suppresses the spurious feedback loop between RS velocity and SFS acceleration, and greatly improves model predictions of the anisotropic structure of SFS acceleration and resolved velocity fields. Unlike algebraic models, the RSFS model accurately captures SFS acceleration intensity and RS-SFS energy flux, even during the nonequilibrium transient, and properly partitions SFS acceleration between SFS stress divergence and SFS pressure force.

Modeling the pressure-strain correlation of turbulence: An invariant dynamical systems approach

NASA Technical Reports Server (NTRS)

Speziale, Charles G.; Sarkar, Sutanu; Gatski, Thomas B.

1990-01-01

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Modelling the pressure-strain correlation of turbulence - An invariant dynamical systems approach

NASA Technical Reports Server (NTRS)

Speziale, Charles G.; Sarkar, Sutanu; Gatski, Thomas B.

1991-01-01

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Prediction of Airfoil Characteristics With Higher Order Turbulence Models

NASA Technical Reports Server (NTRS)

Gatski, Thomas B.

1996-01-01

This study focuses on the prediction of airfoil characteristics, including lift and drag over a range of Reynolds numbers. Two different turbulence models, which represent two different types of models, are tested. The first is a standard isotropic eddy-viscosity two-equation model, and the second is an explicit algebraic stress model (EASM). The turbulent flow field over a general-aviation airfoil (GA(W)-2) at three Reynolds numbers is studied. At each Reynolds number, predicted lift and drag values at different angles of attack are compared with experimental results, and predicted variations of stall locations with Reynolds number are compared with experimental data. Finally, the size of the separation zone predicted by each model is analyzed, and correlated with the behavior of the lift coefficient near stall. In summary, the EASM model is able to predict the lift and drag coefficients over a wider range of angles of attack than the two-equation model for the three Reynolds numbers studied. However, both models are unable to predict the correct lift and drag behavior near the stall angle, and for the lowest Reynolds number case, the two-equation model did not predict separation on the airfoil near stall.

On turbulent flows dominated by curvature effects

NASA Technical Reports Server (NTRS)

Cheng, G. C.; Farokhi, S.

1992-01-01

A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows and the flow in a curved duct are examples of flowfields where streamline curvature plays a dominant role. New algebraic formulations for the eddy viscosity incorporating the k-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature.

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.

Calculation of turbulence-driven secondary motion in ducts with arbitrary cross section

NASA Technical Reports Server (NTRS)

Demuren, A. O.

1989-01-01

Calculation methods for turbulent duct flows are generalized for ducts with arbitrary cross-sections. The irregular physical geometry is transformed into a regular one in computational space, and the flow equations are solved with a finite-volume numerical procedure. The turbulent stresses are calculated with an algebraic stress model derived by simplifying model transport equations for the individual Reynolds stresses. Two variants of such a model are considered. These procedures enable the prediction of both the turbulence-driven secondary flow and the anisotropy of the Reynolds stresses, in contrast to some of the earlier calculation methods. Model predictions are compared to experimental data for developed flow in triangular duct, trapezoidal duct and a rod-bundle geometry. The correct trends are predicted, and the quantitative agreement is mostly fair. The simpler variant of the algebraic stress model procured better agreement with the measured data.

A Modified Mixing Length Turbulence Model for Zero and Adverse Pressure Gradients. M.S. Thesis - Akron Univ., 1993

NASA Technical Reports Server (NTRS)

Conley, Julianne M.; Leonard, B. P.

1994-01-01

The modified mixing length (MML) turbulence model was installed in the Proteus Navier-Stokes code, then modified to make it applicable to a wider range of flows typical of aerospace propulsion applications. The modifications are based on experimental data for three flat-plate flows having zero, mild adverse, and strong adverse pressure gradients. Three transonic diffuser test cases were run with the new version of the model in order to evaluate its performance. All results are compared with experimental data and show improvements over calculations made using the Baldwin-Lomax turbulence model, the standard algebraic model in Proteus.

A numerical study of confined turbulent jets

NASA Technical Reports Server (NTRS)

Zhu, J.; Shih, T.-H.

1993-01-01

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Modelling thermal radiation in buoyant turbulent diffusion flames

NASA Astrophysics Data System (ADS)

Consalvi, J. L.; Demarco, R.; Fuentes, A.

2012-10-01

This work focuses on the numerical modelling of radiative heat transfer in laboratory-scale buoyant turbulent diffusion flames. Spectral gas and soot radiation is modelled by using the Full-Spectrum Correlated-k (FSCK) method. Turbulence-Radiation Interactions (TRI) are taken into account by considering the Optically-Thin Fluctuation Approximation (OTFA), the resulting time-averaged Radiative Transfer Equation (RTE) being solved by the Finite Volume Method (FVM). Emission TRIs and the mean absorption coefficient are then closed by using a presumed probability density function (pdf) of the mixture fraction. The mean gas flow field is modelled by the Favre-averaged Navier-Stokes (FANS) equation set closed by a buoyancy-modified k-ɛ model with algebraic stress/flux models (ASM/AFM), the Steady Laminar Flamelet (SLF) model coupled with a presumed pdf approach to account for Turbulence-Chemistry Interactions, and an acetylene-based semi-empirical two-equation soot model. Two sets of experimental pool fire data are used for validation: propane pool fires 0.3 m in diameter with Heat Release Rates (HRR) of 15, 22 and 37 kW and methane pool fires 0.38 m in diameter with HRRs of 34 and 176 kW. Predicted flame structures, radiant fractions, and radiative heat fluxes on surrounding surfaces are found in satisfactory agreement with available experimental data across all the flames. In addition further computations indicate that, for the present flames, the gray approximation can be applied for soot with a minor influence on the results, resulting in a substantial gain in Computer Processing Unit (CPU) time when the FSCK is used to treat gas radiation.

Isolating Curvature Effects in Computing Wall-Bounded Turbulent Flows

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Gatski, Thomas B.

2001-01-01

The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

An Initial Investigation of the Effects of Turbulence Models on the Convergence of the RK/Implicit Scheme

NASA Technical Reports Server (NTRS)

Swanson, R. C.; Rossow, C.-C.

2008-01-01

A three-stage Runge-Kutta (RK) scheme with multigrid and an implicit preconditioner has been shown to be an effective solver for the fluid dynamic equations. This scheme has been applied to both the compressible and essentially incompressible Reynolds-averaged Navier-Stokes (RANS) equations using the algebraic turbulence model of Baldwin and Lomax (BL). In this paper we focus on the convergence of the RK/implicit scheme when the effects of turbulence are represented by either the Spalart-Allmaras model or the Wilcox k-! model, which are frequently used models in practical fluid dynamic applications. Convergence behavior of the scheme with these turbulence models and the BL model are directly compared. For this initial investigation we solve the flow equations and the partial differential equations of the turbulence models indirectly coupled. With this approach we examine the convergence behavior of each system. Both point and line symmetric Gauss-Seidel are considered for approximating the inverse of the implicit operator of the flow solver. To solve the turbulence equations we use a diagonally dominant alternating direction implicit (DDADI) scheme. Computational results are presented for three airfoil flow cases and comparisons are made with experimental data. We demonstrate that the two-dimensional RANS equations and transport-type equations for turbulence modeling can be efficiently solved with an indirectly coupled algorithm that uses the RK/implicit scheme for the flow equations.

Prediction of free turbulent mixing using a turbulent kinetic energy method

NASA Technical Reports Server (NTRS)

Harsha, P. T.

1973-01-01

Free turbulent mixing of two-dimensional and axisymmetric one- and two-stream flows is analyzed by a relatively simple turbulent kinetic energy method. This method incorporates a linear relationship between the turbulent shear and the turbulent kinetic energy and an algebraic relationship for the length scale appearing in the turbulent kinetic energy equation. Good results are obtained for a wide variety of flows. The technique is shown to be especially applicable to flows with heat and mass transfer, for which nonunity Prandtl and Schmidt numbers may be assumed.

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.

Four-component numerical simulation model of radiative convective interactions in large-scale oxygen-hydrogen turbulent fire balls

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

Surzhikov, S.T.

1996-12-31

Two-dimensional radiative gas dynamics model for numerical simulation of oxygen-hydrogen fire ball which may be generated by an explosion of a launch vehicle with cryogenic (LO{sub 2}-LH{sub 2}) fuel components is presented. The following physical-chemical processes are taken into account in the numerical model: and effective chemical reaction between the gaseous components (O{sub 2}-H{sub 2}) of the propellant, turbulent mixing and diffusion of the components, and radiative heat transfer. The results of numerical investigations of the following problems are presented: The influence of radiative heat transfer on fire ball gas dynamics during the first 13 sec after explosion, the effectmore » of the fuel gaseous components afterburning on fire ball gas dynamics, and the effect of turbulence on fire ball gas dynamics (in a framework of algebraic model of turbulent mixing).« less

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.

Direct numerical simulations and modeling of a spatially-evolving turbulent wake

NASA Technical Reports Server (NTRS)

Cimbala, John M.

1994-01-01

Understanding of turbulent free shear flows (wakes, jets, and mixing layers) is important, not only for scientific interest, but also because of their appearance in numerous practical applications. Turbulent wakes, in particular, have recently received increased attention by researchers at NASA Langley. The turbulent wake generated by a two-dimensional airfoil has been selected as the test-case for detailed high-resolution particle image velocimetry (PIV) experiments. This same wake has also been chosen to enhance NASA's turbulence modeling efforts. Over the past year, the author has completed several wake computations, while visiting NASA through the 1993 and 1994 ASEE summer programs, and also while on sabbatical leave during the 1993-94 academic year. These calculations have included two-equation (K-omega and K-epsilon) models, algebraic stress models (ASM), full Reynolds stress closure models, and direct numerical simulations (DNS). Recently, there has been mutually beneficial collaboration of the experimental and computational efforts. In fact, these projects have been chosen for joint presentation at the NASA Turbulence Peer Review, scheduled for September 1994. DNS calculations are presently underway for a turbulent wake at Re(sub theta) = 1000 and at a Mach number of 0.20. (Theta is the momentum thickness, which remains constant in the wake of a two dimensional body.) These calculations utilize a compressible DNS code written by M. M. Rai of NASA Ames, and modified for the wake by J. Cimbala. The code employs fifth-order accurate upwind-biased finite differencing for the convective terms, fourth-order accurate central differencing for the viscous terms, and an iterative-implicit time-integration scheme. The computational domain for these calculations starts at x/theta = 10, and extends to x/theta = 610. Fully developed turbulent wake profiles, obtained from experimental data from several wake generators, are supplied at the computational inlet, along with

Nonlinear Reynolds stress model for turbulent shear flows

NASA Technical Reports Server (NTRS)

Barton, J. Michael; Rubinstein, R.; Kirtley, K. R.

1991-01-01

A nonlinear algebraic Reynolds stress model, derived using the renormalization group, is applied to equilibrium homogeneous shear flow and fully developed flow in a square duct. The model, which is quadratically nonlinear in the velocity gradients, successfully captures the large-scale inhomogeneity and anisotropy of the flows studied. The ratios of normal stresses, as well as the actual magnitudes of the stresses are correctly predicted for equilibrium homogeneous shear flow. Reynolds normal stress anisotropy and attendant turbulence driven secondary flow are predicted for a square duct. Profiles of mean velocity and normal stresses are in good agreement with measurements. Very close to walls, agreement with measurements diminishes. The model has the benefit of containing no arbitrary constants; all values are determined directly from the theory. It seems that near wall behavior is influenced by more than the large scale anisotropy accommodated in the current model. More accurate near wall calculations may well require a model for anisotropic dissipation.

Evaluation of Turbulence-Model Performance in Jet Flows

NASA Technical Reports Server (NTRS)

Woodruff, S. L.; Seiner, J. M.; Hussaini, M. Y.; Erlebacher, G.

2001-01-01

The importance of reducing jet noise in both commercial and military aircraft applications has made jet acoustics a significant area of research. A technique for jet noise prediction commonly employed in practice is the MGB approach, based on the Lighthill acoustic analogy. This technique requires as aerodynamic input mean flow quantities and turbulence quantities like the kinetic energy and the dissipation. The purpose of the present paper is to assess existing capabilities for predicting these aerodynamic inputs. Two modern Navier-Stokes flow solvers, coupled with several modern turbulence models, are evaluated by comparison with experiment for their ability to predict mean flow properties in a supersonic jet plume. Potential weaknesses are identified for further investigation. Another comparison with similar intent is discussed by Barber et al. The ultimate goal of this research is to develop a reliable flow solver applicable to the low-noise, propulsion-efficient, nozzle exhaust systems being developed in NASA focused programs. These programs address a broad range of complex nozzle geometries operating in high temperature, compressible, flows. Seiner et al. previously discussed the jet configuration examined here. This convergent-divergent nozzle with an exit diameter of 3.6 inches was designed for an exhaust Mach number of 2.0 and a total temperature of 1680 F. The acoustic and aerodynamic data reported by Seiner et al. covered a range of jet total temperatures from 104 F to 2200 F at the fully-expanded nozzle pressure ratio. The aerodynamic data included centerline mean velocity and total temperature profiles. Computations were performed independently with two computational fluid dynamics (CFD) codes, ISAAC and PAB3D. Turbulence models employed include the k-epsilon model, the Gatski-Speziale algebraic-stress model and the Girimaji model, with and without the Sarkar compressibility correction. Centerline values of mean velocity and mean temperature are

Category-theoretic models of algebraic computer systems

NASA Astrophysics Data System (ADS)

Kovalyov, S. P.

2016-01-01

A computer system is said to be algebraic if it contains nodes that implement unconventional computation paradigms based on universal algebra. A category-based approach to modeling such systems that provides a theoretical basis for mapping tasks to these systems' architecture is proposed. The construction of algebraic models of general-purpose computations involving conditional statements and overflow control is formally described by a reflector in an appropriate category of algebras. It is proved that this reflector takes the modulo ring whose operations are implemented in the conventional arithmetic processors to the Łukasiewicz logic matrix. Enrichments of the set of ring operations that form bases in the Łukasiewicz logic matrix are found.

Modeling of turbulence and transition

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing

1992-01-01

The first objective is to evaluate current two-equation and second order closure turbulence models using available direct numerical simulations and experiments, and to identify the models which represent the state of the art in turbulence modeling. The second objective is to study the near-wall behavior of turbulence, and to develop reliable models for an engineering calculation of turbulence and transition. The third objective is to develop a two-scale model for compressible turbulence.

Turbulence modeling: Near-wall turbulence and effects of rotation on turbulence

NASA Technical Reports Server (NTRS)

Shih, T.-H.

1990-01-01

Many Reynolds averaged Navier-Stokes solvers use closure models in conjunction with 'the law of the wall', rather than deal with a thin, viscous sublayer near the wall. This work is motivated by the need for better models to compute near wall turbulent flow. The authors use direct numerical simulation of fully developed channel flow and one of three dimensional turbulent boundary layer flow to develop new models. These direct numerical simulations provide detailed data that experimentalists have not been able to measure directly. Another objective of the work is to examine analytically the effects of rotation on turbulence, using Rapid Distortion Theory (RDT). This work was motivated by the observation that the pressure strain models in all current second order closure models are unable to predict the effects of rotation on turbulence.

Using Students' Interests as Algebraic Models

ERIC Educational Resources Information Center

Whaley, Kenneth A.

2012-01-01

Fostering algebraic thinking is an important goal for middle-grades mathematics teachers. Developing mathematical reasoning requires that teachers cultivate students' habits of mind. Teachers develop students' understanding of algebra by engaging them in tasks that involve modeling and representation. This study was designed to investigate how…

Computational algebraic geometry of epidemic models

NASA Astrophysics Data System (ADS)

Rodríguez Vega, Martín.

2014-06-01

Computational Algebraic Geometry is applied to the analysis of various epidemic models for Schistosomiasis and Dengue, both, for the case without control measures and for the case where control measures are applied. The models were analyzed using the mathematical software Maple. Explicitly the analysis is performed using Groebner basis, Hilbert dimension and Hilbert polynomials. These computational tools are included automatically in Maple. Each of these models is represented by a system of ordinary differential equations, and for each model the basic reproductive number (R0) is calculated. The effects of the control measures are observed by the changes in the algebraic structure of R0, the changes in Groebner basis, the changes in Hilbert dimension, and the changes in Hilbert polynomials. It is hoped that the results obtained in this paper become of importance for designing control measures against the epidemic diseases described. For future researches it is proposed the use of algebraic epidemiology to analyze models for airborne and waterborne diseases.

Compressibility Considerations for kappa-omega Turbulence Models in Hypersonic Boundary Layer Applications

NASA Technical Reports Server (NTRS)

Rumsey, C. L.

2009-01-01

The ability of kappa-omega models to predict compressible turbulent skin friction in hypersonic boundary layers is investigated. Although uncorrected two-equation models can agree well with correlations for hot-wall cases, they tend to perform progressively worse - particularly for cold walls - as the Mach number is increased in the hypersonic regime. Simple algebraic models such as Baldwin-Lomax perform better compared to experiments and correlations in these circumstances. Many of the compressibility corrections described in the literature are summarized here. These include corrections that have only a small influence for kappa-omega models, or that apply only in specific circumstances. The most widely-used general corrections were designed for use with jet or mixing-layer free shear flows. A less well-known dilatation-dissipation correction intended for boundary layer flows is also tested, and is shown to agree reasonably well with the Baldwin-Lomax model at cold-wall conditions. It exhibits a less dramatic influence than the free shear type of correction. There is clearly a need for improved understanding and better overall physical modeling for turbulence models applied to hypersonic boundary layer flows.

Ocean Turbulence I: One-Point Closure Model Momentum and Heat Vertical Diffusivities

NASA Technical Reports Server (NTRS)

Canuto, V. M.; Howard, A.; Cheng, Y.; Dubovikov, M. S.

1999-01-01

Since the early forties, one-point turbulence closure models have been the canonical tools used to describe turbulent flows in many fields. In geophysics, Mellor and Yamada applied such models using the 1980 state-of-the art. Since then, no improvements were introduced to alleviate two major difficulties: 1) closure of the pressure correlations, which affects the correct determination of the critical Richardson number Ri(sub cr) above which turbulent mixing is no longer possible and 2) the need to express the non-local third-order moments (TOM) in terms of lower order moments rather than via the down-gradient approximation as done thus far, since the latter seriously underestimates the TOMs. Since 1) and 2) are still being dealt with adjustable parameters which weaken the credibility of the models, alternative models, not based on turbulence modeling, have been suggested. The aim of this paper is to show that new information, partly derived from the newest 2-point closure model discussed, can be used to solve these shortcomings. The new one-point closure model, which in its simplest form is algebraic and thus simple to implement, is first shown to reproduce a variety of data. Then, it is used in a Ocean-General Circulation Model (O-GCM) where it reproduces well a large variety of ocean data. While phenomenological models are specifically tuned to ocean turbulence, the present model is not. It is first tested against laboratory data on stably stratified flows and then used in an O-GCM. It is more general, more predictive and more resilient, e.g., it can incorporate phenomena like wave-breaking at the surface, salinity diffusivity, non-locality, etc. One important feature that naturally comes out of the new model is that the predicted Richardson critical value Ri(sub cr) is Ri (sub cr approx. = 1) in agreement with both Large Eddy Simulations (LES) and empirical evidence while all previous models predicted Ri (sub cr approx. = 0.2) which led to a considerable

Turbulence Model Comparisons and Reynolds Number Effects Over a High-Speed Aircraft at Transonic Speeds

NASA Technical Reports Server (NTRS)

Rivers, Melissa B.; Wahls, Richard A.

1999-01-01

This paper gives the results of a grid study, a turbulence model study, and a Reynolds number effect study for transonic flows over a high-speed aircraft using the thin-layer, upwind, Navier-Stokes CFL3D code. The four turbulence models evaluated are the algebraic Baldwin-Lomax model with the Degani-Schiff modifications, the one-equation Baldwin-Barth model, the one-equation Spalart-Allmaras model, and Menter's two-equation Shear-Stress-Transport (SST) model. The flow conditions, which correspond to tests performed in the NASA Langley National Transonic Facility (NTF), are a Mach number of 0.90 and a Reynolds number of 30 million based on chord for a range of angle-of-attacks (1 degree to 10 degrees). For the Reynolds number effect study, Reynolds numbers of 10 and 80 million based on chord were also evaluated. Computed forces and surface pressures compare reasonably well with the experimental data for all four of the turbulence models. The Baldwin-Lomax model with the Degani-Schiff modifications and the one-equation Baldwin-Barth model show the best agreement with experiment overall. The Reynolds number effects are evaluated using the Baldwin-Lomax with the Degani-Schiff modifications and the Baldwin-Barth turbulence models. Five angles-of-attack were evaluated for the Reynolds number effect study at three different Reynolds numbers. More work is needed to determine the ability of CFL3D to accurately predict Reynolds number effects.

A mixing timescale model for TPDF simulations of turbulent premixed flames

DOE PAGES

Kuron, Michael; Ren, Zhuyin; Hawkes, Evatt R.; ...

2017-02-06

Transported probability density function (TPDF) methods are an attractive modeling approach for turbulent flames as chemical reactions appear in closed form. However, molecular micro-mixing needs to be modeled and this modeling is considered a primary challenge for TPDF methods. In the present study, a new algebraic mixing rate model for TPDF simulations of turbulent premixed flames is proposed, which is a key ingredient in commonly used molecular mixing models. The new model aims to properly account for the transition in reactive scalar mixing rate behavior from the limit of turbulence-dominated mixing to molecular mixing behavior in flamelets. An a priorimore » assessment of the new model is performed using direct numerical simulation (DNS) data of a lean premixed hydrogen–air jet flame. The new model accurately captures the mixing timescale behavior in the DNS and is found to be a significant improvement over the commonly used constant mechanical-to-scalar mixing timescale ratio model. An a posteriori TPDF study is then performed using the same DNS data as a numerical test bed. The DNS provides the initial conditions and time-varying input quantities, including the mean velocity, turbulent diffusion coefficient, and modeled scalar mixing rate for the TPDF simulations, thus allowing an exclusive focus on the mixing model. Here, the new mixing timescale model is compared with the constant mechanical-to-scalar mixing timescale ratio coupled with the Euclidean Minimum Spanning Tree (EMST) mixing model, as well as a laminar flamelet closure. It is found that the laminar flamelet closure is unable to properly capture the mixing behavior in the thin reaction zones regime while the constant mechanical-to-scalar mixing timescale model under-predicts the flame speed. Furthermore, the EMST model coupled with the new mixing timescale model provides the best prediction of the flame structure and flame propagation among the models tested, as the dynamics of reactive

A mixing timescale model for TPDF simulations of turbulent premixed flames

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

Kuron, Michael; Ren, Zhuyin; Hawkes, Evatt R.

Transported probability density function (TPDF) methods are an attractive modeling approach for turbulent flames as chemical reactions appear in closed form. However, molecular micro-mixing needs to be modeled and this modeling is considered a primary challenge for TPDF methods. In the present study, a new algebraic mixing rate model for TPDF simulations of turbulent premixed flames is proposed, which is a key ingredient in commonly used molecular mixing models. The new model aims to properly account for the transition in reactive scalar mixing rate behavior from the limit of turbulence-dominated mixing to molecular mixing behavior in flamelets. An a priorimore » assessment of the new model is performed using direct numerical simulation (DNS) data of a lean premixed hydrogen–air jet flame. The new model accurately captures the mixing timescale behavior in the DNS and is found to be a significant improvement over the commonly used constant mechanical-to-scalar mixing timescale ratio model. An a posteriori TPDF study is then performed using the same DNS data as a numerical test bed. The DNS provides the initial conditions and time-varying input quantities, including the mean velocity, turbulent diffusion coefficient, and modeled scalar mixing rate for the TPDF simulations, thus allowing an exclusive focus on the mixing model. Here, the new mixing timescale model is compared with the constant mechanical-to-scalar mixing timescale ratio coupled with the Euclidean Minimum Spanning Tree (EMST) mixing model, as well as a laminar flamelet closure. It is found that the laminar flamelet closure is unable to properly capture the mixing behavior in the thin reaction zones regime while the constant mechanical-to-scalar mixing timescale model under-predicts the flame speed. Furthermore, the EMST model coupled with the new mixing timescale model provides the best prediction of the flame structure and flame propagation among the models tested, as the dynamics of reactive

Workshop on Computational Turbulence Modeling

NASA Technical Reports Server (NTRS)

1993-01-01

This document contains presentations given at Workshop on Computational Turbulence Modeling held 15-16 Sep. 1993. The purpose of the meeting was to discuss the current status and future development of turbulence modeling in computational fluid dynamics for aerospace propulsion systems. Papers cover the following topics: turbulence modeling activities at the Center for Modeling of Turbulence and Transition (CMOTT); heat transfer and turbomachinery flow physics; aerothermochemistry and computational methods for space systems; computational fluid dynamics and the k-epsilon turbulence model; propulsion systems; and inlet, duct, and nozzle flow.

Comparative Study of Advanced Turbulence Models for Turbomachinery

NASA Technical Reports Server (NTRS)

Hadid, Ali H.; Sindir, Munir M.

1996-01-01

A computational study has been undertaken to study the performance of advanced phenomenological turbulence models coded in a modular form to describe incompressible turbulent flow behavior in two dimensional/axisymmetric and three dimensional complex geometry. The models include a variety of two equation models (single and multi-scale k-epsilon models with different near wall treatments) and second moment algebraic and full Reynolds stress closure models. These models were systematically assessed to evaluate their performance in complex flows with rotation, curvature and separation. The models are coded as self contained modules that can be interfaced with a number of flow solvers. These modules are stand alone satellite programs that come with their own formulation, finite-volume discretization scheme, solver and boundary condition implementation. They will take as input (from any generic Navier-Stokes solver) the velocity field, grid (structured H-type grid) and computational domain specification (boundary conditions), and will deliver, depending on the model used, turbulent viscosity, or the components of the Reynolds stress tensor. There are separate 2D/axisymmetric and/or 3D decks for each module considered. The modules are tested using Rocketdyn's proprietary code REACT. The code utilizes an efficient solution procedure to solve Navier-Stokes equations in a non-orthogonal body-fitted coordinate system. The differential equations are discretized over a finite-volume grid using a non-staggered variable arrangement and an efficient solution procedure based on the SIMPLE algorithm for the velocity-pressure coupling is used. The modules developed have been interfaced and tested using finite-volume, pressure-correction CFD solvers which are widely used in the CFD community. Other solvers can also be used to test these modules since they are independently structured with their own discretization scheme and solver methodology. Many of these modules have been

Numerical Simulation of High-Speed Turbulent Reacting Flows

NASA Technical Reports Server (NTRS)

Givi, P.; Taulbee, D. B.; Madnia, C. K.; Jaberi, F. A.; Colucci, P. J.; Gicquel, L. Y. M.; Adumitroaie, V.; James, S.

1999-01-01

The objectives of this research are: (1) to develop and implement a new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. (2) To develop algebraic turbulence closures for statistical description of chemically reacting turbulent flows.

Aircraft Dynamic Modeling in Turbulence

NASA Technical Reports Server (NTRS)

Morelli, Eugene A.; Cunninham, Kevin

2012-01-01

A method for accurately identifying aircraft dynamic models in turbulence was developed and demonstrated. The method uses orthogonal optimized multisine excitation inputs and an analytic method for enhancing signal-to-noise ratio for dynamic modeling in turbulence. A turbulence metric was developed to accurately characterize the turbulence level using flight measurements. The modeling technique was demonstrated in simulation, then applied to a subscale twin-engine jet transport aircraft in flight. Comparisons of modeling results obtained in turbulent air to results obtained in smooth air were used to demonstrate the effectiveness of the approach.

Application of a range of turbulence energy models to the determination of M4 tidal current profiles

NASA Astrophysics Data System (ADS)

Xing, Jiuxing; Davies, Alan M.

1996-04-01

A fully nonlinear, three-dimensional hydrodynamic model of the Irish Sea, using a range of turbulence energy sub-models, is used to examine the influence of the turbulence closure method upon the vertical variation of the current profile of the fundamental and higher harmonics of the tide in the region. Computed tidal current profiles are compared with previous calculations using a spectral model with eddy viscosity related to the flow field. The model has a sufficiently fine grid to resolve the advection terms, in particular the advection of turbulence and momentum. Calculations show that the advection of turbulence energy does not have a significant influence upon the current profile of either the fundamental or higher harmonic of the tide, although the advection of momentum is important in the region of headlands. The simplification of the advective terms by only including them in their vertically integrated form does not appear to make a significant difference to current profiles, but does reduce the computational effort by a significant amount. Computed current profiles both for the fundamental and the higher harmonic determined with a prognostic equation for turbulence and an algebraic mixing length formula, are as accurate as those determined with a two prognostic equation model (the so called q2- q2l model), provided the mixing length is specified correctly. A simple, flow-dependent eddy viscosity with a parabolic variation of viscosity also performs equally well.

Turbulence kinetic energy equation for dilute suspensions

NASA Technical Reports Server (NTRS)

Abou-Arab, T. W.; Roco, M. C.

1989-01-01

A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

Chiral algebras in Landau-Ginzburg models

NASA Astrophysics Data System (ADS)

Dedushenko, Mykola

2018-03-01

Chiral algebras in the cohomology of the {\\overline{Q}}+ supercharge of two-dimensional N=(0,2) theories on flat spacetime are discussed. Using the supercurrent multiplet, we show that the answer is renormalization group invariant for theories with an R-symmetry. For N=(0,2) Landau-Ginzburg models, the chiral algebra is determined by the operator equations of motion, which preserve their classical form, and quantum renormalization of composite operators. We study these theories and then specialize to the N=(2,2) models and consider some examples.

Tensor models, Kronecker coefficients and permutation centralizer algebras

NASA Astrophysics Data System (ADS)

Geloun, Joseph Ben; Ramgoolam, Sanjaye

2017-11-01

We show that the counting of observables and correlators for a 3-index tensor model are organized by the structure of a family of permutation centralizer algebras. These algebras are shown to be semi-simple and their Wedderburn-Artin decompositions into matrix blocks are given in terms of Clebsch-Gordan coefficients of symmetric groups. The matrix basis for the algebras also gives an orthogonal basis for the tensor observables which diagonalizes the Gaussian two-point functions. The centres of the algebras are associated with correlators which are expressible in terms of Kronecker coefficients (Clebsch-Gordan multiplicities of symmetric groups). The color-exchange symmetry present in the Gaussian model, as well as a large class of interacting models, is used to refine the description of the permutation centralizer algebras. This discussion is extended to a general number of colors d: it is used to prove the integrality of an infinite family of number sequences related to color-symmetrizations of colored graphs, and expressible in terms of symmetric group representation theory data. Generalizing a connection between matrix models and Belyi maps, correlators in Gaussian tensor models are interpreted in terms of covers of singular 2-complexes. There is an intriguing difference, between matrix and higher rank tensor models, in the computational complexity of superficially comparable correlators of observables parametrized by Young diagrams.

Workshop on Computational Turbulence Modeling

NASA Technical Reports Server (NTRS)

Shabbir, A. (Compiler); Shih, T.-H. (Compiler); Povinelli, L. A. (Compiler)

1994-01-01

The purpose of this meeting was to discuss the current status and future development of turbulence modeling in computational fluid dynamics for aerospace propulsion systems. Various turbulence models have been developed and applied to different turbulent flows over the past several decades and it is becoming more and more urgent to assess their performance in various complex situations. In order to help users in selecting and implementing appropriate models in their engineering calculations, it is important to identify the capabilities as well as the deficiencies of these models. This also benefits turbulence modelers by permitting them to further improve upon the existing models. This workshop was designed for exchanging ideas and enhancing collaboration between different groups in the Lewis community who are using turbulence models in propulsion related CFD. In this respect this workshop will help the Lewis goal of excelling in propulsion related research. This meeting had seven sessions for presentations and one panel discussion over a period of two days. Each presentation session was assigned to one or two branches (or groups) to present their turbulence related research work. Each group was asked to address at least the following points: current status of turbulence model applications and developments in the research; progress and existing problems; and requests about turbulence modeling. The panel discussion session was designed for organizing committee members to answer management and technical questions from the audience and to make concluding remarks.

General Algebraic Modeling System Tutorial | High-Performance Computing |

Science.gov Websites

power generation from two different fuels. The goal is to minimize the cost for one of the fuels while Here's a basic tutorial for modeling optimization problems with the General Algebraic Modeling System (GAMS). Overview The GAMS (General Algebraic Modeling System) package is essentially a compiler for a

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

Evolution of phase singularities of vortex beams propagating in atmospheric turbulence.

PubMed

Ge, Xiao-Lu; Wang, Ben-Yi; Guo, Cheng-Shan

2015-05-01

Optical vortex beams propagating through atmospheric turbulence are studied by numerical modeling, and the phase singularities of the vortices existing in the turbulence-distorted beams are calculated. It is found that the algebraic sum of topological charges (TCs) of all the phase singularities existing in test aperture is approximately equal to the TC of the input vortex beam. This property provides us a possible approach for determining the TC of the vortex beam propagating through the atmospheric turbulence, which could have potential application in optical communication using optical vortices.

Algebra for Enterprise Ontology: towards analysis and synthesis of enterprise models

NASA Astrophysics Data System (ADS)

Suga, Tetsuya; Iijima, Junichi

2018-03-01

Enterprise modeling methodologies have made enterprises more likely to be the object of systems engineering rather than craftsmanship. However, the current state of research in enterprise modeling methodologies lacks investigations of the mathematical background embedded in these methodologies. Abstract algebra, a broad subfield of mathematics, and the study of algebraic structures may provide interesting implications in both theory and practice. Therefore, this research gives an empirical challenge to establish an algebraic structure for one aspect model proposed in Design & Engineering Methodology for Organizations (DEMO), which is a major enterprise modeling methodology in the spotlight as a modeling principle to capture the skeleton of enterprises for developing enterprise information systems. The results show that the aspect model behaves well in the sense of algebraic operations and indeed constructs a Boolean algebra. This article also discusses comparisons with other modeling languages and suggests future work.

Multifractal vector fields and stochastic Clifford algebra.

PubMed

Schertzer, Daniel; Tchiguirinskaia, Ioulia

2015-12-01

In the mid 1980s, the development of multifractal concepts and techniques was an important breakthrough for complex system analysis and simulation, in particular, in turbulence and hydrology. Multifractals indeed aimed to track and simulate the scaling singularities of the underlying equations instead of relying on numerical, scale truncated simulations or on simplified conceptual models. However, this development has been rather limited to deal with scalar fields, whereas most of the fields of interest are vector-valued or even manifold-valued. We show in this paper that the combination of stable Lévy processes with Clifford algebra is a good candidate to bridge up the present gap between theory and applications. We show that it indeed defines a convenient framework to generate multifractal vector fields, possibly multifractal manifold-valued fields, based on a few fundamental and complementary properties of Lévy processes and Clifford algebra. In particular, the vector structure of these algebra is much more tractable than the manifold structure of symmetry groups while the Lévy stability grants a given statistical universality.

Computation of Turbulent Wake Flows in Variable Pressure Gradient

NASA Technical Reports Server (NTRS)

Duquesne, N.; Carlson, J. R.; Rumsey, C. L.; Gatski, T. B.

1999-01-01

Transport aircraft performance is strongly influenced by the effectiveness of high-lift systems. Developing wakes generated by the airfoil elements are subjected to strong pressure gradients and can thicken very rapidly, limiting maximum lift. This paper focuses on the effects of various pressure gradients on developing symmetric wakes and on the ability of a linear eddy viscosity model and a non-linear explicit algebraic stress model to accurately predict their downstream evolution. In order to reduce the uncertainties arising from numerical issues when assessing the performance of turbulence models, three different numerical codes with the same turbulence models are used. Results are compared to available experimental data to assess the accuracy of the computational results.

Computation of flows in a turn-around duct and a turbine cascade using advanced turbulence models

NASA Technical Reports Server (NTRS)

Lakshminarayana, B.; Luo, J.

1993-01-01

Numerical investigation has been carried out to evaluate the capability of the Algebraic Reynolds Stress Model (ARSM) and the Nonlinear Stress Model (NLSM) to predict strongly curved turbulent flow in a turn-around duct (TAD). The ARSM includes the near-wall damping term of pressure-strain correlation phi(sub ij,w), which enables accurate prediction of individual Reynolds stress components in wall flows. The TAD mean flow quantities are reasonably well predicted by various turbulence models. The ARSM yields better predictions for both the mean flow and the turbulence quantities than the NLSM and the k-epsilon (k = turbulent kinetic energy, epsilon = dissipation rate of k) model. The NLSM also shows slight improvement over the k-epsilon model. However, all the models fail to capture the recovery of the flow from strong curvature effects. The formulation for phi(sub ij,w) appears to be incorrect near the concave surface. The hybrid k-epsilon/ARSM, Chien's k-epsilon model, and Coakley's q-omega (q = the square root of k, omega = epsilon/k) model have also been employed to compute the aerodynamics and heat transfer of a transonic turbine cascade. The surface pressure distributions and the wake profiles are predicted well by all the models. The k-epsilon model and the k-epsilon/ARSM model provide better predictions of heat transfer than the q-omega model. The k-epsilon/ARSM solutions show significant differences in the predicted skin friction coefficients, heat transfer rates and the cascade performance parameters, as compared to the k-epsilon model. The k-epsilon/ARSM model appears to capture, qualitatively, the anisotropy associated with by-pass transition.

A rational approach to the use of Prandtl's mixing length model in free turbulent shear flow calculations

NASA Technical Reports Server (NTRS)

Rudy, D. H.; Bushnell, D. M.

1973-01-01

Prandtl's basic mixing length model was used to compute 22 test cases on free turbulent shear flows. The calculations employed appropriate algebraic length scale equations and single values of mixing length constant for planar and axisymmetric flows, respectively. Good agreement with data was obtained except for flows, such as supersonic free shear layers, where large sustained sensitivity changes occur. The inability to predict the more gradual mixing in these flows is tentatively ascribed to the presence of a significant turbulence-induced transverse static pressure gradient which is neglected in conventional solution procedures. Some type of an equation for length scale development was found to be necessary for successful computation of highly nonsimilar flow regions such as jet or wake development from thick wall flows.

Modeling of near-wall turbulence

NASA Technical Reports Server (NTRS)

Shih, T. H.; Mansour, N. N.

1990-01-01

An improved k-epsilon model and a second order closure model is presented for low Reynolds number turbulence near a wall. For the k-epsilon model, a modified form of the eddy viscosity having correct asymptotic near wall behavior is suggested, and a model for the pressure diffusion term in the turbulent kinetic energy equation is proposed. For the second order closure model, the existing models are modified for the Reynolds stress equations to have proper near wall behavior. A dissipation rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The calculations are compared with direct numerical simulations. It is shown that the present models, both the k-epsilon model and the second order closure model, perform well in predicting the behavior of the near wall turbulence. Significant improvements over previous models are obtained.

Computation of turbulent reacting flow in a solid-propellant ducted rocket

NASA Astrophysics Data System (ADS)

Chao, Yei-Chin; Chou, Wen-Fuh; Liu, Sheng-Shyang

1995-05-01

A mathematical model for computation of turbulent reacting flows is developed under general curvilinear coordinate systems. An adaptive, streamline grid system is generated to deal with the complex flow structures in a multiple-inlet solid-propellant ducted rocket (SDR) combustor. General tensor representations of the k-epsilon and algebraic stress (ASM) turbulence models are derived in terms of contravariant velocity components, and modification caused by the effects of compressible turbulence is also included in the modeling. The clipped Gaussian probability density function is incorporated in the combustion model to account for fluctuations of properties. Validation of the above modeling is first examined by studying mixing and reacting characteristics in a confined coaxial-jet problem. This is followed by study of nonreacting and reacting SDR combustor flows. The results show that Gibson and Launder's ASM incorporated with Sarkar's modification for compressible turbulence effects based on the general curvilinear coordinate systems yields the most satisfactory prediction for this complicated SDR flowfield.

Efficiency and Accuracy of Time-Accurate Turbulent Navier-Stokes Computations

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Sanetrik, Mark D.; Biedron, Robert T.; Melson, N. Duane; Parlette, Edward B.

1995-01-01

The accuracy and efficiency of two types of subiterations in both explicit and implicit Navier-Stokes codes are explored for unsteady laminar circular-cylinder flow and unsteady turbulent flow over an 18-percent-thick circular-arc (biconvex) airfoil. Grid and time-step studies are used to assess the numerical accuracy of the methods. Nonsubiterative time-stepping schemes and schemes with physical time subiterations are subject to time-step limitations in practice that are removed by pseudo time sub-iterations. Computations for the circular-arc airfoil indicate that a one-equation turbulence model predicts the unsteady separated flow better than an algebraic turbulence model; also, the hysteresis with Mach number of the self-excited unsteadiness due to shock and boundary-layer separation is well predicted.

A Realizable Reynolds Stress Algebraic Equation Model

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Zhu, Jiang; Lumley, John L.

1993-01-01

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

A turbulence model for iced airfoils and its validation

NASA Technical Reports Server (NTRS)

Shin, Jaiwon; Chen, Hsun H.; Cebeci, Tuncer

1992-01-01

A turbulence model based on the extension of the algebraic eddy viscosity formulation of Cebeci and Smith developed for two dimensional flows over smooth and rough surfaces is described for iced airfoils and validated for computed ice shapes obtained for a range of total temperatures varying from 28 to -15 F. The validation is made with an interactive boundary layer method which uses a panel method to compute the inviscid flow and an inverse finite difference boundary layer method to compute the viscous flow. The interaction between inviscid and viscous flows is established by the use of the Hilbert integral. The calculated drag coefficients compare well with recent experimental data taken at the NASA-Lewis Icing Research Tunnel (IRT) and show that, in general, the drag increase due to ice accretion can be predicted well and efficiently.

The effects of streamline curvature and swirl on turbulent flows in curved ducts

NASA Technical Reports Server (NTRS)

Cheng, Chih-Hsiung; Farokhi, Saeed

1990-01-01

A technique for improving the numerical predictions of turbulent flows with the effect of streamline curvature is developed. Separated flows, the flow in a curved duct, and swirling flows are examples of flow fields where streamline curvature plays a dominant role. A comprehensive literature review on the effect of streamline curvature was conducted. New algebraic formulations for the eddy viscosity incorporating the kappa-epsilon turbulence model are proposed to account for various effects of streamline curvature. The loci of flow reversal of the separated flows over various backward-facing steps are employed to test the capability of the proposed turbulence model in capturing the effect of local curvature. The inclusion of the effect of longitudinal curvature in the proposed turbulence model is validated by predicting the distributions of the static pressure coefficients in an S-bend duct and in 180 degree turn-around ducts. The proposed turbulence model embedded with transverse curvature modification is substantiated by predicting the decay of the axial velocities in the confined swirling flows. The numerical predictions of different curvature effects by the proposed turbulence models are also reported.

An Algebraic Formulation of Level One Wess-Zumino Models

NASA Astrophysics Data System (ADS)

Böckenhauer, Jens

The highest weight modules of the chiral algebra of orthogonal WZW models at level one possess a realization in fermionic representation spaces; the Kac-Moody and Virasoro generators are represented as unbounded limits of even CAR algebras. It is shown that the representation theory of the underlying even CAR algebras reproduces precisely the sectors of the chiral algebra. This fact allows to develop a theory of local von Neumann algebras on the punctured circle, fitting nicely in the Doplicher-Haag-Roberts framework. The relevant localized endomorphisms which generate the charged sectors are explicitly constructed by means of Bogoliubov transformations. Using CAR theory, the fusion rules in terms of sector equivalence classes are proven.

Shell models of magnetohydrodynamic turbulence

NASA Astrophysics Data System (ADS)

Plunian, Franck; Stepanov, Rodion; Frick, Peter

2013-02-01

Shell models of hydrodynamic turbulence originated in the seventies. Their main aim was to describe the statistics of homogeneous and isotropic turbulence in spectral space, using a simple set of ordinary differential equations. In the eighties, shell models of magnetohydrodynamic (MHD) turbulence emerged based on the same principles as their hydrodynamic counter-part but also incorporating interactions between magnetic and velocity fields. In recent years, significant improvements have been made such as the inclusion of non-local interactions and appropriate definitions for helicities. Though shell models cannot account for the spatial complexity of MHD turbulence, their dynamics are not over simplified and do reflect those of real MHD turbulence including intermittency or chaotic reversals of large-scale modes. Furthermore, these models use realistic values for dimensionless parameters (high kinetic and magnetic Reynolds numbers, low or high magnetic Prandtl number) allowing extended inertial range and accurate dissipation rate. Using modern computers it is difficult to attain an inertial range of three decades with direct numerical simulations, whereas eight are possible using shell models. In this review we set up a general mathematical framework allowing the description of any MHD shell model. The variety of the latter, with their advantages and weaknesses, is introduced. Finally we consider a number of applications, dealing with free-decaying MHD turbulence, dynamo action, Alfvén waves and the Hall effect.

Modeling the turbulent kinetic energy equation for compressible, homogeneous turbulence

NASA Technical Reports Server (NTRS)

Aupoix, B.; Blaisdell, G. A.; Reynolds, William C.; Zeman, Otto

1990-01-01

The turbulent kinetic energy transport equation, which is the basis of turbulence models, is investigated for homogeneous, compressible turbulence using direct numerical simulations performed at CTR. It is shown that the partition between dilatational and solenoidal modes is very sensitive to initial conditions for isotropic decaying turbulence but not for sheared flows. The importance of the dilatational dissipation and of the pressure-dilatation term is evidenced from simulations and a transport equation is proposed to evaluate the pressure-dilatation term evolution. This transport equation seems to work well for sheared flows but does not account for initial condition sensitivity in isotropic decay. An improved model is proposed.

Advancements in engineering turbulence modeling

NASA Technical Reports Server (NTRS)

Shih, T.-H.

1991-01-01

Some new developments in two-equation models and second order closure models are presented. Two-equation models (k-epsilon models) have been widely used in computational fluid dynamics (CFD) for engineering problems. Most of low-Reynolds number two-equation models contain some wall-distance damping functions to account for the effect of wall on turbulence. However, this often causes the confusion and difficulties in computing flows with complex geometry and also needs an ad hoc treatment near the separation and reattachment points. A set of modified two-equation models is proposed to remove the aforementioned shortcomings. The calculations using various two-equation models are compared with direct numerical simulations of channel flow and flat boundary layers. Development of a second order closure model is also discussed with emphasis on the modeling of pressure related correlation terms and dissipation rates in the second moment equations. All the existing models poorly predict the normal stresses near the wall and fail to predict the 3-D effect of mean flow on the turbulence (e.g. decrease in the shear stress caused by the cross flow in the boundary layer). The newly developed second order near-wall turbulence model is described and is capable of capturing the near-wall behavior of turbulence as well as the effect of 3-D mean flow on the turbulence.

Comparison of numerical predictions of horizontal nonisothermal jet in a room with three turbulence models -- {kappa}-{epsilon} EVM, ASM, and DSM

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

Murakami, Shuzo; Kato, Shinsuke; Ooka, Ryozo

1994-12-31

A three-dimensional nonisothermal jet in a room is analyzed numerically by the standard {kappa}-{epsilon} eddy viscosity model (EVM) and two second-moment closure models-the algebraic stress model (ASM) (Hossain and Rodi 1982) and the differential stress model (DSM) (Launder et al. 1975). Numerical results given by these turbulence models are compared with experimental results, and the prediction errors existing in the results are examined, thus clarifying the relative structural differences between the {kappa}-{epsilon} EVM and the second-moment closure models. Since the second moment closure models clearly manifest the turbulence structures of the flow field, they are more accurate than the {kappa}-{epsilon}more » EVM. A small difference between the DSM and the ASM -- one based on an inappropriate approximation of the convection and diffusion terms in the Reynolds stress transport equations in the ASM -- is also observed.« less

Numerical analysis of hypersonic turbulent film cooling flows

NASA Technical Reports Server (NTRS)

Chen, Y. S.; Chen, C. P.; Wei, H.

1992-01-01

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

On Thermodynamic Constraints upon Turbulence Modeling

NASA Astrophysics Data System (ADS)

Huang, Yu-Ning; Durst, Franz

2000-11-01

Turbulence is a continuum phenomenon which can be described within the framework of continuum mechanics. Such foundation has the potential for improving turbulence modeling, making it less heuristic and more rational. In the present research, we consider the compatibility of turbulence modeling with the second law of thermodynamics. We show that the Clausius-Planck inequality, as an expression of the principle of entropy growth, places a thermodynamic restriction upon the turbulence modeling of an incompressible Navier-Stokes fluid in an isothermal temperature field. This thermodynamic restriction is given in the form of an inequality, which ensures non-negativeness of the mean internal dissipation. As an illustration, we show the thermodynamic constraints on the modeling of a few typical homogeneous turbulent flows.

Reynolds-Averaged Turbulence Model Assessment for a Highly Back-Pressured Isolator Flowfield

NASA Technical Reports Server (NTRS)

Baurle, Robert A.; Middleton, Troy F.; Wilson, L. G.

2012-01-01

The use of computational fluid dynamics in scramjet engine component development is widespread in the existing literature. Unfortunately, the quantification of model-form uncertainties is rarely addressed with anything other than sensitivity studies, requiring that the computational results be intimately tied to and calibrated against existing test data. This practice must be replaced with a formal uncertainty quantification process for computational fluid dynamics to play an expanded role in the system design, development, and flight certification process. Due to ground test facility limitations, this expanded role is believed to be a requirement by some in the test and evaluation community if scramjet engines are to be given serious consideration as a viable propulsion device. An effort has been initiated at the NASA Langley Research Center to validate several turbulence closure models used for Reynolds-averaged simulations of scramjet isolator flows. The turbulence models considered were the Menter BSL, Menter SST, Wilcox 1998, Wilcox 2006, and the Gatski-Speziale explicit algebraic Reynolds stress models. The simulations were carried out using the VULCAN computational fluid dynamics package developed at the NASA Langley Research Center. A procedure to quantify the numerical errors was developed to account for discretization errors in the validation process. This procedure utilized the grid convergence index defined by Roache as a bounding estimate for the numerical error. The validation data was collected from a mechanically back-pressured constant area (1 2 inch) isolator model with an isolator entrance Mach number of 2.5. As expected, the model-form uncertainty was substantial for the shock-dominated, massively separated flowfield within the isolator as evidenced by a 6 duct height variation in shock train length depending on the turbulence model employed. Generally speaking, the turbulence models that did not include an explicit stress limiter more closely

Modification of the MML turbulence model for adverse pressure gradient flows. M.S. Thesis - Akron Univ., 1993

NASA Technical Reports Server (NTRS)

Conley, Julianne M.

1994-01-01

Computational fluid dynamics is being used increasingly to predict flows for aerospace propulsion applications, yet there is still a need for an easy to use, computationally inexpensive turbulence model capable of accurately predicting a wide range of turbulent flows. The Baldwin-Lomax model is the most widely used algebraic model, even though it has known difficulties calculating flows with strong adverse pressure gradients and large regions of separation. The modified mixing length model (MML) was developed specifically to handle the separation which occurs on airfoils and has given significantly better results than the Baldwin-Lomax model. The success of these calculations warrants further evaluation and development of MML. The objective of this work was to evaluate the performance of MML for zero and adverse pressure gradient flows, and modify it as needed. The Proteus Navier-Stokes code was used for this study and all results were compared with experimental data and with calculations made using the Baldwin-Lomax algebraic model, which is currently available in Proteus. The MML model was first evaluated for zero pressure gradient flow over a flat plate, then modified to produce the proper boundary layer growth. Additional modifications, based on experimental data for three adverse pressure gradient flows, were also implemented. The adapted model, called MMLPG (modified mixing length model for pressure gradient flows), was then evaluated for a typical propulsion flow problem, flow through a transonic diffuser. Three cases were examined: flow with no shock, a weak shock and a strong shock. The results of these calculations indicate that the objectives of this study have been met. Overall, MMLPG is capable of accurately predicting the adverse pressure gradient flows examined in this study, giving generally better agreement with experimental data than the Baldwin-Lomax model.

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.

Modeling of near wall turbulence and modeling of bypass transition

NASA Technical Reports Server (NTRS)

Yang, Z.

1992-01-01

The objectives for this project are as follows: (1) Modeling of the near wall turbulence: We aim to develop a second order closure for the near wall turbulence. As a first step of this project, we try to develop a kappa-epsilon model for near wall turbulence. We require the resulting model to be able to handle both near wall turbulence and turbulent flows away from the wall, computationally robust, and applicable for complex flow situations, flow with separation, for example, and (2) Modeling of the bypass transition: We aim to develop a bypass transition model which contains the effect of intermittency. Thus, the model can be used for both the transitional boundary layers and the turbulent boundary layers. We require the resulting model to give a good prediction of momentum and heat transfer within the transitional boundary and a good prediction of the effect of freestream turbulence on transitional boundary layers.

A second-order modelling of a stably stratified sheared turbulence submitted to a non-vertical shear

NASA Astrophysics Data System (ADS)

Bouzaiane, Mounir; Ben Abdallah, Hichem; Lili, Taieb

2004-09-01

In this work, the evolution of homogeneous stably stratified turbulence submitted to a non-vertical shear is studied using second-order closure models. Two cases of turbulent flows are considered. Firstly, the case of a purely horizontal shear is considered. In this case, the evolution of the turbulence is studied according to the Richardson number Ri which is varied from 0.2 to 2.0 when other parameters are kept constant. In the second case, two components of shear are present. The turbulence is submitted to a vertical component Sv = partU1/partx3 = S cos(thgr) and a horizontal component Sh = partU1/partx2 = S sin(thgr). In this case, we study the influence of shear inclination angle thgr on the evolution of turbulence. In both cases, we are referred respectively to the recent direct numerical simulations of Jacobitz (2002 J. Turbulence 3 055) and Jacobitz and Sarkar (1998 Phys. Fluids 10 1158-68) which are, to our knowledge, the most recent results of the above-mentioned flows. Transport equations of second-order moments \\overline{u_{i} u_{j}} , \\overline{u_{i} \\rho } , \\overline{\\rho^{2}} are derived. The Shih-Lumley (SL) (Shih T H 1996 Turbulence Transition and Modeling ed H D S Henningson, A V Johansson and P H Alfredsson (Dordrecht: Kluwer); Shih and Lumley J L 1989 27th Aerospace Meeting 9-12 January, Center of Turbulent Research, Nevada) and the Craft-Launder (CL) (Craft T J and Launder B E 1989 Turbulent Shear Flow Stanford University, USA, pp 12-1-12-6 Launder B E 1996 Turbulence Transition and Modeling ed H D S Henningson, A V Johansson and P H Alfredsson (Dordrecht: Kluwer)) second-order models are retained for the pressure-strain correlation phgrij and the pressure-scalar gradient correlation phgrirgr. The corresponding models are also retained for the dissipation egr of the turbulent kinetic energy and an algebraic model is retained for the dissipation egrrgrrgr of the variance of the scalar. A fourth-order Runge-Kutta method is used for the

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

Modeling of Turbulent Swirling Flows

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Zhu, Jiang; Liou, William; Chen, Kuo-Huey; Liu, Nan-Suey; Lumley, John L.

1997-01-01

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Turbulence Modelling in Wind Turbine Wakes =

NASA Astrophysics Data System (ADS)

Olivares Espinosa, Hugo

With the expansion of the wind energy industry, wind parks have become a common appearance in our landscapes. Owing to restrictions of space or to economic reasons, wind turbines are located close to each other in wind farms. This causes interference problems which reduce the efficiency of the array. In particular, the wind turbine wakes increase the level of turbulence and cause a momentum defect that may lead to an increase of mechanical loads and to a reduction of power output. Thus, it is important for the wind energy industry to predict the characteristics of the turbulence field in the wakes with the purpose of increasing the efficiency of the power extraction. Since this is a phenomenon of intrinsically non-linear nature, it can only be accurately described by the full set of the Navier-Stokes equations. Furthermore, a proper characterization of turbulence cannot be made without resolving the turbulent motions, so neither linearized models nor the widely used Reynolds-Averaged Navier-Stokes model can be employed. Instead, Large-Eddy Simulations (LES) provide a feasible alternative, where the energy containing fluctuations of the velocity field are resolved and the effects of the smaller eddies are modelled through a sub-grid scale component. The objective of this work is the modelling of turbulence in wind turbine wakes in a homogeneous turbulence inflow. A methodology has been developed to fulfill this objective. Firstly, a synthetic turbulence field is introduced into a computational domain where LES are performed to simulate a decaying turbulence flow. Secondly, the Actuator Disk (AD) technique is employed to simulate the effect of a rotor in the incoming flow and produce a turbulent wake. The implementation is carried out in OpenFOAM, an open-source CFD platform, resembling a well documented procedure previously used for wake flow simulations. Results obtained with the proposed methodology are validated by comparing with values obtained from wind tunnel

A spatial operator algebra for manipulator modeling and control

NASA Technical Reports Server (NTRS)

Rodriguez, G.; Jain, A.; Kreutz-Delgado, K.

1991-01-01

A recently developed spatial operator algebra for manipulator modeling, control, and trajectory design is discussed. The elements of this algebra are linear operators whose domain and range spaces consist of forces, moments, velocities, and accelerations. The effect of these operators is equivalent to a spatial recursion along the span of a manipulator. Inversion of operators can be efficiently obtained via techniques of recursive filtering and smoothing. The operator algebra provides a high-level framework for describing the dynamic and kinematic behavior of a manipulator and for control and trajectory design algorithms. The interpretation of expressions within the algebraic framework leads to enhanced conceptual and physical understanding of manipulator dynamics and kinematics.

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.

Algebraic approach to small-world network models

NASA Astrophysics Data System (ADS)

Rudolph-Lilith, Michelle; Muller, Lyle E.

2014-01-01

We introduce an analytic model for directed Watts-Strogatz small-world graphs and deduce an algebraic expression of its defining adjacency matrix. The latter is then used to calculate the small-world digraph's asymmetry index and clustering coefficient in an analytically exact fashion, valid nonasymptotically for all graph sizes. The proposed approach is general and can be applied to all algebraically well-defined graph-theoretical measures, thus allowing for an analytical investigation of finite-size small-world graphs.

Computation of turbulent reacting flow in a solid-propellant ducted rocket

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

Chao, Y.; Chou, W.; Liu, S.

1995-05-01

A mathematical model for computation of turbulent reacting flows is developed under general curvilinear coordinate systems. An adaptive, streamline grid system is generated to deal with the complex flow structures in a multiple-inlet solid-propellant ducted rocket (SDR) combustor. General tensor representations of the k-epsilon and algebraic stress (ASM) turbulence models are derived in terms of contravariant velocity components, and modification caused by the effects of compressible turbulence is also included in the modeling. The clipped Gaussian probability density function is incorporated in the combustion model to account for fluctuations of properties. Validation of the above modeling is first examined bymore » studying mixing and reacting characteristics in a confined coaxial-jet problem. This is followed by study of nonreacting and reacting SDR combustor flows. The results show that Gibson and Launder`s ASM incorporated with Sarkar`s modification for compressible turbulence effects based on the general curvilinear coordinate systems yields the most satisfactory prediction for this complicated SDR flowfield. 36 refs.« less

Shapes and stability of algebraic nuclear models

NASA Technical Reports Server (NTRS)

Lopez-Moreno, Enrique; Castanos, Octavio

1995-01-01

A generalization of the procedure to study shapes and stability of algebraic nuclear models introduced by Gilmore is presented. One calculates the expectation value of the Hamiltonian with respect to the coherent states of the algebraic structure of the system. Then equilibrium configurations of the resulting energy surface, which depends in general on state variables and a set of parameters, are classified through the Catastrophe theory. For one- and two-body interactions in the Hamiltonian of the interacting Boson model-1, the critical points are organized through the Cusp catastrophe. As an example, we apply this Separatrix to describe the energy surfaces associated to the Rutenium and Samarium isotopes.

Numerical Simulation of High-Speed Turbulent Reacting Flows

NASA Technical Reports Server (NTRS)

Givi, P.; Taulbee, D. B.; Madnia, C. K.; Jaberi, F. A.; Colucci, P. J.; Gicquel, L. Y. M.; Adumitroaie, V.; James, S.

1999-01-01

The objectives of this research are: (1) to develop and implement a new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. (2) To develop algebraic turbulence closures for statistical description of chemically reacting turbulent flows. We have just completed the third year of Phase III of this research. This is the Final Report of our activities on this research sponsored by the NASA LaRC.

The lagRST Model: A Turbulence Model for Non-Equilibrium Flows

NASA Technical Reports Server (NTRS)

Lillard, Randolph P.; Oliver, A. Brandon; Olsen, Michael E.; Blaisdell, Gregory A.; Lyrintzis, Anastasios S.

2011-01-01

This study presents a new class of turbulence model designed for wall bounded, high Reynolds number flows with separation. The model addresses deficiencies seen in the modeling of nonequilibrium turbulent flows. These flows generally have variable adverse pressure gradients which cause the turbulent quantities to react at a finite rate to changes in the mean flow quantities. This "lag" in the response of the turbulent quantities can t be modeled by most standard turbulence models, which are designed to model equilibrium turbulent boundary layers. The model presented uses a standard 2-equation model as the baseline for turbulent equilibrium calculations, but adds transport equations to account directly for non-equilibrium effects in the Reynolds Stress Tensor (RST) that are seen in large pressure gradients involving shock waves and separation. Comparisons are made to several standard turbulence modeling validation cases, including an incompressible boundary layer (both neutral and adverse pressure gradients), an incompressible mixing layer and a transonic bump flow. In addition, a hypersonic Shock Wave Turbulent Boundary Layer Interaction with separation is assessed along with a transonic capsule flow. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWTBLI flows 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.

Modeling near-wall turbulent flows

NASA Astrophysics Data System (ADS)

Marusic, Ivan; Mathis, Romain; Hutchins, Nicholas

2010-11-01

The near-wall region of turbulent boundary layers is a crucial region for turbulence production, but it is also a region that becomes increasing difficult to access and make measurements in as the Reynolds number becomes very high. Consequently, it is desirable to model the turbulence in this region. Recent studies have shown that the classical description, with inner (wall) scaling alone, is insufficient to explain the behaviour of the streamwise turbulence intensities with increasing Reynolds number. Here we will review our recent near-wall model (Marusic et al., Science 329, 2010), where the near-wall turbulence is predicted given information from only the large-scale signature at a single measurement point in the logarithmic layer, considerably far from the wall. The model is consistent with the Townsend attached eddy hypothesis in that the large-scale structures associated with the log-region are felt all the way down to the wall, but also includes a non-linear amplitude modulation effect of the large structures on the near-wall turbulence. Detailed predicted spectra across the entire near- wall region will be presented, together with other higher order statistics over a large range of Reynolds numbers varying from laboratory to atmospheric flows.

Current algebra, statistical mechanics and quantum models

NASA Astrophysics Data System (ADS)

Vilela Mendes, R.

2017-11-01

Results obtained in the past for free boson systems at zero and nonzero temperatures are revisited to clarify the physical meaning of current algebra reducible functionals which are associated to systems with density fluctuations, leading to observable effects on phase transitions. To use current algebra as a tool for the formulation of quantum statistical mechanics amounts to the construction of unitary representations of diffeomorphism groups. Two mathematical equivalent procedures exist for this purpose. One searches for quasi-invariant measures on configuration spaces, the other for a cyclic vector in Hilbert space. Here, one argues that the second approach is closer to the physical intuition when modelling complex systems. An example of application of the current algebra methodology to the pairing phenomenon in two-dimensional fermion systems is discussed.

Approximate Model for Turbulent Stagnation Point Flow.

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

Dechant, Lawrence

2016-01-01

Here we derive an approximate turbulent self-similar model for a class of favorable pressure gradient wedge-like flows, focusing on the stagnation point limit. While the self-similar model provides a useful gross flow field estimate this approach must be combined with a near wall model is to determine skin friction and by Reynolds analogy the heat transfer coefficient. The combined approach is developed in detail for the stagnation point flow problem where turbulent skin friction and Nusselt number results are obtained. Comparison to the classical Van Driest (1958) result suggests overall reasonable agreement. Though the model is only valid near themore » stagnation region of cylinders and spheres it nonetheless provides a reasonable model for overall cylinder and sphere heat transfer. The enhancement effect of free stream turbulence upon the laminar flow is used to derive a similar expression which is valid for turbulent flow. Examination of free stream enhanced laminar flow suggests that the rather than enhancement of a laminar flow behavior free stream disturbance results in early transition to turbulent stagnation point behavior. Excellent agreement is shown between enhanced laminar flow and turbulent flow behavior for high levels, e.g. 5% of free stream turbulence. Finally the blunt body turbulent stagnation results are shown to provide realistic heat transfer results for turbulent jet impingement problems.« less

The dynamics of turbulent premixed flames: Mechanisms and models for turbulence-flame interaction

NASA Astrophysics Data System (ADS)

Steinberg, Adam M.

The use of turbulent premixed combustion in engines has been garnering renewed interest due to its potential to reduce NOx emissions. However there are many aspects of turbulence-flame interaction that must be better understood before such flames can be accurately modeled. The focus of this dissertation is to develop an improved understanding for the manner in which turbulence interacts with a premixed flame in the 'thin flamelet regime'. To do so, two new diagnostics were developed and employed in a turbulent slot Bunsen flame. These diagnostics, Cinema-Stereoscopic Particle Image Velocimetry and Orthogonal-Plane Cinema-Stereoscopic Particle Image Velocimetry, provided temporally resolved velocity and flame surface measurements in two- and three-dimensions with rates of up to 3 kHz and spatial resolutions as low as 280 mum. Using these measurements, the mechanisms with which turbulence generates flame surface area were studied. It was found that the previous concept that flame stretch is characterized by counter-rotating vortex pairs does not accurately describe real turbulence-flame interactions. Analysis of the experimental data showed that the straining of the flame surface is determined by coherent structures of fluid dynamic strain rate, while the wrinkling is caused by vortical structures. Furthermore, it was shown that the canonical vortex pair configuration is not an accurate reflection of the real interaction geometry. Hence, models developed based on this geometry are unlikely to be accurate. Previous models for the strain rate, curvature stretch rate, and turbulent burning velocity were evaluated. It was found that the previous models did not accurately predict the measured data for a variety of reasons: the assumed interaction geometries did not encompass enough possibilities to describe the possible effects of real turbulence, the turbulence was not properly characterized, and the transport of flame surface area was not always considered. New models

PDF turbulence modeling and DNS

NASA Technical Reports Server (NTRS)

Hsu, A. T.

1992-01-01

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

The Standard Model Algebra - a summary

NASA Astrophysics Data System (ADS)

Cristinel Stoica, Ovidiu

2017-08-01

A generation of leptons and quarks and the gauge symmetries of the Standard Model can be obtained from the Clifford algebra ℂℓ 6. An instance of ℂℓ 6 is implicitly generated by the Dirac algebra combined with the electroweak symmetry, while the color symmetry gives another instance of ℂℓ 6 with a Witt decomposition. The minimal mathematical model proposed here results by identifying the two instances of ℂℓ 6. The left ideal decomposition generated by the Witt decomposition represents the leptons and quarks, and their antiparticles. The SU(3)c and U(1)em symmetries of the SM are the symmetries of this ideal decomposition. The patterns of electric charges, colors, chirality, weak isospins, and hypercharges, follow from this, without predicting additional particles or forces, or proton decay. The electroweak symmetry is present in its broken form, due to the geometry. The predicted Weinberg angle is given by sin2 W = 0.25. The model shares common features with previously known models, particularly with Chisholm and Farwell, 1996, Trayling and Baylis, 2004, and Furey, 2016.

An Algebraic Approach to the Eigenstates of the Calogero Model

NASA Astrophysics Data System (ADS)

Ujino, Hideaki

2002-11-01

An algebraic treatment of the eigenstates of the (AN-1-) Calogero model is presented, which provides an algebraic construction of the nonsymmetric orthogonal eigenvectors, symmetrization, antisymmetrization and calculation of square norms in a unified way.

Two-fluid models of turbulence

NASA Technical Reports Server (NTRS)

Spalding, D. B.

1985-01-01

The defects of turbulence models are summarized and the importance of so-called nongradient diffusion in turbulent fluxes is discussed. The mathematical theory of the flow of two interpenetrating continua is reviewed, and the mathematical formulation of the two fluid model is outlined. Results from plane wake, axisymmetric jet, and combustion studies are shown.

Laminar flamelet modeling of turbulent diffusion flames

NASA Technical Reports Server (NTRS)

Mell, W. E.; Kosaly, G.; Planche, O.; Poinsot, T.; Ferziger, J. H.

1990-01-01

In modeling turbulent combustion, decoupling the chemistry from the turbulence is of great practical significance. In cases in which the equilibrium chemistry model breaks down, laminar flamelet modeling (LFM) is a promising approach to decoupling. Here, the validity of this approach is investigated using direct numerical simulation of a simple chemical reaction in two-dimensional turbulence.

Kolmogorov Behavior of Near-Wall Turbulence and Its Application in Turbulence Modeling

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Lumley, John L.

1992-01-01

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder and followers. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterized by Kolmogorov microscales. According to this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become Kolmogorov eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As an example, the concept is incorporated in the standard k-epsilon model which is then applied to channel and boundary flows. Using appropriate boundary conditions (based on Kolmogorov behavior of near-wall turbulence), there is no need for any wall-modification to the k-epsilon equations (including model constants). Results compare very well with the DNS and experimental data.

A spatial operator algebra for manipulator modeling and control

NASA Technical Reports Server (NTRS)

Rodriguez, G.; Kreutz, Kenneth; Jain, Abhinandan

1989-01-01

A recently developed spatial operator algebra, useful for modeling, control, and trajectory design of manipulators is discussed. The elements of this algebra are linear operators whose domain and range spaces consist of forces, moments, velocities, and accelerations. The effect of these operators is equivalent to a spatial recursion along the span of a manipulator. Inversion of operators can be efficiently obtained via techniques of recursive filtering and smoothing. The operator algebra provides a high level framework for describing the dynamic and kinematic behavior of a manipulator and control and trajectory design algorithms. The interpretation of expressions within the algebraic framework leads to enhanced conceptual and physical understanding of manipulator dynamics and kinematics. Furthermore, implementable recursive algorithms can be immediately derived from the abstract operator expressions by inspection. Thus, the transition from an abstract problem formulation and solution to the detailed mechanizaton of specific algorithms is greatly simplified. The analytical formulation of the operator algebra, as well as its implementation in the Ada programming language are discussed.

A k-epsilon modeling of near wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1991-01-01

A k-epsilon model is proposed for turbulent bounded flows. In this model, the turbulent velocity scale and turbulent time scale are used to define the eddy viscosity. The time scale is shown to be bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using the time scale, removing the need to introduce the pseudo-dissipation. A damping function is chosen such that the shear stress satisfies the near wall asymptotic behavior. The model constants used are the same as the model constants in the commonly used high turbulent Reynolds number k-epsilon model. Fully developed turbulent channel flows and turbulent boundary layer flows over a flat plate at various Reynolds numbers are used to validate the model. The model predictions were found to be in good agreement with the direct numerical simulation data.

Collaborative testing of turbulence models

NASA Technical Reports Server (NTRS)

Bradshaw, Peter; Launder, Brian E.; Lumley, John L.

1991-01-01

A review is given of an ongoing international project, in which data from experiments on, and simulations of, turbulent flows are distributed to developers of (time-averaged) engineering turbulence models. The predictions of each model are sent to the organizers and redistributed to all the modelers, plus some experimentalists and other experts (total approx. 120), for comment. The 'reaction time' of modelers has proved to be much longer than anticipated, partly because the comparisons with data have prompted many modelers to improve their models or numerics.

Turbulence modeling

NASA Technical Reports Server (NTRS)

Bardina, Jorge E.

1995-01-01

The objective of this work is to develop, verify, and incorporate the baseline two-equation turbulence models which account for the effects of compressibility into the three-dimensional Reynolds averaged Navier-Stokes (RANS) code and to provide documented descriptions of the models and their numerical procedures so that they can be implemented into 3-D CFD codes for engineering applications.

Evaluation and application of the Baldwin-Lomax turbulence model in two-dimensional, unsteady, compressible boundary layers with and without separation in engine inlets

NASA Technical Reports Server (NTRS)

Sakowski, Barbara; Darling, Douglas; Roach, Robert L.; Vandewall, Allan

1992-01-01

There is a practical need to model high speed flows that exist in jet engine inlets. The boundary layers that form in these inlets may be turbulent or laminar and either separated or attached. Also, unsteady supersonic inlets may be subject to frequent changes in operating conditions. Some changes in the operating conditions of the inlets may include varying the inlet geometry, bleeds and bypasses, and rotating or translating the centerbody. In addition, the inlet may be either started or unstarted. Therefore, a CFD code, used to model these inlets, may have to run for several different cases. Also, since the flow conditions through an unsteady inlet may be continually fluctuating, the CFD code which models these flows may have to be run over many time steps. Therefore, it would be beneficial that the code run quickly. Many turbulence models, however, are cumbersome to implement and require a lot of computer time to run, since they add to the number of differential equations to be solved to model a flow. The Baldwin-Lomax turbulence model is a popular model. It is an algebraic, eddy viscosity model. The Baldwin-Lomax model is used in many CFD codes because it is quick and easy to implement. In this paper, we will discuss implementing the Baldwin-Lomax turbulence model for both steady and unsteady compressible flows. In addition, these flows may be either separated or attached. In order to apply this turbulence model to flows which may be subjected to these conditions, certain modifications should be made to the original Baldwin-Lomax model. We will discuss these modifications and determine whether the Baldwin-Lomax model is a viable turbulence model that produces reasonably accurate results for high speed flows that can be found in engine inlets.

A representation for the turbulent mass flux contribution to Reynolds-stress and two-equation closures for compressible turbulence

NASA Technical Reports Server (NTRS)

Ristorcelli, J. R.

1993-01-01

The turbulent mass flux, or equivalently the fluctuating Favre velocity mean, appears in the first and second moment equations of compressible kappa-epsilon and Reynolds stress closures. Mathematically it is the difference between the unweighted and density-weighted averages of the velocity field and is therefore a measure of the effects of compressibility through variations in density. It appears to be fundamental to an inhomogeneous compressible turbulence, in which it characterizes the effects of the mean density gradients, in the same way the anisotropy tensor characterizes the effects of the mean velocity gradients. An evolution equation for the turbulent mass flux is derived. A truncation of this equation produces an algebraic expression for the mass flux. The mass flux is found to be proportional to the mean density gradients with a tensor eddy-viscosity that depends on both the mean deformation and the Reynolds stresses. The model is tested in a wall bounded DNS at Mach 4.5 with notable results.

High Reynolds number analysis of flat plate and separated afterbody flow using non-linear turbulence models

NASA Technical Reports Server (NTRS)

Carlson, John R.

1996-01-01

The ability of the three-dimensional Navier-Stokes method, PAB3D, to simulate the effect of Reynolds number variation using non-linear explicit algebraic Reynolds stress turbulence modeling was assessed. Subsonic flat plate boundary-layer flow parameters such as normalized velocity distributions, local and average skin friction, and shape factor were compared with DNS calculations and classical theory at various local Reynolds numbers up to 180 million. Additionally, surface pressure coefficient distributions and integrated drag predictions on an axisymmetric nozzle afterbody were compared with experimental data from 10 to 130 million Reynolds number. The high Reynolds data was obtained from the NASA Langley 0.3m Transonic Cryogenic Tunnel. There was generally good agreement of surface static pressure coefficients between the CFD and measurement. The change in pressure coefficient distributions with varying Reynolds number was similar to the experimental data trends, though slightly over-predicting the effect. The computational sensitivity of viscous modeling and turbulence modeling are shown. Integrated afterbody pressure drag was typically slightly lower than the experimental data. The change in afterbody pressure drag with Reynolds number was small both experimentally and computationally, even though the shape of the distribution was somewhat modified with Reynolds number.

Development of a One-Equation Eddy Viscosity Turbulence Model for Application to Complex Turbulent Flows

NASA Astrophysics Data System (ADS)

Wray, Timothy J.

Computational fluid dynamics (CFD) is routinely used in performance prediction and design of aircraft, turbomachinery, automobiles, and in many other industrial applications. Despite its wide range of use, deficiencies in its prediction accuracy still exist. One critical weakness is the accurate simulation of complex turbulent flows using the Reynolds-Averaged Navier-Stokes equations in conjunction with a turbulence model. The goal of this research has been to develop an eddy viscosity type turbulence model to increase the accuracy of flow simulations for mildly separated flows, flows with rotation and curvature effects, and flows with surface roughness. It is accomplished by developing a new zonal one-equation turbulence model which relies heavily on the flow physics; it is now known in the literature as the Wray-Agarwal one-equation turbulence model. The effectiveness of the new model is demonstrated by comparing its results with those obtained by the industry standard one-equation Spalart-Allmaras model and two-equation Shear-Stress-Transport k - o model and experimental data. Results for subsonic, transonic, and supersonic flows in and about complex geometries are presented. It is demonstrated that the Wray-Agarwal model can provide the industry and CFD researchers an accurate, efficient, and reliable turbulence model for the computation of a large class of complex turbulent flows.

Structure and modeling of turbulence

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

Novikov, E.A.

The {open_quotes}vortex strings{close_quotes} scale l{sub s} {approximately} LRe{sup -3/10} (L-external scale, Re - Reynolds number) is suggested as a grid scale for the large-eddy simulation. Various aspects of the structure of turbulence and subgrid modeling are described in terms of conditional averaging, Markov processes with dependent increments and infinitely divisible distributions. The major request from the energy, naval, aerospace and environmental engineering communities to the theory of turbulence is to reduce the enormous number of degrees of freedom in turbulent flows to a level manageable by computer simulations. The vast majority of these degrees of freedom is in the small-scalemore » motion. The study of the structure of turbulence provides a basis for subgrid-scale (SGS) models, which are necessary for the large-eddy simulations (LES).« less

Polynomial algebra of discrete models in systems biology.

PubMed

Veliz-Cuba, Alan; Jarrah, Abdul Salam; Laubenbacher, Reinhard

2010-07-01

An increasing number of discrete mathematical models are being published in Systems Biology, ranging from Boolean network models to logical models and Petri nets. They are used to model a variety of biochemical networks, such as metabolic networks, gene regulatory networks and signal transduction networks. There is increasing evidence that such models can capture key dynamic features of biological networks and can be used successfully for hypothesis generation. This article provides a unified framework that can aid the mathematical analysis of Boolean network models, logical models and Petri nets. They can be represented as polynomial dynamical systems, which allows the use of a variety of mathematical tools from computer algebra for their analysis. Algorithms are presented for the translation into polynomial dynamical systems. Examples are given of how polynomial algebra can be used for the model analysis. alanavc@vt.edu Supplementary data are available at Bioinformatics online.

Philosophies and fallacies in turbulence modeling

NASA Astrophysics Data System (ADS)

Spalart, Philippe R.

2015-04-01

We present a set of positions, likely to be controversial, on turbulence modeling for the Reynolds-Averaged Navier Stokes (RANS) equations. The paper has three themes. First is what we call the "fundamental paradox" of turbulence modeling, between the local character of the Partial Differential Equations strongly favored by CFD methods and the nonlocal physical nature of turbulence. Second, we oppose two philosophies. The "Systematic" philosophy attempts to model the exact transport equations for the Reynolds stresses or possibly higher moments term by term, gradually relegating the Closure Problem to higher moments and invoking the "Principle of Receding Influence" (although rarely formulating it). In contrast, the "Openly Empirical" philosophy produces models which satisfy strict constraints such as Galilean invariance, but lack an explicit connection with terms in the exact turbulence equations. The prime example is the eddy-viscosity assumption. Third, we explain a series of what we perceive as fallacies, many of them widely held and by senior observers, in turbulence knowledge, leading to turbulence models. We divide them into "hard" fallacies for which a short mathematical argument demonstrates that a particular statement is wrong or meaningless, and "soft" fallacies for which approximate physical arguments can be opposed, but we contend that a clear debate is overdue and wishful thinking has been involved. Some fallacies appear to be "intermediate." An example in the hard class is the supposed isotropy of the diagonal Reynolds stresses. Examples in the soft class are the need to match the decay rate of isotropic turbulence, and the value of realizability in a model. Our hope is to help the direct effort in this field away from simplistic and hopeless lines of work, and to foster debates.

A multiple-scale turbulence model for incompressible flow

NASA Technical Reports Server (NTRS)

Duncan, B. S.; Liou, W. W.; Shih, T. H.

1993-01-01

A multiple-scale eddy viscosity model is described. This model splits the energy spectrum into a high wave number regime and a low wave number regime. Dividing the energy spectrum into multiple regimes simplistically emulates the cascade of energy through the turbulence spectrum. The constraints on the model coefficients are determined by examining decaying turbulence and homogeneous turbulence. A direct link between the partitioned energies and the energy transfer process is established through the coefficients. This new model was calibrated and tested for boundary-free turbulent shear flows. Calculations of mean and turbulent properties show good agreement with experimental data for two mixing layers, a plane jet and a round jet.

Modeling of Turbulence Effect on Liquid Jet Atomization

NASA Technical Reports Server (NTRS)

Trinh, H. P.

2007-01-01

Recent studies indicate that turbulence behaviors within a liquid jet have considerable effect on the atomization process. Such turbulent flow phenomena are encountered in most practical applications of common liquid spray devices. This research aims to model the effects of turbulence occurring inside a cylindrical liquid jet to its atomization process. The two widely used atomization models Kelvin-Helmholtz (KH) instability of Reitz and the Taylor analogy breakup (TAB) of O'Rourke and Amsden portraying primary liquid jet disintegration and secondary droplet breakup, respectively, are examined. Additional terms are formulated and appropriately implemented into these two models to account for the turbulence effect. Results for the flow conditions examined in this study indicate that the turbulence terms are significant in comparison with other terms in the models. In the primary breakup regime, the turbulent liquid jet tends to break up into large drops while its intact core is slightly shorter than those without turbulence. In contrast, the secondary droplet breakup with the inside liquid turbulence consideration produces smaller drops. Computational results indicate that the proposed models provide predictions that agree reasonably well with available measured data.

A spatial operator algebra for manipulator modeling and control

NASA Technical Reports Server (NTRS)

Rodriguez, G.; Kreutz, K.; Milman, M.

1988-01-01

A powerful new spatial operator algebra for modeling, control, and trajectory design of manipulators is discussed along with its implementation in the Ada programming language. Applications of this algebra to robotics include an operator representation of the manipulator Jacobian matrix; the robot dynamical equations formulated in terms of the spatial algebra, showing the complete equivalence between the recursive Newton-Euler formulations to robot dynamics; the operator factorization and inversion of the manipulator mass matrix which immediately results in O(N) recursive forward dynamics algorithms; the joint accelerations of a manipulator due to a tip contact force; the recursive computation of the equivalent mass matrix as seen at the tip of a manipulator; and recursive forward dynamics of a closed chain system. Finally, additional applications and current research involving the use of the spatial operator algebra are discussed in general terms.

Uncertainty Assessments of 2D and Axisymmetric Hypersonic Shock Wave - Turbulent Boundary Layer Interaction Simulations at Compression Corners

NASA Technical Reports Server (NTRS)

Gnoffo, Peter A.; Berry, Scott A.; VanNorman, John W.

2011-01-01

This paper is one of a series of five papers in a special session organized by the NASA Fundamental Aeronautics Program that addresses uncertainty assessments for CFD simulations in hypersonic flow. Simulations of a shock emanating from a compression corner and interacting with a fully developed turbulent boundary layer are evaluated herein. Mission relevant conditions at Mach 7 and Mach 14 are defined for a pre-compression ramp of a scramjet powered vehicle. Three compression angles are defined, the smallest to avoid separation losses and the largest to force a separated flow engaging more complicated flow physics. The Baldwin-Lomax and the Cebeci-Smith algebraic models, the one-equation Spalart-Allmaras model with the Catrix-Aupoix compressibility modification and two-equation models including Menter SST, Wilcox k-omega 98, and Wilcox k-omega 06 turbulence models are evaluated. Each model is fully defined herein to preclude any ambiguity regarding model implementation. Comparisons are made to existing experimental data and Van Driest theory to provide preliminary assessment of model form uncertainty. A set of coarse grained uncertainty metrics are defined to capture essential differences among turbulence models. Except for the inability of algebraic models to converge for some separated flows there is no clearly superior model as judged by these metrics. A preliminary metric for the numerical component of uncertainty in shock-turbulent-boundary-layer interactions at compression corners sufficiently steep to cause separation is defined as 55%. This value is a median of differences with experimental data averaged for peak pressure and heating and for extent of separation captured in new, grid-converged solutions presented here. This value is consistent with existing results in a literature review of hypersonic shock-turbulent-boundary-layer interactions by Roy and Blottner and with more recent computations of MacLean.

An experimental and numerical investigation of shock-wave induced turbulent boundary-layer separation at hypersonic speeds

NASA Technical Reports Server (NTRS)

Marvin, J. G.; Horstman, C. C.; Rubesin, M. W.; Coakley, T. J.; Kussoy, M. I.

1975-01-01

An experiment designed to test and guide computations of the interaction of an impinging shock wave with a turbulent boundary layer is described. Detailed mean flow-field and surface data are presented for two shock strengths which resulted in attached and separated flows, respectively. Numerical computations, employing the complete time-averaged Navier-Stokes equations along with algebraic eddy-viscosity and turbulent Prandtl number models to describe shear stress and heat flux, are used to illustrate the dependence of the computations on the particulars of the turbulence models. Models appropriate for zero-pressure-gradient flows predicted the overall features of the flow fields, but were deficient in predicting many of the details of the interaction regions. Improvements to the turbulence model parameters were sought through a combination of detailed data analysis and computer simulations which tested the sensitivity of the solutions to model parameter changes. Computer simulations using these improvements are presented and discussed.

A Structural Model of Algebra Achievement: Computational Fluency and Spatial Visualisation as Mediators of the Effect of Working Memory on Algebra Achievement

ERIC Educational Resources Information Center

Tolar, Tammy Daun; Lederberg, Amy R.; Fletcher, Jack M.

2009-01-01

The goal of this study was to develop and evaluate a structural model of the relations among cognitive abilities and arithmetic skills and college students' algebra achievement. The model of algebra achievement was compared to a model of performance on the Scholastic Assessment in Mathematics (SAT-M) to determine whether the pattern of relations…

Wall-resolved spectral cascade-transport turbulence model

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

Brown, C. S.; Shaver, D. R.; Lahey, R. T.

A spectral cascade-transport model has been developed and applied to turbulent channel flows (Reτ= 550, 950, and 2000 based on friction velocity, uτ ; or ReδΜ= 8,500; 14,800 and 31,000, based on the mean velocity and channel half-width). This model is an extension of a spectral model previously developed for homogeneous single and two-phase decay of isotropic turbulence and uniform shear flows; and a spectral turbulence model for wall-bounded flows without resolving the boundary layer. Data from direct numerical simulation (DNS) of turbulent channel flow was used to help develop this model and to assess its performance in the 1Dmore » direction across the channel width. The resultant spectral model is capable of predicting the mean velocity, turbulent kinetic energy and energy spectrum distributions for single-phase wall-bounded flows all the way to the wall, where the model source terms have been developed to account for the wall influence. We implemented the model into the 3D multiphase CFD code NPHASE-CMFD and the latest results are within reasonable error of the 1D predictions.« less

Wall-resolved spectral cascade-transport turbulence model

DOE PAGES

Brown, C. S.; Shaver, D. R.; Lahey, R. T.; ...

2017-07-08

A spectral cascade-transport model has been developed and applied to turbulent channel flows (Reτ= 550, 950, and 2000 based on friction velocity, uτ ; or ReδΜ= 8,500; 14,800 and 31,000, based on the mean velocity and channel half-width). This model is an extension of a spectral model previously developed for homogeneous single and two-phase decay of isotropic turbulence and uniform shear flows; and a spectral turbulence model for wall-bounded flows without resolving the boundary layer. Data from direct numerical simulation (DNS) of turbulent channel flow was used to help develop this model and to assess its performance in the 1Dmore » direction across the channel width. The resultant spectral model is capable of predicting the mean velocity, turbulent kinetic energy and energy spectrum distributions for single-phase wall-bounded flows all the way to the wall, where the model source terms have been developed to account for the wall influence. We implemented the model into the 3D multiphase CFD code NPHASE-CMFD and the latest results are within reasonable error of the 1D predictions.« less

Turbulent motion of mass flows. Mathematical modeling

NASA Astrophysics Data System (ADS)

Eglit, Margarita; Yakubenko, Alexander; Yakubenko, Tatiana

2016-04-01

New mathematical models for unsteady turbulent mass flows, e.g., dense snow avalanches and landslides, are presented. Such models are important since most of large scale flows are turbulent. In addition to turbulence, the two other important points are taken into account: the entrainment of the underlying material by the flow and the nonlinear rheology of moving material. The majority of existing models are based on the depth-averaged equations and the turbulent character of the flow is accounted by inclusion of drag proportional to the velocity squared. In this paper full (not depth-averaged) equations are used. It is assumed that basal entrainment takes place if the bed friction equals the shear strength of the underlying layer (Issler D, M. Pastor Peréz. 2011). The turbulent characteristics of the flow are calculated using a three-parameter differential model (Lushchik et al., 1978). The rheological properties of moving material are modeled by one of the three types of equations: 1) Newtonian fluid with high viscosity, 2) power-law fluid and 3) Bingham fluid. Unsteady turbulent flows down long homogeneous slope are considered. The flow dynamical parameters and entrainment rate behavior in time as well as their dependence on properties of moving and underlying materials are studied numerically. REFERENCES M.E. Eglit and A.E. Yakubenko, 2014. Numerical modeling of slope flows entraining bottom material. Cold Reg. Sci. Technol., 108, 139-148 Margarita E. Eglit and Alexander E. Yakubenko, 2016. The effect of bed material entrainment and non-Newtonian rheology on dynamics of turbulent slope flows. Fluid Dynamics, 51(3) Issler D, M. Pastor Peréz. 2011. Interplay of entrainment and rheology in snow avalanches; a numerical study. Annals of Glaciology, 52(58), 143-147 Lushchik, V.G., Paveliev, A.A. , and Yakubenko, A.E., 1978. Three-parameter model of shear turbulence. Fluid Dynamics, 13, (3), 350-362

Algebraic K-theory, K-regularity, and -duality of -stable C ∗-algebras

NASA Astrophysics Data System (ADS)

Mahanta, Snigdhayan

2015-12-01

We develop an algebraic formalism for topological -duality. More precisely, we show that topological -duality actually induces an isomorphism between noncommutative motives that in turn implements the well-known isomorphism between twisted K-theories (up to a shift). In order to establish this result we model topological K-theory by algebraic K-theory. We also construct an E ∞ -operad starting from any strongly self-absorbing C ∗-algebra . Then we show that there is a functorial topological K-theory symmetric spectrum construction on the category of separable C ∗-algebras, such that is an algebra over this operad; moreover, is a module over this algebra. Along the way we obtain a new symmetric spectra valued functorial model for the (connective) topological K-theory of C ∗-algebras. We also show that -stable C ∗-algebras are K-regular providing evidence for a conjecture of Rosenberg. We conclude with an explicit description of the algebraic K-theory of a x+ b-semigroup C ∗-algebras coming from number theory and that of -stabilized noncommutative tori.

A multiple-scale turbulence model for incompressible flow

NASA Technical Reports Server (NTRS)

Duncan, B. S.; Liou, W. W.; Shih, T. H.

1993-01-01

A multiple-scale eddy viscosity model is described in this paper. This model splits the energy spectrum into a high wave number regime and a low wave number regime. Dividing the energy spectrum into multiple regimes simplistically emulates the cascade of energy through the turbulence spectrum. The constraints on the model coefficients are determined by examining decaying turbulence and homogeneous turbulence. A direct link between the partitioned energies and the energy transfer process is established through the coefficients. This new model has been calibrated and tested for boundary-free turbulent shear flows. Calculations of mean and turbulent properties show good agreement with experimental data for two mixing layers, a plane jet and a round jet.

An algebraic cluster model based on the harmonic oscillator basis

NASA Technical Reports Server (NTRS)

Levai, Geza; Cseh, J.

1995-01-01

We discuss the semimicroscopic algebraic cluster model introduced recently, in which the internal structure of the nuclear clusters is described by the harmonic oscillator shell model, while their relative motion is accounted for by the Vibron model. The algebraic formulation of the model makes extensive use of techniques associated with harmonic oscillators and their symmetry group, SU(3). The model is applied to some cluster systems and is found to reproduce important characteristics of nuclei in the sd-shell region. An approximate SU(3) dynamical symmetry is also found to hold for the C-12 + C-12 system.

Gyrofluid Modeling of Turbulent, Kinetic Physics

NASA Astrophysics Data System (ADS)

Despain, Kate Marie

2011-12-01

Gyrofluid models to describe plasma turbulence combine the advantages of fluid models, such as lower dimensionality and well-developed intuition, with those of gyrokinetics models, such as finite Larmor radius (FLR) effects. This allows gyrofluid models to be more tractable computationally while still capturing much of the physics related to the FLR of the particles. We present a gyrofluid model derived to capture the behavior of slow solar wind turbulence and describe the computer code developed to implement the model. In addition, we describe the modifications we made to a gyrofluid model and code that simulate plasma turbulence in tokamak geometries. Specifically, we describe a nonlinear phase mixing phenomenon, part of the E x B term, that was previously missing from the model. An inherently FLR effect, it plays an important role in predicting turbulent heat flux and diffusivity levels for the plasma. We demonstrate this importance by comparing results from the updated code to studies done previously by gyrofluid and gyrokinetic codes. We further explain what would be necessary to couple the updated gyrofluid code, gryffin, to a turbulent transport code, thus allowing gryffin to play a role in predicting profiles for fusion devices such as ITER and to explore novel fusion configurations. Such a coupling would require the use of Graphical Processing Units (GPUs) to make the modeling process fast enough to be viable. Consequently, we also describe our experience with GPU computing and demonstrate that we are poised to complete a gryffin port to this innovative architecture.

Turbulence modeling in three-dimensional stenosed arterial bifurcations.

PubMed

Banks, J; Bressloff, N W

2007-02-01

Under normal healthy conditions, blood flow in the carotid artery bifurcation is laminar. However, in the presence of a stenosis, the flow can become turbulent at the higher Reynolds numbers during systole. There is growing consensus that the transitional k-omega model is the best suited Reynolds averaged turbulence model for such flows. Further confirmation of this opinion is presented here by a comparison with the RNG k-epsilon model for the flow through a straight, nonbifurcating tube. Unlike similar validation studies elsewhere, no assumptions are made about the inlet profile since the full length of the experimental tube is simulated. Additionally, variations in the inflow turbulence quantities are shown to have no noticeable affect on downstream turbulence intensity, turbulent viscosity, or velocity in the k-epsilon model, whereas the velocity profiles in the transitional k-omega model show some differences due to large variations in the downstream turbulence quantities. Following this validation study, the transitional k-omega model is applied in a three-dimensional parametrically defined computer model of the carotid artery bifurcation in which the sinus bulb is manipulated to produce mild, moderate, and severe stenosis. The parametric geometry definition facilitates a powerful means for investigating the effect of local shape variation while keeping the global shape fixed. While turbulence levels are generally low in all cases considered, the mild stenosis model produces higher levels of turbulent viscosity and this is linked to relatively high values of turbulent kinetic energy and low values of the specific dissipation rate. The severe stenosis model displays stronger recirculation in the flow field with higher values of vorticity, helicity, and negative wall shear stress. The mild and moderate stenosis configurations produce similar lower levels of vorticity and helicity.

Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures

NASA Technical Reports Server (NTRS)

Ameri, Ali A.; Arnone, Andrea

1992-01-01

Navier-Stokes calculations were carried out in order to predict the heat-transfer rates on turbine blades. The calculations were performed using TRAF2D which is a k-epsilon, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using Coakley's q-omega and Chien's k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a nonperiodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the SSME high-pressure fuel turbine. The comparison serves to highlight the weaknesses of the turbulence models for use in turbomachinery heat-transfer calculations.

Testing an algebraic model of self-reflexion.

PubMed

Grice, James W; McDaniel, Brenda L; Thompsen, Dana

2005-06-01

Self-reflexion is the conscious process of taking the position of an observer in relation to one's own thoughts, feelings, and experiences. Building on the work of Lefebvre, Lefebvre, and Adams-Webber, we used a formal algebraic model of self-reflexion to derive several predictions regarding the frequencies with which individuals would rate themselves and others positively on bipolar scales anchored by adjective terms. The current results from 108 participants (41 men, 67 women; M age= 20.2 yr.) confirmed two predictions derived from the model. Three other predictions, however, were not supported even though the observed frequencies were close to the predicted values. Although not as promising as results reported by Lefebvre, et al., these mixed findings were interpreted as encouraging support for the validity of Lefebvre's algebraic model of self-reflexion. Differences between the current methods and those from previous investigations were also examined, and methodological implications for further studies were discussed.

Turbulence modeling for compressible flows

NASA Technical Reports Server (NTRS)

Marvin, J. G.

1977-01-01

Material prepared for a course on Applications and Fundamentals of Turbulence given at the University of Tennessee Space Institute, January 10 and 11, 1977, is presented. A complete concept of turbulence modeling is described, and examples of progess for its use in computational aerodynimics are given. Modeling concepts, experiments, and computations using the concepts are reviewed in a manner that provides an up-to-date statement on the status of this problem for compressible flows.

Second-order closure models for supersonic turbulent flows

NASA Technical Reports Server (NTRS)

Speziale, Charles G.; Sarkar, Sutanu

1991-01-01

Recent work by the authors on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulence closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equation are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

Second-order closure models for supersonic turbulent flows

NASA Technical Reports Server (NTRS)

Speziale, Charles G.; Sarkar, Sutanu

1991-01-01

Recent work on the development of a second-order closure model for high-speed compressible flows is reviewed. This turbulent closure is based on the solution of modeled transport equations for the Favre-averaged Reynolds stress tensor and the solenoidal part of the turbulent dissipation rate. A new model for the compressible dissipation is used along with traditional gradient transport models for the Reynolds heat flux and mass flux terms. Consistent with simple asymptotic analyses, the deviatoric part of the remaining higher-order correlations in the Reynolds stress transport equations are modeled by a variable density extension of the newest incompressible models. The resulting second-order closure model is tested in a variety of compressible turbulent flows which include the decay of isotropic turbulence, homogeneous shear flow, the supersonic mixing layer, and the supersonic flat-plate turbulent boundary layer. Comparisons between the model predictions and the results of physical and numerical experiments are quite encouraging.

A LES-Langevin model for turbulence

NASA Astrophysics Data System (ADS)

Dolganov, Rostislav; Dubrulle, Bérengère; Laval, Jean-Philippe

2006-11-01

The rationale for Large Eddy Simulation is rooted in our inability to handle all degrees of freedom (N˜10^16 for Re˜10^7). ``Deterministic'' models based on eddy-viscosity seek to reproduce the intensification of the energy transport. However, they fail to reproduce backward energy transfer (backscatter) from small to large scale, which is an essentiel feature of the turbulence near wall or in boundary layer. To capture this backscatter, ``stochastic'' strategies have been developed. In the present talk, we shall discuss such a strategy, based on a Rapid Distorsion Theory (RDT). Specifically, we first divide the small scale contribution to the Reynolds Stress Tensor in two parts: a turbulent viscosity and the pseudo-Lamb vector, representing the nonlinear cross terms of resolved and sub-grid scales. We then estimate the dynamics of small-scale motion by the RDT applied to Navier-Stockes equation. We use this to model the cross term evolution by a Langevin equation, in which the random force is provided by sub-grid pressure terms. Our LES model is thus made of a truncated Navier-Stockes equation including the turbulent force and a generalized Langevin equation for the latter, integrated on a twice-finer grid. The backscatter is automatically included in our stochastic model of the pseudo-Lamb vector. We apply this model to the case of homogeneous isotropic turbulence and turbulent channel flow.

Turbulence modeling for hypersonic flight

NASA Technical Reports Server (NTRS)

Bardina, Jorge E.

1992-01-01

The objective of the present work is to develop, verify, and incorporate two equation turbulence models which account for the effect of compressibility at high speeds into a three dimensional Reynolds averaged Navier-Stokes code and to provide documented model descriptions and numerical procedures so that they can be implemented into the National Aerospace Plane (NASP) codes. A summary of accomplishments is listed: (1) Four codes have been tested and evaluated against a flat plate boundary layer flow and an external supersonic flow; (2) a code named RANS was chosen because of its speed, accuracy, and versatility; (3) the code was extended from thin boundary layer to full Navier-Stokes; (4) the K-omega two equation turbulence model has been implemented into the base code; (5) a 24 degree laminar compression corner flow has been simulated and compared to other numerical simulations; and (6) work is in progress in writing the numerical method of the base code including the turbulence model.

Turbulence and modeling in transonic flow

NASA Technical Reports Server (NTRS)

Rubesin, Morris W.; Viegas, John R.

1989-01-01

A review is made of the performance of a variety of turbulence models in the evaluation of a particular well documented transonic flow. This is done to supplement a previous attempt to calibrate and verify transonic airfoil codes by including many more turbulence models than used in the earlier work and applying the calculations to an experiment that did not suffer from uncertainties in angle of attack and was free of wind tunnel interference. It is found from this work, as well as in the earlier study, that the Johnson-King turbulence model is superior for transonic flows over simple aerodynamic surfaces, including moderate separation. It is also shown that some field equation models with wall function boundary conditions can be competitive with it.

Mimicking Natural Laminar to Turbulent Flow Transition: A Systematic CFD Study Using PAB3D

NASA Technical Reports Server (NTRS)

Pao, S. Paul; Abdol-Hamid, Khaled S.

2005-01-01

For applied aerodynamic computations using a general purpose Navier-Stokes code, the common practice of treating laminar to turbulent flow transition over a non-slip surface is somewhat arbitrary by either treating the entire flow as turbulent or forcing the flow to undergo transition at given trip locations in the computational domain. In this study, the possibility of using the PAB3D code, standard k-epsilon turbulence model, and the Girimaji explicit algebraic stresses model to mimic natural laminar to turbulent flow transition was explored. The sensitivity of flow transition with respect to two limiters in the standard k-epsilon turbulence model was examined using a flat plate and a 6:1 aspect ratio prolate spheroid for our computations. For the flat plate, a systematic dependence of transition Reynolds number on background turbulence intensity was found. For the prolate spheroid, the transition patterns in the three-dimensional boundary layer at different flow conditions were sensitive to the free stream turbulence viscosity limit, the reference Reynolds number and the angle of attack, but not to background turbulence intensity below a certain threshold value. The computed results showed encouraging agreements with the experimental measurements at the corresponding geometry and flow conditions.

Statistical modeling of compressible turbulence - Shock-wave/turbulence interactions and buoyancy effects

NASA Astrophysics Data System (ADS)

Yoshizawa, Akira

1991-12-01

A mass-weighted mean compressible turbulence model is presented with the aid of the results from a two-scale DIA. This model aims at dealing with two typical aspects in compressible flows: the interaction of a shock wave with turbulence in high-speed flows and strong buoyancy effects in thermally-driven flows as in stellar convection and conflagration. The former is taken into account through the effect of turbulent dilatation that is related to the density fluctuation and leads to the enhanced kinetic-energy dissipation. The latter is incorporated through the interaction between the gravitational and density-fluctuation effects.

Near-wall k-epsilon turbulence modeling

NASA Technical Reports Server (NTRS)

Mansour, N. N.; Kim, J.; Moin, P.

1987-01-01

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

A turbulence model for pulsatile arterial flows.

PubMed

Younis, B A; Berger, S A

2004-10-01

Difficulties in predicting the behavior of some high Reynolds number flows in the circulatory system stem in part from the severe requirements placed on the turbulence model chosen to close the time-averaged equations of fluid motion. In particular, the successful turbulence model is required to (a) correctly capture the "nonequilibrium" effects wrought by the interactions of the organized mean-flow unsteadiness with the random turbulence, (b) correctly reproduce the effects of the laminar-turbulent transitional behavior that occurs at various phases of the cardiac cycle, and (c) yield good predictions of the near-wall flow behavior in conditions where the universal logarithmic law of the wall is known to be not valid. These requirements are not immediately met by standard models of turbulence that have been developed largely with reference to data from steady, fully turbulent flows in approximate local equilibrium. The purpose of this paper is to report on the development of a turbulence model suited for use in arterial flows. The model is of the two-equation eddy-viscosity variety with dependent variables that are zero-valued at a solid wall and vary linearly with distance from it. The effects of transition are introduced by coupling this model to the local value of the intermittency and obtaining the latter from the solution of a modeled transport equation. Comparisons with measurements obtained in oscillatory transitional flows in circular tubes show that the model produces substantial improvements over existing closures. Further pulsatile-flow predictions, driven by a mean-flow wave form obtained in a diseased human carotid artery, indicate that the intermittency-modified model yields much reduced levels of wall shear stress compared to the original, unmodified model. This result, which is attributed to the rapid growth in the thickness of the viscous sublayer arising from the severe acceleration of systole, argues in favor of the use of the model for the

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.

Hopf algebras of rooted forests, cocyles, and free Rota-Baxter algebras

NASA Astrophysics Data System (ADS)

Zhang, Tianjie; Gao, Xing; Guo, Li

2016-10-01

The Hopf algebra and the Rota-Baxter algebra are the two algebraic structures underlying the algebraic approach of Connes and Kreimer to renormalization of perturbative quantum field theory. In particular, the Hopf algebra of rooted trees serves as the "baby model" of Feynman graphs in their approach and can be characterized by certain universal properties involving a Hochschild 1-cocycle. Decorated rooted trees have also been applied to study Feynman graphs. We will continue the study of universal properties of various spaces of decorated rooted trees with such a 1-cocycle, leading to the concept of a cocycle Hopf algebra. We further apply the universal properties to equip a free Rota-Baxter algebra with the structure of a cocycle Hopf algebra.

Review and assessment of turbulence models for hypersonic flows

NASA Astrophysics Data System (ADS)

Roy, Christopher J.; Blottner, Frederick G.

2006-10-01

Accurate aerodynamic prediction is critical for the design and optimization of hypersonic vehicles. Turbulence modeling remains a major source of uncertainty in the computational prediction of aerodynamic forces and heating for these systems. The first goal of this article is to update the previous comprehensive review of hypersonic shock/turbulent boundary-layer interaction experiments published in 1991 by Settles and Dodson (Hypersonic shock/boundary-layer interaction database. NASA CR 177577, 1991). In their review, Settles and Dodson developed a methodology for assessing experiments appropriate for turbulence model validation and critically surveyed the existing hypersonic experiments. We limit the scope of our current effort by considering only two-dimensional (2D)/axisymmetric flows in the hypersonic flow regime where calorically perfect gas models are appropriate. We extend the prior database of recommended hypersonic experiments (on four 2D and two 3D shock-interaction geometries) by adding three new geometries. The first two geometries, the flat plate/cylinder and the sharp cone, are canonical, zero-pressure gradient flows which are amenable to theory-based correlations, and these correlations are discussed in detail. The third geometry added is the 2D shock impinging on a turbulent flat plate boundary layer. The current 2D hypersonic database for shock-interaction flows thus consists of nine experiments on five different geometries. The second goal of this study is to review and assess the validation usage of various turbulence models on the existing experimental database. Here we limit the scope to one- and two-equation turbulence models where integration to the wall is used (i.e., we omit studies involving wall functions). A methodology for validating turbulence models is given, followed by an extensive evaluation of the turbulence models on the current hypersonic experimental database. A total of 18 one- and two-equation turbulence models are reviewed

Directed Abelian algebras and their application to stochastic models.

PubMed

Alcaraz, F C; Rittenberg, V

2008-10-01

With each directed acyclic graph (this includes some D-dimensional lattices) one can associate some Abelian algebras that we call directed Abelian algebras (DAAs). On each site of the graph one attaches a generator of the algebra. These algebras depend on several parameters and are semisimple. Using any DAA, one can define a family of Hamiltonians which give the continuous time evolution of a stochastic process. The calculation of the spectra and ground-state wave functions (stationary state probability distributions) is an easy algebraic exercise. If one considers D-dimensional lattices and chooses Hamiltonians linear in the generators, in finite-size scaling the Hamiltonian spectrum is gapless with a critical dynamic exponent z=D. One possible application of the DAA is to sandpile models. In the paper we present this application, considering one- and two-dimensional lattices. In the one-dimensional case, when the DAA conserves the number of particles, the avalanches belong to the random walker universality class (critical exponent sigma_(tau)=32 ). We study the local density of particles inside large avalanches, showing a depletion of particles at the source of the avalanche and an enrichment at its end. In two dimensions we did extensive Monte-Carlo simulations and found sigma_(tau)=1.780+/-0.005 .

On the Subgrid-Scale Modeling of Compressible Turbulence

NASA Technical Reports Server (NTRS)

Squires, Kyle; Zeman, Otto

1990-01-01

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented.

Horizontal atmospheric turbulence, beam propagation, and modeling

NASA Astrophysics Data System (ADS)

Wilcox, Christopher C.; Santiago, Freddie; Martinez, Ty; Judd, K. Peter; Restaino, Sergio R.

2017-05-01

The turbulent effect from the Earth's atmosphere degrades the performance of an optical imaging system. Many studies have been conducted in the study of beam propagation in a turbulent medium. Horizontal beam propagation and correction presents many challenges when compared to vertical due to the far harsher turbulent conditions and increased complexity it induces. We investigate the collection of beam propagation data, analysis, and use for building a mathematical model of the horizontal turbulent path and the plans for an adaptive optical system to use this information to correct for horizontal path atmospheric turbulence.

Predicting Turbulent Convective Heat Transfer in Three-Dimensional Duct Flows

NASA Technical Reports Server (NTRS)

Rokni, M.; Gatski, T. B.

1999-01-01

The performance of an explicit algebraic stress model is assessed in predicting the turbulent flow and forced heat transfer in straight ducts, with square, rectangular, trapezoidal and triangular cross-sections, under fully developed conditions over a range of Reynolds numbers. Iso-thermal conditions are imposed on the duct walls and the turbulent heat fluxes are modeled by gradient-diffusion type models. At high Reynolds numbers (>/= 10(exp 5)), wall functions are used for the velocity and temperature fields; while at low Reynolds numbers damping functions are introduced into the models. Hydraulic parameters such as friction factor and Nusselt number are well predicted even when damping functions are used, and the present formulation imposes minimal demand on the number of grid points without any convergence or stability problems. Comparison between the models is presented in terms of the hydraulic parameters, friction factor and Nusselt number, as well as in terms of the secondary flow patterns occurring within the ducts.

A note on probabilistic models over strings: the linear algebra approach.

PubMed

Bouchard-Côté, Alexandre

2013-12-01

Probabilistic models over strings have played a key role in developing methods that take into consideration indels as phylogenetically informative events. There is an extensive literature on using automata and transducers on phylogenies to do inference on these probabilistic models, in which an important theoretical question is the complexity of computing the normalization of a class of string-valued graphical models. This question has been investigated using tools from combinatorics, dynamic programming, and graph theory, and has practical applications in Bayesian phylogenetics. In this work, we revisit this theoretical question from a different point of view, based on linear algebra. The main contribution is a set of results based on this linear algebra view that facilitate the analysis and design of inference algorithms on string-valued graphical models. As an illustration, we use this method to give a new elementary proof of a known result on the complexity of inference on the "TKF91" model, a well-known probabilistic model over strings. Compared to previous work, our proving method is easier to extend to other models, since it relies on a novel weak condition, triangular transducers, which is easy to establish in practice. The linear algebra view provides a concise way of describing transducer algorithms and their compositions, opens the possibility of transferring fast linear algebra libraries (for example, based on GPUs), as well as low rank matrix approximation methods, to string-valued inference problems.

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.

Virasoro algebra in the KN algebra; Bosonic string with fermionic ghosts on Riemann surfaces

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

Koibuchi, H.

1991-10-10

In this paper the bosonic string model with fermionic ghosts is considered in the framework of the KN algebra. The authors' attentions are paid to representations of KN algebra and a Clifford algebra of the ghosts. The authors show that a Virasoro-like algebra is obtained from KN algebra when KN algebra has certain antilinear anti-involution, and that it is isomorphic to the usual Virasoro algebra. The authors show that there is an expected relation between a central charge of this Virasoro-like algebra and an anomaly of the combined system.

A set for relational reasoning: Facilitation of algebraic modeling by a fraction task.

PubMed

DeWolf, Melissa; Bassok, Miriam; Holyoak, Keith J

2016-12-01

Recent work has identified correlations between early mastery of fractions and later math achievement, especially in algebra. However, causal connections between aspects of reasoning with fractions and improved algebra performance have yet to be established. The current study investigated whether relational reasoning with fractions facilitates subsequent algebraic reasoning using both pre-algebra students and adult college students. Participants were first given either a relational reasoning fractions task or a fraction algebra procedures control task. Then, all participants solved word problems and constructed algebraic equations in either multiplication or division format. The word problems and the equation construction tasks involved simple multiplicative comparison statements such as "There are 4 times as many students as teachers in a classroom." Performance on the algebraic equation construction task was enhanced for participants who had previously completed the relational fractions task compared with those who completed the fraction algebra procedures task. This finding suggests that relational reasoning with fractions can establish a relational set that promotes students' tendency to model relations using algebraic expressions. Copyright © 2016 Elsevier Inc. All rights reserved.

Flame-conditioned turbulence modeling for reacting flows

NASA Astrophysics Data System (ADS)

Macart, Jonathan F.; Mueller, Michael E.

2017-11-01

Conventional approaches to turbulence modeling in reacting flows rely on unconditional averaging or filtering, that is, consideration of the momentum equations only in physical space, implicitly assuming that the flame only weakly affects the turbulence, aside from a variation in density. Conversely, for scalars, which are strongly coupled to the flame structure, their evolution equations are often projected onto a reduced-order manifold, that is, conditionally averaged or filtered, on a flame variable such as a mixture fraction or progress variable. Such approaches include Conditional Moment Closure (CMC) and related variants. However, recent observations from Direct Numerical Simulation (DNS) have indicated that the flame can strongly affect turbulence in premixed combustion at low Karlovitz number. In this work, a new approach to turbulence modeling for reacting flows is investigated in which conditionally averaged or filtered equations are evolved for the momentum. The conditionally-averaged equations for the velocity and its covariances are derived, and budgets are evaluated from DNS databases of turbulent premixed planar jet flames. The most important terms in these equations are identified, and preliminary closure models are proposed.

Apparent Transition Behavior of Widely-Used Turbulence Models

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.

2006-01-01

The Spalart-Allmaras and the Menter SST kappa-omega turbulence models are shown to have the undesirable characteristic that, for fully turbulent computations, a transition region can occur whose extent varies with grid density. Extremely fine two-dimensional grids over the front portion of an airfoil are used to demonstrate the effect. As the grid density is increased, the laminar region near the nose becomes larger. In the Spalart-Allmaras model this behavior is due to convergence to a laminar-behavior fixed point that occurs in practice when freestream turbulence is below some threshold. It is the result of a feature purposefully added to the original model in conjunction with a special trip function. This degenerate fixed point can also cause nonuniqueness regarding where transition initiates on a given grid. Consistent fully turbulent results can easily be achieved by either using a higher freestream turbulence level or by making a simple change to one of the model constants. Two-equation kappa-omega models, including the SST model, exhibit strong sensitivity to numerical resolution near the area where turbulence initiates. Thus, inconsistent apparent transition behavior with grid refinement in this case does not appear to stem from the presence of a degenerate fixed point. Rather, it is a fundamental property of the kappa-omega model itself, and is not easily remedied.

Apparent Transition Behavior of Widely-Used Turbulence Models

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.

2007-01-01

The Spalart-Allmaras and the Menter SST k-omega turbulence models are shown to have the undesirable characteristic that, for fully turbulent computations, a transition region can occur whose extent varies with grid density. Extremely fine two-dimensional grids over the front portion of an airfoil are used to demonstrate the effect. As the grid density is increased, the laminar region near the nose becomes larger. In the Spalart-Allmaras model this behavior is due to convergence to a laminar-behavior fixed point that occurs in practice when freestream turbulence is below some threshold. It is the result of a feature purposefully added to the original model in conjunction with a special trip function. This degenerate fixed point can also cause non-uniqueness regarding where transition initiates on a given grid. Consistent fully turbulent results can easily be achieved by either using a higher freestream turbulence level or by making a simple change to one of the model constants. Two-equation k-omega models, including the SST model, exhibit strong sensitivity to numerical resolution near the area where turbulence initiates. Thus, inconsistent apparent transition behavior with grid refinement in this case does not appear to stem from the presence of a degenerate fixed point. Rather, it is a fundamental property of the k-omega model itself, and is not easily remedied.

Finite-element numerical modeling of atmospheric turbulent boundary layer

NASA Technical Reports Server (NTRS)

Lee, H. N.; Kao, S. K.

1979-01-01

A dynamic turbulent boundary-layer model in the neutral atmosphere is constructed, using a dynamic turbulent equation of the eddy viscosity coefficient for momentum derived from the relationship among the turbulent dissipation rate, the turbulent kinetic energy and the eddy viscosity coefficient, with aid of the turbulent second-order closure scheme. A finite-element technique was used for the numerical integration. In preliminary results, the behavior of the neutral planetary boundary layer agrees well with the available data and with the existing elaborate turbulent models, using a finite-difference scheme. The proposed dynamic formulation of the eddy viscosity coefficient for momentum is particularly attractive and can provide a viable alternative approach to study atmospheric turbulence, diffusion and air pollution.

Modeling of Turbulence Effects on Liquid Jet Atomization and Breakup

NASA Technical Reports Server (NTRS)

Trinh, Huu; Chen, C. P.

2004-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. For certain flow regimes, it has been observed that the liquid jet surface is highly turbulent. This turbulence characteristic plays a key role on the breakup of the liquid jet near to the injector exit. Other experiments also showed that the breakup length of the liquid core is sharply shortened as the liquid jet is changed from the laminar to the turbulent flow conditions. In the numerical and physical modeling arena, most of commonly used atomization models do not include the turbulence effect. Limited attempts have been made in modeling the turbulence phenomena on the liquid jet disintegration. The subject correlation and models treat the turbulence either as an only source or a primary driver in the breakup process. This study aims to model the turbulence effect in the atomization process of a cylindrical liquid jet. In the course of this study, two widely used models, Reitz's primary atomization (blob) and Taylor-Analogy-Break (TAB) secondary droplet breakup by O Rourke et al. are examined. Additional terms are derived and implemented appropriately into these two models to account for the turbulence effect on the atomization process. Since this enhancement effort is based on a framework of the two existing atomization models, it is appropriate to denote the two present models as T-blob and T-TAB for the primary and secondary atomization predictions, respectively. In the primary breakup model, the level of the turbulence effect on the liquid breakup depends on the characteristic time scales and the initial flow conditions. This treatment offers a balance of contributions of individual physical phenomena on the liquid breakup process. For the secondary breakup, an addition turbulence force acted on parent drops is modeled and integrated into the TAB governing equation. The drop size

Streamwise Vorticity Generation in Laminar and Turbulent Jets

NASA Technical Reports Server (NTRS)

Demuren, Aodeji O.; Wilson, Robert V.

1999-01-01

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

Implementation of algebraic stress models in a general 3-D Navier-Stokes method (PAB3D)

NASA Technical Reports Server (NTRS)

Abdol-Hamid, Khaled S.

1995-01-01

A three-dimensional multiblock Navier-Stokes code, PAB3D, which was developed for propulsion integration and general aerodynamic analysis, has been used extensively by NASA Langley and other organizations to perform both internal (exhaust) and external flow analysis of complex aircraft configurations. This code was designed to solve the simplified Reynolds Averaged Navier-Stokes equations. A two-equation k-epsilon turbulence model has been used with considerable success, especially for attached flows. Accurate predicting of transonic shock wave location and pressure recovery in separated flow regions has been more difficult. Two algebraic Reynolds stress models (ASM) have been recently implemented in the code that greatly improved the code's ability to predict these difficult flow conditions. Good agreement with Direct Numerical Simulation (DNS) for a subsonic flat plate was achieved with ASM's developed by Shih, Zhu, and Lumley and Gatski and Speziale. Good predictions were also achieved at subsonic and transonic Mach numbers for shock location and trailing edge boattail pressure recovery on a single-engine afterbody/nozzle model.

Inverse Modelling Problems in Linear Algebra Undergraduate Courses

ERIC Educational Resources Information Center

Martinez-Luaces, Victor E.

2013-01-01

This paper will offer an analysis from a theoretical point of view of mathematical modelling, applications and inverse problems of both causation and specification types. Inverse modelling problems give the opportunity to establish connections between theory and practice and to show this fact, a simple linear algebra example in two different…

Modeling of Turbulent Free Shear Flows

NASA Technical Reports Server (NTRS)

Yoder, Dennis A.; DeBonis, James R.; Georgiadis, Nicolas J.

2013-01-01

The modeling of turbulent free shear flows is crucial to the simulation of many aerospace applications, yet often receives less attention than the modeling of wall boundary layers. Thus, while turbulence model development in general has proceeded very slowly in the past twenty years, progress for free shear flows has been even more so. This paper highlights some of the fundamental issues in modeling free shear flows for propulsion applications, presents a review of past modeling efforts, and identifies areas where further research is needed. Among the topics discussed are differences between planar and axisymmetric flows, development versus self-similar regions, the effect of compressibility and the evolution of compressibility corrections, the effect of temperature on jets, and the significance of turbulent Prandtl and Schmidt numbers for reacting shear flows. Large eddy simulation greatly reduces the amount of empiricism in the physical modeling, but is sensitive to a number of numerical issues. This paper includes an overview of the importance of numerical scheme, mesh resolution, boundary treatment, sub-grid modeling, and filtering in conducting a successful simulation.

Applied Algebra: The Modeling Technique of Least Squares

ERIC Educational Resources Information Center

Zelkowski, Jeremy; Mayes, Robert

2008-01-01

The article focuses on engaging students in algebra through modeling real-world problems. The technique of least squares is explored, encouraging students to develop a deeper understanding of the method. (Contains 2 figures and a bibliography.)

Dynamic and Thermal Turbulent Time Scale Modelling for Homogeneous Shear Flows

NASA Technical Reports Server (NTRS)

Schwab, John R.; Lakshminarayana, Budugur

1994-01-01

A new turbulence model, based upon dynamic and thermal turbulent time scale transport equations, is developed and applied to homogeneous shear flows with constant velocity and temperature gradients. The new model comprises transport equations for k, the turbulent kinetic energy; tau, the dynamic time scale; k(sub theta), the fluctuating temperature variance; and tau(sub theta), the thermal time scale. It offers conceptually parallel modeling of the dynamic and thermal turbulence at the two equation level, and eliminates the customary prescription of an empirical turbulent Prandtl number, Pr(sub t), thus permitting a more generalized prediction capability for turbulent heat transfer in complex flows and geometries. The new model also incorporates constitutive relations, based upon invariant theory, that allow the effects of nonequilibrium to modify the primary coefficients for the turbulent shear stress and heat flux. Predictions of the new model, along with those from two other similar models, are compared with experimental data for decaying homogeneous dynamic and thermal turbulence, homogeneous turbulence with constant temperature gradient, and homogeneous turbulence with constant temperature gradient and constant velocity gradient. The new model offers improvement in agreement with the data for most cases considered in this work, although it was no better than the other models for several cases where all the models performed poorly.

A minimal model of self-sustaining turbulence

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

Thomas, Vaughan L.; Gayme, Dennice F.; Farrell, Brian F.

2015-10-15

In this work, we examine the turbulence maintained in a Restricted Nonlinear (RNL) model of plane Couette flow. This model is a computationally efficient approximation of the second order statistical state dynamics obtained by partitioning the flow into a streamwise averaged mean flow and perturbations about that mean, a closure referred to herein as the RNL{sub ∞} model. The RNL model investigated here employs a single member of the infinite ensemble that comprises the covariance of the RNL{sub ∞} dynamics. The RNL system has previously been shown to support self-sustaining turbulence with a mean flow and structural features that aremore » consistent with direct numerical simulations (DNS). Regardless of the number of streamwise Fourier components used in the simulation, the RNL system’s self-sustaining turbulent state is supported by a small number of streamwise varying modes. Remarkably, further truncation of the RNL system’s support to as few as one streamwise varying mode can suffice to sustain the turbulent state. The close correspondence between RNL simulations and DNS that has been previously demonstrated along with the results presented here suggest that the fundamental mechanisms underlying wall-turbulence can be analyzed using these highly simplified RNL systems.« less

Mathematical modelling in engineering: an alternative way to teach Linear Algebra

NASA Astrophysics Data System (ADS)

Domínguez-García, S.; García-Planas, M. I.; Taberna, J.

2016-10-01

Technological advances require that basic science courses for engineering, including Linear Algebra, emphasize the development of mathematical strengths associated with modelling and interpretation of results, which are not limited only to calculus abilities. Based on this consideration, we have proposed a project-based learning, giving a dynamic classroom approach in which students modelled real-world problems and turn gain a deeper knowledge of the Linear Algebra subject. Considering that most students are digital natives, we use the e-portfolio as a tool of communication between students and teachers, besides being a good place making the work visible. In this article, we present an overview of the design and implementation of a project-based learning for a Linear Algebra course taught during the 2014-2015 at the 'ETSEIB'of Universitat Politècnica de Catalunya (UPC).

Turbulence Modeling: Progress and Future Outlook

NASA Technical Reports Server (NTRS)

Marvin, Joseph G.; Huang, George P.

1996-01-01

Progress in the development of the hierarchy of turbulence models for Reynolds-averaged Navier-Stokes codes used in aerodynamic applications is reviewed. Steady progress is demonstrated, but transfer of the modeling technology has not kept pace with the development and demands of the computational fluid dynamics (CFD) tools. An examination of the process of model development leads to recommendations for a mid-course correction involving close coordination between modelers, CFD developers, and application engineers. In instances where the old process is changed and cooperation enhanced, timely transfer is realized. A turbulence modeling information database is proposed to refine the process and open it to greater participation among modeling and CFD practitioners.

Assessment of Turbulent CFD Against STS-128 Hypersonic Flight Data

NASA Technical Reports Server (NTRS)

Wood, William A.; Kleb, William L.; Hyatt, Andrew J.

2010-01-01

Turbulent CFD simulations are compared against surface temperature measurements of the space shuttle orbiter windward tiles at reentry flight conditions. Algebraic turbulence models are used within both the LAURA and DPLR CFD codes. The flight data are from temperature measurements obtained by seven thermocouples during the STS-128 mission (September 2009). The flight data indicate boundary layer transition onset over the Mach number range 13.5{15.5, depending upon the location on the vehicle. But the boundary layer flow appeared to be transitional down through Mach 12, based upon the flight data and CFD trends. At Mach 9 the simulations match the flight data on average within 20 F/11 C, where typical surface temperatures were approximately 1600 F/870 C.

Algebraic model checking for Boolean gene regulatory networks.

PubMed

Tran, Quoc-Nam

2011-01-01

We present a computational method in which modular and Groebner bases (GB) computation in Boolean rings are used for solving problems in Boolean gene regulatory networks (BN). In contrast to other known algebraic approaches, the degree of intermediate polynomials during the calculation of Groebner bases using our method will never grow resulting in a significant improvement in running time and memory space consumption. We also show how calculation in temporal logic for model checking can be done by means of our direct and efficient Groebner basis computation in Boolean rings. We present our experimental results in finding attractors and control strategies of Boolean networks to illustrate our theoretical arguments. The results are promising. Our algebraic approach is more efficient than the state-of-the-art model checker NuSMV on BNs. More importantly, our approach finds all solutions for the BN problems.

Gasdynamic model of turbulent combustion in an explosion

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

Kuhl, A.L.; Ferguson, R.E.; Chien, K.Y.

1994-08-31

Proposed here is a gasdynamic model of turbulent combustion in explosions. It is used to investigate turbulent mixing aspects of afterburning found in TNT charges detonated in air. Evolution of the turbulent velocity field was calculated by a high-order Godunov solution of the gasdynamic equations. Adaptive Mesh Refinement (AMR) was used to follow convective-mixing processes on the computational grid. Combustion was then taken into account by a simplified sub-grid model, demonstrating that it was controlled by turbulent mixing. The rate of fuel consumption decayed inversely with time, and was shown to be insensitive to grid resolution.

Computations of turbulent lean premixed combustion using conditional moment closure

NASA Astrophysics Data System (ADS)

Amzin, Shokri; Swaminathan, Nedunchezhian

2013-12-01

Conditional Moment Closure (CMC) is a suitable method for predicting scalars such as carbon monoxide with slow chemical time scales in turbulent combustion. Although this method has been successfully applied to non-premixed combustion, its application to lean premixed combustion is rare. In this study the CMC method is used to compute piloted lean premixed combustion in a distributed combustion regime. The conditional scalar dissipation rate of the conditioning scalar, the progress variable, is closed using an algebraic model and turbulence is modelled using the standard k-ɛ model. The conditional mean reaction rate is closed using a first order CMC closure with the GRI-3.0 chemical mechanism to represent the chemical kinetics of methane oxidation. The PDF of the progress variable is obtained using a presumed shape with the Beta function. The computed results are compared with the experimental measurements and earlier computations using the transported PDF approach. The results show reasonable agreement with the experimental measurements and are consistent with the transported PDF computations. When the compounded effects of shear-turbulence and flame are strong, second order closures may be required for the CMC.

Conference on Complex Turbulent Flows: Comparison of Computation and Experiment, Stanford University, Stanford, CA, September 14-18, 1981, Proceedings. Volume 2 - Taxonomies, reporters' summaries, evaluation, and conclusions

NASA Technical Reports Server (NTRS)

Kline, S. J. (Editor); Cantwell, B. J. (Editor); Lilley, G. M.

1982-01-01

Computational techniques for simulating turbulent flows were explored, together with the results of experimental investigations. Particular attention was devoted to the possibility of defining a universal closure model, applicable for all turbulence situations; however, conclusions were drawn that zonal models, describing localized structures, were the most promising techniques to date. The taxonomy of turbulent flows was summarized, as were algebraic, differential, integral, and partial differential methods for numerical depiction of turbulent flows. Numerous comparisons of theoretically predicted and experimentally obtained data for wall pressure distributions, velocity profiles, turbulent kinetic energy profiles, Reynolds shear stress profiles, and flows around transonic airfoils were presented. Simplifying techniques for reducing the necessary computational time for modeling complex flowfields were surveyed, together with the industrial requirements and applications of computational fluid dynamics techniques.

A non-isotropic multiple-scale turbulence model

NASA Technical Reports Server (NTRS)

Chen, C. P.

1990-01-01

A newly developed non-isotropic multiple scale turbulence model (MS/ASM) is described for complex flow calculations. This model focuses on the direct modeling of Reynolds stresses and utilizes split-spectrum concepts for modeling multiple scale effects in turbulence. Validation studies on free shear flows, rotating flows and recirculating flows show that the current model perform significantly better than the single scale k-epsilon model. The present model is relatively inexpensive in terms of CPU time which makes it suitable for broad engineering flow applications.

Modeling turbulent energy behavior and sudden viscous dissipation in compressing plasma turbulence

DOE PAGES

Davidovits, Seth; Fisch, Nathaniel J.

2017-12-21

Here, we present a simple model for the turbulent kinetic energy behavior of subsonic plasma turbulence undergoing isotropic three-dimensional compression, which may exist in various inertial confinement fusion experiments or astrophysical settings. The plasma viscosity depends on both the temperature and the ionization state, for which many possible scalings with compression are possible. For example, in an adiabatic compression the temperature scales as 1/L 2, with L the linear compression ratio, but if thermal energy loss mechanisms are accounted for, the temperature scaling may be weaker. As such, the viscosity has a wide range of net dependencies on the compression.more » The model presented here, with no parameter changes, agrees well with numerical simulations for a range of these dependencies. This model permits the prediction of the partition of injected energy between thermal and turbulent energy in a compressing plasma.« less

Modeling turbulent energy behavior and sudden viscous dissipation in compressing plasma turbulence

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

Davidovits, Seth; Fisch, Nathaniel J.

Here, we present a simple model for the turbulent kinetic energy behavior of subsonic plasma turbulence undergoing isotropic three-dimensional compression, which may exist in various inertial confinement fusion experiments or astrophysical settings. The plasma viscosity depends on both the temperature and the ionization state, for which many possible scalings with compression are possible. For example, in an adiabatic compression the temperature scales as 1/L 2, with L the linear compression ratio, but if thermal energy loss mechanisms are accounted for, the temperature scaling may be weaker. As such, the viscosity has a wide range of net dependencies on the compression.more » The model presented here, with no parameter changes, agrees well with numerical simulations for a range of these dependencies. This model permits the prediction of the partition of injected energy between thermal and turbulent energy in a compressing plasma.« less

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.

Scaling and modeling of turbulent suspension flows

NASA Technical Reports Server (NTRS)

Chen, C. P.

1989-01-01

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.

Turbulence modeling and experiments

NASA Technical Reports Server (NTRS)

Shabbir, Aamir

1992-01-01

The best way of verifying turbulence is to do a direct comparison between the various terms and their models. The success of this approach depends upon the availability of the data for the exact correlations (both experimental and DNS). The other approach involves numerically solving the differential equations and then comparing the results with the data. The results of such a computation will depend upon the accuracy of all the modeled terms and constants. Because of this it is sometimes difficult to find the cause of a poor performance by a model. However, such a calculation is still meaningful in other ways as it shows how a complete Reynolds stress model performs. Thirteen homogeneous flows are numerically computed using the second order closure models. We concentrate only on those models which use a linear (or quasi-linear) model for the rapid term. This, therefore, includes the Launder, Reece and Rodi (LRR) model; the isotropization of production (IP) model; and the Speziale, Sarkar, and Gatski (SSG) model. Which of the three models performs better is examined along with what are their weaknesses, if any. The other work reported deal with the experimental balances of the second moment equations for a buoyant plume. Despite the tremendous amount of activity toward the second order closure modeling of turbulence, very little experimental information is available about the budgets of the second moment equations. Part of the problem stems from our inability to measure the pressure correlations. However, if everything else appearing in these equations is known from the experiment, pressure correlations can be obtained as the closing terms. This is the closest we can come to in obtaining these terms from experiment, and despite the measurement errors which might be present in such balances, the resulting information will be extremely useful for the turbulence modelers. The purpose of this part of the work was to provide such balances of the Reynolds stress and heat

Global identifiability of linear compartmental models--a computer algebra algorithm.

PubMed

Audoly, S; D'Angiò, L; Saccomani, M P; Cobelli, C

1998-01-01

A priori global identifiability deals with the uniqueness of the solution for the unknown parameters of a model and is, thus, a prerequisite for parameter estimation of biological dynamic models. Global identifiability is however difficult to test, since it requires solving a system of algebraic nonlinear equations which increases both in nonlinearity degree and number of terms and unknowns with increasing model order. In this paper, a computer algebra tool, GLOBI (GLOBal Identifiability) is presented, which combines the topological transfer function method with the Buchberger algorithm, to test global identifiability of linear compartmental models. GLOBI allows for the automatic testing of a priori global identifiability of general structure compartmental models from general multi input-multi output experiments. Examples of usage of GLOBI to analyze a priori global identifiability of some complex biological compartmental models are provided.

Industry-Wide Workshop on Computational Turbulence Modeling

NASA Technical Reports Server (NTRS)

Shabbir, Aamir (Compiler)

1995-01-01

This publication contains the presentations made at the Industry-Wide Workshop on Computational Turbulence Modeling which took place on October 6-7, 1994. The purpose of the workshop was to initiate the transfer of technology developed at Lewis Research Center to industry and to discuss the current status and the future needs of turbulence models in industrial CFD.

Geometric model of topological insulators from the Maxwell algebra

NASA Astrophysics Data System (ADS)

Palumbo, Giandomenico

2017-11-01

We propose a novel geometric model of time-reversal-invariant topological insulators in three dimensions in presence of an external electromagnetic field. Their gapped boundary supports relativistic quantum Hall states and is described by a Chern-Simons theory, where the gauge connection takes values in the Maxwell algebra. This represents a non-central extension of the Poincaré algebra and takes into account both the Lorentz and magnetic-translation symmetries of the surface states. In this way, we derive a relativistic version of the Wen-Zee term and we show that the non-minimal coupling between the background geometry and the electromagnetic field in the model is in agreement with the main properties of the relativistic quantum Hall states in the flat space.

Piloted Evaluation of a UH-60 Mixer Equivalent Turbulence Simulation Model

NASA Technical Reports Server (NTRS)

Lusardi, Jeff A.; Blanken, Chris L.; Tischeler, Mark B.

2002-01-01

A simulation study of a recently developed hover/low speed Mixer Equivalent Turbulence Simulation (METS) model for the UH-60 Black Hawk helicopter was conducted in the NASA Ames Research Center Vertical Motion Simulator (VMS). The experiment was a continuation of previous work to develop a simple, but validated, turbulence model for hovering rotorcraft. To validate the METS model, two experienced test pilots replicated precision hover tasks that had been conducted in an instrumented UH-60 helicopter in turbulence. Objective simulation data were collected for comparison with flight test data, and subjective data were collected that included handling qualities ratings and pilot comments for increasing levels of turbulence. Analyses of the simulation results show good analytic agreement between the METS model and flight test data, with favorable pilot perception of the simulated turbulence. Precision hover tasks were also repeated using the more complex rotating-frame SORBET (Simulation Of Rotor Blade Element Turbulence) model to generate turbulence. Comparisons of the empirically derived METS model with the theoretical SORBET model show good agreement providing validation of the more complex blade element method of simulating turbulence.

Analysis of two-equation turbulence models for recirculating flows

NASA Technical Reports Server (NTRS)

Thangam, S.

1991-01-01

The two-equation kappa-epsilon model is used to analyze turbulent separated flow past a backward-facing step. It is shown that if the model constraints are modified to be consistent with the accepted energy decay rate for isotropic turbulence, the dominant features of the flow field, namely the size of the separation bubble and the streamwise component of the mean velocity, can be accurately predicted. In addition, except in the vicinity of the step, very good predictions for the turbulent shear stress, the wall pressure, and the wall shear stress are obtained. The model is also shown to provide good predictions for the turbulence intensity in the region downstream of the reattachment point. Estimated long time growth rates for the turbulent kinetic energy and dissipation rate of homogeneous shear flow are utilized to develop an optimal set of constants for the two equation kappa-epsilon model. The physical implications of the model performance are also discussed.

Collaborative testing of turbulence models

NASA Astrophysics Data System (ADS)

Bradshaw, P.

1992-12-01

This project, funded by AFOSR, ARO, NASA, and ONR, was run by the writer with Profs. Brian E. Launder, University of Manchester, England, and John L. Lumley, Cornell University. Statistical data on turbulent flows, from lab. experiments and simulations, were circulated to modelers throughout the world. This is the first large-scale project of its kind to use simulation data. The modelers returned their predictions to Stanford, for distribution to all modelers and to additional participants ('experimenters')--over 100 in all. The object was to obtain a consensus on the capabilities of present-day turbulence models and identify which types most deserve future support. This was not completely achieved, mainly because not enough modelers could produce results for enough test cases within the duration of the project. However, a clear picture of the capabilities of various modeling groups has appeared, and the interaction has been helpful to the modelers. The results support the view that Reynolds-stress transport models are the most accurate.

On the quasi-conical flowfield structure of the swept shock wave-turbulent boundary layer interaction

NASA Technical Reports Server (NTRS)

Knight, Doyle D.; Badekas, Dias

1991-01-01

The swept oblique shock-wave/turbulent-boundary-layer interaction generated by a 20-deg sharp fin at Mach 4 and Reynolds number 21,000 is investigated via a series of computations using both conical and three-dimensional Reynolds-averaged Navier-Stokes equations with turbulence incorporated through the algebraic turbulent eddy viscosity model of Baldwin-Lomax. Results are compared with known experimental data, and it is concluded that the computed three-dimensional flowfield is quasi-conical (in agreement with the experimental data), the computed three-dimensional and conical surface pressure and surface flow direction are in good agreement with the experiment, and the three-dimensional and conical flows significantly underpredict the peak experimental skin friction. It is pointed out that most of the features of the conical flowfield model in the experiment are observed in the conical computation which also describes the complete conical streamline pattern not included in the model of the experiment.

ADAM: analysis of discrete models of biological systems using computer algebra.

PubMed

Hinkelmann, Franziska; Brandon, Madison; Guang, Bonny; McNeill, Rustin; Blekherman, Grigoriy; Veliz-Cuba, Alan; Laubenbacher, Reinhard

2011-07-20

Many biological systems are modeled qualitatively with discrete models, such as probabilistic Boolean networks, logical models, Petri nets, and agent-based models, to gain a better understanding of them. The computational complexity to analyze the complete dynamics of these models grows exponentially in the number of variables, which impedes working with complex models. There exist software tools to analyze discrete models, but they either lack the algorithmic functionality to analyze complex models deterministically or they are inaccessible to many users as they require understanding the underlying algorithm and implementation, do not have a graphical user interface, or are hard to install. Efficient analysis methods that are accessible to modelers and easy to use are needed. We propose a method for efficiently identifying attractors and introduce the web-based tool Analysis of Dynamic Algebraic Models (ADAM), which provides this and other analysis methods for discrete models. ADAM converts several discrete model types automatically into polynomial dynamical systems and analyzes their dynamics using tools from computer algebra. Specifically, we propose a method to identify attractors of a discrete model that is equivalent to solving a system of polynomial equations, a long-studied problem in computer algebra. Based on extensive experimentation with both discrete models arising in systems biology and randomly generated networks, we found that the algebraic algorithms presented in this manuscript are fast for systems with the structure maintained by most biological systems, namely sparseness and robustness. For a large set of published complex discrete models, ADAM identified the attractors in less than one second. Discrete modeling techniques are a useful tool for analyzing complex biological systems and there is a need in the biological community for accessible efficient analysis tools. ADAM provides analysis methods based on mathematical algorithms as a web

ADAM: Analysis of Discrete Models of Biological Systems Using Computer Algebra

PubMed Central

2011-01-01

Background Many biological systems are modeled qualitatively with discrete models, such as probabilistic Boolean networks, logical models, Petri nets, and agent-based models, to gain a better understanding of them. The computational complexity to analyze the complete dynamics of these models grows exponentially in the number of variables, which impedes working with complex models. There exist software tools to analyze discrete models, but they either lack the algorithmic functionality to analyze complex models deterministically or they are inaccessible to many users as they require understanding the underlying algorithm and implementation, do not have a graphical user interface, or are hard to install. Efficient analysis methods that are accessible to modelers and easy to use are needed. Results We propose a method for efficiently identifying attractors and introduce the web-based tool Analysis of Dynamic Algebraic Models (ADAM), which provides this and other analysis methods for discrete models. ADAM converts several discrete model types automatically into polynomial dynamical systems and analyzes their dynamics using tools from computer algebra. Specifically, we propose a method to identify attractors of a discrete model that is equivalent to solving a system of polynomial equations, a long-studied problem in computer algebra. Based on extensive experimentation with both discrete models arising in systems biology and randomly generated networks, we found that the algebraic algorithms presented in this manuscript are fast for systems with the structure maintained by most biological systems, namely sparseness and robustness. For a large set of published complex discrete models, ADAM identified the attractors in less than one second. Conclusions Discrete modeling techniques are a useful tool for analyzing complex biological systems and there is a need in the biological community for accessible efficient analysis tools. ADAM provides analysis methods based on

Modeling the pressure-dilatation correlation

NASA Technical Reports Server (NTRS)

Sarkar, S.

1991-01-01

It is generally accepted that pressure dilatation, which is an additional compressibility term in turbulence transport equations, may be important for high speed flows. Recent direct simulations of homogeneous shear turbulence have given concrete evidence that the pressure dilatation is important insofar that it contributes to the reduced growth of turbulent kinetic energy due to compressibility effects. The problem of modeling pressure dilatation is addressed. A component of the pressure dilatation is isolated which exhibits temporal oscillations and, using direct numerical simulations of homogeneous shear turbulence and isotropic turbulence, show that it has a negligible contribution to the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous turbulence is performed to obtain a model for the nonoscillatory pressure dilatation. This model algebraically relates the pressure dilatation to quantities traditionally obtained in incompressible turbulence closures. The model is validated by direct comparison with the pressure dilatation data obtained from the simulations.

A new turbulence-based model for sand transport

NASA Astrophysics Data System (ADS)

Mayaud, Jerome; Wiggs, Giles; Bailey, Richard

2016-04-01

Knowledge of the changing rate of sediment flux in space and time is essential for quantifying surface erosion and deposition in desert landscapes. While many aeolian studies have relied on time-averaged parameters such as wind velocity (U) and wind shear velocity (u*) to determine sediment flux, there is increasing evidence that high-frequency turbulence is an important driving force behind the entrainment and transport of sand. However, turbulence has yet to be incorporated into a functional sand transport model that can be used for predictive purposes. In this study we present a new transport model (the 'turbulence model') that accounts for high-frequency variations in the horizontal (u) and vertical (w) components of wind flow. The turbulence model is fitted to wind velocity and sediment transport data from a field experiment undertaken in Namibia's Skeleton Coast National Park, and its performance at three temporal resolutions (10 Hz, 1 Hz, 1 min) is compared to two existing models that rely on time-averaged wind velocity data (Radok, 1977; Dong et al., 2003). The validity of the three models is analysed under a variety of saltation conditions, using a 2-hour (1 Hz measurement resolution) dataset from the Skeleton Coast and a 5-hour (1 min measurement resolution) dataset from the southwestern Kalahari Desert. The turbulence model is shown to outperform the Radok and Dong models when predicting total saltation count over the three experimental periods. For all temporal resolutions presented in this study (10 Hz-10 min), the turbulence model predicted total saltation count to within at least 0.34%, whereas the Radok and Dong models over- or underestimated total count by up to 5.50% and 20.53% respectively. The strong performance of the turbulence model can be attributed to a lag in mass flux response built into its formulation, which can be adapted depending on the temporal resolution of investigation. This accounts for the inherent lag within the physical

Model-free simulations of turbulent reactive flows

NASA Technical Reports Server (NTRS)

Givi, Peyman

1989-01-01

The current computational methods for solving transport equations of turbulent reacting single-phase flows are critically reviewed, with primary attention given to those methods that lead to model-free simulations. In particular, consideration is given to direct numerical simulations using spectral (Galerkin) and pseudospectral (collocation) methods, spectral element methods, and Lagrangian methods. The discussion also covers large eddy simulations and turbulence modeling.

A geometrical optics approach for modeling atmospheric turbulence

NASA Astrophysics Data System (ADS)

Yuksel, Heba; Atia, Walid; Davis, Christopher C.

2005-08-01

Atmospheric turbulence has a significant impact on the quality of a laser beam propagating through the atmosphere over long distances. Turbulence causes the optical phasefront to become distorted from propagation through turbulent eddies of varying sizes and refractive index. Turbulence also results in intensity scintillation and beam wander, which can severely impair the operation of target designation and free space optical (FSO) communications systems. We have developed a new model to assess the effects of turbulence on laser beam propagation in such applications. We model the atmosphere along the laser beam propagation path as a spatial distribution of spherical bubbles or curved interfaces. The size and refractive index discontinuity represented by each bubble are statistically distributed according to various models. For each statistical representation of the atmosphere, the path of a single ray, or a bundle of rays, is analyzed using geometrical optics. These Monte Carlo techniques allow us to assess beam wander, beam spread, and phase shifts along the path. An effective Cn2 can be determined by correlating beam wander behavior with the path length. This model has already proved capable of assessing beam wander, in particular the (Range)3 dependence of mean-squared beam wander, and in estimating lateral phase decorrelations that develop across the laser phasefront as it propagates through turbulence. In addition, we have developed efficient computational techniques for various correlation functions that are important in assessing the effects of turbulence. The Monte Carlo simulations are compared and show good agreement with the predictions of wave theory.

A Modeling-Based College Algebra Course and Its Effect on Student Achievement

ERIC Educational Resources Information Center

Ellington, Aimee J.

2005-01-01

In Fall 2004, Virginia Commonwealth University (VCU) piloted a modeling-based approach to college algebra. This paper describes the course and an assessment that was conducted to determine the effect of this approach on student achievement in comparison to a traditional approach to college algebra. The results show that compared with their…

Full-coverage film cooling: 3-dimensional measurements of turbulence structure and prediction of recovery region hydrodynamics

NASA Technical Reports Server (NTRS)

Yavuzkurt, S.; Moffat, R. J.; Kays, W. M.

1979-01-01

Hydrodynamic measurements were made with a triaxial hot-wire in the full-coverage region and the recovery region following an array of injection holes inclined downstream, at 30 degrees to the surface. The data were taken under isothermal conditions at ambient temperature and pressure for two blowing ratios: M = 0.9 and M = 0.4. Profiles of the three main velocity components and the six Reynolds stresses were obtained at several spanwise positions at each of the five locations down the test plate. A one-equation model of turbulence (using turbulent kinetic energy with an algebraic mixing length) was used in a two-dimensional computer program to predict the mean velocity and turbulent kinetic energy profiles in the recovery region. A new real-time hotwire scheme was developed to make measurements in the three-dimensional turbulent boundary layer over the full-coverage surface.

An algebraic method for constructing stable and consistent autoregressive filters

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

Harlim, John, E-mail: jharlim@psu.edu; Department of Meteorology, the Pennsylvania State University, University Park, PA 16802; Hong, Hoon, E-mail: hong@ncsu.edu

2015-02-15

In this paper, we introduce an algebraic method to construct stable and consistent univariate autoregressive (AR) models of low order for filtering and predicting nonlinear turbulent signals with memory depth. By stable, we refer to the classical stability condition for the AR model. By consistent, we refer to the classical consistency constraints of Adams–Bashforth methods of order-two. One attractive feature of this algebraic method is that the model parameters can be obtained without directly knowing any training data set as opposed to many standard, regression-based parameterization methods. It takes only long-time average statistics as inputs. The proposed method provides amore » discretization time step interval which guarantees the existence of stable and consistent AR model and simultaneously produces the parameters for the AR models. In our numerical examples with two chaotic time series with different characteristics of decaying time scales, we find that the proposed AR models produce significantly more accurate short-term predictive skill and comparable filtering skill relative to the linear regression-based AR models. These encouraging results are robust across wide ranges of discretization times, observation times, and observation noise variances. Finally, we also find that the proposed model produces an improved short-time prediction relative to the linear regression-based AR-models in forecasting a data set that characterizes the variability of the Madden–Julian Oscillation, a dominant tropical atmospheric wave pattern.« 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.

An improved k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Shih, T. H.; Hsu, Andrew T.

1991-01-01

An improved k-epsilon model for low Reynolds number turbulence near a wall is presented. In the first part of this work, the near-wall asymptotic behavior of the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation are analyzed. Based on these analyses, a modified eddy viscosity model with the correct near-wall behavior is suggested, and a model for the pressure transport term in the k-equation is proposed. In addition, a modeled dissipation rate equation is reformulated, and a boundary condition for the dissipation rate is suggested. In the second part of the work, one of the deficiencies of the existing k-epsilon models, namely, the wall distance dependency of the equations and the damping functions, is examined. An improved model that does not depend on any wall distance is introduced. Fully developed turbulent channel flows and turbulent boundary layers over a flat plate are studied as validations for the proposed new models. Numerical results obtained from the present and other previous k-epsilon models are compared with data from direct numerical simulation. The results show that the present k-epsilon model, with added robustness, performs as well as or better than other existing models in predicting the behavior of near-wall turbulence.

Modeling molecular mixing in a spatially inhomogeneous turbulent flow

NASA Astrophysics Data System (ADS)

Meyer, Daniel W.; Deb, Rajdeep

2012-02-01

Simulations of spatially inhomogeneous turbulent mixing in decaying grid turbulence with a joint velocity-concentration probability density function (PDF) method were conducted. The inert mixing scenario involves three streams with different compositions. The mixing model of Meyer ["A new particle interaction mixing model for turbulent dispersion and turbulent reactive flows," Phys. Fluids 22(3), 035103 (2010)], the interaction by exchange with the mean (IEM) model and its velocity-conditional variant, i.e., the IECM model, were applied. For reference, the direct numerical simulation data provided by Sawford and de Bruyn Kops ["Direct numerical simulation and lagrangian modeling of joint scalar statistics in ternary mixing," Phys. Fluids 20(9), 095106 (2008)] was used. It was found that velocity conditioning is essential to obtain accurate concentration PDF predictions. Moreover, the model of Meyer provides significantly better results compared to the IECM model at comparable computational expense.

Turbulence modeling of gas-solid suspension flows

NASA Technical Reports Server (NTRS)

Chen, C. P.

1988-01-01

The purpose here is to discuss and review advances in two-phase turbulent modeling techniques and their applications in various gas-solid suspension flow situations. In addition to the turbulence closures, heat transfer effect, particle dispersion and wall effects are partially covered.

Turbulence Modeling Validation, Testing, and Development

NASA Technical Reports Server (NTRS)

Bardina, J. E.; Huang, P. G.; Coakley, T. J.

1997-01-01

The primary objective of this work is to provide accurate numerical solutions for selected flow fields and to compare and evaluate the performance of selected turbulence models with experimental results. Four popular turbulence models have been tested and validated against experimental data often turbulent flows. The models are: (1) the two-equation k-epsilon model of Wilcox, (2) the two-equation k-epsilon model of Launder and Sharma, (3) the two-equation k-omega/k-epsilon SST model of Menter, and (4) the one-equation model of Spalart and Allmaras. The flows investigated are five free shear flows consisting of a mixing layer, a round jet, a plane jet, a plane wake, and a compressible mixing layer; and five boundary layer flows consisting of an incompressible flat plate, a Mach 5 adiabatic flat plate, a separated boundary layer, an axisymmetric shock-wave/boundary layer interaction, and an RAE 2822 transonic airfoil. The experimental data for these flows are well established and have been extensively used in model developments. The results are shown in the following four sections: Part A describes the equations of motion and boundary conditions; Part B describes the model equations, constants, parameters, boundary conditions, and numerical implementation; and Parts C and D describe the experimental data and the performance of the models in the free-shear flows and the boundary layer flows, respectively.

Chaotic Lagrangian models for turbulent relative dispersion.

PubMed

Lacorata, Guglielmo; Vulpiani, Angelo

2017-04-01

A deterministic multiscale dynamical system is introduced and discussed as a prototype model for relative dispersion in stationary, homogeneous, and isotropic turbulence. Unlike stochastic diffusion models, here trajectory transport and mixing properties are entirely controlled by Lagrangian chaos. The anomalous "sweeping effect," a known drawback common to kinematic simulations, is removed through the use of quasi-Lagrangian coordinates. Lagrangian dispersion statistics of the model are accurately analyzed by computing the finite-scale Lyapunov exponent (FSLE), which is the optimal measure of the scaling properties of dispersion. FSLE scaling exponents provide a severe test to decide whether model simulations are in agreement with theoretical expectations and/or observation. The results of our numerical experiments cover a wide range of "Reynolds numbers" and show that chaotic deterministic flows can be very efficient, and numerically low-cost, models of turbulent trajectories in stationary, homogeneous, and isotropic conditions. The mathematics of the model is relatively simple, and, in a geophysical context, potential applications may regard small-scale parametrization issues in general circulation models, mixed layer, and/or boundary layer turbulence models as well as Lagrangian predictability studies.

Chaotic Lagrangian models for turbulent relative dispersion

NASA Astrophysics Data System (ADS)

Lacorata, Guglielmo; Vulpiani, Angelo

2017-04-01

A deterministic multiscale dynamical system is introduced and discussed as a prototype model for relative dispersion in stationary, homogeneous, and isotropic turbulence. Unlike stochastic diffusion models, here trajectory transport and mixing properties are entirely controlled by Lagrangian chaos. The anomalous "sweeping effect," a known drawback common to kinematic simulations, is removed through the use of quasi-Lagrangian coordinates. Lagrangian dispersion statistics of the model are accurately analyzed by computing the finite-scale Lyapunov exponent (FSLE), which is the optimal measure of the scaling properties of dispersion. FSLE scaling exponents provide a severe test to decide whether model simulations are in agreement with theoretical expectations and/or observation. The results of our numerical experiments cover a wide range of "Reynolds numbers" and show that chaotic deterministic flows can be very efficient, and numerically low-cost, models of turbulent trajectories in stationary, homogeneous, and isotropic conditions. The mathematics of the model is relatively simple, and, in a geophysical context, potential applications may regard small-scale parametrization issues in general circulation models, mixed layer, and/or boundary layer turbulence models as well as Lagrangian predictability studies.

Optical linear algebra processors: noise and error-source modeling.

PubMed

Casasent, D; Ghosh, A

1985-06-01

The modeling of system and component noise and error sources in optical linear algebra processors (OLAP's) are considered, with attention to the frequency-multiplexed OLAP. General expressions are obtained for the output produced as a function of various component errors and noise. A digital simulator for this model is discussed.

Optical linear algebra processors - Noise and error-source modeling

NASA Technical Reports Server (NTRS)

Casasent, D.; Ghosh, A.

1985-01-01

The modeling of system and component noise and error sources in optical linear algebra processors (OLAPs) are considered, with attention to the frequency-multiplexed OLAP. General expressions are obtained for the output produced as a function of various component errors and noise. A digital simulator for this model is discussed.

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Spalart, Philippe R.

2008-01-01

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave transitionally in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Spalart, Philippe R.

2008-01-01

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave 'transitionally' in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Examination of various turbulence models for application in liquid rocket thrust chambers

NASA Technical Reports Server (NTRS)

Hung, R. J.

1991-01-01

There is a large variety of turbulence models available. These models include direct numerical simulation, large eddy simulation, Reynolds stress/flux model, zero equation model, one equation model, two equation k-epsilon model, multiple-scale model, etc. Each turbulence model contains different physical assumptions and requirements. The natures of turbulence are randomness, irregularity, diffusivity and dissipation. The capabilities of the turbulence models, including physical strength, weakness, limitations, as well as numerical and computational considerations, are reviewed. Recommendations are made for the potential application of a turbulence model in thrust chamber and performance prediction programs. The full Reynolds stress model is recommended. In a workshop, specifically called for the assessment of turbulence models for applications in liquid rocket thrust chambers, most of the experts present were also in favor of the recommendation of the Reynolds stress model.

Fractional Order Modeling of Atmospheric Turbulence - A More Accurate Modeling Methodology for Aero Vehicles

NASA Technical Reports Server (NTRS)

Kopasakis, George

2014-01-01

The presentation covers a recently developed methodology to model atmospheric turbulence as disturbances for aero vehicle gust loads and for controls development like flutter and inlet shock position. The approach models atmospheric turbulence in their natural fractional order form, which provides for more accuracy compared to traditional methods like the Dryden model, especially for high speed vehicle. The presentation provides a historical background on atmospheric turbulence modeling and the approaches utilized for air vehicles. This is followed by the motivation and the methodology utilized to develop the atmospheric turbulence fractional order modeling approach. Some examples covering the application of this method are also provided, followed by concluding remarks.

Turbulent transport modeling of shear flows around an aerodynamic wing. Development of turbulent near-wall model and its application to recirculating flows

NASA Technical Reports Server (NTRS)

Amano, R. S.

1982-01-01

Progress in implementing and refining two near-wall turbulence models in which the near-wall region is divided into either two or three zones is outlined. These models were successfully applied to the computation of recirculating flows. The research was further extended to obtaining experimental results of two different recirculating flow conditions in order to check the validity of the present models. Two different experimental apparatuses were set up: axisymmetric turbulent impinging jets on a flat plate, and turbulent flows in a circular pipe with a abrupt pipe expansion. It is shown that generally better results are obtained by using the present near-wall models, and among the models the three-zone model is superior to the two-zone model.

A Galilean and tensorial invariant k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1993-01-01

A k-epsilon model is proposed for wall bounded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation rate equation is reformulated using this time scale and no singularity exists at the wall. A new parameter R = k/S(nu) is introduced to characterize the damping function in the eddy viscosity. This parameter is determined by local properties of both the mean and the turbulent flow fields and is free from any geometry parameter. The proposed model is then Galilean and tensorial invariant. The model constants used are the same as in the high Reynolds number Standard k-epsilon Model. Thus, the proposed model will also be suitable for flows far from the wall. Turbulent channel flows and turbulent boundary layer flows with and without pressure gradients are calculated. Comparisons with the data from direct numerical simulations and experiments show that the model predictions are excellent for turbulent channel flows and turbulent boundary layers with favorable pressure gradients, good for turbulent boundary layers with zero pressure gradients, and fair for turbulent boundary layer with adverse pressure gradients.

The analysis and modelling of dilatational terms in compressible turbulence

NASA Technical Reports Server (NTRS)

Sarkar, S.; Erlebacher, G.; Hussaini, M. Y.; Kreiss, H. O.

1991-01-01

It is shown that the dilatational terms that need to be modeled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of these dilatational terms in homogeneous turbulence is explored by asymptotic analysis of the compressible Navier-Stokes equations. A non-dimensional parameter which characterizes some compressible effects in moderate Mach number, homogeneous turbulence is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

The analysis and modeling of dilatational terms in compressible turbulence

NASA Technical Reports Server (NTRS)

Sarkar, S.; Erlebacher, G.; Hussaini, M. Y.; Kreiss, H. O.

1989-01-01

It is shown that the dilatational terms that need to be modeled in compressible turbulence include not only the pressure-dilatation term but also another term - the compressible dissipation. The nature of these dilatational terms in homogeneous turbulence is explored by asymptotic analysis of the compressible Navier-Stokes equations. A non-dimensional parameter which characterizes some compressible effects in moderate Mach number, homogeneous turbulence is identified. Direct numerical simulations (DNS) of isotropic, compressible turbulence are performed, and their results are found to be in agreement with the theoretical analysis. A model for the compressible dissipation is proposed; the model is based on the asymptotic analysis and the direct numerical simulations. This model is calibrated with reference to the DNS results regarding the influence of compressibility on the decay rate of isotropic turbulence. An application of the proposed model to the compressible mixing layer has shown that the model is able to predict the dramatically reduced growth rate of the compressible mixing layer.

Particle dispersion in homogeneous turbulence using the one-dimensional turbulence model

DOE PAGES

Sun, Guangyuan; Lignell, David O.; Hewson, John C.; ...

2014-10-09

Lagrangian particle dispersion is studied using the one-dimensional turbulence (ODT) model in homogeneous decaying turbulence configurations. The ODT model has been widely and successfully applied to a number of reacting and nonreacting flow configurations, but only limited application has been made to multiphase flows. We present a version of the particle implementation and interaction with the stochastic and instantaneous ODT eddy events. The model is characterized by comparison to experimental data of particle dispersion for a range of intrinsic particle time scales and body forces. Particle dispersion, velocity, and integral time scale results are presented. Moreover, the particle implementation introducesmore » a single model parameter β p , and sensitivity to this parameter and behavior of the model are discussed. Good agreement is found with experimental data and the ODT model is able to capture the particle inertial and trajectory crossing effects. Our results serve as a validation case of the multiphase implementations of ODT for extensions to other flow configurations.« less

Modeling of confined turbulent fluid-particle flows using Eulerian and Lagrangian schemes

NASA Technical Reports Server (NTRS)

Adeniji-Fashola, A.; Chen, C. P.

1990-01-01

Two important aspects of fluid-particulate interaction in dilute gas-particle turbulent flows (the turbulent particle dispersion and the turbulence modulation effects) are addressed, using the Eulerian and Lagrangian modeling approaches to describe the particulate phase. Gradient-diffusion approximations are employed in the Eulerian formulation, while a stochastic procedure is utilized to simulate turbulent dispersion in the Lagrangina formulation. The k-epsilon turbulence model is used to characterize the time and length scales of the continuous phase turbulence. Models proposed for both schemes are used to predict turbulent fully-developed gas-solid vertical pipe flow with reasonable accuracy.

A Model for Jet-Surface Interaction Noise Using Physically Realizable Upstream Turbulence Conditions

NASA Technical Reports Server (NTRS)

Afsar, Mohammed Z.; Leib, Stewart J.; Bozak, Richard F.

2016-01-01

This paper is a continuation of previous work in which a generalized Rapid Distortion Theory (RDT) formulation was used to model low-frequency trailing-edge noise. The research was motivated by proposed next-generation aircraft configurations where the exhaust system is tightly integrated with the airframe. Data from recent experiments at NASA on the interaction between high-Reynolds-number subsonic jet flows and an external flat plate showed that the power spectral density (PSD) of the far-field pressure underwent considerable amplification at low frequencies. For example, at the 90deg observation angle, the low-frequency noise could be as much as 10 dB greater than the jet noise itself. In this paper, we present predictions of the noise generated by the interaction of a rectangular jet with the trailing edge of a semi-infinite flat plate. The calculations are based on a formula for the acoustic spectrum of this noise source derived from an exact formal solution of the linearized Euler equations involving (in this case) one arbitrary convected scalar quantity and a Rayleigh equation Green's function. A low-frequency asymptotic approximation for the Green's function based on a two-dimensional mean flow is used in the calculations along with a physically realizable upstream turbulence spectrum, which includes a finite decorrelation region. Numerical predictions of the sound field, based on three-dimensional RANS solutions to determine the mean flow, turbulent kinetic energy and turbulence length and time scales, for a range of subsonic acoustic Mach number jets and nozzle aspect ratios are compared with experimental data. Comparisons of the RANS results with flow data are also presented for selected cases. We find that a finite decorrelation region in the turbulence spectrum increases the low-frequency algebraic decay (the low frequency "roll-off") of the acoustic spectrum with angular frequency thereby producing much closer agreement with noise data for Strouhal

A review of Reynolds stress models for turbulent shear flows

NASA Technical Reports Server (NTRS)

Speziale, Charles G.

1995-01-01

A detailed review of recent developments in Reynolds stress modeling for incompressible turbulent shear flows is provided. The mathematical foundations of both two-equation models and full second-order closures are explored in depth. It is shown how these models can be systematically derived for two-dimensional mean turbulent flows that are close to equilibrium. A variety of examples are provided to demonstrate how well properly calibrated versions of these models perform for such flows. However, substantial problems remain for the description of more complex turbulent flows where there are large departures from equilibrium. Recent efforts to extend Reynolds stress models to nonequilibrium turbulent flows are discussed briefly along with the major modeling issues relevant to practical naval hydrodynamics applications.

Coupled turbulence and aerosol dynamics modeling of vehicle exhaust plumes using the CTAG model

NASA Astrophysics Data System (ADS)

Wang, Yan Jason; Zhang, K. Max

2012-11-01

This paper presents the development and evaluation of an environmental turbulent reacting flow model, the Comprehensive Turbulent Aerosol Dynamics and Gas Chemistry (CTAG) model. CTAG is designed to simulate transport and transformation of multiple air pollutants, e.g., from emission sources to ambient background. For the on-road and near-road applications, CTAG explicitly couples the major turbulent mixing processes, i.e., vehicle-induced turbulence (VIT), road-induced turbulence (RIT) and atmospheric boundary layer turbulence with gas-phase chemistry and aerosol dynamics. CTAG's transport model is referred to as CFD-VIT-RIT. This paper presents the evaluation of the CTAG model in simulating the dynamics of individual plumes in the “tailpipe-to-road” stage, i.e., VIT behind a moving van and aerosol dynamics in the wake of a diesel car by comparing the modeling results against the respective field measurements. Combined with sensitivity studies, we analyze the relative roles of VIT, sulfuric acid induced nucleation, condensation of organic compounds and presence of soot-mode particles in capturing the dynamics of exhaust plumes as well as their implications in vehicle emission controls.

Helicity statistics in homogeneous and isotropic turbulence and turbulence models

NASA Astrophysics Data System (ADS)

Sahoo, Ganapati; De Pietro, Massimo; Biferale, Luca

2017-02-01

We study the statistical properties of helicity in direct numerical simulations of fully developed homogeneous and isotropic turbulence and in a class of turbulence shell models. We consider correlation functions based on combinations of vorticity and velocity increments that are not invariant under mirror symmetry. We also study the scaling properties of high-order structure functions based on the moments of the velocity increments projected on a subset of modes with either positive or negative helicity (chirality). We show that mirror symmetry is recovered at small scales, i.e., chiral terms are subleading and they are well captured by a dimensional argument plus anomalous corrections. These findings are also supported by a high Reynolds numbers study of helical shell models with the same chiral symmetry of Navier-Stokes equations.

Lattice Virasoro algebra and corner transfer matrices in the Baxter eight-vertex model

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

Itoyama, H.; Thacker, H.B.

1987-04-06

A lattice Virasoro algebra is constructed for the Baxter eight-vertex model. The operator L/sub 0/ is obtained from the logarithm of the corner transfer matrix and is given by the first moment of the XYZ spin-chain Hamiltonian. The algebra is valid even when the Hamiltonian includes a mass term, in which case it represents lattice coordinate transformations which distinguish between even and odd sublattices. We apply the quantum inverse scattering method to demonstrate that the Virasoro algebra follows from the Yang-Baxter relations.

The Lag Model, a Turbulence Model for Wall Bounded Flows Including Separation

NASA Technical Reports Server (NTRS)

Olsen, Michael E.; Coakley, Thomas J.; Kwak, Dochan (Technical Monitor)

2001-01-01

A new class of turbulence model is described for wall bounded, high Reynolds number flows. A specific turbulence model is demonstrated, with results for favorable and adverse pressure gradient flowfields. Separation predictions are as good or better than either Spalart Almaras or SST models, do not require specification of wall distance, and have similar or reduced computational effort compared with these models.

A Two-length Scale Turbulence Model for Single-phase Multi-fluid Mixing

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

Schwarzkopf, J. D.; Livescu, D.; Baltzer, J. R.

2015-09-08

A two-length scale, second moment turbulence model (Reynolds averaged Navier-Stokes, RANS) is proposed to capture a wide variety of single-phase flows, spanning from incompressible flows with single fluids and mixtures of different density fluids (variable density flows) to flows over shock waves. The two-length scale model was developed to address an inconsistency present in the single-length scale models, e.g. the inability to match both variable density homogeneous Rayleigh-Taylor turbulence and Rayleigh-Taylor induced turbulence, as well as the inability to match both homogeneous shear and free shear flows. The two-length scale model focuses on separating the decay and transport length scales,more » as the two physical processes are generally different in inhomogeneous turbulence. This allows reasonable comparisons with statistics and spreading rates over such a wide range of turbulent flows using a common set of model coefficients. The specific canonical flows considered for calibrating the model include homogeneous shear, single-phase incompressible shear driven turbulence, variable density homogeneous Rayleigh-Taylor turbulence, Rayleigh-Taylor induced turbulence, and shocked isotropic turbulence. The second moment model shows to compare reasonably well with direct numerical simulations (DNS), experiments, and theory in most cases. The model was then applied to variable density shear layer and shock tube data and shows to be in reasonable agreement with DNS and experiments. Additionally, the importance of using DNS to calibrate and assess RANS type turbulence models is highlighted.« less

Compressible Turbulent Channel Flows: DNS Results and Modeling

NASA Technical Reports Server (NTRS)

Huang, P. G.; Coleman, G. N.; Bradshaw, P.; Rai, Man Mohan (Technical Monitor)

1994-01-01

The present paper addresses some topical issues in modeling compressible turbulent shear flows. The work is based on direct numerical simulation of two supersonic fully developed channel flows between very cold isothermal walls. Detailed decomposition and analysis of terms appearing in the momentum and energy equations are presented. The simulation results are used to provide insights into differences between conventional time-and Favre-averaging of the mean-flow and turbulent quantities. Study of the turbulence energy budget for the two cases shows that the compressibility effects due to turbulent density and pressure fluctuations are insignificant. In particular, the dilatational dissipation and the mean product of the pressure and dilatation fluctuations are very small, contrary to the results of simulations for sheared homogeneous compressible turbulence and to recent proposals for models for general compressible turbulent flows. This provides a possible explanation of why the Van Driest density-weighted transformation is so successful in correlating compressible boundary layer data. Finally, it is found that the DNS data do not support the strong Reynolds analogy. A more general representation of the analogy is analysed and shown to match the DNS data very well.

Modeling of turbulent transport as a volume process

NASA Technical Reports Server (NTRS)

Jennings, Mark J.; Morel, Thomas

1987-01-01

An alternative type of modeling was proposed for the turbulent transport terms in Reynolds-averaged equations. One particular implementation of the model was considered, based on the two-point velocity correlations. The model was found to reproduce the trends but not the magnitude of the nonisotropic behavior of the turbulent transport. Some interesting insights were developed concerning the shape of the contracted two-point correlation volume. This volume is strongly deformed by mean shear from the spherical shape found in unstrained flows. Of particular interest is the finding that the shape is sharply waisted, indicating preferential lines of communication, which should have a direct effect on turbulent transfer and on other processes.

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.

A Particle Representation Model for the Deformation of Homogeneous Turbulence

NASA Technical Reports Server (NTRS)

Kassinos, S. C.; Reynolds, W. C.

1996-01-01

In simple flows, where the mean deformation rates are mild and the turbulence has time to come to equilibrium with the mean flow, the Reynolds stresses are determined by the applied strain rate. Hence in these flows, it is often adequate to use an eddy-viscosity representation. The modern family of kappa-epsilon models has been very useful in predicting near equilibrium turbulent flows, where the rms deformation rate S is small compared to the reciprocal time scale of the turbulence (epsilon/kappa). In modern engineering applications, turbulence models are quite often required to predict flows with very rapid deformations (large S kappa/epsilon). In these flows, the structure takes some time to respond and eddy viscosity models are inadequate. The response of turbulence to rapid deformations is given by rapid distortion theory (RDT). Under RDT the nonlinear effects due to turbulence-turbulence interactions are neglected in the governing equations, but even when linearized in this fashion, the governing equations are unclosed at the one-point level due to the non-locality of the pressure fluctuations.

Consequences of Symmetries on the Analysis and Construction of Turbulence Models

NASA Astrophysics Data System (ADS)

Razafindralandy, Dina; Hamdouni, Aziz

2006-05-01

Since they represent fundamental physical properties in turbulence (conservation laws, wall laws, Kolmogorov energy spectrum, ...), symmetries are used to analyse common turbulence models. A class of symmetry preserving turbulence models is proposed. This class is refined such that the models respect the second law of thermodynamics. Finally, an example of model belonging to the class is numerically tested.

Analysis of Unsteady Simulations to Inform Turbulence Modeling

NASA Technical Reports Server (NTRS)

Vyas, Manan; Waindim, Mbu; Gaitonde, Datta

2016-01-01

In this work, budgets of the turbulent kinetic energy are presented for a two-dimensional shock wave boundary-layer interaction (SBLI). The work should be of interest to the SBLI research and turbulence modeling community.

Modelling of nanoscale quantum tunnelling structures using algebraic topology method

NASA Astrophysics Data System (ADS)

Sankaran, Krishnaswamy; Sairam, B.

2018-05-01

We have modelled nanoscale quantum tunnelling structures using Algebraic Topology Method (ATM). The accuracy of ATM is compared to the analytical solution derived based on the wave nature of tunnelling electrons. ATM provides a versatile, fast, and simple model to simulate complex structures. We are currently expanding the method for modelling electrodynamic systems.

Modelling Unsteady Wall Pressures Beneath Turbulent Boundary Layers

NASA Technical Reports Server (NTRS)

Ahn, B-K.; Graham, W. R.; Rizzi, S. A.

2004-01-01

As a structural entity of turbulence, hairpin vortices are believed to play a major role in developing and sustaining the turbulence process in the near wall region of turbulent boundary layers and may be regarded as the simplest conceptual model that can account for the essential features of the wall pressure fluctuations. In this work we focus on fully developed typical hairpin vortices and estimate the associated surface pressure distributions and their corresponding spectra. On the basis of the attached eddy model, we develop a representation of the overall surface pressure spectra in terms of the eddy size distribution. Instantaneous wavenumber spectra and spatial correlations are readily derivable from this representation. The model is validated by comparison of predicted wavenumber spectra and cross-correlations with existing emperical models and experimental data.

Highest weight representation for Sklyanin algebra sl(3)(u) with application to the Gaudin model

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

Burdik, C., E-mail: burdik@kmlinux.fjfi.cvut.cz; Navratil, O.

2011-06-15

We study the infinite-dimensional Sklyanin algebra sl(3)(u). Specifically we construct the highest weight representation for this algebra in an explicit form. Its application to the Gaudin model is mentioned.

Stratified turbulent Bunsen flames: flame surface analysis and flame surface density modelling

NASA Astrophysics Data System (ADS)

Ramaekers, W. J. S.; van Oijen, J. A.; de Goey, L. P. H.

2012-12-01

In this paper it is investigated whether the Flame Surface Density (FSD) model, developed for turbulent premixed combustion, is also applicable to stratified flames. Direct Numerical Simulations (DNS) of turbulent stratified Bunsen flames have been carried out, using the Flamelet Generated Manifold (FGM) reduction method for reaction kinetics. Before examining the suitability of the FSD model, flame surfaces are characterized in terms of thickness, curvature and stratification. All flames are in the Thin Reaction Zones regime, and the maximum equivalence ratio range covers 0.1⩽φ⩽1.3. For all flames, local flame thicknesses correspond very well to those observed in stretchless, steady premixed flamelets. Extracted curvature radii and mixing length scales are significantly larger than the flame thickness, implying that the stratified flames all burn in a premixed mode. The remaining challenge is accounting for the large variation in (subfilter) mass burning rate. In this contribution, the FSD model is proven to be applicable for Large Eddy Simulations (LES) of stratified flames for the equivalence ratio range 0.1⩽φ⩽1.3. Subfilter mass burning rate variations are taken into account by a subfilter Probability Density Function (PDF) for the mixture fraction, on which the mass burning rate directly depends. A priori analysis point out that for small stratifications (0.4⩽φ⩽1.0), the replacement of the subfilter PDF (obtained from DNS data) by the corresponding Dirac function is appropriate. Integration of the Dirac function with the mass burning rate m=m(φ), can then adequately model the filtered mass burning rate obtained from filtered DNS data. For a larger stratification (0.1⩽φ⩽1.3), and filter widths up to ten flame thicknesses, a β-function for the subfilter PDF yields substantially better predictions than a Dirac function. Finally, inclusion of a simple algebraic model for the FSD resulted only in small additional deviations from DNS data

Recent Developments on the Turbulence Modeling Resource Website (Invited)

NASA Technical Reports Server (NTRS)

Rumssey, Christopher L.

2015-01-01

The NASA Langley Turbulence Model Resource (TMR) website has been active for over five years. Its main goal of providing a one-stop, easily accessible internet site for up-to-date information on Reynolds-averaged Navier-Stokes turbulence models remains unchanged. In particular, the site strives to provide an easy way for users to verify their own implementations of widely-used turbulence models, and to compare the results from different models for a variety of simple unit problems covering a range of flow physics. Some new features have been recently added to the website. This paper documents the site's features, including recent developments, future plans, and open questions.

The significance of turbulent flow representation in single-continuum models

USGS Publications Warehouse

Reimann, T.; Rehrl, C.; Shoemaker, W.B.; Geyer, T.; Birk, S.

2011-01-01

Karst aquifers exhibit highly conductive features caused from rock dissolution processes. Flow within these structures can become turbulent and therefore can be expressed by nonlinear gradient functions. One way to account for these effects is by coupling a continuum model with a conduit network. Alternatively, turbulent flow can be considered by adapting the hydraulic conductivity within the continuum model. Consequently, the significance of turbulent flow on the dynamic behavior of karst springs is investigated by an enhanced single-continuum model that results in conduit-type flow in continuum cells (CTFC). The single-continuum approach CTFC represents laminar and turbulent flow as well as more complex hybrid models that require additional programming and numerical efforts. A parameter study is conducted to investigate the effects of turbulent flow on the response of karst springs to recharge events using the new CTFC approach, existing hybrid models, and MODFLOW-2005. Results reflect the importance of representing (1) turbulent flow in karst conduits and (2) the exchange between conduits and continuum cells. More specifically, laminar models overestimate maximum spring discharge and underestimate hydraulic gradients within the conduit. It follows that aquifer properties inferred from spring hydrographs are potentially impaired by ignoring flow effects due to turbulence. The exchange factor used for hybrid models is necessary to account for the scale dependency between hydraulic properties of the matrix continuum and conduits. This functionality, which is not included in CTFC, can be mimicked by appropriate use of the Horizontal Flow Barrier package for MODFLOW. Copyright 2011 by the American Geophysical Union.

Assessment of CFD Hypersonic Turbulent Heating Rates for Space Shuttle Orbiter

NASA Technical Reports Server (NTRS)

Wood, William A.; Oliver, A. Brandon

2011-01-01

Turbulent CFD codes are assessed for the prediction of convective heat transfer rates at turbulent, hypersonic conditions. Algebraic turbulence models are used within the DPLR and LAURA CFD codes. The benchmark heat transfer rates are derived from thermocouple measurements of the Space Shuttle orbiter Discovery windward tiles during the STS-119 and STS-128 entries. The thermocouples were located underneath the reaction-cured glass coating on the thermal protection tiles. Boundary layer transition flight experiments conducted during both of those entries promoted turbulent flow at unusually high Mach numbers, with the present analysis considering Mach 10{15. Similar prior comparisons of CFD predictions directly to the flight temperature measurements were unsatisfactory, showing diverging trends between prediction and measurement for Mach numbers greater than 11. In the prior work, surface temperatures and convective heat transfer rates had been assumed to be in radiative equilibrium. The present work employs a one-dimensional time-accurate conduction analysis to relate measured temperatures to surface heat transfer rates, removing heat soak lag from the flight data, in order to better assess the predictive accuracy of the numerical models. The turbulent CFD shows good agreement for turbulent fuselage flow up to Mach 13. But on the wing in the wake of the boundary layer trip, the inclusion of tile conduction effects does not explain the prior observed discrepancy in trends between simulation and experiment; the flight heat transfer measurements are roughly constant over Mach 11-15, versus an increasing trend with Mach number from the CFD.

Potential capabilities of Reynolds stress turbulence model in the COMMIX-RSM code

NASA Technical Reports Server (NTRS)

Chang, F. C.; Bottoni, M.

1994-01-01

A Reynolds stress turbulence model has been implemented in the COMMIX code, together with transport equations describing turbulent heat fluxes, variance of temperature fluctuations, and dissipation of turbulence kinetic energy. The model has been verified partially by simulating homogeneous turbulent shear flow, and stable and unstable stratified shear flows with strong buoyancy-suppressing or enhancing turbulence. This article outlines the model, explains the verifications performed thus far, and discusses potential applications of the COMMIX-RSM code in several domains, including, but not limited to, analysis of thermal striping in engineering systems, simulation of turbulence in combustors, and predictions of bubbly and particulate flows.

Report on PDF Models for Turbulence Chemistry Interaction

DTIC Science & Technology

2014-03-01

significantly within the flowfield (like rocket plumes or scramjet combustors). For multi-species flows turbulence can increase the apparent mass...Variable Turbulent Schmidt-Number Formulation for Scramjet Applications, AIAA Journal, 44(3), 593–599. [12] Xiao, X., Hassan, H.A., and Baurle, R.A...2006), Modeling Scramjet Flows with Variable Turbulent Prandtl and Schmidt Numbers. AIAA Paper 2006-128. [13] Xiao, X., Hassan, H.A., and Baurle, R.A

Fan noise caused by the ingestion of anisotropic turbulence - A model based on axisymmetric turbulence theory

NASA Technical Reports Server (NTRS)

Kerschen, E. J.; Gliebe, P. R.

1980-01-01

An analytical model of fan noise caused by inflow turbulence, a generalization of earlier work by Mani, is presented. Axisymmetric turbulence theory is used to develop a statistical representation of the inflow turbulence valid for a wide range of turbulence properties. Both the dipole source due to rotor blade unsteady forces and the quadrupole source resulting from the interaction of the turbulence with the rotor potential field are considered. The effects of variations in turbulence properties and fan operating conditions are evaluated. For turbulence axial integral length scales much larger than the blade spacing, the spectrum is shown to consist of sharp peaks at the blade passing frequency and its harmonics, with negligible broadband content. The analysis can then be simplified considerably and the total sound power contained within each spectrum peak becomes independent of axial length scale, while the width of the peak is inversely proportional to this parameter. Large axial length scales are characteristic of static fan test facilities, where the transverse contraction of the inlet flow produces highly anisotropic turbulence. In this situation, the rotor/turbulence interaction noise is mainly caused by the transverse component of turbulent velocity.

Effect of Turbulence Modeling on Hovering Rotor Flows

NASA Technical Reports Server (NTRS)

Yoon, Seokkwan; Chaderjian, Neal M.; Pulliam, Thomas H.; Holst, Terry L.

2015-01-01

The effect of turbulence models in the off-body grids on the accuracy of solutions for rotor flows in hover has been investigated. Results from the Reynolds-Averaged Navier-Stokes and Laminar Off-Body models are compared. Advection of turbulent eddy viscosity has been studied to find the mechanism leading to inaccurate solutions. A coaxial rotor result is also included.

Turbulence and transition modeling for high-speed flows

NASA Technical Reports Server (NTRS)

Wilcox, David C.

1993-01-01

Research conducted during the past three and a half years aimed at developing and testing a turbulence/transition model applicable to high-speed turbulent flows is summarized. The first two years of the project focused on fully turbulent flows, while emphasis shifted to boundary-layer development in the transition region during the final year and a half. A brief summary of research accomplished during the first three years is included and publications that describe research results in greater detail are cited. Research conducted during the final six months of the period of performance is summarized. The primary results of the last six months of the project are elimination of the k-omega model's sensitivity to the freestream value of omega and development of a method for triggering transition at a specified location, independent of the freestream turbulence level.

Enlarged symmetry algebras of spin chains, loop models, and S-matrices

NASA Astrophysics Data System (ADS)

Read, N.; Saleur, H.

2007-08-01

The symmetry algebras of certain families of quantum spin chains are considered in detail. The simplest examples possess m states per site ( m⩾2), with nearest-neighbor interactions with U(m) symmetry, under which the sites transform alternately along the chain in the fundamental m and its conjugate representation m¯. We find that these spin chains, even with arbitrary coefficients of these interactions, have a symmetry algebra A much larger than U(m), which implies that the energy eigenstates fall into sectors that for open chains (i.e., free boundary conditions) can be labeled by j=0,1,…,L, for the 2 L-site chain such that the degeneracies of all eigenvalues in the jth sector are generically the same and increase rapidly with j. For large j, these degeneracies are much larger than those that would be expected from the U(m) symmetry alone. The enlarged symmetry algebra A(2L) consists of operators that commute in this space of states with the Temperley-Lieb algebra that is generated by the set of nearest-neighbor interaction terms; A(2L) is not a Yangian. There are similar results for supersymmetric chains with gl(m+n|n) symmetry of nearest-neighbor interactions, and a richer representation structure for closed chains (i.e., periodic boundary conditions). The symmetries also apply to the loop models that can be obtained from the spin chains in a spacetime or transfer matrix picture. In the loop language, the symmetries arise because the loops cannot cross. We further define tensor products of representations (for the open chains) by joining chains end to end. The fusion rules for decomposing the tensor product of representations labeled j and j take the same form as the Clebsch-Gordan series for SU(2). This and other structures turn the symmetry algebra A into a ribbon Hopf algebra, and we show that this is "Morita equivalent" to the quantum group U(sl) for m=q+q. The open-chain results are extended to the cases |m|<2 for which the algebras are no longer

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.

The use of direct numerical simulation data in turbulence modeling

NASA Technical Reports Server (NTRS)

Mansour, N. N.

1991-01-01

Direct numerical simulations (DNS) of turbulent flows provide a complete data base to develop and to test turbulence models. In this article, the progress made in developing models for the dissipation rate equation is reviewed. New scaling arguments for the various terms in the dissipation rate equation were tested using data from DNS of homogeneous shear flows. Modifications to the epsilon-equation model that take into account near-wall effects were developed using DNS of turbulent channel flows. Testing of new models for flows under mean compression was carried out using data from DNS of isotropically compressed turbulence. In all of these studies the data from the simulations was essential in guiding the model development. The next generation of DNS will be at higher Reynolds numbers, and will undoubtedly lead to improved models for computations of flows of practical interest.

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

NASA Technical Reports Server (NTRS)

Canuto, V. M.

1994-01-01

latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

NASA Astrophysics Data System (ADS)

Canuto, V. M.

1994-06-01

, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Progress in modeling hypersonic turbulent boundary layers

NASA Technical Reports Server (NTRS)

Zeman, Otto

1993-01-01

A good knowledge of the turbulence structure, wall heat transfer, and friction in turbulent boundary layers (TBL) at high speeds is required for the design of hypersonic air breathing airplanes and reentry space vehicles. This work reports on recent progress in the modeling of high speed TBL flows. The specific research goal described here is the development of a second order closure model for zero pressure gradient TBL's for the range of Mach numbers up to hypersonic speeds with arbitrary wall cooling requirements.

The continuous adjoint approach to the k-ε turbulence model for shape optimization and optimal active control of turbulent flows

NASA Astrophysics Data System (ADS)

Papoutsis-Kiachagias, E. M.; Zymaris, A. S.; Kavvadias, I. S.; Papadimitriou, D. I.; Giannakoglou, K. C.

2015-03-01

The continuous adjoint to the incompressible Reynolds-averaged Navier-Stokes equations coupled with the low Reynolds number Launder-Sharma k-ε turbulence model is presented. Both shape and active flow control optimization problems in fluid mechanics are considered, aiming at minimum viscous losses. In contrast to the frequently used assumption of frozen turbulence, the adjoint to the turbulence model equations together with appropriate boundary conditions are derived, discretized and solved. This is the first time that the adjoint equations to the Launder-Sharma k-ε model have been derived. Compared to the formulation that neglects turbulence variations, the impact of additional terms and equations is evaluated. Sensitivities computed using direct differentiation and/or finite differences are used for comparative purposes. To demonstrate the need for formulating and solving the adjoint to the turbulence model equations, instead of merely relying upon the 'frozen turbulence assumption', the gain in the optimization turnaround time offered by the proposed method is quantified.

One-dimensional wave bottom boundary layer model comparison: specific eddy viscosity and turbulence closure models

USGS Publications Warehouse

Puleo, J.A.; Mouraenko, O.; Hanes, D.M.

2004-01-01

Six one-dimensional-vertical wave bottom boundary layer models are analyzed based on different methods for estimating the turbulent eddy viscosity: Laminar, linear, parabolic, k—one equation turbulence closure, k−ε—two equation turbulence closure, and k−ω—two equation turbulence closure. Resultant velocity profiles, bed shear stresses, and turbulent kinetic energy are compared to laboratory data of oscillatory flow over smooth and rough beds. Bed shear stress estimates for the smooth bed case were most closely predicted by the k−ω model. Normalized errors between model predictions and measurements of velocity profiles over the entire computational domain collected at 15° intervals for one-half a wave cycle show that overall the linear model was most accurate. The least accurate were the laminar and k−ε models. Normalized errors between model predictions and turbulence kinetic energy profiles showed that the k−ω model was most accurate. Based on these findings, when the smallest overall velocity profile prediction error is required, the processing requirements and error analysis suggest that the linear eddy viscosity model is adequate. However, if accurate estimates of bed shear stress and TKE are required then, of the models tested, the k−ω model should be used.

One-dimensional Turbulence Models of Type I X-ray Bursts

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

Hou, Chen

Type I X-ray bursts are caused by thermonuclear explosions occurring on the surface of an accreting neutron star in a binary star system. Observations and simulations of these phenomena are of great importance for understanding the fundamental properties of neutron stars and dense matter because the equation of state for cold dense matter can be constrained by the mass-radius relationship of neutron stars. During the bursts, turbulence plays a key role in mixing the fuels and driving the unstable nuclear burning process. This dissertation presents one-dimensional models of photospheric radius expansion bursts with a new approach to simulate turbulent advection.more » Compared with the traditional mixing length theory, the one-dimensional turbulence (ODT) model represents turbulent motions by a sequence of maps that are generated according to a stochastic process. The light curves I obtained with the ODT models are in good agreement with those of the KEPLER model in which the mixing length theory and various diffusive processes are applied. The abundance comparison, however, indicates that the differences in turbulent regions and turbulent diffusivities result in more 12C survival during the bursts in the ODT models, which can make a difference in the superbursts phenomena triggered by unstable carbon burning.« less

New time scale based k-epsilon model for near-wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1993-01-01

A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

Dynamical systems defined on infinite dimensional lie algebras of the ''current algebra'' or ''Kac-Moody'' type

NASA Astrophysics Data System (ADS)

Hermann, Robert

1982-07-01

Recent work by Morrison, Marsden, and Weinstein has drawn attention to the possibility of utilizing the cosymplectic structure of the dual of the Lie algebra of certain infinite dimensional Lie groups to study hydrodynamical and plasma systems. This paper treats certain models arising in elementary particle physics, considered by Lee, Weinberg, and Zumino; Sugawara; Bardacki, Halpern, and Frishman; Hermann; and Dolan. The lie algebras involved are associated with the ''current algebras'' of Gell-Mann. This class of Lie algebras contains certain of the algebras that are called ''Kac-Moody algebras'' in the recent mathematics and mathematical physics literature.

Modeling of Fine-Particle Formation in Turbulent Flames

NASA Astrophysics Data System (ADS)

Raman, Venkat; Fox, Rodney O.

2016-01-01

The generation of nanostructured particles in high-temperature flames is important both for the control of emissions from combustion devices and for the synthesis of high-value chemicals for a variety of applications. The physiochemical processes that lead to the production of fine particles in turbulent flames are highly sensitive to the flow physics and, in particular, the history of thermochemical compositions and turbulent features they encounter. Consequently, it is possible to change the characteristic size, structure, composition, and yield of the fine particles by altering the flow configuration. This review describes the complex multiscale interactions among turbulent fluid flow, gas-phase chemical reactions, and solid-phase particle evolution. The focus is on modeling the generation of soot particles, an unwanted pollutant from automobile and aircraft engines, as well as metal oxides, a class of high-value chemicals sought for specialized applications, including emissions control. Issues arising due to the numerical methods used to approximate the particle number density function, the modeling of turbulence-chemistry interactions, and model validation are also discussed.

On the modeling of wave-enhanced turbulence nearshore

NASA Astrophysics Data System (ADS)

Moghimi, Saeed; Thomson, Jim; Özkan-Haller, Tuba; Umlauf, Lars; Zippel, Seth

2016-07-01

A high resolution k-ω two-equation turbulence closure model, including surface wave forcing was employed to fully resolve turbulence dissipation rate profiles close to the ocean surface. Model results were compared with observations from Surface Wave Instrument Floats with Tracking (SWIFTs) in the nearshore region at New River Inlet, North Carolina USA, in June 2012. A sensitivity analysis for different physical parameters and wave and turbulence formulations was performed. The flux of turbulent kinetic energy (TKE) prescribed by wave dissipation from a numerical wave model was compared with the conventional prescription using the wind friction velocity. A surface roughness length of 0.6 times the significant wave height was proposed, and the flux of TKE was applied at a distance below the mean sea surface that is half of this roughness length. The wave enhanced layer had a total depth that is almost three times the significant wave height. In this layer the non-dimensionalized Terray scaling with power of - 1.8 (instead of - 2) was applicable.

Transition Heat Transfer Modeling Based on the Characteristics of Turbulent Spots

NASA Technical Reports Server (NTRS)

Simon, Fred; Boyle, Robert

1998-01-01

While turbulence models are being developed which show promise for simulating the transition region on a turbine blade or vane, it is believed that the best approach with the greatest potential for practical use is the use of models which incorporate the physics of turbulent spots present in the transition region. This type of modeling results in the prediction of transition region intermittency which when incorporated in turbulence models give a good to excellent prediction of the transition region heat transfer. Some models are presented which show how turbulent spot characteristics and behavior can be employed to predict the effect of pressure gradient and Mach number on the transition region. The models predict the spot formation rate which is needed, in addition to the transition onset location, in the Narasimha concentrated breakdown intermittency equation. A simplified approach is taken for modeling turbulent spot growth and interaction in the transition region which utilizes the turbulent spot variables governing transition length and spot generation rate. The models are expressed in terms of spot spreading angle, dimensionless spot velocity, dimensionless spot area, disturbance frequency and Mach number. The models are used in conjunction with a computer code to predict the effects of pressure gradient and Mach number on the transition region and compared with VKI experimental turbine data.

Using computer algebra and SMT-solvers to analyze a mathematical model of cholera propagation

NASA Astrophysics Data System (ADS)

Trujillo Arredondo, Mariana

2014-06-01

We analyze a mathematical model for the transmission of cholera. The model is already defined and involves variables such as the pathogen agent, which in this case is the bacterium Vibrio cholera, and the human population. The human population is divided into three classes: susceptible, infectious and removed. Using Computer Algebra, specifically Maple we obtain two equilibrium states: the disease free state and the endemic state. Using Maple it is possible to prove that the disease free state is locally asymptotically stable if and only if R0 < 1. Using Maple it is possible to prove that the endemic equilibrium state is locally stable when it exists, it is to say when R0 > 1. Using the package Red-Log of the Computer algebra system Reduce and the SMT-Solver Z3Py it is possible to obtain numerical conditions for the model. The formula for the basic reproductive number makes a synthesis with all epidemic parameters in the model. Also it is possible to make numerical simulations which are very illustrative about the epidemic patters that are expected to be observed in real situations. We claim that these kinds of software are very useful in the analysis of epidemic models given that the symbolic computation provides algebraic formulas for the basic reproductive number and such algebraic formulas are very useful to derive control measures. For other side, computer algebra software is a powerful tool to make the stability analysis for epidemic models given that the all steps in the stability analysis can be made automatically: finding the equilibrium points, computing the jacobian, computing the characteristic polynomial for the jacobian, and applying the Routh-Hurwitz theorem to the characteristic polynomial. Finally, using SMT-Solvers is possible to make automatically checks of satisfiability, validity and quantifiers elimination being these computations very useful to analyse complicated epidemic models.

Phase relations in a forced turbulent boundary layer: implications for modelling of high Reynolds number wall turbulence.

PubMed

Duvvuri, Subrahmanyam; McKeon, Beverley

2017-03-13

Phase relations between specific scales in a turbulent boundary layer are studied here by highlighting the associated nonlinear scale interactions in the flow. This is achieved through an experimental technique that allows for targeted forcing of the flow through the use of a dynamic wall perturbation. Two distinct large-scale modes with well-defined spatial and temporal wavenumbers were simultaneously forced in the boundary layer, and the resulting nonlinear response from their direct interactions was isolated from the turbulence signal for the study. This approach advances the traditional studies of large- and small-scale interactions in wall turbulence by focusing on the direct interactions between scales with triadic wavenumber consistency. The results are discussed in the context of modelling high Reynolds number wall turbulence.This article is part of the themed issue 'Toward the development of high-fidelity models of wall turbulence at large Reynolds number'. © 2017 The Author(s).

Phase relations in a forced turbulent boundary layer: implications for modelling of high Reynolds number wall turbulence

PubMed Central

2017-01-01

Phase relations between specific scales in a turbulent boundary layer are studied here by highlighting the associated nonlinear scale interactions in the flow. This is achieved through an experimental technique that allows for targeted forcing of the flow through the use of a dynamic wall perturbation. Two distinct large-scale modes with well-defined spatial and temporal wavenumbers were simultaneously forced in the boundary layer, and the resulting nonlinear response from their direct interactions was isolated from the turbulence signal for the study. This approach advances the traditional studies of large- and small-scale interactions in wall turbulence by focusing on the direct interactions between scales with triadic wavenumber consistency. The results are discussed in the context of modelling high Reynolds number wall turbulence. This article is part of the themed issue ‘Toward the development of high-fidelity models of wall turbulence at large Reynolds number’. PMID:28167576

Turbulence Model Comparisons for a High-Speed Aircraft

NASA Technical Reports Server (NTRS)

Rivers, Melissa B.; Wahls, Richard A.

1999-01-01

Four turbulence models are described and evaluated for transonic flows over the High-Speed Research/industry baseline configuration known as Reference H by using the thin-layer, upwind, Navier-Stokes solver known as CFL3D. The turbulence models studied are the equilibrium model of Baldwin-Lomax (B-L) with the Degani-Schiff (D-S) modifications, the one-equation Baldwin-Barth (B-B) model, the one-equation Spalart-Allmaras (S-A) model, and Menter's two-equation Shear Stress Transport (SST) model. The flow conditions, which correspond to tests performed in the National Transonic Facility (NTF) at Langley Research Center, are a Mach number of 0.90 and a Reynolds number of 30 x 10 (exp. 6) based on mean aerodynamic chord for angles of attack of 1 deg., 5 deg., and 10 deg. The effects of grid topology and the representation of the actual wind tunnel model geometry are also investigated. Computed forces and surface pressures compare reasonably well with the experimental data for all four turbulence models.

Evaluation of Turbulence-Model Performance as Applied to Jet-Noise Prediction

NASA Technical Reports Server (NTRS)

Woodruff, S. L.; Seiner, J. M.; Hussaini, M. Y.; Erlebacher, G.

1998-01-01

The accurate prediction of jet noise is possible only if the jet flow field can be predicted accurately. Predictions for the mean velocity and turbulence quantities in the jet flowfield are typically the product of a Reynolds-averaged Navier-Stokes solver coupled with a turbulence model. To evaluate the effectiveness of solvers and turbulence models in predicting those quantities most important to jet noise prediction, two CFD codes and several turbulence models were applied to a jet configuration over a range of jet temperatures for which experimental data is available.

Geometric Model of Topological Insulators from the Maxwell Algebra

NASA Astrophysics Data System (ADS)

Palumbo, Giandomenico

I propose a novel geometric model of time-reversal-invariant topological insulators in three dimensions in presence of an external electromagnetic field. Their gapped boundary supports relativistic quantum Hall states and is described by a Chern-Simons theory, where the gauge connection takes values in the Maxwell algebra. This represents a non-central extension of the Poincare' algebra and takes into account both the Lorentz and magnetic-translation symmetries of the surface states. In this way, I derive a relativistic version of the Wen-Zee term and I show that the non-minimal coupling between the background geometry and the electromagnetic field in the model is in agreement with the main properties of the relativistic quantum Hall states in the flat space. This work is part of the DITP consortium, a program of the Netherlands Organisation for Scientific Research (NWO) that is funded by the Dutch Ministry of Education, Culture and Science (OCW).

Entanglement in a model for Hawking radiation: An application of quadratic algebras

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

Bambah, Bindu A., E-mail: bbsp@uohyd.ernet.in; Mukku, C., E-mail: mukku@iiit.ac.in; Shreecharan, T., E-mail: shreecharan@gmail.com

2013-03-15

Quadratic polynomially deformed su(1,1) and su(2) algebras are utilized in model Hamiltonians to show how the gravitational system consisting of a black hole, infalling radiation and outgoing (Hawking) radiation can be solved exactly. The models allow us to study the long-time behaviour of the black hole and its outgoing modes. In particular, we calculate the bipartite entanglement entropies of subsystems consisting of (a) infalling plus outgoing modes and (b) black hole modes plus the infalling modes, using the Janus-faced nature of the model. The long-time behaviour also gives us glimpses of modifications in the character of Hawking radiation. Finally, wemore » study the phenomenon of superradiance in our model in analogy with atomic Dicke superradiance. - Highlights: Black-Right-Pointing-Pointer We examine a toy model for Hawking radiation with quantized black hole modes. Black-Right-Pointing-Pointer We use quadratic polynomially deformed su(1,1) algebras to study its entanglement properties. Black-Right-Pointing-Pointer We study the 'Dicke Superradiance' in black hole radiation using quadratically deformed su(2) algebras. Black-Right-Pointing-Pointer We study the modification of the thermal character of Hawking radiation due to quantized black hole modes.« less

On the Yakhot-Orszag renormalization group method for deriving turbulence statistics and models

NASA Technical Reports Server (NTRS)

Smith, L. M.; Reynolds, W. C.

1992-01-01

An independent, comprehensive, critical review of the 'renormalization group' (RNG) theory of turbulence developed by Yakhot and Orszag (1986) is provided. Their basic theory for the Navier-Stokes equations is confirmed, and approximations in the scale removal procedure are discussed. The YO derivations of the velocity-derivative skewness and the transport equation for the energy dissipation rate are examined. An algebraic error in the derivation of the skewness is corrected. The corrected RNG skewness value of -0.59 is in agreement with experiments at moderate Reynolds numbers. Several problems are identified in the derivation of the energy dissipation rate equations which suggest that the derivation should be reformulated.

Algebraic Structure of tt * Equations for Calabi-Yau Sigma Models

NASA Astrophysics Data System (ADS)

Alim, Murad

2017-08-01

The tt * equations define a flat connection on the moduli spaces of {2d, \\mathcal{N}=2} quantum field theories. For conformal theories with c = 3 d, which can be realized as nonlinear sigma models into Calabi-Yau d-folds, this flat connection is equivalent to special geometry for threefolds and to its analogs in other dimensions. We show that the non-holomorphic content of the tt * equations, restricted to the conformal directions, in the cases d = 1, 2, 3 is captured in terms of finitely many generators of special functions, which close under derivatives. The generators are understood as coordinates on a larger moduli space. This space parameterizes a freedom in choosing representatives of the chiral ring while preserving a constant topological metric. Geometrically, the freedom corresponds to a choice of forms on the target space respecting the Hodge filtration and having a constant pairing. Linear combinations of vector fields on that space are identified with the generators of a Lie algebra. This Lie algebra replaces the non-holomorphic derivatives of tt * and provides these with a finer and algebraic meaning. For sigma models into lattice polarized K3 manifolds, the differential ring of special functions on the moduli space is constructed, extending known structures for d = 1 and 3. The generators of the differential rings of special functions are given by quasi-modular forms for d = 1 and their generalizations in d = 2, 3. Some explicit examples are worked out including the case of the mirror of the quartic in {\\mathbbm{P}^3}, where due to further algebraic constraints, the differential ring coincides with quasi modular forms.

Constructing and Modeling Algebraic Statements in the Multiplicative Domain: Investigating Fourth-Grade Student and Teacher Learning

ERIC Educational Resources Information Center

Grandau, Laura

2013-01-01

This study of fourth-grade students and teachers explores mathematics teaching and learning that focuses on discovering and modeling algebraic relationships. The study has two parts: an investigation of how students learn to construct algebraic statements and models for comparisons and measurement situations in the multiplicative domain, and an…

Implementation of Dryden Continuous Turbulence Model into Simulink for LSA-02 Flight Test Simulation

NASA Astrophysics Data System (ADS)

Ichwanul Hakim, Teuku Mohd; Arifianto, Ony

2018-04-01

Turbulence is a movement of air on small scale in the atmosphere that caused by instabilities of pressure and temperature distribution. Turbulence model is integrated into flight mechanical model as an atmospheric disturbance. Common turbulence model used in flight mechanical model are Dryden and Von Karman model. In this minor research, only Dryden continuous turbulence model were made. Dryden continuous turbulence model has been implemented, it refers to the military specification MIL-HDBK-1797. The model was implemented into Matlab Simulink. The model will be integrated with flight mechanical model to observe response of the aircraft when it is flight through turbulence field. The turbulence model is characterized by multiplying the filter which are generated from power spectral density with band-limited Gaussian white noise input. In order to ensure that the model provide a good result, model verification has been done by comparing the implemented model with the similar model that is provided in aerospace blockset. The result shows that there are some difference for 2 linear velocities (vg and wg), and 3 angular rate (pg, qg and rg). The difference is instantly caused by different determination of turbulence scale length which is used in aerospace blockset. With the adjustment of turbulence length in the implemented model, both model result the similar output.

Near-wall modeling of compressible turbulent flow

NASA Technical Reports Server (NTRS)

So, Ronald M. C.

1991-01-01

A near-wall two-equation model for compressible flows is proposed. The model is formulated by relaxing the assumption of dynamic field similarity between compressible and incompressible flows. A postulate is made to justify the extension of incompressible models to ammount for compressibility effects. This requires formulation the turbulent kinetic energy equation in a form similar to its incompressible counterpart. As a result, the compressible dissipation function has to be split into a solenoidal part, which is not sensitive to changes of compressibility indicators, and a dilatational part, which is directly affected by these changes. A model with an explicit dependence on the turbulent Mach number is proposed for the dilatational dissipation rate.

On the Connection Between One-and Two-Equation Models of Turbulence

NASA Technical Reports Server (NTRS)

Menter, F. R.; Rai, Man Mohan (Technical Monitor)

1994-01-01

A formalism will be presented that allows the transformation of two-equation eddy viscosity turbulence models into one-equation models. The transformation is based on an assumption that is widely accepted over a large range of boundary layer flows and that has been shown to actually improve predictions when incorporated into two-equation models of turbulence. Based on that assumption, a new one-equation turbulence model will be derived. The new model will be tested in great detail against a previously introduced one-equation model and against its parent two-equation model.

Multifractal Modeling of Turbulent Mixing

NASA Astrophysics Data System (ADS)

Samiee, Mehdi; Zayernouri, Mohsen; Meerschaert, Mark M.

2017-11-01

Stochastic processes in random media are emerging as interesting tools for modeling anomalous transport phenomena. Applications include intermittent passive scalar transport with background noise in turbulent flows, which are observed in atmospheric boundary layers, turbulent mixing in reactive flows, and long-range dependent flow fields in disordered/fractal environments. In this work, we propose a nonlocal scalar transport equation involving the fractional Laplacian, where the corresponding fractional index is linked to the multifractal structure of the nonlinear passive scalar power spectrum. This work was supported by the AFOSR Young Investigator Program (YIP) award (FA9550-17-1-0150) and partially by MURI/ARO (W911NF-15-1-0562).

A dynamical model of plasma turbulence in the solar wind

PubMed Central

Howes, G. G.

2015-01-01

A dynamical approach, rather than the usual statistical approach, is taken to explore the physical mechanisms underlying the nonlinear transfer of energy, the damping of the turbulent fluctuations, and the development of coherent structures in kinetic plasma turbulence. It is argued that the linear and nonlinear dynamics of Alfvén waves are responsible, at a very fundamental level, for some of the key qualitative features of plasma turbulence that distinguish it from hydrodynamic turbulence, including the anisotropic cascade of energy and the development of current sheets at small scales. The first dynamical model of kinetic turbulence in the weakly collisional solar wind plasma that combines self-consistently the physics of Alfvén waves with the development of small-scale current sheets is presented and its physical implications are discussed. This model leads to a simplified perspective on the nature of turbulence in a weakly collisional plasma: the nonlinear interactions responsible for the turbulent cascade of energy and the formation of current sheets are essentially fluid in nature, while the collisionless damping of the turbulent fluctuations and the energy injection by kinetic instabilities are essentially kinetic in nature. PMID:25848075

VNAP2: A Computer Program for Computation of Two-dimensional, Time-dependent, Compressible, Turbulent Flow

NASA Technical Reports Server (NTRS)

Cline, M. C.

1981-01-01

A computer program, VNAP2, for calculating turbulent (as well as laminar and inviscid), steady, and unsteady flow is presented. It solves the two dimensional, time dependent, compressible Navier-Stokes equations. The turbulence is modeled with either an algebraic mixing length model, a one equation model, or the Jones-Launder two equation model. The geometry may be a single or a dual flowing stream. The interior grid points are computed using the unsplit MacCormack scheme. Two options to speed up the calculations for high Reynolds number flows are included. The boundary grid points are computed using a reference plane characteristic scheme with the viscous terms treated as source functions. An explicit artificial viscosity is included for shock computations. The fluid is assumed to be a perfect gas. The flow boundaries may be arbitrary curved solid walls, inflow/outflow boundaries, or free jet envelopes. Typical problems that can be solved concern nozzles, inlets, jet powered afterbodies, airfoils, and free jet expansions. The accuracy and efficiency of the program are shown by calculations of several inviscid and turbulent flows. The program and its use are described completely, and six sample cases and a code listing are included.

Stochastic modeling of turbulent reacting flows

NASA Technical Reports Server (NTRS)

Fox, R. O.; Hill, J. C.; Gao, F.; Moser, R. D.; Rogers, M. M.

1992-01-01

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

The Dissipation Rate Transport Equation and Subgrid-Scale Models in Rotating Turbulence

NASA Technical Reports Server (NTRS)

Rubinstein, Robert; Ye, Zhou

1997-01-01

The dissipation rate transport equation remains the most uncertain part of turbulence modeling. The difficulties arc increased when external agencies like rotation prevent straightforward dimensional analysis from determining the correct form of the modelled equation. In this work, the dissipation rate transport equation and subgrid scale models for rotating turbulence are derived from an analytical statistical theory of rotating turbulence. In the strong rotation limit, the theory predicts a turbulent steady state in which the inertial range energy spectrum scales as k(sup -2) and the turbulent time scale is the inverse rotation rate. This scaling has been derived previously by heuristic arguments.

A Study of Two-Equation Turbulence Models on the Elliptic Streamline Flow

NASA Technical Reports Server (NTRS)

Blaisdell, Gregory A.; Qin, Jim H.; Shariff, Karim; Rai, Man Mohan (Technical Monitor)

1995-01-01

Several two-equation turbulence models are compared to data from direct numerical simulations (DNS) of the homogeneous elliptic streamline flow, which combines rotation and strain. The models considered include standard two-equation models and models with corrections for rotational effects. Most of the rotational corrections modify the dissipation rate equation to account for the reduced dissipation rate in rotating turbulent flows, however, the DNS data shows that the production term in the turbulent kinetic energy equation is not modeled correctly by these models. Nonlinear relations for the Reynolds stresses are considered as a means of modifying the production term. Implications for the modeling of turbulent vortices will be discussed.

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.

Designing Tasks for Math Modeling in College Algebra: A Critical Review

ERIC Educational Resources Information Center

Staats, Susan; Robertson, Douglas

2014-01-01

Over the last decade, the pedagogical approach known as mathematical modeling has received increased interest in college algebra classes in the United States. Math modeling assignments ask students to develop their own problem-solving tools to address non-routine, realistic scenarios. The open-ended quality of modeling activities creates dilemmas…

Modeling Scramjet Flows with Variable Turbulent Prandtl and Schmidt Numbers

NASA Technical Reports Server (NTRS)

Xiao, X.; Hassan, H. A.; Baurle, R. A.

2006-01-01

A complete turbulence model, where the turbulent Prandtl and Schmidt numbers are calculated as part of the solution and where averages involving chemical source terms are modeled, is presented. The ability of avoiding the use of assumed or evolution Probability Distribution Functions (PDF's) results in a highly efficient algorithm for reacting flows. The predictions of the model are compared with two sets of experiments involving supersonic mixing and one involving supersonic combustion. The results demonstrate the need for consideration of turbulence/chemistry interactions in supersonic combustion. In general, good agreement with experiment is indicated.

Numerical methods in laminar and turbulent flow; Proceedings of the 7th International Conference, Stanford Univ., CA, July 15-19, 1991. Vol. 7, pts. 1 & 2

NASA Technical Reports Server (NTRS)

Taylor, C. (Editor); Chin, J. H. (Editor); Homsy, G. M. (Editor)

1991-01-01

Consideration is given to the impulse response of a laminar boundary layer and receptivity; numerical transition to turbulence in plane Poiseuille flow; large eddy simulation of turbulent wake flow; a viscous model and loss calculation of a multisplitter cascade; vortex initiation during dynamic stall of an airfoil; a numerical analysis of isothermal flow in a combustion chamber; and compressible flow calculations with a two-equation turbulence model and unstructured grids. Attention is also given to a 2D calculation of a buoyant flow around a burning sphere, a fast multigrid method for 3D turbulent incompressible flows, a streaming flow induced by an oscillating cascade of circular cylinders, an algebraic multigrid scheme for solving the Navier-Stokes equations on unstructured meshes; and nonlinear coupled multigrid solutions to thermal problems employing different nodal grid arrangements and convective transport approximations.

Quantum cluster algebras and quantum nilpotent algebras.

PubMed

Goodearl, Kenneth R; Yakimov, Milen T

2014-07-08

A major direction in the theory of cluster algebras is to construct (quantum) cluster algebra structures on the (quantized) coordinate rings of various families of varieties arising in Lie theory. We prove that all algebras in a very large axiomatically defined class of noncommutative algebras possess canonical quantum cluster algebra structures. Furthermore, they coincide with the corresponding upper quantum cluster algebras. We also establish analogs of these results for a large class of Poisson nilpotent algebras. Many important families of coordinate rings are subsumed in the class we are covering, which leads to a broad range of applications of the general results to the above-mentioned types of problems. As a consequence, we prove the Berenstein-Zelevinsky conjecture [Berenstein A, Zelevinsky A (2005) Adv Math 195:405-455] for the quantized coordinate rings of double Bruhat cells and construct quantum cluster algebra structures on all quantum unipotent groups, extending the theorem of Geiß et al. [Geiß C, et al. (2013) Selecta Math 19:337-397] for the case of symmetric Kac-Moody groups. Moreover, we prove that the upper cluster algebras of Berenstein et al. [Berenstein A, et al. (2005) Duke Math J 126:1-52] associated with double Bruhat cells coincide with the corresponding cluster algebras.

Quantum cluster algebras and quantum nilpotent algebras

PubMed Central

Goodearl, Kenneth R.; Yakimov, Milen T.

2014-01-01

A major direction in the theory of cluster algebras is to construct (quantum) cluster algebra structures on the (quantized) coordinate rings of various families of varieties arising in Lie theory. We prove that all algebras in a very large axiomatically defined class of noncommutative algebras possess canonical quantum cluster algebra structures. Furthermore, they coincide with the corresponding upper quantum cluster algebras. We also establish analogs of these results for a large class of Poisson nilpotent algebras. Many important families of coordinate rings are subsumed in the class we are covering, which leads to a broad range of applications of the general results to the above-mentioned types of problems. As a consequence, we prove the Berenstein–Zelevinsky conjecture [Berenstein A, Zelevinsky A (2005) Adv Math 195:405–455] for the quantized coordinate rings of double Bruhat cells and construct quantum cluster algebra structures on all quantum unipotent groups, extending the theorem of Geiß et al. [Geiß C, et al. (2013) Selecta Math 19:337–397] for the case of symmetric Kac–Moody groups. Moreover, we prove that the upper cluster algebras of Berenstein et al. [Berenstein A, et al. (2005) Duke Math J 126:1–52] associated with double Bruhat cells coincide with the corresponding cluster algebras. PMID:24982197

Mathematical Modelling in Engineering: A Proposal to Introduce Linear Algebra Concepts

ERIC Educational Resources Information Center

Cárcamo Bahamonde, Andrea; Gómez Urgelles, Joan; Fortuny Aymemí, Josep

2016-01-01

The modern dynamic world requires that basic science courses for engineering, including linear algebra, emphasise the development of mathematical abilities primarily associated with modelling and interpreting, which are not exclusively calculus abilities. Considering this, an instructional design was created based on mathematical modelling and…

Characterizing the Severe Turbulence Environments Associated With Commercial Aviation Accidents: A Real-Time Turbulence Model (RTTM) Designed for the Operational Prediction of Hazardous Aviation Turbulence Environments

NASA Technical Reports Server (NTRS)

Kaplan, Michael L.; Lux, Kevin M.; Cetola, Jeffrey D.; Huffman, Allan W.; Riordan, Allen J.; Slusser, Sarah W.; Lin, Yuh-Lang; Charney, Joseph J.; Waight, Kenneth T.

2004-01-01

Real-time prediction of environments predisposed to producing moderate-severe aviation turbulence is studied. We describe the numerical model and its postprocessing system designed for said prediction of environments predisposed to severe aviation turbulence as well as presenting numerous examples of its utility. The numerical model is MASS version 5.13, which is integrated over three different grid matrices in real time on a university work station in support of NASA Langley Research Center s B-757 turbulence research flight missions. The postprocessing system includes several turbulence-related products, including four turbulence forecasting indices, winds, streamlines, turbulence kinetic energy, and Richardson numbers. Additionally, there are convective products including precipitation, cloud height, cloud mass fluxes, lifted index, and K-index. Furthermore, soundings, sounding parameters, and Froude number plots are also provided. The horizontal cross-section plot products are provided from 16 000 to 46 000 ft in 2000-ft intervals. Products are available every 3 hours at the 60- and 30-km grid interval and every 1.5 hours at the 15-km grid interval. The model is initialized from the NWS ETA analyses and integrated two times a day.

Intelligently deciphering unintelligible designs: algorithmic algebraic model checking in systems biology

PubMed Central

Mishra, Bud

2009-01-01

Systems biology, as a subject, has captured the imagination of both biologists and systems scientists alike. But what is it? This review provides one researcher's somewhat idiosyncratic view of the subject, but also aims to persuade young scientists to examine the possible evolution of this subject in a rich historical context. In particular, one may wish to read this review to envision a subject built out of a consilience of many interesting concepts from systems sciences, logic and model theory, and algebra, culminating in novel tools, techniques and theories that can reveal deep principles in biology—seen beyond mere observations. A particular focus in this review is on approaches embedded in an embryonic program, dubbed ‘algorithmic algebraic model checking’, and its powers and limitations. PMID:19364723

Turbulent flow separation in three-dimensional asymmetric diffusers

NASA Astrophysics Data System (ADS)

Jeyapaul, Elbert

2011-12-01

separation is due to effects of lateral straining, streamline curvature, secondary flow of second kind, transverse pressure gradient on turbulence. Resolving these effects is possible with anisotropy turbulence models as the Explicit Algebraic Reynolds stress model (EARSM). This model has provided accurate prediction of streamwise and transverse velocity, however the wall pressure is under predicted. An improved EARSM model is developed by correcting the coefficients, which predicts a more accurate wall pressure. There exists scope for improvement of this model, by including convective effects and dynamics of velocity gradient invariants.

Computations of Flow over a Hump Model Using Higher Order Method with Turbulence Modeling

NASA Technical Reports Server (NTRS)

Balakumar, P.

2005-01-01

Turbulent separated flow over a two-dimensional hump is computed by solving the RANS equations with k - omega (SST) turbulence model for the baseline, steady suction and oscillatory blowing/suction flow control cases. The flow equations and the turbulent model equations are solved using a fifth-order accurate weighted essentially. nonoscillatory (WENO) scheme for space discretization and a third order, total variation diminishing (TVD) Runge-Kutta scheme for time integration. Qualitatively the computed pressure distributions exhibit the same behavior as those observed in the experiments. The computed separation regions are much longer than those observed experimentally. However, the percentage reduction in the separation region in the steady suction case is closer to what was measured in the experiment. The computations did not predict the expected reduction in the separation length in the oscillatory case. The predicted turbulent quantities are two to three times smaller than the measured values pointing towards the deficiencies in the existing turbulent models when they are applied to strong steady/unsteady separated flows.

Numerical Modeling of Turbulent Combustion

NASA Technical Reports Server (NTRS)

Ghoneim, A. F.; Chorin, A. J.; Oppenheim, A. K.

1983-01-01

The work in numerical modeling is focused on the use of the random vortex method to treat turbulent flow fields associated with combustion while flame fronts are considered as interfaces between reactants and products, propagating with the flow and at the same time advancing in the direction normal to themselves at a prescribed burning speed. The latter is associated with the generation of specific volume (the flame front acting, in effect, as the locus of volumetric sources) to account for the expansion of the flow field due to the exothermicity of the combustion process. The model was applied to the flow in a channel equipped with a rearward facing step. The results obtained revealed the mechanism of the formation of large scale turbulent structure in the wake of the step, while it showed the flame to stabilize on the outer edges of these eddies.

A spatial operator algebra for manipulator modeling and control

NASA Technical Reports Server (NTRS)

Rodriguez, G.; Kreutz, K.; Jain, A.

1989-01-01

A spatial operator algebra for modeling the control and trajectory design of manipulation is discussed, with emphasis on its analytical formulation and implementation in the Ada programming language. The elements of this algebra are linear operators whose domain and range spaces consist of forces, moments, velocities, and accelerations. The effect of these operators is equivalent to a spatial recursion along the span of the manipulator. Inversion is obtained using techniques of recursive filtering and smoothing. The operator alegbra provides a high-level framework for describing the dynamic and kinematic behavior of a manipulator and control and trajectory design algorithms. Implementable recursive algorithms can be immediately derived from the abstract operator expressions by inspection, thus greatly simplifying the transition from an abstract problem formulation and solution to the detailed mechanization of a specific algorithm.

High Reynolds number turbulence model of rotating shear flows

NASA Astrophysics Data System (ADS)

Masuda, S.; Ariga, I.; Koyama, H. S.

1983-09-01

A Reynolds stress closure model for rotating turbulent shear flows is developed. Special attention is paid to keeping the model constants independent of rotation. First, general forms of the model of a Reynolds stress equation and a dissipation rate equation are derived, the only restrictions of which are high Reynolds number and incompressibility. The model equations are then applied to two-dimensional equilibrium boundary layers and the effects of Coriolis acceleration on turbulence structures are discussed. Comparisons with the experimental data and with previous results in other external force fields show that there exists a very close analogy between centrifugal, buoyancy and Coriolis force fields. Finally, the model is applied to predict the two-dimensional boundary layers on rotating plane walls. Comparisons with existing data confirmed its capability of predicting mean and turbulent quantities without employing any empirical relations in rotating fields.

Spatial-Operator Algebra For Robotic Manipulators

NASA Technical Reports Server (NTRS)

Rodriguez, Guillermo; Kreutz, Kenneth K.; Milman, Mark H.

1991-01-01

Report discusses spatial-operator algebra developed in recent studies of mathematical modeling, control, and design of trajectories of robotic manipulators. Provides succinct representation of mathematically complicated interactions among multiple joints and links of manipulator, thereby relieving analyst of most of tedium of detailed algebraic manipulations. Presents analytical formulation of spatial-operator algebra, describes some specific applications, summarizes current research, and discusses implementation of spatial-operator algebra in the Ada programming language.

Boundary-layer computational model for predicting the flow and heat transfer in sudden expansions

NASA Technical Reports Server (NTRS)

Lewis, J. P.; Pletcher, R. H.

1986-01-01

Fully developed turbulent and laminar flows through symmetric planar and axisymmetric expansions with heat transfer were modeled using a finite-difference discretization of the boundary-layer equations. By using the boundary-layer equations to model separated flow in place of the Navier-Stokes equations, computational effort was reduced permitting turbulence modelling studies to be economically carried out. For laminar flow, the reattachment length was well predicted for Reynolds numbers as low as 20 and the details of the trapped eddy were well predicted for Reynolds numbers above 200. For turbulent flows, the Boussinesq assumption was used to express the Reynolds stresses in terms of a turbulent viscosity. Near-wall algebraic turbulence models based on Prandtl's-mixing-length model and the maximum Reynolds shear stress were compared.

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 new energy transfer model for turbulent free shear flow

NASA Technical Reports Server (NTRS)

Liou, William W.-W.

1992-01-01

A new model for the energy transfer mechanism in the large-scale turbulent kinetic energy equation is proposed. An estimate of the characteristic length scale of the energy containing large structures is obtained from the wavelength associated with the structures predicted by a weakly nonlinear analysis for turbulent free shear flows. With the inclusion of the proposed energy transfer model, the weakly nonlinear wave models for the turbulent large-scale structures are self-contained and are likely to be independent flow geometries. The model is tested against a plane mixing layer. Reasonably good agreement is achieved. Finally, it is shown by using the Liapunov function method, the balance between the production and the drainage of the kinetic energy of the turbulent large-scale structures is asymptotically stable as their amplitude saturates. The saturation of the wave amplitude provides an alternative indicator for flow self-similarity.

Comparative Study on the Prediction of Aerodynamic Characteristics of Aircraft with Turbulence Models

NASA Astrophysics Data System (ADS)

Jang, Yujin; Huh, Jinbum; Lee, Namhun; Lee, Seungsoo; Park, Youngmin

2018-04-01

The RANS equations are widely used to analyze complex flows over aircraft. The equations require a turbulence model for turbulent flow analyses. A suitable turbulence must be selected for accurate predictions of aircraft aerodynamic characteristics. In this study, numerical analyses of three-dimensional aircraft are performed to compare the results of various turbulence models for the prediction of aircraft aerodynamic characteristics. A 3-D RANS solver, MSAPv, is used for the aerodynamic analysis. The four turbulence models compared are the Sparlart-Allmaras (SA) model, Coakley's q-ω model, Huang and Coakley's k-ɛ model, and Menter's k-ω SST model. Four aircrafts are considered: an ARA-M100, DLR-F6 wing-body, DLR-F6 wing-body-nacelle-pylon from the second drag prediction workshop, and a high wing aircraft with nacelles. The CFD results are compared with experimental data and other published computational results. The details of separation patterns, shock positions, and Cp distributions are discussed to find the characteristics of the turbulence models.

Low-Level Turbulence Forecasts From Fine-Scale Models

DTIC Science & Technology

2014-02-01

aircraft which often fly above 6097 m for much of the flight. Sharman et al. (35) have developed the Graphical Turbulence Guidance ( GTG ), which is a...original GTG was used above 6097-m AGL and is for MOD or greater CAT. The Rapid Refresh Model, which is an hourly updated assimilation model operational...applications and tactical decision aids. The recent version of GTG (V2.5) forecasts turbulence as low as 3354 m. Silberberg (36) notes that the Aviation

Continual Lie algebras and noncommutative counterparts of exactly solvable models

NASA Astrophysics Data System (ADS)

Zuevsky, A.

2004-01-01

Noncommutative counterparts of exactly solvable models are introduced on the basis of a generalization of Saveliev-Vershik continual Lie algebras. Examples of noncommutative Liouville and sin/h-Gordon equations are given. The simplest soliton solution to the noncommutative sine-Gordon equation is found.

A near-wall two-equation model for compressible turbulent flows

NASA Technical Reports Server (NTRS)

Zhang, H. S.; So, R. M. C.; Speziale, C. G.; Lai, Y. G.

1991-01-01

A near-wall two-equation turbulence model of the K - epsilon type is developed for the description of high-speed compressible flows. The Favre-averaged equations of motion are solved in conjunction with modeled transport equations for the turbulent kinetic energy and solenoidal dissipation wherein a variable density extension of the asymptotically consistent near-wall model of So and co-workers is supplemented with new dilatational models. The resulting compressible two-equation model is tested in the supersonic flat plate boundary layer - with an adiabatic wall and with wall cooling - for Mach numbers as large as 10. Direct comparisons of the predictions of the new model with raw experimental data and with results from the K - omega model indicate that it performs well for a wide range of Mach numbers. The surprising finding is that the Morkovin hypothesis, where turbulent dilatational terms are neglected, works well at high Mach numbers, provided that the near wall model is asymptotically consistent. Instances where the model predictions deviate from the experiments appear to be attributable to the assumption of constant turbulent Prandtl number - a deficiency that will be addressed in a future paper.

Yang-Baxter algebras, integrable theories and Bethe Ansatz

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

De Vega, H.J.

1990-03-10

This paper presents the Yang-Baxter algebras (YBA) in a general framework stressing their power to exactly solve the lattice models associated to them. The algebraic Behe Ansatz is developed as an eigenvector construction based on the YBA. The six-vertex model solution is given explicitly. The generalization of YB algebras to face language is considered. The algebraic BA for the SOS model of Andrews, Baxter and Forrester is described using these face YB algebras. It is explained how these lattice models yield both solvable massive QFT and conformal models in appropriated scaling (continuous) limits within the lattice light-cone approach. This approachmore » permit to define and solve rigorously massive QFT as an appropriate continuum limit of gapless vertex models. The deep links between the YBA and Lie algebras are analyzed including the quantum groups that underlay the trigonometric/hyperbolic YBA. Braid and quantum groups are derived from trigonometric/hyperbolic YBA in the limit of infinite spectral parameter. To conclude, some recent developments in the domain of integrable theories are summarized.« less

Turbulence Model Selection for Low Reynolds Number Flows

PubMed Central

2016-01-01

One of the major flow phenomena associated with low Reynolds number flow is the formation of separation bubbles on an airfoil’s surface. NACA4415 airfoil is commonly used in wind turbines and UAV applications. The stall characteristics are gradual compared to thin airfoils. The primary criterion set for this work is the capture of laminar separation bubble. Flow is simulated for a Reynolds number of 120,000. The numerical analysis carried out shows the advantages and disadvantages of a few turbulence models. The turbulence models tested were: one equation Spallart Allmars (S-A), two equation SST K-ω, three equation Intermittency (γ) SST, k-kl-ω and finally, the four equation transition γ-Reθ SST. However, the variation in flow physics differs between these turbulence models. Procedure to establish the accuracy of the simulation, in accord with previous experimental results, has been discussed in detail. PMID:27104354

Turbulence Model Selection for Low Reynolds Number Flows.

PubMed

Aftab, S M A; Mohd Rafie, A S; Razak, N A; Ahmad, K A

2016-01-01

One of the major flow phenomena associated with low Reynolds number flow is the formation of separation bubbles on an airfoil's surface. NACA4415 airfoil is commonly used in wind turbines and UAV applications. The stall characteristics are gradual compared to thin airfoils. The primary criterion set for this work is the capture of laminar separation bubble. Flow is simulated for a Reynolds number of 120,000. The numerical analysis carried out shows the advantages and disadvantages of a few turbulence models. The turbulence models tested were: one equation Spallart Allmars (S-A), two equation SST K-ω, three equation Intermittency (γ) SST, k-kl-ω and finally, the four equation transition γ-Reθ SST. However, the variation in flow physics differs between these turbulence models. Procedure to establish the accuracy of the simulation, in accord with previous experimental results, has been discussed in detail.

New Models for Velocity/Pressure-Gradient Correlations in Turbulent Boundary Layers

NASA Astrophysics Data System (ADS)

Poroseva, Svetlana; Murman, Scott

2014-11-01

To improve the performance of Reynolds-Averaged Navier-Stokes (RANS) turbulence models, one has to improve the accuracy of models for three physical processes: turbulent diffusion, interaction of turbulent pressure and velocity fluctuation fields, and dissipative processes. The accuracy of modeling the turbulent diffusion depends on the order of a statistical closure chosen as a basis for a RANS model. When the Gram-Charlier series expansions for the velocity correlations are used to close the set of RANS equations, no assumption on Gaussian turbulence is invoked and no unknown model coefficients are introduced into the modeled equations. In such a way, this closure procedure reduces the modeling uncertainty of fourth-order RANS (FORANS) closures. Experimental and direct numerical simulation data confirmed the validity of using the Gram-Charlier series expansions in various flows including boundary layers. We will address modeling the velocity/pressure-gradient correlations. New linear models will be introduced for the second- and higher-order correlations applicable to two-dimensional incompressible wall-bounded flows. Results of models' validation with DNS data in a channel flow and in a zero-pressure gradient boundary layer over a flat plate will be demonstrated. A part of the material is based upon work supported by NASA under award NNX12AJ61A.

Improved turbulence models based on large eddy simulation of homogeneous, incompressible turbulent flows

NASA Technical Reports Server (NTRS)

Bardino, J.; Ferziger, J. H.; Reynolds, W. C.

1983-01-01

The physical bases of large eddy simulation and subgrid modeling are studied. A subgrid scale similarity model is developed that can account for system rotation. Large eddy simulations of homogeneous shear flows with system rotation were carried out. Apparently contradictory experimental results were explained. The main effect of rotation is to increase the transverse length scales in the rotation direction, and thereby decrease the rates of dissipation. Experimental results are shown to be affected by conditions at the turbulence producing grid, which make the initial states a function of the rotation rate. A two equation model is proposed that accounts for effects of rotation and shows good agreement with experimental results. In addition, a Reynolds stress model is developed that represents the turbulence structure of homogeneous shear flows very well and can account also for the effects of system rotation.

Turbulence model development and application at Lockheed Fort Worth Company

NASA Technical Reports Server (NTRS)

Smith, Brian R.

1995-01-01

This viewgraph presentation demonstrates that computationally efficient k-l and k-kl turbulence models have been developed and implemented at Lockheed Fort Worth Company. Many years of experience have been gained applying two equation turbulence models to complex three-dimensional flows for design and analysis.

Status, Emerging Ideas and Future Directions of Turbulence Modeling Research in Aeronautics

NASA Technical Reports Server (NTRS)

Duraisamy, Karthik; Spalart, Philippe R.; Rumsey, Christopher L.

2017-01-01

In July 2017, a three-day Turbulence Modeling Symposium sponsored by the University of Michigan and NASA was held in Ann Arbor, Michigan. This meeting brought together nearly 90 experts from academia, government and industry, with good international participation, to discuss the state of the art in turbulence modeling, emerging ideas, and to wrestle with questions surrounding its future. Emphasis was placed on turbulence modeling in a predictive context in complex problems, rather than on turbulence theory or descriptive modeling. This report summarizes many of the questions, discussions, and conclusions from the symposium, and suggests immediate next steps.

Turbulent flow in a 180 deg bend: Modeling and computations

NASA Technical Reports Server (NTRS)

Kaul, Upender K.

1989-01-01

A low Reynolds number k-epsilon turbulence model was presented which yields accurate predictions of the kinetic energy near the wall. The model is validated with the experimental channel flow data of Kreplin and Eckelmann. The predictions are also compared with earlier results from direct simulation of turbulent channel flow. The model is especially useful for internal flows where the inflow boundary condition of epsilon is not easily prescribed. The model partly derives from some observations based on earlier direct simulation results of near-wall turbulence. The low Reynolds number turbulence model together with an existing curvature correction appropriate to spinning cylinder flows was used to simulate the flow in a U-bend with the same radius of curvature as the Space Shuttle Main Engine (SSME) Turn-Around Duct (TAD). The present computations indicate a space varying curvature correction parameter as opposed to a constant parameter as used in the spinning cylinder flows. Comparison with limited available experimental data is made. The comparison is favorable, but detailed experimental data is needed to further improve the curvature model.

Models of inertial range spectra of interplanetary magnetohydrodynamic turbulence

NASA Technical Reports Server (NTRS)

Zhou, YE; Matthaeus, William H.

1990-01-01

A framework based on turbulence theory is presented to develop approximations for the local turbulence effects that are required in transport models. An approach based on Kolmogoroff-style dimensional analysis is presented as well as one based on a wave-number diffusion picture. Particular attention is given to the case of MHD turbulence with arbitrary cross helicity and with arbitrary ratios of the Alfven time scale and the nonlinear time scale.

A new time scale based k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)