ALBERT

Description: A theoretical proof of the optimal rate of convergence for the least-squares method is developed for the Stokes problem based on the velocity-pressure-vorticity formula. The 2D Stokes problem is analyzed to define the product space and its inner product, and the a priori estimates are derived to give the finite-element approximation. The least-squares method is found to converge at the optimal rate for equal-order interpolation.

Description: Aluminum plates of various length/width ratios loaded in compression are examined analytically in terms of buckling and postbuckling for comparison with classical theory. The plates are considered to be subjected to longitudinal compressive displacement with no edge stress and are assumed to be simply supported at the side edges. The average longitudinal direct stresses are computed to plot characteristic curves by means of the following theories: 3D flexibility, higher-order shear deformation, first-order shear deformation, and classical von Karman. The 3D flexibility approach yields the lowest results with more unknowns than the higher-order shear-deformation method. The 3D flexibility approach is considered to be the most accurate, and it is shown that the calculated resultants and displacements tend to vary when consideration is given to the effects of transverse shear.

Description: Present interest in hypersonic vehicles has resulted in a renewed interest in thermal stress analysis of airframe structures. While there are numerous texts and papers on thermal stress analysis, practical examples and experience on light gage aircraft structures are fairly limited. A research program has been undertaken at General Dynamics to demonstrate the present state of the art, verify methods of analysis, gain experience in their use, and develop engineering judgement in thermal stress analysis. The approach for this project has been to conduct a series of analyses of this sample problem and compare analysis results with test data. This comparison will give an idea of how to use our present methods of thermal stress analysis, and how accurate we can expect them to be.

Description: This paper presents the fuzzy dynamical reliability and failure probability as well as the basic principles and the analytical method of loss assessment for nonlinear seismic steel structures. Also presented is the optimization formulation and a numerical example for double objectives, initial construction cost and expected failure loss, and dynamical reliability constraints. The earthquake ground motion is based on a stationary filtered non-white noise and the fuzzy damage grade is described by damage index.

Description: Effective methods of approximate eigensolution reanalysis of modified nonclassically damped structures are developed. For structures with passive or active discrete damping devices or with damping treatment, the system becomes nonproportionally damped and the computation of its dynamic responses may require the use of complex modes. For larger systems, the computation of complex modes is very expensive. Thus it is desirable to have approximate reanalysis techniques for the efficient evaluation of the effect of design changes. In recent years, the assumed mode reanalysis method was successfully applied to minimum weight design of undamped structures with natural frequency constraints. The accuracy of the assumed mode reanalysis method can be improved dramatically if the global approximation function includes the normal modes of the original system and their derivatives. This approach was demonstrated to be effective even for a system with shape changes. The approach used by Noor et al. for eigensolution reanalysis of undamped structures is extended to treat a nonclassically damped system.

Description: With structural design in mind, a new unified variational model has been developed which represents the mechanics of deformation elasto-plasticity with unilateral contact conditions. For a design problem formulated as maximization of the load carrying capacity of a structure under certain constraints, the unified model allows for a simultaneous analysis and design synthesis for a whole range of mechanical behavior.

Description: Personnel from NASA Ames Research Center presented a paper on establishing a benchmark experimental data base for generic hypersonic vehicle shape for validation and/or calibration of advanced computational fluid dynamics computer codes. The need for this capability is based on a requirement for extensive hypersonic data to fully validate CFD codes to be used for NASP and other hypersonic vehicles. The use of wind tunnel models in the Ames 3.5-ft Hypersonic Wind Tunnel to obtain pertinent surface and flow-field data over a broad range of test conditions is described.

Description: Personnel from NASA Lewis Research Center spoke on the subject of inlet duct and nozzle high speed validation experiments to include crossing shocks and boundary layer interaction, unsteady shock/boundary layer interactions, and vortex generators. Attention was also paid to the subjects of high speed mixing and transition ducts. Specific application was made to the NASP hypermixing concepts.

Description: NASA Lewis Research Center personnel presented a summary of LeRC validation experiments. The subjects discussed included: iced wing, linear transonic cascade, transition duct aerodynamics, transition duct heat transfer, low speed centrifugal compressor, turbomachinery blade row interactions, and three dimensional fluid mechanics.

Description: A presentation by Langley Research Center covered subjects of: LaRC approach to CFD code validation, experimental CFD perceptions, CFD code validation program experiment, and highlights of the experiment. The objective of the validation program and the approach taken are discussed.

Description: The real benefit of structural optimization techniques is in the application of these techniques to large structures such as full vehicles or full aircraft. For these structures, however, the sequential computer's time and memory requirements prohibit the solutions. With the rapid development of parallel computers, parallel processing of large scale structural optimization problems is achievable. In this paper we discuss the parallel processing of structural optimization problems with parallel structural analysis. Two different types of interface between the optimization and analysis routines are developed and tested.

Description: In most structural optimization problems the implicit behavior constraints are evaluated for successive modifications in the design. For each trial design, the analysis equations must be solved and the multiple repeated analyses usually involve extensive computational effort. This difficulty motivated several studies on explicit approximations of the structural behavior in terms of the design variables. The latter approach can considerably reduce the amount of computations, but the quality of the approximations might not be sufficient. Many of the approximate behavior models proposed in the past are valid only for relatively small changes in the design variables. The accuracy of the results is often insufficient for large changes in the design. The object of this study is to present efficient and high quality approximations of the structural behavior. It will be shown that the quality of the approximations can greatly be improved by combining scaling of the initial design, using intervening variables, and scaling a set of fictitious loads. Integrating these means, a powerful solution procedure can be introduced. In addition, the errors in satisfying the analysis equations can readily be evaluated. A numerical example illustrates the solution methodology and the effectiveness of the proposed approach.

Description: The weight savings due to usage of composite materials in aircraft structural applications is well known. Significant weight and cost benefits are achievable by developing structurally tailored concepts and efficiently integrating them with suitable material and manufacturing technologies. The proposed paper will describe such an efficient concept for application to primary aircraft structures.

Description: This paper first presents a brief review on the application of optimization and active control of seismic structures along with some of the author's recent work. It then assesses the practicality and future development of seismic structural optimization, and some practical problems associated with active control.

Description: The flow characteristics of a low speed centrifugal compressor were examined at NASA Lewis Research Center to improve understanding of the flow in centrifugal compressors, to provide models of various flow phenomena, and to acquire benchmark data for three dimensional viscous flow code validation. The paper describes the objectives, test facilities' instrumentation, and experiment preliminary comparisons.

