Showing papers by "Langley Research Center published in 2000"

Theory of the lattice boltzmann method: dispersion, dissipation, isotropy, galilean invariance, and stability

Abstract: The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Electrohydrodynamic Flow Control with a Glow-Discharge Surface Plasma

Abstract: Multiplee ow diagnostics have been applied to planarpanels covered by strips of glow-discharge surface plasma in atmospheric pressure air generated by the one atmosphere uniform glow discharge plasma. Direct drag measurements, smoke wire and titanium tetrachloride e ow visualization, and boundary-layer velocity proe les were obtained. The plasma generated along streamwise-oriented, symmetric strip electrodes is shown to cause a large increase in drag, whereas the plasma along spanwise-oriented, asymmetric strip electrodes can generate a signie cant thrust. Flow visualization and mean velocity measurements show the primary cause of the phenomena to be a combination of mass transport and vortical structures induced by strong electrohydrodynamic body forces on the e ow, known as paraelectric forcing.

Low-cost piezocomposite actuator for structural control applications

Abstract: The design, manufacture, and testing of a low-cost, flexible, planar composite piezoceramic actuator device will be presented. The actuator uses interdigitated electrodes for poling and subsequent actuation of an internal layer of machined piezoceramic fibers. The fiber sheets are formed from monolithic piezoceramic wafers and conventional computer controlled wafer-dicing methods. The fabrication and use of fiber sheets allows precise handing and alignment of piezoceramic fibers during subsequent phases of actuator assembly. Test show that the actuator is capable of producing large, directional in-plane strains; on order of 2000 parts-per-million under a 4000 V peak-to-peak applied voltage cycle. Preliminary endurance testing indicates that the device is relatively durable, with no reductions in free-strain performance up to 90 million electrical cycles.

Distribution and fate of selected oxygenated organic species in the troposphere and lower stratosphere over the Atlantic

Abstract: A large number of oxygenated organic chemicals (peroxyacyl nitrates, alkyl nitrates, acetone, formaldehyde, methanol, methylhydroperoxide, acetic acid and formic acid) were measured during the 1997 Subsonic Assessment (SASS) Ozone and Nitrogen Oxide Experiment (SONEX) airborne field campaign over the Atlantic. In this paper, we present a first picture of the distribution of these oxygenated organic chemicals (Ox-organic) in the troposphere and the lower stratosphere, and assess their source and sink relationships. In both the troposphere and the lower stratosphere, the total atmospheric abundance of these oxygenated species (ΣOx-organic) nearly equals that of total nonmethane hydrocarbons (ΣNMHC), which have been traditionally measured. A sizable fraction of the reactive nitrogen (10–30%) is present in its oxygenated organic form. The organic reactive nitrogen fraction is dominated by peroxyacetyl nitrate (PAN), with alkyl nitrates and peroxypropionyl nitrate (PPN) accounting for <5% of total NOy. Comparison of observations with the predictions of the Harvard three-dimensional global model suggests that in many key areas (e.g., formaldehyde and peroxides) substantial differences between measurements and theory are present and must be resolved. In the case of CH3OH, there appears to be a large mismatch between atmospheric concentrations and estimated sources, indicating the presence of major unknown removal processes. Instrument intercomparisons as well as disagreements between observations and model predictions are used to identify needed improvements in key areas. The atmospheric chemistry and sources of this group of chemicals is poorly understood even though their fate is intricately linked with upper tropospheric NOx and HOx cycles.

The NASA Tetrahedral Unstructured Software System (TETRUSS)

Abstract: The NASA Tetrahedral Unstructured Software System (TetrUSS) was developed during the 1990's to provide a rapid aerodynamic analysis and design capability to applied aerodynamicists. The system is comprised of loosely integrated, user-friendly software that enables the application of advanced Euler and Navier-Stokes tetrahedral finite volume technology to complex aerodynamic problems. TetrUSS has matured well because of the generous feedback from many willing users representing a broad cross-section of background and skill levels. This paper presents an overview of the current capabilities of the TetrUSS system along with some representative results from selected applications.

A closed-loop system for examining psychophysiological measures for adaptive task allocation.

Abstract: A closed-loop system was evaluated for its efficacy in using psychophysiological indexes to moderate workload. Participants were asked to perform either 1 or 3 tasks from the Multiattribute Task Battery and complete the NASA Task Load Index after each trial. An electroencephalogram (EEG) was sampled continuously while they performed the tasks, and an EEG index (beta/alpha plus theta) was derived. The system made allocation decisions as a function of the level of operator engagement based on the value of the EEG index. The results of the study demonstrated that it was possible to moderate an operator's level of engagement through a closed-loop system driven by the operator's own EEG. In addition, the system had a significant impact on behavioral, subjective, and psychophysiological correlates of workload as task load increased. The theoretical and practical implications of these results for adaptive automation are discussed.

