Results for “Closure problem” — Page 4

NASA NTRS1994

LES, DNS and RANS for the analysis of high-speed turbulent reacting flows

The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES), direct numerical simulation (DNS), and Reynolds averaged Navier Stokes (RANS) methods for the computational analysis of high-speed reacting turbulent flows. In the second phase of this work, covering the period 1 Sep. 1993 - 1 Sep. 1994, we have focused our efforts on two research problems: (1) developments of 'algebraic' moment closures for statistical descriptions of nonpremixed reacting systems, and (2) assessments of the Dirichlet frequency in presumed scalar probability density function (PDF) methods in stochastic description of turbulent reacting flows. This report provides a complete description of our efforts during this past year as supported by the NASA Langley Research Center under Grant NAG1-1122.

Givi, Peyman↗

NASA NTRS1998

Transport Coefficients in Rotating Weakly Compressible Turbulence

Analytical studies of compressible turbulence have found that compressible velocity fluctuations create both effective fluid transport properties and an effective equation of state. This paper investigates the effects of rotation on compressible turbulence. It is shown that rotation modifies the transport properties of compressible turbulence by replacing the turbulence time scale by a rotational time scale, much as rotation modifies the transport properties of incompressible turbulence. But thermal equilibrium properties are modified in a more complex manner. Two regimes are possible: one dominated by incompressible fluctuations, in which the sound speed is modified as it is in non-rotating turbulence, and a rotation dominated regime in which the sound speed enhancement is rotation dependent. The dimensionless parameter which discriminates between regimes is identified. In general, rotation is found to suppress the effects of compressibility. A novel feature of the present analysis is the use of a non-Kolmogorov steady state as the reference state of turbulence. introduction of such steady states expands the power and utility of analytical turbulence closures to a wider range of problems.

Rubinstein, Robert↗

NASA NTRS1993

Discrete ordinates solutions of nongray radiative transfer with diffusely reflecting walls

Nongray gas radiation in a plane parallel slab bounded by gray, diffusely reflecting walls is studied using the discrete ordinates method. The spectral equation of transfer is averaged over a narrow wavenumber interval preserving the spectral correlation effect. The governing equations are derived by considering the history of multiple reflections between two reflecting wails. A closure approximation is applied so that only a finite number of reflections have to be explicitly included. The closure solutions express the physics of the problem to a very high degree and show relatively little error. Numerical solutions are obtained by applying a statistical narrow-band model for gas properties and a discrete ordinates code. The net radiative wail heat fluxes and the radiative source distributions are obtained for different temperature profiles. A zeroth-degree formulation, where no wall reflection is handled explicitly, is sufficient to predict the radiative transfer accurately for most cases considered, when compared with increasingly accurate solutions based on explicitly tracing a larger number of wail reflections without any closure approximation applied.

Menart, J. A.↗

NASA NTRS1976

Calculation methods for compressible turbulent boundary layers

Calculation procedures for non-reacting compressible two- and three-dimensional turbulent boundary layers were reviewed. Integral, transformation and correlation methods, as well as finite difference solutions of the complete boundary layer equations summarized. Alternative numerical solution procedures were examined, and both mean field and mean turbulence field closure models were considered. Physics and related calculation problems peculiar to compressible turbulent boundary layers are described. A catalog of available solution procedures of the finite difference, finite element, and method of weighted residuals genre is included. Influence of compressibility, low Reynolds number, wall blowing, and pressure gradient upon mean field closure constants are reported.

Bushnell, D. M.↗

NASA NTRS1973

A linear stochastic model of the human operator

A linear stochastic model of the human operator is developed and applied to the problem of piloted control of an aircraft. The pilot and aircraft are modeled as linear time-invariant systems containing both process and measurement noise. The loop closure by the pilot is determined by formulating the problem as an optimal stochastic control problem. The solution to the optimal control problem yields not only the pilot's optimal control output which he uses to control the vehicle, but also the optimal combination of his observations of the vehicle states upon which the pilot bases his control. A method is presented so that, using experimental pilot vehicle data, the cost functional which is minimized in the optimal control problem will be numerically equal to the pilot rating that the pilot would associate with the given vehicle and task.

Durrett, J. C.↗

NASA NTRS1990

Mathematical theory of a relaxed design problem in structural optimization