Description: The characteristics of flow over a rearward facing step is discussed in a paper presented by NASA Lewis Research Center personnel. The objective is to obtain data which will be used for validation of Direct Numerical Simulations being developed at NASA Langley and Ames. Two phases of the work are identified as a study of large scale structures in the flow using hot-wire/cold wire techniques and the development of a statistical data base for using three-component laser velocimetry.

Description: NASA Ames Center personnel presented data on stages of code development and corresponding experiments in the application of computational fluid dynamics for aeronautical investigations. Specific subjects included algorithms, grid generation, facilities, instrumentation, and data acquisition. Numerical simulation and flow modelling were described to show the procedure for calibration and validation.

Description: The paper presents the theoretical bases and implementation techniques of sensitivity analyses for efficient structural optimization of large structures, based on finite element static and dynamic analysis methods. The sensitivity analyses have been implemented in conjunction with two methods for optimization, namely, the Mathematical Programming and Optimality Criteria methods. The paper discusses the implementation of the sensitivity analysis method into our in-house software package, AutoDesign.

Description: Ill-conditioned systems arising in analysis and optimization can display a high sensitivity to numerical precision for changes and errors in data input. Such data may be in the form of system parameter input or desired system response. The ill-conditioning referred to generally arises from the lack of sufficient independent data to define a complex system or the weak sensitivity of response to source input parameters. It is shown how small errors in data and assumed fixed and known parameters can lead to highly erroneous results in ill-conditioned linear algebraic equations. A simplified detection and correction of critical input data arising in the coefficient matrix and desired response (i.e., right hand side) is proposed.

Description: An identification procedure proposed by Shen and Taylor to determine the crack characteristics (location and size of the crack) from dynamic measurements is tested. This procedure was based on minimization of either the 'mean-square' measure of difference between measurement data (natural frequencies and mode shapes) and the corresponding predictions obtained from the computational model. The procedure is tested for simulated damage in the form of symmetric cracks in a simply-supported Bernoulli-Euler beam. The sensitivity of the solution of damage identification problems to the values of parameters that characterize damage is discussed. A sensitivity formula is derived.

Description: The issue of the utility of multilevel decomposition and optimization remains controversial. To date, only the structural optimization community has actively developed and promoted multilevel optimization techniques. However, even this community acknowledges that multilevel optimization is ideally suited for a rather limited set of problems. It is warned that decomposition typically requires eliminating local variables by using global variables and that this in turn causes ill-conditioning of the multilevel optimization by adding equality constraints. The purpose is to suggest a new multilevel optimization technique. This technique uses behavior variables, in addition to design variables and constraints, to decompose the problem. The new technique removes the need for equality constraints, simplifies the decomposition of the design problem, simplifies the programming task, and improves the convergence speed of multilevel optimization compared to conventional optimization.

Description: In recent years there have been several hierarchic multilevel optimization algorithms proposed and implemented in design studies. Equality constraints are often imposed between levels in these multilevel optimizations to maintain system and subsystem variable continuity. Equality constraints of this nature will be referred to as coupling equality constraints. In many implementation studies these coupling equality constraints have been handled indirectly. This indirect handling has been accomplished using the coupling equality constraints' explicit functional relations to eliminate design variables (generally at the subsystem level), with the resulting optimization taking place in a reduced design space. In one multilevel optimization study where the coupling equality constraints were handled directly, the researchers encountered numerical difficulties which prevented their multilevel optimization from reaching the same minimum found in conventional single level solutions. The researchers did not explain the exact nature of the numerical difficulties other than to associate them with the direct handling of the coupling equality constraints. The coupling equality constraints are handled directly, by employing the Generalized Reduced Gradient (GRG) method as the optimizer within a multilevel linear decomposition scheme based on the Sobieski hierarchic algorithm. Two engineering design examples are solved using this approach. The results show that the direct handling of coupling equality constraints in a multilevel optimization does not introduce any problems when the GRG method is employed as the internal optimizer. The optimums achieved are comparable to those achieved in single level solutions and in multilevel studies where the equality constraints have been handled indirectly.

Description: Reliability was determined for two degrading dynamic systems subject to random load processes. Damage is caused by loss of components for Daniels systems and crack extension for plates with cracks. The analysis accounted for the coupling between response and current damage state of the system. It is based on mean crossing rates of conditional processes and properties of diffusion models. Simple systems are used to illustrate proposed methods for estimating reliability.

Description: A passive vibration damping technique that is referred to as 'Non-Obstructive Particle Damping (NOPD)' is presented. The NOPD technique consists of making small diameter holes (or cavities) at appropriate locations inside vibrating structures and filling these holes to appropriate levels with particles which yield the maximum damping effectiveness for the desired mode (or modes). Powders, spherical shaped, metallic, non-metallic or liquid particles (or mixtures) with different densities, viscosities and adhesive or cohesive characteristics can be used.

Description: Recent results from sensitivity analysis for strain energy with anisotropic elasticity are applied to thickness and orientational design of laminated membranes. Primarily, the first order gradients of the total elastic energy are used in an optimality criteria based method. This traditional method is shown to give slow convergence with respect to design parameters, although the convergence of strain energy is very good. To get a deeper insight into this rather general characteristic, second order derivatives are included and it is shown how they can be obtained by first order sensitivity analysis. Examples of only thickness design, only orientational design, and combined thickness--orientational design are presented.

Description: NASA Marshall Space Flight Center personnel presented a paper on the status of MSFC computational fluid dynamics application and validation activities. Subjects discussed included the Space Shuttle Main Engine studies, unsteady multistage turbine loads, fuel pump discharge volutes, and injector LOX inlet results based on fundamental flows, subcomponents, and interactive components/systems.

Description: The objective is to set up and analyze benchmark problems on multibody dynamics and to verify the predictions of two multibody computer simulation codes. TREETOPS and DISCOS have been used to run three example problems - one degree-of-freedom spring mass dashpot system, an inverted pendulum system, and a triple pendulum. To study the dynamics and control interaction, an inverted planar pendulum with an external body force and a torsional control spring was modeled as a hinge connected two-rigid body system. TREETOPS and DISCOS affected the time history simulation of this problem. System state space variables and their time derivatives from two simulation codes were compared.