A case-study application of RTCA DO-254: design assurance guidance for airborne electronic hardware

Abstract: In a joint project with the FAA, NASA Langley is developing a hardware design in accordance with RTCA DO-254: Design Assurance Guidance for Airborne Electronic Hardware. The purpose of the case study is to gain understanding of the new guidance document and generate an example suitable for use in training. For the case study, we have selected a core subsystem of the Scalable Processor-Independent Design for Electromagnetic Resilience (SPIDER). SPIDER is a new fault-tolerant architecture under development at NASA Langley Research Center.

The GCM–Reality Intercomparison Project for SPARC (GRIPS): Scientific Issues and Initial Results

Abstract: To investigate the effects of the middle atmosphere on climate, the World Climate Research Programme is supporting the project "Stratospheric Processes and their Role in Climate" (SPARC). A central theme of SPARC, to examine model simulations of the coupled troposphere—middle atmosphere system, is being performed through the initiative called GRIPS (GCM—Reality Intercomparison Project for SPARC). In this paper, an overview of the objectives of GRIPS is given. Initial activities include an assessment of the performance of middle atmosphere climate models, and preliminary results from this evaluation are presented here. It is shown that although all 13 models evaluated represent most major features of the mean atmospheric state, there are deficiencies in the magnitude and location of the features, which cannot easily be traced to the formulation (resolution or the parameterizations included) of the models. Most models show a cold bias in all locations, apart from the tropical tropopause region wher...

Photochemistry of HOx in the upper troposphere at northern midlatitudes

Abstract: The factors controlling the concentrations of HOx radicals (= OH + peroxy) in the upper troposphere (8–12 km) are examined using concurrent aircraft observations of OH, HO2, H2O2, CH3OOH, and CH2O made during the Subsonic Assessment Ozone and Nitrogen Oxide Experiment (SONEX) at northern midlatitudes in the fall. These observations, complemented by concurrent measurements of O3, H2O, NO, peroxyacetyl nitrate (PAN), HNO3, CH4, CO, acetone, hydrocarbons, actinic fluxes, and aerosols, allow a highly constrained mass balance analysis of HOx and of the larger chemical family HOy (= HOx + 2 H2O2 + 2 CH3OOH + HNO2 + HNO4). Observations of OH and HO2 are successfully simulated to within 40% by a diel steady state model constrained with observed H2O2 and CH3OOH. The model captures 85% of the observed HOx variance, which is driven mainly by the concentrations of NOx (= NO + NO2) and by the strength of the HOx primary sources. Exceptions to the good agreement between modeled and observed HOx are at sunrise and sunset, where the model is too low by factors of 2–5, and inside cirrus clouds, where the model is too high by factors of 1.2–2. Heterogeneous conversion of NO2 to HONO on aerosols (γNO2 = 10−3) during the night followed by photolysis of HONO could explain part of the discrepancy at sunrise. Heterogeneous loss of HO2 on ice crystals (γice_HO2 = 0.025) could explain the discrepancy in cirrus. Primary sources of HOx from O(1D)+H2O and acetone photolysis were of comparable magnitude during SONEX. The dominant sinks of HOy were OH+HO2 (NOx 50 pptv). Observed H2O2 concentrations are reproduced by model calculations to within 50% if one allows in the model for heterogeneous conversion of HO2 to H2O2 on aerosols (γHO2 = 0.2). Observed CH3OOH concentrations are underestimated by a factor of 2 on average. Observed CH2O concentrations were usually below the 50 pptv detection limit, consistent with model results; however, frequent occurrences of high values in the observations (up to 350 pptv) are not captured by the model. These high values are correlated with high CH3OH and with cirrus clouds. Heterogeneous oxidation of CH3OH to CH2O on aerosols or ice crystals might provide an explanation (γice_CH3OH ∼ 0.01 would be needed).

Software Fault Tolerance: A Tutorial

Abstract: Because of our present inability to produce error-free software, software fault tolerance is and will continue to be an important consideration in software systems. The root cause of software design errors is the complexity of the systems. Compounding the problems in building correct software is the difficulty in assessing the correctness of software for highly complex systems. After a brief overview of the software development processes, we note how hard-to-detect design faults are likely to be introduced during development and how software faults tend to be state-dependent and activated by particular input sequences. Although component reliability is an important quality measure for system level analysis, software reliability is hard to characterize and the use of post-verification reliability estimates remains a controversial issue. For some applications software safety is more important than reliability, and fault tolerance techniques used in those applications are aimed at preventing catastrophes. Single version software fault tolerance techniques discussed include system structuring and closure, atomic actions, inline fault detection, exception handling, and others. Multiversion techniques are based on the assumption that software built differently should fail differently and thus, if one of the redundant versions fails, it is expected that at least one of the other versions will provide an acceptable output. Recovery blocks, N-version programming, and other multiversion techniques are reviewed.

