Search NASASearch

Engineering topics

Hussaini, M. Yousuff

Publications and source records attributed to Hussaini, M. Yousuff.

At least 19 records

A Stochastic Collocation Algorithm for Uncertainty Analysis

This report describes a stochastic collocation method to adequately handle a physically intrinsic uncertainty in the variables of a numerical simulation. For instance, while the standard Galerkin approach to Polynomial Chaos requires multi-dimensional summations over the stochastic basis functions, the stochastic collocation method enables to collapse those summations to a one-dimensional summation only. This report furnishes the essential algorithmic details of the new stochastic collocation method and provides as a numerical example the solution of the Riemann problem with the stochastic collocation method used for the discretization of the stochastic parameters.

Mathelin, Lionel

Uncertainty Propagation for Turbulent, Compressible Flow in a Quasi-1D Nozzle Using Stochastic Methods

This paper describes a fully spectral, Polynomial Chaos method for the propagation of uncertainty in numerical simulations of compressible, turbulent flow, as well as a novel stochastic collocation algorithm for the same application. The stochastic collocation method is key to the efficient use of stochastic methods on problems with complex nonlinearities, such as those associated with the turbulence model equations in compressible flow and for CFD schemes requiring solution of a Riemann problem. Both methods are applied to compressible flow in a quasi-one-dimensional nozzle. The stochastic collocation method is roughly an order of magnitude faster than the fully Galerkin Polynomial Chaos method on the inviscid problem.

Zang, Thomas A.

On the Exploitation of Sensitivity Derivatives for Improving Sampling Methods

Many application codes, such as finite-element structural analyses and computational fluid dynamics codes, are capable of producing many sensitivity derivatives at a small fraction of the cost of the underlying analysis. This paper describes a simple variance reduction method that exploits such inexpensive sensitivity derivatives to increase the accuracy of sampling methods. Three examples, including a finite-element structural analysis of an aircraft wing, are provided that illustrate an order of magnitude improvement in accuracy for both Monte Carlo and stratified sampling schemes.

Cao, Yanzhao

Towards Perfectly Absorbing Boundary Conditions for Euler Equations

In this paper, we examine the effectiveness of absorbing layers as non-reflecting computational boundaries for the Euler equations. The absorbing-layer equations are simply obtained by splitting the governing equations in the coordinate directions and introducing absorption coefficients in each split equation. This methodology is similar to that used by Berenger for the numerical solutions of Maxwell's equations. Specifically, we apply this methodology to three physical problems shock-vortex interactions, a plane free shear flow and an axisymmetric jet- with emphasis on acoustic wave propagation. Our numerical results indicate that the use of absorbing layers effectively minimizes numerical reflection in all three problems considered.

Hayder, M. Ehtesham

Self-Contained Automated Methodology for Optimal Flow Control

This paper describes a self-contained, automated methodology for active flow control which couples the time-dependent Navier-Stokes system with an adjoint Navier-Stokes system and optimality conditions from which optimal states, i.e., unsteady flow fields and controls (e.g., actuators), may be determined. The problem of boundary layer instability suppression through wave cancellation is used as the initial validation case to test the methodology. Here, the objective of control is to match the stress vector along a portion of the boundary to a given vector; instability suppression is achieved by choosing the given vector to be that of a steady base flow. Control is effected through the injection or suction of fluid through a single orifice on the boundary. The results demonstrate that instability suppression can be achieved without any a priori knowledge of the disturbance, which is significant because other control techniques have required some knowledge of the flow unsteadiness such as frequencies, instability type, etc. The present methodology has been extended to three dimensions and may potentially be applied to separation control, re-laminarization, and turbulence control applications using one to many sensors and actuators.

Joslin, Ronald D.

A self-contained, automated methodology for optimal flow control validated for transition delay