Description: The non-linear dynamic analysis of large structures is always very time, effort and CPU consuming. Whenever possible the reduction of the size of the mathematical model involved is of main importance to speed up the computational procedures. Such reduction can be performed for the part of the structure which perform linearly. Most of the time, the classical Guyan reduction process is used. For non-linear dynamic process where the non-linearity is present at interfaces between different structures, Craig-Bampton models can provide a very rich information, and allow easy selection of the relevant modes with respect to the phenomenon driving the non-linearity. The paper presents the employment of Craig-Bampton models combined with Newmark direct integration for solving non-linear friction problems appearing at the interface between the Hubble Space Telescope and its solar arrays during in-orbit maneuvers. Theory, implementation in the FEM code ASKA, and practical results are shown.

Description: Structural optimization problems are mostly solved under constraints from statics, such as stresses, strains, or displacements under static loads. But in the design process, dynamic quantities like eigenfrequencies or accelerations under dynamic loads become more and more important. Therefore, it is obvious that constraints from dynamics must be considered in structural optimization packages. This paper addresses the dynamics branch in MBB-LAGRANGE. It will concentrate on two topics, namely on the different formulations for eigenfrequency constraints and on frequency response constraints. For the latter the necessity of a system reduction is emphasized. The methods implemented in LAGRANGE are presented and examples are given.

Description: This paper addresses the weight minimization of a circular plate-like structure which resulted in a 26 percent weight reduction. The optimization was performed numerically with the COPES/ADS program using the modified method of feasible directions. Design parameters were the inner thickness and outer thickness of the plate with constraints on maximum yield stress and maximum transverse displacement. Also, constraints were specified for the upper and lower bounds of the fundamental frequency and plate thicknesses. The MSC/NASTRAN finite element program was used for the evaluation of response variables. Original and final designs of the plate were tested using an Instron tension-compression machine to compare finite element results to measured strain data. The difference between finite element strain components and measured strain data was within engineering accuracy.

Description: Any optimization of structures for maximum stability or for maximum dynamic stiffness deals with an eigenvalue problem. The goal of this optimization is to raise the lowest eigenvalue (or eigenvalues) of the problem to its highest (optimal) level at a constant volume of the structure. Likely the lowest eigenvalue may be either inherently multi-modal or it can become multi-modal as a result of the optimization process. The multimodeness introduces some ambiguity to the eigenvalue problem and make the optimization difficult to handle. Thus far, only the simplest cases of multi-modal structures have been effectively optimized using rather elaborate analytical methods. Numerous publications report design of a minimum volume structure with different eigenvalues constraints, in which, however, the modality of the problem is assumed a priori. The method presented here utilizes a multi-modal optimality criteria and allows for inclusion of an arbitrary number of buckling or vibrations modes which might influence the optimization process. The real multi-modality of the problem, that is the number of modes participating in the final optimal design is determined iteratively. Because of a natural use of the FEM technique the method is easy to program and might be helpful in design of large flexible space structures.

Description: The structure-control system optimization problem is formulated with constraints on the closed-loop eigenvalues and the efficiency of the reduced order system. The feasibility of the approach is illustrated by designing the ACOSS-FOUR structure with a reduced order system and improving the efficiency characteristics of the structures-control system.

Description: In this paper an optimal design procedure is introduced to improve the overall performance of nonlinear framed structures. The design methodology presented here is a multiple-objective optimization procedure whose objective functions involve the buckling eigenvalues and eigenvectors of the structure. A constant volume with bounds on the design variables is used in conjunction with an optimality criterion approach. The method provides a general tool for solving complex design problems and generally leads to structures with better limit strength and stability. Many algorithms have been developed to improve the limit strength of structures. In most applications geometrically linear analysis is employed with the consequence that overall strength of the design is overestimated. Directly optimizing the limit load of the structure would require a full nonlinear analysis at each iteration which would be prohibitively expensive. The objective of this paper is to develop an algorithm that can improve the limit-load of geometrically nonlinear framed structures while avoiding the nonlinear analysis. One of the novelties of the new design methodology is its ability to efficiently model and design structures under multiple loading conditions. These loading conditions can be different factored loads or any kind of loads that can be applied to the structure simultaneously or independently. Attention is focused on optimal design of space framed structures. Three-dimensional design problems are more complicated to carry out, but they yield insight into real behavior of the structure and can help avoiding some of the problems that might appear in planar design procedure such as the need for out-of-plane buckling constraint. Although researchers in the field of structural engineering generally agree that optimum design of three-dimension building frames especially in the seismic regions would be beneficial, methods have been slow to emerge. Most of the research in this area has dealt with the optimization of truss and plane frame structures.

Description: Fluid management impacts strongly on the optimization of space construction. Large quantities of liquids are needed for propellants and life support. The mass of propellant liquids is comparable to that required for the structures. There may be a strong dynamic interaction between the stored liquids and the space structure unless the design minimizes the interaction. The constraints of cost and time required optimization of the supply/resupply strategy. The proper selection and design of the fluid management methods for: slosh control; stratification control; acquisition; transfer; gauging; venting; dumping; contamination control; selection of tank configuration and size; the storage state and the control system can improve the entire system performance substantially. Our effort consists of building mathematical/computer models of the various fluid management methods and testing them against the available experimental data. The results of the models are used as inputs to the system operations studies. During the past year, the emphasis has been on modeling: the transfer of cryogens; sloshing and the storage configuration. The work has been intermeshed with ongoing NASA design and development studies to leverage the funds provided by the Center.

Description: In the design of complex built-up structures that are made of truss, beam, membrane, shell, and solid, there are five different kinds of design variables: material property, sizing, shape, configuration, and topological variables. Previous research has shown that the improvement in performances obtained by altering the configuration of structural components can be much more significant than those obtained when the geometry is assumed to be fixed. Using the variational approach, a unified design sensitivity was developed for the first three kinds of design variables, and was further extended recently in many structural analysis problems such as nonlinear, structural dynamics, and frequency response analysis. A continuum design sensitivity analysis method is developed for the configuration design variable of built-up structures.

Description: A method is described for generating unstructured meshes of triangles or tetrahedra for computational domains of complex geometrical shape. To illustrate the power of the approach, it is applied to the solution of flows past several complete aircraft configurations. The advocated approach allows for the natural incorporation of mesh adaptivity and this is demonstrated for both inviscid and viscous computations in two and three dimensions.

