SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
We study the three dimensional instability and nonlinear growth of the two dimensional flow described by the stream function Psi = (A sin b1 x sin b2 y)/(b1(exp 2) + b2(exp 2)) where b1 = pi/L1, b2 = pi/L2. This is a swirling flow in a box which is bounded by 0 less than x less than L1, 0 less than y less than L2 and is infinite in the z direction. This flow is a solution of the Navier-Stokes equation with A = exp(-v(b1(exp 2) + b2(exp 2))t) which slowly decays. We seek a viscous solution which starts near this one and slips along but does not penetrate the bounding walls. The vorticity of the basic flow is w(sub z) = A sin b1 x sin b2 y which has maximum value A at the center of the box and drops to zero at the boundaries. We can think of the resulting flow as that of a captive vortex.
Stability of the elliptic vortex attracted interest in the past decade. Cambon (1982), and Cambon, Teissedre and Jeandel (1985) have studied the stability of such flows with spatially uniform velocity gradient, and have provided RDT solutions for a wide range of the parameter S/Omega (where the strain rate S and the vorticity 2 Omega define the velocity gradient matrix). The range studied included those of hyperbolic streamlines (strain dominated, S/Omega is greater than 1), linear streamlines (simple shear, S/Omega = 1), and elliptical streamlines (vorticity dominated, S/Omega is less than 1). The latter class has more recently attracted interest and several studies appeared (Pierrehumbert 1986, Bayly 1986, Craik and Criminale 1986). These studies will be collectively referred to as PBCC.
The linearized attenuation theory of NACA Technical Note 3375 is modified in the following manner: (a) an unsteady compressible local skin-friction coefficient is employed rather than the equivalent steady-flow incompressible coefficient; (b) a nonlinear approach is used to permit application of the theory to large attenuations; and (c) transition effects are considered. Curves are presented for predicting attenuation for a shock pressure ratio up to 20 and a range of shock-tube Reynolds numbers. Comparison of theory and experimental data for shock-wave strengths between 1.5 and 10 over a wide range of Reynolds numbers shows good agreement with the nonlinear theory evaluated for a transition Reynolds number of 2.5 X 10(exp 5).
When an alternating magnetic field interacts with induced eddy currents in a conducting body, there will be a repulsive force between the body and the driving coil system generating the field. This repulsive force is the basis of electromagnetic levitation, which allows containerless processing of different materials. The eddy currents in the conducting body also generate Joule heating. Axial rotation of electromagnetically levitated objects is a common observation in levitation systems and often an undesirable side effect of such experiments on 1-g and -g. There have been recent efforts to use magnetic damping and suppress this tendency of body rotation. The first report of rotation in EML drops was attributed to a slight asymmetry of the shape and location of the levitation coils could change the axis and speed of rotation. Other theories of sample rotation include a frequency difference in the traveling electromagnetic waves and a phase difference in two different applied fields of the same frequency. All of these different mechanisms share the following characteristics: the torque is small, constant for constant field strength, and very weakly dependent on the sample's temperature and phase (solid or liquid). During experiments on the MSL-1 (First Microgravity Science Laboratory) mission of the Space Shuttle (STS-83 and STS-94, April and July 1997), a droplet of palladium-silicon alloy was electromagnetically levitated for viscosity measurements. For the non-deforming droplet, the resultant MHD flow inside the drop is inferred from motion of impurities on the surface. These observations indicate formation of a pair of co-rotating toroidal flow structures inside the spheroidal levitated drop that undergo secondary flow instabilities. As rise in the fluid temperature rises, the viscosity falls and the internal flow accelerates and becomes oscillatory; and beyond a point in the experiments, the surface impurities exhibit non-coherent chaotic motion signifying emergence of turbulence inside the drop. In this work, a background of these set of observations will be given followed by a presentation of our results on the digital particle tracking analysis that has been performed on a number of available videos. The analysis indicates that the levitated drop attains a constant rotational speed during the melting phase and formation of the co-rotating axi-symmetric laminar toroidal structures. However, the rate of axial rotation increases dramatically during the deformation of the toroidal structures anti their breakup into chaotic entities. This new data suggests an interaction between the flow inside the levitated molten drop and the driving coils in the experiments. Possible mechanisms for this interaction will be reviewed. The data will also be used to make an assessment of existing theories on droplet rotation.
The objectives of this paper are two-fold. First, a numerical scheme for the simulation of a buoyant, reacting jet is presented with special attention given to boundary conditions. In the absence of coflow, a jet flame is particularly sensitive to boundary conditions enforced upon the computational domain. However, careful consideration of proper boundary conditions can minimize their effect upon the overall simulation. Second, results of some preliminary simulations are presented over a range of Froude and Damkohler numbers. This range was chosen so as to produce lifted flames in both normal gravity and microgravity environments.
Explore the source record for details and available documents.
No abstract available
Software assurance is the planned and systematic set of activities that ensures that software processes and products conform to requirements, standards, and procedures. Examples of such activities are the following: code inspections, unit tests, design reviews, performance analyses, construction of traceability matrices, etc. In practice, software development projects have only limited resources (e.g., schedule, budget, and availability of personnel) to cover the entire development effort, of which assurance is but a part. Projects must therefore select judiciously from among the possible assurance activities. At its heart, this can be viewed as an optimization problem; namely, to determine the allocation of limited resources (time, money, and personnel) to minimize risk or, alternatively, to minimize the resources needed to reduce risk to an acceptable level. The end result of the work reported here is a means to optimize quality-assurance processes used in developing software. This is achieved by combining two prior programs in an innovative manner
No abstract available
The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.
Optimal initial conditions for transient growth in a two-dimensional boundary layer flow correspond to stationary, counter-rotating vortices that subsequently develop into streamwise elongated streaks, which are characterized by an alternating pattern of low and high streamwise velocity. For incompressible flows, previous studies have shown that boundary layer modulation due to streaks below a threshold amplitude level can stabilize the Tollmien-Schlichting instability waves, resulting in a delay in the onset of laminar-turbulent transition. In the supersonic regime, the linearly, most-amplified waves become three-dimensional, corresponding to oblique, first-mode waves. This change in the character of dominant instabilities leads to an important change in the transition process, which is now dominated by oblique breakdown via nonlinear interactions between pairs of first-mode waves that propagate at equal but opposite angles with respect to the free stream. Because the oblique breakdown process is characterized by a rapid amplification of stationary streamwise streaks, artificial excitation of such streaks may be expected to promote transition in a supersonic boundary layer. Indeed, suppression of those streaks has been shown to delay the onset of transition in prior literature. Consistent with those findings, the present study shows that optimally growing stationary streaks indeed destabilize the first-mode waves, but only when the spanwise wavelength of the instability waves is equal to or smaller than twice the streak spacing. Transition in a benign disturbance environment typically involves first-mode waves with significantly longer spanwise wavelengths, and hence, these waves are stabilized by the optimal growth streaks. Thus, as long as the amplification factors for the destabilized, short wavelength instability waves remain below the threshold level for transition, a significant net stabilization is achieved, yielding a transition delay that is comparable to the length of the laminar region in the uncontrolled case.
Explore the source record for details and available documents.
The flow field over an airfoil contains several distinct flow characteristics, e.g. laminar, transitional, turbulent boundary layer flow, flow separation, unstable free shear layers, and a wake. This diversity of flow regimes taxes the presently available Reynolds averaged turbulence models. Such models are generally tuned to predict a particular flow regime, and adjustments are necessary for the prediction of a different flow regime. Similar difficulties are likely to emerge when the large eddy simulation technique is applied with the widely used Smagorinsky model. This model has not been successful in correctly representing different turbulent flow fields with a single universal constant and has an incorrect near-wall behavior. Germano et al. (1991) and Ghosal, Lund & Moin have developed a new subgrid-scale model, the dynamic model, which is very promising in alleviating many of the persistent inadequacies of the Smagorinsky model: the model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model has been remarkably successful in prediction of several turbulent and transitional flows. We plan to simulate turbulent flow over a '2D' airfoil using the large eddy simulation technique. Our primary objective is to assess the performance of the newly developed dynamic subgrid-scale model for computation of complex flows about aircraft components and to compare the results with those obtained using the Reynolds average approach and experiments. The present computation represents the first application of large eddy simulation to a flow of aeronautical interest and a key demonstration of the capabilities of the large eddy simulation technique.
Accurate, robust, and efficient prediction of transition in viscous flows is a significant challenge in computational fluid dynamics. We present a coupled, high-fidelity, iterative approach that leverages the FUN3D flow solver and the LASTRAC stability code to predict transition in low-disturbance environments, initiated by the linear growth of boundary-layer instability modes. Our method integrates the ability of FUN3D to compute mixed laminar-transitional-turbulent mean flows via transition-sensitized Reynolds-Averaged Navier-Stokes equations with the ability of LASTRAC to perform linear stability analysis, all within an automated framework that requires no intermediate user involvement. Unlike conventional frameworks that rely on classical stability theory or reduced-order metamodels, our approach employs the parabolized stability equations to provide more accurate and reliable estimates of disturbance growth for multiple instability mechanisms, including Tollmien-Schlichting, Kelvin-Helmholtz, and crossflow modes. By accounting for the effects of mean-flow nonparallelism as well as the surface curvature, this approach lays the foundation for improved N-factor correlations for transition onset prediction in a broad class of flows. We apply this method to three distinct flow configurations: 1) flow over a zero-pressure-gradient flat plate, 2) the NLF-0416 airfoil with both natural and separation-induced transitions, and 3) a 6:1 prolate spheroid, where transition is primarily driven by crossflow instability. For the two-dimensional cases, a formulated intermittency distribution is used to model the transition zone in between the laminar and fully turbulent flows. The results include comparisons with experimental measurements, similar numerical approaches, and transport-equations-based models, demonstrating good agreement in surface pressure coefficients, transition onset locations, and skin-friction coefficients for all three configurations. Besides contributing a couple of new insights into boundary-layer transition in these canonical cases, this study provides a powerful tool for transition modeling in both research and design applications in aerodynamics.
This paper describes recent improvements to the LAURA and HARA codes. LAURA is a CFD code for aerothermodynamics, and HARA evaluates the shock-layer radiation that provides the radiative source term for the flowfield energy equations and radiative heating to a surface. The next release of LAURA and HARA includes a variety of new capabilities. These new capabilities include an automated uncertainty quantification workflow for radiative heat transfer, options for specifying surface roughness and turbulent transition location in the algebraic turbulence models, and improved grid and solution interpolation techniques. Additionally, the computational efficiency of both LAURA and HARA have been improved. Optimization of the MPI communication routines in LAURA are shown to improve the parallel efficiency of the primary flow when running with multiple processes per block, and recent optimization of HARA leverage graphics processing unit (GPU) acceleration in the radiation calculations. Using GPU acceleration of HARA is shown to decrease the cost of the radiation line-of-sight calculation by approximately one order of magnitude for a 10.5 km/s Earth entry simulation.
The calculation of engineering flows undergoing laminar-turbulent transition presents special problems. Mean-flow quantities obey neither the fully laminar nor the fully turbulent correlations. In addition, local maxima in skin friction, wall temperature, and heat transfer often occur near the end of the transition region. Traditionally, modeling this region has been important for the design of turbine blades, where the transition region is long in relation to the chord length of the blade. More recently, the need for better transition-region models has been recognized by designers of hypersonic vehicles where the high Mach number, the low Reynolds number, and the low-disturbance flight environment emphasize the importance of the transition region. Needless to say, a model that might work well for the transitional flows typically found in gas turbines will not necessarily work well for the external surface of a hypersonic vehicle. In Section 2 of this report, some of the important flow features that control the transition region will be discussed. In Section 3, different approaches to the modeling problem will be summarized and cataloged. Fully turbulent flow models will be discussed in detail in Section 4; models specifically designed for transitional flow, in Section 5; and the evaluation of models, in Section 6.
In the near future, RANS computations will continue to play a significant role in the design of hypersonic vehicles with complex geometries. Therefore, it is imperative to continue testing, benchmarking, and refining the RANS models. In the present study, we evaluate RANS-like, transport equations-based models for predicting laminar-turbulent transition over a full-scale scale model of the BOLT flight configuration that was tested in the CUBRC LENS-II wind tunnel test facility. Based on the availability of the experimental results, comparisons are made between computations from several pre-existing transition models and Computational Fluid Dynamics (CFD) codes, with an emphasis on using the same computational meshes and flow conditions for all computations. The analysis covers the sensitivity of the transition predictions to the input parameters for five different transition models (four mainstream models designed for low-speed flows and a uniquely high-speed model), grid resolution, and the details of model implementation across three different flow solvers. The results show that the phenomenological models can describe significant aspects of the measured transition front. However, a number of additional improvements are required before these models can offer more reliable estimates of transition in high-speed flows.