Various attempts have been made to construct a rigorous mathematical theory of optimization for size, shape, and topology (i.e. layout) of an elastic structure. If these are represented by a finite number of parametric functions, as Armand described, it is possible to construct an existence theory of the optimum design using compactness argument in a finite dimensional design space or a closed admissible set of a finite dimensional design space. However, if the admissible design set is a subset of non-reflexive Banach space such as L(sup infinity)(Omega), construction of the existence theory of the optimum design becomes suddenly difficult and requires to extend (i.e. generalize) the design problem to much more wider class of design that is compatible to mechanics of structures in the sense of variational principle. Starting from the study by Cheng and Olhoff, Lurie, Cherkaev, and Fedorov introduced a new concept of convergence of design variables in a generalized sense and construct the 'G-Closure' theory of an extended (relaxed) optimum design problem. A similar attempt, but independent in large extent, can also be found in Kohn and Strang in which the shape and topology optimization problem is relaxed to allow to use of perforated composites rather than restricting it to usual solid structures. An identical idea is also stated in Murat and Tartar using the notion of the homogenization theory. That is, introducing possibility of micro-scale perforation together with the theory of homogenization, the optimum design problem is relaxed to construct its mathematical theory. It is also noted that this type of relaxed design problem is perfectly matched to the variational principle in structural mechanics.

Kikuchi, Noboru↗

NASA NTRS1980

Self-replicating systems: A systems engineering approach

A first approach to conceptualize self-replicating systems was developed from past and present abstract theories. The engineering elements of self-replicating systems are defined in terms of a basic reference system. A number of options are investigated. The growth characteristics and their problems are analyzed, the mathematics of various exponential growth options are outlined, and the problems of universal parts production and systems closure are discussed. Selected areas of further study are defined and a 20 year development and demonstration program is presented.

Vontiesenhausen, G.↗

NASA NTRS2002

Fluid Mechanics of Capillary-Elastic Instabilities in Microgravity Environment

The aim of this project is to investigate the closure and reopening of lung airways due to surface tension forces, coupled with airway elasticity. Airways are liquid-lined, flexible tubes and closure of airways can occur by a Rayleigh instability of the liquid lining, or an instability of the elastic support for the airway as the surface tension of the air-liquid interface pulls the tube shut, or both. Regardless of the mechanism, the airway is closed because the liquid lining has created a plug that prevents axial gas exchange. In the microgravity environment, surface tension forces dominate lung mechanics and would lead to more prevalent, and more uniformly distributed air-way closure, thereby creating a potential for respiratory problems for astronauts. Once closed the primary option for reopening an airway is by deep inspiration. This maneuver will pull the flexible airways open and force the liquid plug to flow distally by the incoming air stream. Airway reopening depends to a large extent on this plug flow and how it may lead to plug rupture to regain the continuity of gas between the environment and the alveoli. In addition to mathematical modeling of plug flows in liquid-lined, flexible tubes, this work has involved benchtop studies of propagating liquid plugs down tube networks that mimic the human airway tree. We have extended the work to involve animal models of liquid plug propagation in rat lungs. The liquid is radio-opaque and x-ray video imaging is used to ascertain the movement and distribution of the liquid plugs so that comparisons to theory may be made. This research has other uses, such as the delivery of liquids or drugs into the lung that may be used for surfactant replacement therapy or for liquid ventilation.

Grotberg, James B.↗

NASA NTRS2012

Study of Shortwave Spectra in Fully 3D Environment: Synergy Between Scanning Radars and Spectral Radiation Measurements

The main theme for our research is the understanding and closure of the surface spectral shortwave radiation problem in fully 3D cloud situations by combining the new ARM scanning radars, shortwave spectrometers, and microwave radiometers with the arsenal of radiative transfer tools developed by our group. In particular, we define first a large number of cloudy test cases spanning all 3D possibilities not just the customary uniform-overcast ones. Second, for each case, we define a "Best Estimate of Clouds That Affect Shortwave Radiation" using all relevant ARM instruments, notably the new scanning radars, and contribute this to the ARM Archive. Third, we test the ASR-signature radiative transfer model RRTMG_SW for those cases, focusing on the near-IR because of long-standing problems in this spectral region, and work with the developers to improve RRTMG_SW in order to increase its penetration into the modeling community.

Wiscombe, Warren J.↗