Description: Techniques and applications of algebraic grid generation are described. The techniques are univariate interpolations and transfinite assemblies of univariate interpolations. Because algebraic grid generation is computationally efficient, the use of interactive graphics in conjunction with the techniques is advocated. A flexible approach, which works extremely well in an interactive environment, called the control point form of algebraic grid generation is described. The applications discussed are three-dimensional grids constructed about airplane and submarine configurations.

Description: Progress to date on the development of a method for turbulent, wall-bounded flow which uses the defect stream function formulation in the outer layer and an analytic law of the wall and wake formulation in the inner region is reviewed. This two-formulation approach avoids the need to computationally resolve the high-gradient inner layer. One of the most appealing recent developments is the transformation of the compressible governing equation for the defect stream function into a linear, second-order differential equation which has analytic solutions for many problems of practical interest. Numerical and analytic results for incompressible and compressible flows are shown to be in excellent agreement with experimental results. In this paper the two-formulation approach is applied to primitive-variable computations. Excellent comparisons with experiment are presented for two compressible flat plate flows.

Description: Computer generated databases containing velocity and pressure fields in three-dimensional space at a sequence of time-steps, were used for the investigation of near-wall turbulence structures, their space-time evolution, and their associated pressure fields. The main body of the results were obtained from simulation data for turbulent channel flow at a Reynolds number of 180 (based on half-channel height and friction velocity) with a grid of 128 x 129 x and 128 points. The flow was followed over a total time of 141 viscous time units. Spanwise centering of the detected structures was found to be essential in order to obtain a correct magnitude of the associated Reynolds stress contribution. A positive wall-pressure peak is found immediately beneath the center of the structure. The maximum amplitude of the pressure pattern was, however, found in the buffer region at the center of the shear-layer. It was also found that these flow structures often reach a maximum strength in connection with an asymmetric spanwise motion, which motivated the construction of a conditional sampling scheme that preserved this asymmetry.

Description: A review is presented of the approach, motivation and early results of a reevaluation of the knowledge collected by the research community during thirty years of research on the structure of turbulent boundary layers. Four distinctions or criteria concerning the need to improve the practices in this field are discussed: data versus inferences, possible versus actual versus significant events, one structure versus many, and one or two flows versus the totality of data. Attention given to the known quasi-coherent structures show that three factors significantly reduce the information available from laboratory data as a basis for forming a complete model of the quasi-coherent structure in the turbulent boundary layer: (1) the necessity to ensemble average probe output with resulting loss of phase information, (2) the inability to see more than one or two of the various types of structures simultaneously, and (3) the inability to see the spatial relations between the various quasi-coherent structures resulting from (1) and (2).

Description: The two-point correlation tensor is used to analyze near-wall structures in turbulent channel flow. Simulation results have been utilized to obtain the correlation tensor with sufficient spatial resolution to permit accurate differentiation for computing the vorticity correlation, and the resolution of sharp features such as the shear layer and its corresponding region of spanwise vorticity. Proper orthogonal decomposition and stochastic estimation were observed to yield similar results when appropriate conditions were employed for the estimation. Resulting structures were such that they would be detected by standard conditional sampling methods. The results of decomposing the vorticity field and the velocity field were significantly different, indicating that if a coherent structure is found that dominates the velocity fields, the curl of that velocity structure will not dominate the vorticity field.

Description: Recent analysis of databases generated by direct numerical simulations of homogeneous turbulent shear flows have revealed the presence of coherent structures similar to those in turbulent boundary layers. In this paper these findings and tentative conclusions on their significance are discussed.

Description: This paper presents a simple new algorithm for constructing second-order models for any given flexible structure, with these models satisfying the constraint that each leading principal submatrix of the mass matrix must have a specified condition number. Now, the effects of rounding errors in practical computation generally increase with increasing condition number, i.e. as the matrix considered approaches singularity. Thus, the new models should prove useful for the systematic testing of the numerical properties of many algorithms in the dynamics and control of flexible structures.

Description: The hopscotch scheme is examined to see whether it can compete effectively with implicit schemes for the integration of the Navier-Stokes equations on a vector processing machine. This scheme is stable as long as the Courant number is less than or equal to one and it does not underestimate viscous effects. The accuracy of the scheme is tested on one- and two-dimensional problems whose exact solutions are known. The scheme is then used to simulate flows around a circular cylinder with Reynolds numbers 200, 500 and 1000.

Description: The stress-ratio algorithm associated with fully stressed design philosophy has been used as a convenient tool to achieve minimum weight design of strength-limited structures. The algorithm is effective and converges quickly for many cases. However, it presents extremely slow oscillatory iteration histories for plate-thickness design problems that involve transverse bending loads. Modification of the basic algorithm presented in this paper provides an effective remedy to this problem when both membrane and bending loads are present. The modified resizing algorithm requires numerical solutions of a fourth-order algebraic equation. No additional data, beyond the ordinary static analysis results, are required.

Description: A method for performing a global/local stress analysis is described and its capabilities are demonstrated. The method employs spline interpolation functions which satisfy the linear plate bending equation to determine displacements and rotations from a global model which are used as 'boundary conditions' for the local model. Then, the local model is analyzed independently of the global model of the structure. This approach can be used to determine local, detailed stress states for specific structural regions using independent, refined local models which exploit information from less-refined global models. The method presented is not restricted to having a priori knowledge of the location of the regions requiring local detailed stress analysis. This approach also reduces the computational effort necessary to obtain the detailed stress state. Criteria for applying the method are developed. The effectiveness of the method is demonstrated using a classical stress concentration problem and a graphite-epoxy blade-stiffened panel with a discontinuous stiffener.

Description: An analytical investigation is conducted to determine the shape of a growing delamination and the distribution of the energy release rate along the delamination front in a laminated composite double cantilever beam specimen. Distributions of the energy release rate for specimens with straight delamination fronts and delamination front contours for delaminations whose growth is governed by the fracture criterion that G = Gc at all points are predicted as a function of material properties and delamination length. The predicted delamination front contours are utilized to ascertain the effect of the changing shape of the delamination front on the value of the critical strain energy release rate as computed from double cantilever beam fracture toughness test data.

