Engineering topics
Publications and source records attributed to Zang, T. A..
Shielding design on large human-rated systems allows minimization of radiation impact on electronic systems. Shielding design tools require adequate methods for evaluation of design layouts, guiding qualification testing, and adequate follow-up on final design evaluation.
An overview of multidisciplinary optimization (MDO) methodology and two applications of this methodology to the preliminary design phase are presented. These applications are being undertaken to improve, develop, validate and demonstrate MDO methods. Each is presented to illustrate different aspects of this methodology. The first application is an MDO preliminary design problem for defining the geometry and structure of an aerospike nozzle of a linear aerospike rocket engine. The second application demonstrates the use of the Framework for Interdisciplinary Design Optimization (FIDO), which is a computational environment system, by solving a preliminary design problem for a High-Speed Civil Transport (HSCT). The two sample problems illustrate the advantages to performing preliminary design with an MDO process.
As a tool for large-scale computations in fluid dynamics, spectral methods were prophesized in 1944, born in 1954, virtually buried in the mid-1960's, resurrected in 1969, evangalized in the 1970's, and catholicized in the 1980's. The use of spectral methods for meteorological problems was proposed by Blinova in 1944 and the first numerical computations were conducted by Silberman (1954). By the early 1960's computers had achieved sufficient power to permit calculations with hundreds of degrees of freedom. For problems of this size the traditional way of computing the nonlinear terms in spectral methods was expensive compared with finite-difference methods. Consequently, spectral methods fell out of favor. The expense of computing nonlinear terms remained a severe drawback until Orszag (1969) and Eliasen, Machenauer, and Rasmussen (1970) developed the transform methods that still form the backbone of many large-scale spectral computations. The original proselytes of spectral methods were meteorologists involved in global weather modeling and fluid dynamicists investigating isotropic turbulence. The converts who were inspired by the successes of these pioneers remained, for the most part, confined to these and closely related fields throughout the 1970's. During that decade spectral methods appeared to be well-suited only for problems governed by ordinary diSerential eqllations or by partial differential equations with periodic boundary conditions. And, of course, the solution itself needed to be smooth. Some of the obstacles to wider application of spectral methods were: (1) poor resolution of discontinuous solutions; (2) inefficient implementation of implicit methods; and (3) drastic geometric constraints. All of these barriers have undergone some erosion during the 1980's, particularly the latter two. As a result, the applicability and appeal of spectral methods for computational fluid dynamics has broadened considerably. The motivation for the use of spectral methods in numerical calculations stems from the attractive approximation properties of orthogonal polynomial expansions.
We consider the transition to turbulence in three-dimensional reconnection of a magnetic neutral sheet. We find that the transition can occur via a three-step process. First, the sheet undergoes the usual tearing instability. Second, the tearing mode saturates to form a two-dimensional quasi-steady state. Third, this secondary equilibrium is itself unstable when it is perturbed by three-dimensional disturbances. Most of this paper is devoted to the analysis and simulation of the three-dimensional linear stability properties of the two-dimensional saturated tearing layer. The numerical simulations are performed with a semi-implicit, pseudospectral-Fourier collocation algorithm. We identify a three-dimensional secondary linear stability which grows on the ideal timescale. An examination of the modal energetics reveals that the largest energy transfer is from the mean field to the three-dimensional field, with the two-dimensional field acting as a catalyst.
The ongoing effort at NASA Langley Research Center aimed at developing transitional zone models for Reynolds-averaged Navier-Stokes calculations is described. Results from three models that we have examined - two zero-equation models and a two-equation model of the k-omega type - are presented for high-speed flows. The use of analytical and computational tools to provide guidance and to provide the transitional flow-field data bases necessary to develop such models is discussed. Some results from validation of these tools are given; qualitative comparisons between modeled and theoretical/computational representations of the transitional zone are provided.
New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.
This paper reports results from the numerical implementation and testing of the compressible large eddy-simulation (LES). Relevant quantities from 32-cubed coarse grid LES solutions are compared with results generated from direct numerical simulations (DNS) of three-dimensional compressible turbulence that have been run both with sufficient resolution, at 96 cubed. The 32 cubed LES results overall agree well with their 96 cubed DNS counterparts. Moreover, the new DNS results confirm several recent conclusions about compressible turbulence that have been based primarily on two-dimensional simulations.
Temporal direct numerical simulation of laminar breakdown via subharmonic secondary instability in high-speed axisymmetric boundary layers has been accomplished using a highly accurate, fully explicit algorithm which combines spectral collocation and high-order compact-difference techniques. Numerical test cases confirm that subharmonic secondary instability is confirmed to be a viable path to transition in high-speed boundary-layer flow. Secondary instability is shown to account for peaks in the Reynolds stresses at or near the critical layer which are not possible from the second-mode primary instability alone. Reynolds stresses spatially reconstructed from the temporal model via the Gaster transformation show a 'spreading angle' of about 12 deg, in qualitative agreement with experimental findings. The rate of broadening of the Reynolds stress peak is a strongly nonlinear phenomenon which cannot be reproduced by secondary instability theory.
The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.
A 3D time-dependent numerical simulation of the foot-point stressing of coronal magnetic field was developed in order to relate coronal activity with the stressing of the coronal magnetic field by foot-point motions at the photosphere. The results of the simulation did not reveal magnetic reconnection, kinking, or the formation of concave-up magnetic field lines suitable for prominence formation. It is concluded that, contrary to many models, photospheric twisting of a single arcade does not lead to the type of processes required to explain solar activity.
A filtering analysis of a turbulent flow was developed which provides details of the path of the kinetic energy of the flow from its creation via thermal production to its dissipation. A low-pass spatial filter is used to split the velocity and the temperature field into a filtered component (composed mainly of scales larger than a specific size, nominally the filter width) and a fluctuation component (scales smaller than a specific size). Variables derived from these fields can fall into one of the above two ranges or be composed of a mixture of scales dominated by scales near the specific size. The filter is used to split the kinetic energy equation into three equations corresponding to the three scale ranges described above. The data from a direct simulation of the Rayleigh-Benard problem for conditions where the flow is turbulent are used to calculate the individual terms in the three kinetic energy equations. This is done for a range of filter widths. These results are used to study the spatial location and the scale range of the thermal energy production, the cascading of kinetic energy, the diffusion of kinetic energy, and the energy dissipation. These results are used to evaluate two subgrid models typically used in large-eddy simulations of turbulence. Subgrid models attempt to model the energy below the filter width that is removed by a low-pass filter.
New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.
This paper summarizes some of the recent progress made at NASA Langley Research Center in the understanding, prediction and modeling of high speed boundary-layer transition. Linear and nonlinear theories together with large-eddy and direct numerical simulations have been used to understand various aspects of the transition problem while low disturbance 'quiet' tunnels provide means for validating the theoretical results.
Many aspects of solar activity are believed to be due to the stressing of the coronal magnetic field by footpoint motions at the photosphere. The results are presented of a fully spectral numerical simulation which is the first 3-D time dependent simulation of footpoint stressing in a geometry appropriate for the corona. An arcade is considered that is initially current-free and impose a smooth footpoint motion that produces a twist in the field of approx 2 pi. The footprints were fixed and the evolution was followed until the field relaxes to another current-free state. No evidence was seen for any instability, either ideal or resistive and no evidence for current sheet formation. The most striking feature of the evolution is that in response to photospheric motions, the field expands rapidly upward to minimize the stress. The expansion has two important effects. First, it suppresses the development of dips in the field that could support dense, cool material. For the motions assumed, the magnetic field does not develop a geometry suitable for prominence formation. Second, the expansion inhibits ideal instabilities such as kinking. The results indicate that simple stearing of a single arcade is unlikely to lead to solar activity such as flares or prominences. Effects are discussed that might possibly lead to such activity.
Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.
A subgrid-scale model recently derived for use in the large-eddy simulation of compressible turbulent flows is examined from a fundamental theoretical and computational standpoint. It is demonstrated that this model, which is applicable only to compressible turbulent flows in the limit of small density fluctuations, correlates somewhat poorly with the results of direct numerical simulations of compressible isotropic turbulence at low Mach numbers. An alternative model, based on Favre-filtered fields, is suggested which appears to reduce these limitations.
New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.
The overall architecture of the multipurpose parallel-processing Navier-Stokes Computer (NSC) being developed by Princeton and NASA Langley (Nosenchuck et al., 1986) is described and illustrated with extensive diagrams, and the NSC implementation of an elementary multigrid algorithm for simulating isotropic turbulence (based on solution of the incompressible time-dependent Navier-Stokes equations with constant viscosity) is characterized in detail. The present NSC design concept calls for 64 nodes, each with the performance of a class VI supercomputer, linked together by a fiber-optic hypercube network and joined to a front-end computer by a global bus. In this configuration, the NSC would have a storage capacity of over 32 Gword and a peak speed of over 40 Gflops. The multigrid Navier-Stokes code discussed would give sustained operation rates of about 25 Gflops.