A Shell/3D Modeling Technique for the Analysis of Delaminated Composite Laminates

Abstract: A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a shell finite element model. Multi-point constraints provided a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with shell finite elements. Double Cantilever Beam, End Notched Flexure, and Single Leg Bending specimens were analyzed first using full 3D finite element models to obtain reference solutions. Mixed mode strain energy release rate distributions were computed using the virtual crack closure technique. The analyses were repeated using the shell/3D technique to study the feasibility for pure mode I, mode II and mixed mode I/II cases. Specimens with a unidirectional layup and with a multidirectional layup were simulated. For a local 3D model, extending to a minimum of about three specimen thicknesses on either side of the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures the shell/3D modeling technique offers a great potential for reducing the model size, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.

Report on the Loss of the Mars Polar Lander and Deep Space 2 Missions

Abstract: NASA's Mars Surveyor Program (MSP) began in 1994 with plans to send spacecraft to Mars every 26 months. Mars Global Surveyor (MGS), a global mapping mission, was launched in 1996 and is currently orbiting Mars. Mars Surveyor '98 consisted of Mars Climate Orbiter (MCO) and Mars Polar Lander (MPL). Lockheed Martin Astronautics (LMA) was the prime contractor for Mars Surveyor '98. The Jet Propulsion Laboratory (JPL), California Institute of Technology, manages the Mars Surveyor Program for NASA's Office of Space Science. MPL was developed under very tight funding constraints. The combined development cost of MPL and MCO, including the cost of the two launch vehicles, was approximately the same as the development cost of the Mars Pathfinder mission, including the cost of its single launch vehicle. The MPL project accepted the challenge to develop effective implementation methodologies consistent with programmatic requirements.

Structures Technology for Future Aerospace Systems

Abstract: An overview of structures technology for future aerospace systems is given. Discussion focuses on developments in component technologies that will improve the vehicle performance, advance the technology exploitation process, and reduce system life-cycle costs. The component technologies described are smart materials and structures, multifunctional materials and structures, affordable composite structures, extreme environment structures, flexible load bearing structures, and computational methods and simulation-based design. The trends in each of the component technologies are discussed and the applicability of these technologies to future aerospace vehicles is described.

Recent Turbulence Model Advances Applied to Multi-Element Airfoil Computations

Abstract: A one-equation linear turbulence model and two-equation nonlinear explicit algebraic stress model (EASM) are applied to the flow over a multielement airfoil. The effect of the K-epsilon and K-omega forms of the two-equation model are explored, and the K-epsilon form is shown to be deficient in the wall-bounded regions of adverse pressure gradient flows. A new K-omega form of EASM is introduced. Nonlinear terms present in EASM are shown to improve predictions of turbulent shear stress behind the trailing edge of the main element and near midflap. Curvature corrections are applied to both the one- and two-equation turbulence models and yield only relatively small local differences in the flap region, where the flow field undergoes the greatest curvature. Predictions of maximum lifts are essentially unaffected by the turbulence model variations studied.

Transonic Flutter Suppression Control Law Design and Wind-Tunnel Test Results

Abstract: Fluttersuppressioncontrollawdesignandwind-tunneltestresultsintransonice owforaNACA0012wingmodel, under the benchmark active control technology program at NASA Langley Research Center, will be presented. Two control law design processes using classical and minimax techniques are described. Design considerations for improving the multivariable system robustness are outlined. The classical control law was digitally implemented and tested in the NASA Langley Transonic Dynamics Tunnel. In wind-tunnel tests in air and heavy gas medium, the closed-loop e utter dynamic pressure was increased by over 50% up to the wind-tunnel upper limit. The active e utter suppression system also provided signie cant robustness, relative to gain and phase perturbations, even in the presence of transonic shocks and e ow separation.

Bilevel Integrated System Synthesis with Response Surfaces

Abstract: Two variants of the BLISS method are proposed for multidisciplinary design optimization (MDO). These variants utilize experimental design methods and response surfaces to reduce the number of expensive system analysis required for solution of the MDO problem. BLISS is an MDO method for decomposition-based optimization of engineering systems that involves system optimization with a relatively small number of design variables and a number of subsystem optimizations that could each have a large number of local variables. In BLISS, the optimum sensitivity analysis data are used to relate the subsystem optimization solutions with the system optimizations. Instead, with the proposed variants, polynomial response surface approximations using either the system analysis or the subsystem optimization results are used. Additionally, the response surface construction process is well suited for computing in a concurrent processing environment. The proposed variants are implemented and evaluated on a conceptual-level aircraft and ship design problem.