Description: An adaptive unstructured remeshing technique is applied to transient thermal-structural analysis. The effectiveness of the technique, together with the finite element method and an error estimation technique, is evaluated by two applications which have exact solutions: (1) the steady-state thermal analysis of a plate subjected to a highly localized surface heating, and (2) the transient thermal-structural analysis of a simulated convectively cooled leading edge subjected to a translating heat source. These applications demonstrate that the remeshing technique significantly reduces the problem size as well as the analysis solution error as compared to the results produced using standard structured meshes.

Description: The evolution of two miscible liquids meeting at an initially sharp interface inside a cavity under microgravity g-jitter conditions is studied numerically. The response of the interface and kinematics of the flowfield to various g-jitter accelerations and aspect ratio variations is shown. The interface region acts like a vortex source sheet, and it can be unstable to Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The vortices produced along the interface can serve as a stirring mechanism to promote local mixing. Below the critical Stokes-Reynolds number, the destabilization of the interface results in deformation into wavy structures. In some parameter regions, these structures oscillate in time; in others, they are quasi-steady. Above the critical Stokes-Reynolds number, 'chaotic' instability results, and the interface breaks into concentration pockets. The morphology of the initial breakup is similar to observed Rayleigh-Taylor instability. Subsequent mixing of the two fluids after the breakup of the interface is then very rapid.

Description: Development of probabilistic structural analysis methods for hot structures is a major activity at NASA-Lewis, and consists of five program elements: (1) probabilistic loads, (2) probabilistic finite element analysis, (3) probabilistic material behavior, (4) assessment of reliability and risk, and (5) probabilistic structural performance evaluation. Attention is given to quantification of the effects of uncertainties for several variables on High Pressure Fuel Turbopump blade temperature, pressure, and torque of the Space Shuttle Main Engine; the evaluation of the cumulative distribution function for various structural response variables based on assumed uncertainties in primitive structural variables; evaluation of the failure probability; reliability and risk-cost assessment; and an outline of an emerging approach for eventual hot structures certification. Collectively, the results demonstrate that the structural durability/reliability of hot structural components can be effectively evaluated in a formal probabilistic framework. In addition, the approach can be readily extended to computationally simulate certification of hot structures for aerospace environments.

Description: This paper focuses on the design, analysis, and test of a high temperature structural panel, constructed of Ti-6-4, subjected to a variety of thermal and mechanical load conditions. A follow-on panel, constructed of TMC, is also discussed in less detail. The design constraints and test set-up are discussed, as well as the test methods that were used: the grid shadow moire method and a single gage force stiffness method. The agreement between the test data and analysis for this test program provides confidence in the methods that are currently being used to design structures for hypersonic vehicles. The agreement also suggests that postbuckled strength may potentially be used to reduce the vehicle weight.

Description: The time-dependent thermo-viscoplastic response of aerospace structures subjected to intense aerothermal loads is predicted using the finite-element method. The finite-element analysis uses the Bodner-Partom unified viscoplastic constitutive relations to determine rate-dependent nonlinear material behavior. The methodology is verified by comparison with experimental data and other numerical results for a uniaxially-loaded bar. The method is then used (1) to predict the structural response of a rectangular plate subjected to line heating along a centerline, and (2) to predict the thermal-structural response of a convectively-cooled engine cowl leading edge subjected to aerodynamic shock-shock interference heating. Compared to linear elastic analysis, the viscoplastic analysis results in lower peak stresses and regions of plastic deformations.

Description: Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Description: An efficient procedure for obtaining the compatibility conditions of finite-element models involves the generation of both field and compatibility conditions from deformation-displacement relations, using (1) the compatibility bandwidth, and (2) the node-determinacy concept. A computer program thus structured will generate sparse and banded compatibility conditions for a structure that is idealized by the finite elements.

Description: In addition to characterizing the various concepts reported in the literature on longitudinally-ribbed surface for aerodynamic surface drag reduction, the present development status evaluation of this technology correlates all available experimental data. An analysis of these data is then conducted to ascertain the parameters most directly involved in drag reduction, and to evaluate the effects which have thus far been exerted on turbulent boundary layer structures. Such advanced riblet techniques as compound and three-dimensional riblets, riblets in combination with large-eddy breakup devices, and riblets with suction/blowing, are also discussed.

Description: The physical processes involved in the boundary layer instability-generation and transition process seem to pose basic restrictions on the implementation of active, wave-based transition-control methods. While suppression short of elimination is possible by these means for primary disturbances, wave-interaction instabilities demand that the control be implemented almost immediately after the appearance of primary disturbances; even slight delays can negate the intended beneficial effects, and this basic problem is exacerbated at the higher Reynolds numbers typical of aircraft in cruising flight. Three-dimensional disturbances are noted to be important in this context, together with the continuous regeneration of TS waves.

Description: The present paper describes recent advances and trends in finite element developments and applications for solidification problems. In particular, in comparison to traditional methods of approach, new enthalpy-based architectures based on a generalized trapezoidal family of representations are presented which provide different perspectives, physical interpretation and solution architectures for effective numerical simulation of phase change processes encountered in solidification problems. Various numerical test models are presented and the results support the proposition for employing such formulations for general phase change applications.

Description: Formulations of inviscid flux splitting algorithms for chemical nonequilibrium gases are presented. A chemical system for air dissociation and recombination is described. Numerical results for one-dimensional shock tube and nozzle flows of air in chemical nonequilibrium are examined.

Description: A robust self-starting explicit architecture for computational structural dynamics is described. The proposed methodology involves expressing the governing equations of motion in conservation form and temporal discretization is accomplished in the spirit of the Lax-Wendroff type formulations. The development of the basic methodology is shown. Discretization in space is accomplished by introducing stress-based representations and employing the classical Galerkin scheme. Numerical test model results are presented which validate the architecture.

Description: A systematic analysis is developed which describes the isothermal counterdiffusion of two gases in the presence of a third nondiffusing gas. Four distinct regimes are identified for diffusive slip. The nondiffusing gas reduces the pressure drop in all cases, whether it raises or lowers the mean density. The gas leads to nonzero diffusive slip even in the limit of very disparate masses for the diffusing gases. The implications of these findings for controlling concentration creep in crystal growth are discussed.

Description: The behavior of solitons induced by boundary excitation is investigated at various time-dependent conditions and different unperturbed water depths, using the Korteweg-de Vries (KdV) equation. Then, solitons induced from Boussinesq equations under similar conditions were studied, making it possible to remove the restriction in the KdV equation and to treat soliton head-on collisions (as well as overtaking collisions) and reflections. It is found that the results obtained from the KdV and the Boussinesq equations are in good agreement.