NASA NTRS2011

Spectral Invariance Principles Observed in Spectral Radiation Measurements of the Transition Zone

The main theme for our research is the understanding and closure of the surface spectral shortwave radiation problem in fully 3D cloud situations by combining the new ARM scanning radars, shortwave spectrometers, and microwave radiometers with the arsenal of radiative transfer tools developed by our group. In particular, we define first a large number of cloudy test cases spanning all 3D possibilities not just the customary uniform-overcast ones. Second, for each case, we define a "Best Estimate of Clouds That Affect Shortwave Radiation" using all relevant ARM instruments, notably the new scanning radars, and contribute this to the ARM Archive. Third, we test the ASR-signature radiative transfer model RRTMG_SW for those cases, focusing on the near-IR because of long-standing problems in this spectral region, and work with the developers to improve RRTMG_SW in order to increase its penetration into the modeling community.

Marshak, Alexander↗

NASA NTRS1987

Second-order closure models for rotating turbulent flows

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler k-epsilon and k-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

Speziale, Charles G.↗

NASA NTRS1985

Second-order closure models for rotating turbulent flows

The physical properties of the commonly used second-order closure models are examined theoretically for rotating turbulent flows. Comparisons are made with results which are a rigorous consequence of the Navier-Stokes equations for the problem of a fully-developed turbulent channel flow in a rapidly rotating framework. It is demonstrated that all existing second-order closures yield spurious physical results for this test problem of rotating channel flow. In fact, the results obtained are shown to be substantially more unphysical than those obtained from the simpler K-epsilon and K-l models. Modifications in the basic structure of these second-order closure models are proposed which can alleviate this problem.

Speziale, C. G.↗

NASA NTRS1987

Reynolds stress closure in jet flows using wave models

Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.

Morris, P. J.↗

NASA NTRS1993

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗

NASA NTRS2017

Uncertainty Quantification of CFD Data Generated for a Model Scramjet Isolator Flowfield

Computational fluid dynamics is now considered to be an indispensable tool for the design and development of scramjet engine components. Unfortunately, the quantification of uncertainties is rarely addressed with anything other than sensitivity studies, so the degree of confidence associated with the numerical results remains exclusively with the subject matter expert that generated them. 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. Given the limitations of current hypersonic ground test facilities, this expanded role is believed to be a requirement by some in the hypersonics community if scramjet engines are to be given serious consideration as a viable propulsion system. The present effort describes a simple, relatively low cost, nonintrusive approach to uncertainty quantification that includes the basic ingredients required to handle both aleatoric (random) and epistemic (lack of knowledge) sources of uncertainty. The nonintrusive nature of the approach allows the computational fluid dynamicist to perform the uncertainty quantification with the flow solver treated as a "black box". Moreover, a large fraction of the process can be automated, allowing the uncertainty assessment to be readily adapted into the engineering design and development workflow. In the present work, the approach is applied to a model scramjet isolator problem where the desire is to validate turbulence closure models in the presence of uncertainty. In this context, the relevant uncertainty sources are determined and accounted for to allow the analyst to delineate turbulence model-form errors from other sources of uncertainty associated with the simulation of the facility flow.

Baurle, R. A.↗

NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of their predictions.

Speziale, Charles G.↗

NASA NTRS1989

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of thier predictions.

Speziale, Charles G.↗

NASA NTRS1990

A critical comparison of turbulence models for homogeneous shear flows in a rotating frame

A variety of turbulence models, including five second-order closure models and four two equation models, are tested for the problem of homogeneous turbulent shear flow in a rotating frame. The model predictions for the time evolution of the turbulent kinetic energy and dissipation rate, as well as those for the equilibrium states, are compared with the results of physical and numerical experiments. Most of the two-equation models predict the same results for all rotation rates (omega/S) in which there is an exponential time growth of the turbulent kinetic energy and dissipation rate. The second-order closures are qualitatively superior since, consistent with physical and numerical experiments, they only predict this type of unstable flow for intermediate rotation rates in the range -0.1 less than or equal to omega/S less than or equal to 1.6. For rotation rates outside this range, there is an exchange of stabilities with a solution whose kinetic energy and dissipation rate decay with time. Although the second-order closures are superior to the two-equation models, there are still problems with the quantitative accuracy of their predictions.

Speziale, Charles G.↗