This paper describes a self-contained, automated methodology for flow control along with a validation of the methodology for the problem of boundary layer instability suppression. The objective of control is to match the stress vector along a portion of the boundary to a given vector; instability suppression is achieved by choosing the given vector to be that of a steady base flow, e.g., Blasius boundary layer. Control is effected through the injection or suction of fluid through a single orifice on the boundary. The present approach couples the time-dependent Navier-Stokes system with an adjoint Navier-Stokes system and optimality conditions from which optimal states, i.e., unsteady flow fields, and control, e.g., actuators, may be determined. The results demonstrate that instability suppression can be achieved without any a priori knowledge of the disturbance, which is significant because other control techniques have required some knowledge of the flow unsteadiness such as frequencies, instability type, etc.

Joslin, Ronald D.

Active control of instabilities in laminar boundary-layer flow. Part 1: An overview

This paper (the first in a series) focuses on using active-control methods to maintain laminar flow in a region of the flow in which the natural instabilities, if left unattended, lead to turbulent flow. The authors review previous studies that examine wave cancellation (currently the most prominent method) and solve the unsteady, nonlinear Navier-Stokes equations to evaluate this method of controlling instabilities. It is definitely shown that instabilities are controlled by the linear summation of waves (i.e., wave cancellation). Although a mathematically complete method for controlling arbitrary instabilities has been developed (but not yet tested), the review, duplication, and physical explanation of previous studies are important steps for providing an independent verification of those studies, for establishing a framework for subsequent work which will involve automated transition control, and for detailing the phenomena by which the automated studies can be used to expand knowledge of flow control.

Joslin, Ronald D.

Active Control of Instabilities in Laminar Boundary-Layer Flow: Use of Sensors and Spectral Controller - Part 2

This study focuses on the suppression of instability growth using an automated active-control technique. The evolution of 2D disturbances that are spatially growing in a flat-plate boundary layer are computed with a spatial DNS code. A controller receives wall sensor information (pressure or shear) as input and provides a signal that controls an actuator response as output. The control law assumes that wave cancellation is valid. The results indicate that a measure of wave cancellation can be obtained for small- and large-amplitude instabilities without feedback; however, feedback is required to optimize the control amplitude and phase for exact wave cancellation.

Joslin, Ronald D.

Linear and nonlinear PSE for compressible boundary layers

Compressible stability of growing boundary layers is studied by numerically solving the partial differential equations under a parabolizing approximation. The resulting parabolized stability equations (PSE) account for nonparallel as well as nonlinear effects. Evolution of disturbances in compressible flat-plate boundary layers are studied for freestream Mach numbers ranging from 0 to 4.5. Results indicate that the effect of boundary-layer growth is important for linear disturbances. Nonlinear calculations are performed for various Mach numbers. Two-dimensional nonlinear results using the PSE approach agree well with those from direct numerical simulations using the full Navier-Stokes equations while the required computational time is less by an order of magnitude. Spatial simulation using PSE were carried out for both the fundamental and subharmonic type breakdown for a Mach 1.6 boundary layer. The promising results obtained show that the PSE method is a powerful tool for studying boundary-layer instabilities and for predicting transition over a wide range of Mach numbers.

Chang, Chau-Lyan

On the linear stability of compressible plane Couette flow

The linear stability of compressible plane Couette flow is investigated. The correct and proper basic velocity and temperature distributions are perturbed by a small amplitude normal mode disturbance. The full small amplitude disturbance equations are solved numerically at finite Reynolds numbers, and the inviscid limit of these equations is then investigated in some detail. It is found that instability can occur, although the stability characteristics of the flow are quite different from unbounded flows. The effects of viscosity are also calculated, asymptotically, and shown to have a stabilizing role in all the cases investigated. Exceptional regimes to the problem occur when the wavespeed of the disturbances approaches the velocity of either of the walls, and these regimes are also analyzed in some detail. Finally, the effect of imposing radiation-type boundary conditions on the upper (moving) wall (in place of impermeability) is investigated, and shown to yield results common to both bounded and unbounded flows.

Duck, Peter W.

Effects of shock on the stability of hypersonic boundary layers