Photogrammetry Applied to Wind Tunnel Testing

Abstract: In image-based measurements, quantitative image data must be mapped to three-dimensional object space. Analytical photogrammetric methods, which may be used to accomplish this task, are discussed from the viewpoint of experimental fluid dynamicists. The Direct Linear Transformation (DLT) for camera calibration, used in pressure sensitive paint, is summarized. An optimization method for camera calibration is developed that can be used to determine the camera calibration parameters, including those describing lens distortion, from a single image. Combined with the DLT method, this method allows a rapid and comprehensive in-situ camera calibration and therefore is particularly useful for quantitative flow visualization and other measurements such as model attitude and deformation in production wind tunnels. The paper also includes a brief description of typical photogrammetric applications to temperature- and pressure-sensitive paint measurements and model deformation measurements in wind tunnels.

Turbulence Model Predictions of Strongly Curved Flow in a U-Duct

Abstract: The ability of three types of turbulence models to accurately predict the effects of curvature on the flow in a U-duct is studied An explicit algebraic stress model performs slightly better than one- or two-equation linear eddy viscosity models, although it is necessary to fully account for the variation of the production-to-dissipation-rate ratio in the algebraic stress model formulation In their original formulations, none of these turbulence models fully captures the suppressed turbulence near the convex wall, whereas a full Reynolds stress model does Some of the underlying assumptions used in the development of algebraic stress models are investigated and compared with the computed flowfield from the full Reynolds stress model Through this analysis, the assumption of Reynolds stress anisotropy equilibrium used in the algebraic stress model formulation is found to be incorrect in regions of strong curvature By the accounting for the local variation of the principal axes of the strain rate tensor, the explicit algebraic stress model correctly predicts the suppressed turbulence in the outer part of the boundary layer near the convex wall

First-Order Model Management With Variable-Fidelity Physics Applied to Multi-Element Airfoil Optimization

Abstract: First-order approximation and model management is a methodology for a systematic use of variable-fidelity models or approximations in optimization. The intent of model management is to attain convergence to high-fidelity solutions with minimal expense in high-fidelity computations. The savings is terms of computationally intensive evaluations depends on the ability of the available lower-fidelity model or a suite of models to predict the improvement trends for the high-fidelity problem. Variable-fidelity models can be represented by data-fitting approximations, variable-resolution models, variable-convergence models, or variable physical fidelity models. The present work considers the use of variable-fidelity physics models. We demonstrate the performance of model management on an aerodynamic optimization of a multi-element airfoil designed to operate in the transonic regime. Reynolds-averaged Navier-Stokes equations represent the high-fidelity model, while the Euler equations represents the low-fidelity model. An unstructured mesh-based analysis code FUN2D evaluates functions and sensitivity derivatives for both models. Model management for the present demonstration problem yields fivefold savings in terms of high-fidelity evaluations compared to optimization done with high-fidelity computations alone.

Sea ice remote sensing using surface reflected GPS signals

Abstract: This paper describes a new research effort to extend the application of Global Positioning System (GPS) signal reflections, received by airborne instruments, to cryospheric remote sensing. Their experimental results indicate that reflected GPS signals have potential to provide information on the presence and condition of sea and fresh-water ice as well as the freeze/thaw state of frozen ground. They show results from aircraft experiments over the ice pack near Barrow, Alaska indicating correlation between forward-scattered GPS returns and RADARSAT backscattered measurements.

Wake Vortex Transport and Decay in Ground Effect: Vortex Linking with the Ground

Abstract: Numerical simulations are carried out with a three-dimensional Large-Eddy Simulation (LES) model to explore the sensitivity of vortex decay and transport in ground effect (IGE). The vortex decay rates are found to be strongly enhanced following maximum descent into ground effect. The nondimensional decay rate is found to be insensitive to the initial values of circulation, height, and vortex separation. The information gained from these simulations is used to construct a simple decay relationship. This relationship compares well with observed data from an IGE case study. Similarly, a relationship for lateral drift due to ground effect is constructed from the LES data. In the second part of this paper, vortex linking with the ground is investigated. Our numerical simulations of wake vortices for IGE show that a vortex may link with its image beneath the ground, if the intensity of the ambient turbulence is moderate to high. This linking with the ground (which is observed in real cases)gives the appearance of a vortex tube that bends to become vertically oriented and which terminates at the ground. From the simulations conducted, the linking time for vortices in the free atmosphere; i.e., a function of ambient turbulence intensity.