Description: The algorithm employed in the present incompressible two-dimensional calculations of an impulsively-started lid-driven cavity has its basis in the time-dependent stream-function equation. While a Crank-Nicholson differencing scheme is used for the diffusion terms, the Adams-Bashforth scheme is used for the convection terms. The periodic asymptotic solutions obtained for Reynolds numbers of 5000 and 10,000 are found to be precisely periodic; it is demonstrated that they have reached asymptotic states. The indicators of that achievement are discussed.

Description: The problem addressed is that of obtaining reduced-order component models for use in simulating the dynamics of a multibody system. In certain cases, nonlinear system models may be constructed using linear dynamic models for each component, but allowing large angle motion between components. Without some form of model reduction, system models constructed in this manner may be too large for use in control system design and simulation trades. This paper analyzes one method of component model reduction that allows systems level requirements (e.g., capturing the effect of body 1 reaction wheel noise on body 2 camera pointing) to aid in the selection of the reduced-order component models. Briefly stated, important modes are selected at the system level and projected onto the components, and reduced-order components are then assembled into a reduced-order system model that retains the projected modes.

Description: A theory proposed by Concus and Finn (1974) and recently developed by Finn (1983 and 1984) yields explicit geometrical criteria for the position of the free surface of a liquid at zero gravity in a cylindrical container of specified cross section. These criteria were applied by Concus and Finn to three container geometries: the bathtub, the trapezoid, and the keyhole. It is possible to find geometrical criteria that promise a liquid interface of finite height, with the base still covered with liquid, or a liquid interface of infinite height, with the liquid wetting a well defined portion of the wall. In the present work, calculations are presented for a fourth geometry, the non-concentric cylinders. In addition, the earlier calculations of Concus and Finn are extended, and a unified graphical presentation of all four geometries is given that can be used directly for the design of containments for liquids at zero gravity.

Description: A simple model, parameterized by the Reynolds stress anisotropy, is proposed for the spectrum of weakly anisotropic turbulence. It contains a model constant that affects its region of realizability. This spectrum model is used to derive a one-point closure to the rapid pressure-strain term. The derived pressure-strain model is linear in the Reynolds stress anisotropy and is of the same form as the closure model of Launder et al. (1975). The spectrum model becomes unrealizable in some regions of wave space for sufficiently large anisotropy of the Reynolds stress, and this is used to infer the region of validity of the linear closure model. It is found that the extent of the valid region is very small when the model constant is set to match rapid distortion theory, and largest for a model-constant set close to the value suggested by LRR. However, even the largest valid domain does not extend very far from isotropy, suggesting inherent weakness in the linear pressure-strain models.

Description: An experimental technique is described for obtaining time-resolved heat flux measurements with high-frequency response (up to 100 kHz) in a steady-flow ambient-temperature facility. The heat transfer test object is preheated and suddenly injected into an established steady flow. Thin-film gages deposited on the test surface detect the unsteady substrate surface temperature. Analog circuitry designed for use in short-duration facilities and based on one-dimensional semiinfinite heat conduction is used to perform the temperature/heat flux transformation. A detailed description of substrate properties, instrumentation, experimental procedure, and data reduction is given, along with representative results obtained in the stagnation region of a circular cylinder subjected to a wake-dominated unsteady flow. An in-depth discussion of related work is also provided.

Description: A procedure for the minimum weight design of helicopter rotor blades with constraints on multiple coupled flap-lag natural frequencies, autorotational inertia, and centrifugal stress is presented. Optimum designs are obtained for blades with both rectangular and tapered planforms and are compared within a reference blade. The effects of higher-frequency constraints and stress constraints on the optimum blade designs are assessed. The results indicate that there is an increase in blade weight and a significant change in the design variable distributions with an increase in the number of frequency constraints. The inclusion of stress constraints has different effects on the wall thickness distributions of rectangular and tapered blades, but tends to increase the magnitude of the nonstructural segment weight distributions for both blade types.

Description: A near-wall turbulence model and its incorporation into a multiple-timescale turbulence model are presented. The near-wall turbulence model is obtained from a k-equation turbulence model and a near-wall analysis. In the method, the equations for the conservation of mass, momentum, and turbulent kinetic energy are integrated up to the wall, and the energy transfer and the dissipation rates inside the near-wall layer are obtained from algebraic equations. Fully developed turbulent channel and pipe flows are solved using a finite element method. The computational results compare favorably with experimental data. It is also shown that the turbulence model can resolve the overshoot phenomena of the turbulent kinetic energy and the dissipation rate in the region very close to the wall.

Description: Current research in the field of advanced adaptive structures for space applications is reviewed. A classification of adaptive structures is proposed whereby such structures are subdivided into adaptive, sensory, controlled, active, and intelligent structures. The definition and properties of each type of adaptive structures are presented, and methods of structure control are discussed.

Description: The simultaneous active control of flexural and extensional vibrations in elastic beams is experimentally investigated. The results demonstrate that using pairs of piezoceramic transducers, whose elements are symmetrically located and independently controlled by a multichannel adaptive controller, enables the high attenuation of both flexural and extensional response. This capability is due to the nature of the piezoceramic element, which when bonded to the surface of the structure and electrically excited, exerts a surface strain on the structure. This strain enables input of both shear forces and moments into the structural system. The results are applicable to many situations where extensional vibrations couple to large flexural vibrations and subsequently radiate significant sound levels.

Description: Results on the transition from soft to hard turbulence in simulations of two-dimensional Boussinesq convection are reported. The computed probability densities for temperature fluctuations are exponential in form in both soft and hard turbulence, unlike what is observed in experiments. In contrast, a change is obtained in the Nusselt number scaling on Rayleigh number in good agreement with the three-dimensional experiments.

Description: 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.

Description: A revised formal solution of the vibrating ribbon problem of hydrodynamic stability is presented. The initial formulation of Gaster (1965) is modified by application of the Briggs method and a careful treatment of the complex double Fourier transform inversions. Expressions are obtained in a natural way for the discrete spectrum as well as for the four branches of the continuous spectra. These correspond to discrete and branch-cut singularities in the complex wavenumber plane. The solutions from the continuous spectra decay both upsteam and downstream of the ribbon, with the decay in the upstream direction being much more rapid than that in the downstream direction. Comments and clarification of related prior work are made.