A set of linearized shock boundary conditions is derived, which is then imposed at the shock to account for the interaction of the shock wave with the boundary/shock layer instability wave; these boundary conditions are used to study the effect of shock on hypersonic boundary layer stability under the assumption of quasi-parallel flow. The result show that the shock has little effect on the boundary layer instability (subsonic first and second mode disturbances) when the shock is located outside the boundary layer edge. When the shock is located near the boundary layer edge, it exerts a stabilizing influence on the first and second modes. The shock also induces unstable supersonic modes with oscillatory structure in the shock layer, but these modes grow slower than the subsonic modes.

Chang, Chau-Lyan

Multiple paths to subharmonic laminar breakdown in a boundary layer

Numerical simulations demonstrate that laminar breakdown in a boundary layer induced by the secondary instability of two-dimensional Tollmien-Schlichting waves to three-dimensional subharmonic disturbancews need not take the conventional lambda vortex/high-shear layer path.

Zang, Thomas A.

On the large-eddy simulation of transitional wall-bounded flows

The structure of the subgrid scale fields in plane channel flow has been studied at various stages of the transition process to turbulence. The residual stress and subgrid scale dissipation calculated using velocity fields generated by direct numerical simulations of the Navier-Stokes equations are significantly different from their counterparts in turbulent flows. The subgrid scale dissipation changes sign over extended areas of the channel, indicating energy flow from the small scales to the large scales. This reversed energy cascade becomes less pronounced at the later stages of transition. Standard residual stress models of the Smagorinsky type are excessively dissipative. Rescaling the model constant improves the prediction of the total (integrated) subgrid scale dissipation, but not that of the local one. Despite the somewhat excessive dissipation of the rescaled Smagorinsky model, the results of a large eddy simulation of transition on a flat-plate boundary layer compare quite well with those of a direct simulation, and require only a small fraction of the computational effort. The inclusion of non-dissipative models, which could lead to further improvements, is proposed.

Piomelli, Ugo

Multiple paths to subharmonic laminar breakdown in a boundary layer

Numerical simulations demonstrate that laminar breakdown in a boundary layer induced by the secondary instability of two-dimensional Tollmien-Schlichting waves to three-dimensional subharmonic disturbances need not take the conventional lambda vortex/high-shear layer path.

Zang, Thomas A.

Spectral methods for CFD

One of the objectives of these notes is to provide a basic introduction to spectral methods with a particular emphasis on applications to computational fluid dynamics. Another objective is to summarize some of the most important developments in spectral methods in the last two years. The fundamentals of spectral methods for simple problems will be covered in depth, and the essential elements of several fluid dynamical applications will be sketched.

Zang, Thomas A.

Resolution requirements for numerical simulations of transition

The resolution requirements for direct numerical simulations of transition to turbulence are investigated. A reliable resolution criterion is determined from the results of several detailed simulations of channel and boundary-layer transition.

Zang, Thomas A.

Multidomain spectral solution of shock-turbulence interactions

The use of a fitted-shock multidomain spectral method for solving the time-dependent Euler equations of gasdynamics is described. The multidomain method allows short spatial scale features near the shock to be resolved throughout the calculation. Examples presented are of a shock-plane wave, shock-hot spot and shock-vortex street interaction.

Kopriva, David A.

A spectral multi-domain technique applied to high-speed chemically reacting flows

The first applications of a spectral multidomain method for viscous compressible flow is presented. The method imposes a global flux balance condition at the interface so that high-order continuity of the solution is preserved. The global flux balance is imposed in terms of a spectral integral of the discrete equations across adjoining domains. Since the discretized equations interior to each domain solved are uncoupled from each other, and since the interface relation has a block structure, the solution scheme can be adapted to the particular requirements of each subdomain. The spectral multidomain technique presented is well-suited for the multiple scales associated with the chemically reacting and transition flows in hypersonic research. A nonstaggered multidomain discretization is used for the chemically reacting flow calculation, and the first implementation of a staggered multidomain mesh is presented for accurately solving the stability equation for a viscous compressible fluid.

Macaraeg, Michele G.