Description: The literature on the mathematical modeling of large space structures is first reviewed, with attention given to continuum models, model order reduction, substructuring, and computational techniques. System identification and mode verification are then discussed with reference to the verification of mathematical models of large space structures. In connection with analysis, the paper surveys recent research on eigensolvers and dynamic response solvers for large-order finite-element-based models.

Description: Numerical solutions of the Navier-Stokes equations using explicit schemes can be obtained at the expense of efficiency. Conventional implicit methods which often achieve fast convergence rates suffer high cost per iteration. A new implicit scheme based on lower-upper factorization and symmetric Gauss-Seidel relaxation offers very low cost per iteration as well as fast convergence. High efficiency is achieved by accomplishing the complete vectorizability of the algorithm on oblique planes of sweep in three dimensions.

Description: Some of the hardware and software tools and techniques in use at NASA's numerical aerodynamic simulation facility for the analysis of computational fluid dynamics are described. The visualization process can be illustrated by video tapes and stereo pictures. Although these visualization tools have dramatically improved the ability to conduct research in fluid dynamics, a comparison of the current environment for analysis with an 'ideal' environment illustrates that there are still major improvements that should be made. The most time-consuming task in future analyses of the increasingly complex computer simulations will be the extraction and clear display of the key features. In addition, the interface between the workstation and the scientist should be improved significantly. Current research on techniques for creating these improvements is described.

Description: This paper discusses the research in adaptive structures conducted at NASA Langley Research Center. The objective in the research program on adaptive truss structures is to develop an integrated approach for the design, test, and evaluation of adaptive trusses. An adaptive structure must include sensors for measuring some of the states of the mechanism, a controller that processes the sensor information and generates command signals, and a device that will cause the struts in the mechanism to change length. Based on the results of the research program, it is concluded that alternative designs for future space applications are offered, and adaptive truss has been demonstrated to be effective for vibration suppression.

Description: The mobility power flow (MPF) approach is used in this paper to describe the flexural behavior of an L-shaped plate structure consisting of thick plates with rotary inertia and shear deformation effects included in the analysis. The introduction of the thick plate effects significantly increases the complexity of the structural mobility functions used in the definitions of the power flow terms; however, because of the substructuring that is used in the MPF approach, the complexity of the problem is significantly reduced as compared to solving for the global structure. Additionally, with the MPF approach the modal behavior is described. The MPF analysis of the L-shaped plate is performed for the case of point force excitation on one plate, with the two plates being identical in both size and thickness. The results of this analysis are compared to results from the finite-element analysis (FEA) and the statistical energy analysis (SEA) and show very good agreement in the low- and high-frequency regimes, respectively.

Description: The study deals with the approximation of the effect of all velocity-increment moments on the spatial evolution of the pdf for maximally correlated structures by the effect of the second moment alone by assuming independence of separated length scales. The dependence of the two-point velocity statistics on length-scale or separation is found to be governed by the gradients of the stress, rather than of the velocity. This may provide a new approach for predicting the domain patterns observed in turbulent flows and a new means of characterizing or classifying different structures, namely, by their nonlinear diffusion coefficient.

Description: A method is presented for tailoring plate and shell composite structures for optimal forced damped dynamic response. The damping of specific vibration modes is optimized with respect to dynamic performance criteria including placement of natural frequencies and minimization of resonance amplitudes. The structural composite damping is synthesized from the properties of the constituent materials, laminate parameters, and structural geometry based on a specialty finite element. Application studies include the optimization of laminated composite beams and composite shells with fiber volume ratios and ply angles as design variables. The results illustrate the significance of damping tailoring to the dynamic performance of composite structures, and the effectiveness of the method in optimizing the structural dynamic response.

Description: Analytical and experimental techniques for the prediction and ground verification of the damped structural dynamics of space structures are developed. The options available for similarity-scaled model testing, including replica and multiple scale approaches, are reviewed. For the case when the distortion of potentially dissipative or nonlinear joints, which would be required in multiple-scale modeling, is impractical, a new type of modeling is introduced, which uses a hybrid of joints at replica scale and connecting elements at a modified multiple scale. The model design requirements for replica, multiple-scale, and hybrid models are developed, and the expected scaling of nonlinear dissipation in joints is derived. A damping prediction scheme is developed that relies on a finite element model of the undamped structure and measurements of the individual joint properties to predict the modal damping of the truss attributable to the joints. A hybrid-scaled model of a segment of the Space Station was built and dynamically tested. The predicted and measured truss damping compared favorably.

Description: Elastic-plastic behavior of fibrous composite laminates is analyzed for coupled in-plane mechanical loads and uniform thermal changes. Constitutive equations of the individual fibrous layers are derived from a vanishing fiber diameter model that represents the essential axial constraint between the phases. This permits derivation of closed form equations for the overall yield condition, stress concentration factors, and instantaneous compliance. Thermoelastic properties of the phases and yield stress of the matrix phase are functions of temperature. The effect of the model assumptions on the predicted behavior of composite laminates is examined by comparing the calculated response under cyclic thermal changes to available theoretical results and experimental measurements.

Description: An elasticity solution has been used to analyze matrix stresses near the fiber/matrix interface in continuous fiber-reinforced metal-matrix composites, modeling the micromechanics in question in terms of a cylindrical fiber and cylindrical matrix sheath which is embedded in an orthotropic medium representing the composite. The model's predictions for lamina thermal and mechanical properties are applied to a laminate analysis determining ply-level stresses due to thermomechanical loading. A comparison is made between these results, which assume cylindrical symmetry, and the predictions yielded by a FEM model in which the fibers are arranged in a square array.

Description: A dual-control test system measuring the strength of graphite/epoxy specimens under tension and compression for a wide range of bearing-bypass load ratios is presented. The bearing-bypass strengths of single-fastener specimens are measured for damage onset and ultimate failure and are plotted against one another. The results indicate that an initiated damage grows to failure in the same mode for most cases. However, if the damage is initiated in the compression-reacted-bearing mode, the specimens fail in the offset-compression mode. It is concluded that this type of transition can occur in multi-fastener joints under compression and complicate the strength predictions.

Description: Acoustic emission in 51-cm diameter graphite/epoxy pressure vessels was monitored during pressurization (hydrotesting). Several vessels were subjected to impact by a blunt impactor, but only after the vessels had been proofed; that is, pressurized to 80 percent of nominal burst pressure as determined from control (unimpacted) vessels. AE activity was then monitored throughout a series of successively higher pressure cycles ranging from 10 to 60 percent of ultimate. Each cycle included a ramp up to pressure followed by a 4-min hold period and then pressure unload. The event rate was high, and especially modified AE analyzers had to be used to acquire the data. This paper presents the AE event count versus pressure history of these tests and demonstrates the ability of the AE technique to monitor the growth of damage and to estimate the effect on ultimate strength. The number of events that occurred during pressure holds proved to be a reasonable estimator of vessel performance.

Description: An experimental investigation was conducted to determine the degree to which the thermal contact conductance at the interface of contacting Aluminum 6061 T6 surfaces could be enhanced through the use of vapor-deposited metallic coatings. Three different coating materials (lead, tin, and indium) were evaluated using four different thicknesses for each coating material. The results verified the existence of an optimum coating thickness, shown to be in the range of 2.0 to 3.0 microns for indium, 1.5 to 2.5 microns for lead, and 0.2 to 0.5 microns for tin. The enhancement factors for thermal contact conductance were found to be on the order of 700, 400, and 50 percent, respectively. Based upon the experimental data, the hardness of the coating materials appears to be the most significant parameter in ranking the substrate and coating material combinations; however, additional experimental data are needed to substantiate this hypothesis. Finally, it was apparent that the thermal contact conductance enhancement effect was greatest at low contact pressures and decreased significantly with increases in the contact pressure.

Description: The spray characteristics of several different simplex pressure-swirl nozzles are examined using water as the working fluid. Measurements of mean drop size, drop-size distribution, effective spray cone angle, and circumferential liquid distribution are carried out over wide ranges of injection pressure. Eight different nozzles are employed in order to achieve a wide variation in the length/diameter ratio of the final discharge orifice. Generally, it is found that an increase in discharge orifice length/diameter ratio, l(o)/d(o), increases the mean drop size in the spray and reduces the spray cone angle. The circumferential liquid distribution is most uniform when l(o)/d(o) = 2. If l(o)/d(o) is raised above or lowered below this optimum value, the circumferential uniformity of the liquid distribution is impaired. The observed effects of l(o)/d(o) on spray characteristics are generally the same regardless of whether the change in l(o)/d(o) is accomplished by varying l(o) or d(o).

Description: The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a decrease in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. Terms in the Reynolds stress transport equations were computed. The balances indicate that the decrease in turbulent kinetic energy is a result of a decrease in turbulence production, along with an increase in turbulent dissipation. Intuitive reasoning and current turbulence models would predict an increase in kinetic energy along with increases in production and dissipation rates as a result of increased mean-flow strain rate. Later in the evolution of the flow, both turbulence production and dissipation increase.

Description: A general boundary condition formalism is developed for all types of boundary conditions to which hyperbolic systems are subject; the formalism makes possible a 'cookbook' approach to boundary conditions, by means of which novel boundary 'recipes' may be derived and previously devised ones may be consulted as required. Numerous useful conditions are derived for such CFD problems as subsonic and supersonic inflows and outflows, nonreflecting boundaries, force-free boundaries, constant pressure boundaries, and constant mass flux. Attention is given to the computation and integration of time derivatives.

Description: Long-wave instabilities in a directionally-solidified binary mixture may occur in several limits. Sivashinsky (1983) identified a small-segregation-coefficient limit and obtained a weakly nonlinear evolution equation governing subcritical two-dimensional bifurcation. Brattkus and Davis (1988) identified a near-absolute-stability limit and obtained a strongly nonlinear evolution equation governing supercritical two-dimensional bifurcation. The present investigation identifies a third strongly nonlinear evolution equation, arising in the small-segregation-coefficient, large-surface-energy limit. This equation links both of the former and describes the change from the sub- to super-critical bifurcations. This study sets the previous long-wave analyses into a logical framework.

Description: A theory to explain the initial stages of unsteady separation was proposed by Van Dommelen and Cowley (1989). This theory is verified for the separation process that occurs at the equatorial plane of a sphere or a spheroid which is impulsively spun around an axis of symmetry. A Lagrangian numerical scheme is developed which gives results in good agreement with Eulerian computations, but which is significantly more accurate. This increased accuracy, and a simpler structure to the solution, also allows verification of the Eulerian structure, including the presence of logarithmic terms. Further, while the Eulerian computations broke down at the first occurrence of separation, it is found that the Lagrangian computation can be continued. It is argued that this separated solution does provide useful insight into the further evolution of the separated flow. A remarkable conclusion is that an unseparated vorticity layer at the wall, a familiar feature in unsteady separation processes, disappears in finite time.

Description: An experimental study was conducted to investigate the generation process of random small-scale turbulence in an originally laminar mixing layer. The evolutions of the two types of deterministic structures, the spanwise and streamwise vortices, were first clarified in order to determine their roles in the transition process. A scaling rule for the streamwise distance from the trailing edge of the splitter plate to the vortex merging position was found for various velocity ratios. After this stremwise lengthscale was determined, it became clear that the spanwise wavelength of the streamwise vortices doubled after the merging of the spanwise structures which nominally doubled streamwise wavelengths. The most interesting finding was that the random small-scale eddies were produced by the interactions between the merging spanwise structures and the streamwise vortices.

Description: Although unsteady, high-Reynolds number, laminar boundary layers have conventionally been studied in terms of Eulerian coordinates, a Lagrangian approach may have significant analytical and computational advantages. In Lagrangian coordinates the classical boundary layer equations decouple into a momentum equation for the motion parallel to the boundary, and a hyperbolic continuity equation (essentially a conserved Jacobian) for the motion normal to the boundary. The momentum equations, plus the energy equation if the flow is compressible, can be solved independently of the continuity equation. Unsteady separation occurs when the continuity equation becomes singular as a result of touching characteristics, the condition for which can be expressed in terms of the solution of the momentum equations. The solutions to the momentum and energy equations remain regular. Asymptotic structures for a number of unsteady 3-D separating flows follow and depend on the symmetry properties of the flow. In the absence of any symmetry, the singularity structure just prior to separation is found to be quasi 2-D with a displacement thickness in the form of a crescent shaped ridge. Physically the singularities can be understood in terms of the behavior of a fluid element inside the boundary layer which contracts in a direction parallel to the boundary and expands normal to it, thus forcing the fluid above it to be ejected from the boundary layer.