Showing papers by "Langley Research Center published in 1995"

Atmospheric effects of the Mt Pinatubo eruption

Abstract: The eruption of Mt Pinatubo in June 1991 caused the largest perturbation this century to the participate content of the stratosphere. The radiative influence of the injected particles put an end to several years of globally warm surface temperatures. At the same time, the combined effect of volcanic particles and anthropogenic reactive chlorine has led to record low levels of stratospheric ozone.

Absorption of solar radiation by clouds: observations versus models.

Abstract: There has been a long history of unexplained anomalous absorption of solar radiation by clouds. Collocated satellite and surface measurements of solar radiation at five geographically diverse locations showed significant solar absorption by clouds, resulting in about 25 watts per square meter more global-mean absorption by the cloudy atmosphere than predicted by theoretical models. It has often been suggested that tropospheric aerosols could increase cloud absorption. But these aerosols are temporally and spatially heterogeneous, whereas the observed cloud absorption is remarkably invariant with respect to season and location. Although its physical cause is unknown, enhanced cloud absorption substantially alters our understanding of the atmosphere's energy budget.

Robust attitude stabilization of spacecraft using nonlinear quaternion feedback

Abstract: This paper considers the problem of three-axis attitude stabilization of a rigid spacecraft. A nonlinear control law which uses the feedback of the unit quaternion and the measured angular velocities is proposed and is shown to provide global asymptotic stability. The control law does not require the knowledge of the system parameters and is, therefore, robust to modeling errors. The significance of the control law is that it can be used for large-angle maneuvers with guaranteed stability. >

HZETRN: Description of a Free-Space Ion and Nucleon Transport and Shielding Computer Program

Abstract: The high-charge-and-energy (HZE) transport computer program HZETRN is developed to address the problems of free-space radiation transport and shielding. The HZETRN program is intended specifically for the design engineer who is interested in obtaining fast and accurate dosimetric information for the design and construction of space modules and devices. The program is based on a one-dimensional space-marching formulation of the Boltzmann transport equation with a straight-ahead approximation. The effect of the long-range Coulomb force and electron interaction is treated as a continuous slowing-down process. Atomic (electronic) stopping power coefficients with energies above a few A MeV are calculated by using Bethe''s theory including Bragg''s rule, Ziegler''s shell corrections, and effective charge. Nuclear absorption cross sections are obtained from fits to quantum calculations and total cross sections are obtained with a Ramsauer formalism. Nuclear fragmentation cross sections are calculated with a semiempirical abrasion-ablation fragmentation model. The relation of the final computer code to the Boltzmann equation is discussed in the context of simplifying assumptions. A detailed description of the flow of the computer code, input requirements, sample output, and compatibility requirements for non-VAX platforms are provided.

A multispectrum nonlinear least squares fitting technique

Abstract: An extension of the nonlinear least squares spectrum fitting technique has been developed for the simultaneous fitting of multiple spectra. This procedure's smaller number of fitted parameters as compared to fitting one spectrum at a time improves the determination of spectroscopic parameters. A more reliable evaluation of the errors associated with the solution is possible. Correlations among fitted parameters may preclude their determination from a single spectrum fit. If the correlations differ from spectrum to spectrum, however, separation of these parameters is often possible when including the spectra in one solution. Overfilling is also avoided when combining spectra with sufficiently different experimental physical conditions.

Development and Application of the Collaborative Optimization Architecture in a Multidisciplinary Design Environment

Abstract: Collaborative optimization is a design architecture applicable in any multidisciplinary analysis environment but specifically intended for large-scale distributed analysis applications. In this approach, a complex problem is hierarchically decomposed along disciplinary boundaries into a number of subproblems which are brought into multidisciplinary agreement by a system-level coordination process. When applied to problems in a multidisciplinary design environment, this scheme has several advantages over traditional solution strategies. These advantageous features include reducing the amount of information transferred between disciplines, the removal of large iteration-loops, allowing the use of different subspace optimizers among the various analysis groups, an analysis framework which is easily parallelized and can operate on heterogenous equipment, and a structural framework that is well-suited for conventional disciplinary organizations. In this article, the collaborative architecture is developed and its mathematical foundation is presented. An example application is also presented which highlights the potential of this method for use in large-scale design applications.

A predictor-corrector technique for visualizing unsteady flow

Abstract: Presents a method for visualizing unsteady flow by displaying its vortices. The vortices are identified by using a vorticity-predictor pressure-corrector scheme that follows vortex cores. The cross-sections of a vortex at each point along the core can be represented by a Fourier series. A vortex can be faithfully reconstructed from the series as a simple quadrilateral mesh, or its reconstruction can be enhanced to indicate helical motion. The mesh can reduce the representation of the flow features by a factor of 1000 or more compared with the volumetric dataset. With this amount of reduction, it is possible to implement an interactive system on a graphics workstation to permit a viewer to examine, in 3D, the evolution of the vortical structures in a complex, unsteady flow. >

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: a study of the boundary error

Abstract: The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: 1) impose the exact boundary condition only at the end of the complete RK cycle, 2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases , results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Multidisciplinary Design Optimization: An Emerging New Engineering Discipline

Abstract: This paper attempts to define the Multidisciplinary Design Optimization (MDO) as a new field of research endeavor and as an aid in the design of engineering systems. It examines the MDO conceptual components in relation to each other and defines their functions.

Fully-explicit and self-consistent algebraic reynolds stress models

Abstract: A fully-explicit, self-consistent algebraic expression for the Reynolds stress, which is the exact solution to the Reynolds stress transport equation in the `weak equilibrium'' limit for two-dimensional mean flows for all linear and some quasi-linear pressure-strain models, is derived. Current explicit algebraic Reynolds stress models derived by employing the `weak equilibrium'' assumption treat the production-to-dissipation (P varepsilon) ratio implicitly, resulting in an effective viscosity that can be singular away from the equilibrium limit. In the present paper, the set of simultaneous algebraic Reynolds stress equations are solved in the full non-linear form and the eddy viscosity is found to be non-singular. Preliminary tests indicate that the model performs adequately, even for three dimensional mean flow cases. Due to the explicit and non-singular nature of the effective viscosity, this model should mitigate many of the difficulties encountered in computing complex turbulent flows with the algebraic Reynolds stress models.

First global WCRP shortwave surface radiation budget dataset

Abstract: Shortwave radiative fluxes that reach the earth's surface are key factors that influence atmospheric and oceanic circulations as well as surface climate. Yet, information on these fluxes is meager. Surface site data are generally available from only a limited number of observing stations over land. Much less is known about the large-scale variability of the shortwave radiative fluxes over the oceans, which cover most of the globe. Recognizing the need to produce global-scale fields of such fluxes for use in climate research, the World Climate Research Program has initiated activities that led to the establishment of the Surface Radiation Budget Climatology Project with the ultimate goal to determine various components of the surface radiation budget from satellite data. In this paper, the first global products that resulted from this activity are described. Monthly and daily data on a 280-km grid scale are available. Samples of climate parameters obtainable from the dataset are presented. Emphasis is given to validation and limitations of the results. For most of the globe, satellite estimates have bias values between +/- 20 W/sq m and root mean square (rms) values are around 25 W/sq m. There are specific regions with much larger uncertainties however.

A unified multigrid solver for the navier-stokes equations on mixed element meshes

Abstract: A unified multigrid solution technique is presented for solving the Euler and Reynolds-averaged Navier-Stokes equations on unstructured meshes using mixed elements consisting of triangles and quadrilaterals in two dimensions, and of hexahedra, pyramids, prisms and tetrahedra in three dimensions. While the use of mixed elements is by no means a novel idea, the contribution of the paper lies in the formulation of a complete solution technique which can handle structured grids, block structured grids, and unstructured grids of tetrahedra or mixed elements without any modification. This is achieved by discretizing the full Navier-Stokes equations on tetrahedral elements, and the thin layer version of these equations on other types of elements, while using a single edge-based data-structure to construct the discretization over all element types. An agglomeration multigrid algorithm, which naturally handles meshes of any types of elements, is employed to accelerate convergence. An automatic algorithm which reduces the complexity of a given triangular or tetrahedral mesh by merging candidate triangular or tetrahedral elements into quadrilateral or prismatic elements is also described. The gains in computational efficiency afforded by the use of non-simplicial meshes over fully tetrahedral meshes are demonstrated through several examples.

Summary of an Active Flexible Wing program

Abstract: This paper presents a summary of the NASA/ Rockwell Active Flexible Wing program. Major elements of the program are presented. Key program accomplishmerits included singleand multiple-mode flutter suppression, load alleviation and load control during rapid roll maneuvers, and multi-input/multi-output multiple-function active controls tests above the open-loop flutter boundary.

Wing flutter boundary prediction using unsteady Euler aerodynamic method

Abstract: Modifications to an existing three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1) for the aeroelastic analysis of wings are described. These modifications, which were previously added to CFL3D Version 1.0, include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time-integration with the government flow equations. The paper gives a brief description of these modifications and presents unsteady calculations which check the modifications to the code. Euler flutter results for an isolated 45 degree swept-back wing are compared with experimental data for seven freestream Mach numbers which define the flutter boundary over a range of Mach number from 0.499 to 1.14. These comparisons show good agreement in flutter characteristics for freestream Mach numbers below unity. For freestream Mach numbers above unity, the computed aeroelastic results predict a premature rise in the flutter boundary as compared with the experimental boundary. Steady and unsteady contours of surface Mach number and pressure are included to illustrate the basic flow characteristics of the time-marching flutter calculations and to aid in identifying possible causes for the premature rise in the computational flutter boundary.

A perspective on unstructured grid flow solvers

Abstract: This survey paper assesses the status of compressible Euler and Navier-Stokes solvers on unstructured grids. Different spatial and temporal discretization options for steady and unsteady flows are discussed. The integration of these components into an overall framework to solve practical problems is addressed. Issues such as grid adaptation, higher order methods, hybrid discretizations and parallel computing are briefly discussed. Finally, some outstanding issues and future research directions are presented.

THE CORRELATED k-DISTRIBUTION TECHNIQUE AS APPLIED TO THE AVHRR CHANNELS

Abstract: Correlated k-distributions have been created to account for the molecular absorption found in the spectral ranges of the five Advanced Very High Resolution Radiometer (AVHRR) satellite channels. The production of the k-distributions was based upon an exponential-sum fitting of transmissions (ESFT) technique which was applied to reference line-by-line absorptance calculations. To account for the overlap of spectral features from different molecular species, the present routines made use of the multiplication transmissivity property which allows for considerable flexibility, especially when altering relative mixing ratios of the various molecular species. To determine the accuracy of the correlated k-distribution technique as compared to the line-by-line procedure, atmospheric flux and heating rate calculations were run for a wide variety of atmospheric conditions. For the atmospheric conditions taken into consideration, the correlated k-distribution technique has yielded results within about 0.5% for both the cases where the satellite spectral response functions were applied and where they were not. The correlated k-distribution's principal advantages is that it can be incorporated directly into multiple scattering routines that consider scattering as well as absorption by clouds and aerosol particles.

Long-term total solar irradiance variability during sunspot cycle 22

Abstract: Total solar irradiance measurements from the 1984-1993 Earth Radiation Budget Satellite (ERBS) active cavity radiometer and 1978-1993 Nimbus 7 transfer cavity radiometer spacecraft experiments are analyzed to detect the presence of 11-, 22-, and 80-year irradiance variability components. The analyses confirmed the existence of a significant 11-year irradiance variability component, associated with solar magnetic activity and the sunspot cycle. The analyses also suggest the presence of a 22- or 80-year variability component. The earlier Nimbus 7 and Solar Maximum Mission (SMM) spacecraft irradiance measurements decreased approximately 1.2 and 1.3 W/sq m, respectively, between 1980 and 1986. The Nimbus 7 values increased 1.2 W/sq m between 1986 and 1989. The ERBS irradiance measurements increased 1.3 W/sq m during 1986-1989, and then decreased 0.4 W/sq m (at an annual rate of 0.14 W/sq. m/yr) during 1990-1993. Considering the correlations between ERBS, Nimbus 7, and SMM irradiance trends and solar magnetic activity, the total solar irradiance should decrease to minimum levels by 1997 as solar activity decreases to minimum levels, and then increase to maximum levels by the year 2000 as solar activity rises. The ERBS measurements yielded 165.4 +/- 0.7 W/sq m as the mean irradiance value with measurement accuracies and precisions of 0.2% and 0.02%, respectively. The ERBS mean irradiance value is within 0.2% of the 1367.4, 1365.9, and 1366.9 W/sq m mean values for the SMM, Upper Atmosphere Research Satellite (UARS), and Space Shuttle Atmospheric Laboratory for Applications and Science (ATLAS 1) Solar Constant (SOLCON) active cavity radiometer spacecraft experiments, respectively. The Nimbus 7 measurements yielded 1372.1 W/sq m as the mean value with a measurement accuracy of 0.5%. Empirical irradiance model fits, based upon 10.7 -cm solar radio flux (F10) and photometric sunspot index (PSI), were used to assess the quality of the ERBS, Numbus 7, SMM, and the UARS irradiance data sets and to identify irradiance variability trends which may be caused by drifts or shifts in the spacecraft sensor responses. Comparisons among the fits and measured irradiances indicate that the Nimbus 7 radiometer response shifted by a total of 0.8 W/sq m between September 1989 and April 1990 and that the ERBS and UARS radiometers each drifted approximately 0.5 W/sq m during the first 5 months in orbit.

Experimental Investigation of Inlet-Combustor Isolators for a Dual-Mode Scramjet at a Mach Number of 4

Abstract: This report details experimentally derived operational characteristics of numerous two-dimensional planar inlet-combustor isolator configurations at a Mach number of 4. Variations in geometry included (1) inlet cowl length; (2) inlet cowl rotation angle; (3) isolator length; and (4) utilization of a rearward-facing isolator step. To obtain inlet-isolator maximum pressure-rise data relevant to ramjet-engine combustion operation, configurations were mechanically back pressured. Results demonstrated that the combined inlet-isolator maximum back-pressure capability increases as a function of isolator length and contraction ratio, and that the initiation of unstart is nearly independent of inlet cowl length, inlet cowl contraction ratio, and mass capture. Additionally, data are presented quantifying the initiation of inlet unstarts and the corresponding unstart pressure levels.

A coupled analysis method for structures with independently modelled finite element subdomains

Abstract: A new method for analysing plate and shell structures with two or more independently modelled finite element subdomains is presented, assessed, and demonstrated. This method provides a means of coupling local and global finite element models whose nodes do not coincide along their common interface. In general, the method provides a means of coupling structural components (e.g., wing and fuselage) which may have been modelled by different analysts. In both cases, the need for transition modelling, which is often tedious and complicated, is eliminated. The coupling is accomplished through an interface for which three formulations are considered and presented. These formulations are: collocation, discrete least-squares, and hybrid variational. Several benchmark problems are analysed and it is shown that the hybrid variational formulation provides the most accurate solutions.

Issues in Space Radiation Protection: Galactic Cosmic Rays

Abstract: When shielding from cosmic heavy ions, one is faced with limited knowledge about the physical properties and biological responses of these radiations. Herein, the current health is discussed in terms of conventional protection practice and a test biological response model. The impact of biological response on optimum materials selection for cosmic ray shielding is presented in terms of the transmission characteristics of the shield material. Although liquid hydrogen is an optimum shield material, evaluation of the effectiveness of polymeric structural materials must await improvement in our knowledge of both the biological response and the nuclear processes.

Spectral methods on arbitrary grids

Abstract: Stable and spectrally accurate numerical methods are constructed on arbitrary grids for partial differential equations. These new methods are equivalent to conventional spectral methods but do not rely on specific grid distributions. Specifically, we show how to implement Legendre Galerkin, Legendre collocation, and Laguerre Galerkin methodology on arbitrary grids.

ICASE/LaRC Workshop on Benchmark Problems in Computational Aeroacoustics (CAA)

Abstract: The proceedings of the Benchmark Problems in Computational Aeroacoustics Workshop held at NASA Langley Research Center are the subject of this report. The purpose of the Workshop was to assess the utility of a number of numerical schemes in the context of the unusual requirements of aeroacoustical calculations. The schemes were assessed from the viewpoint of dispersion and dissipation -- issues important to long time integration and long distance propagation in aeroacoustics. Also investigated were the effect of implementation of different boundary conditions. The Workshop included a forum in which practical engineering problems related to computational aeroacoustics were discussed. This discussion took the form of a dialogue between an industrial panel and the workshop participants and was an effort to suggest the direction of evolution of this field in the context of current engineering needs.

LIDAR multiple scattering from clouds

Abstract: Multiple-scattering LIDAR return calculations obtained by seven different models for the same specified numerical experiment are compared. This work results from an international joint effort stimulated by the workshop group called MUSCLE for MUltiple SCattering Lidar Experiments. The models include approximations to the radiative-transfer theory, Monte-Carlo calculations, a stochastic model of the process of multiple scattering, and an extension of Mie theory for particles illuminated by direct and scattered light. The model solutions are similar in form but differ by up to a factor of 5 in the strength of the multiple-scattering contributions. Various reasons for the observed differences are explored and their practical significance is discussed.

Spatial direct numerical simulation of high-speed boundary-layer flows part I: Algorithmic considerations and validation

Abstract: A highly accurate algorithm for the direct numerical simulation (DNS) of spatially evolving high-speed boundary-layer flows is described in detail and is carefully validated. To represent the evolution of instability waves faithfully, the fully explicit scheme relies on non-dissipative high-order compact-difference and spectral collocation methods. Several physical, mathematical, and practical issues relevant to the simulation of high-speed transitional flows are discussed. In particular, careful attention is paid to the implementation of inflow, outflow, and far-field boundary conditions. Four validation cases are presented, in which comparisons are made between DNS results and results obtained from either compressible linear stability theory or from the parabolized stability equation (PSE) method, the latter of which is valid for nonparallel flows and moderately nonlinear disturbance amplitudes. The first three test cases consider the propagation of two-dimensional second-mode disturbances in Mach 4.5 flat-plate boundary-layer flows. The final test case considers the evolution of a pair of oblique second-mode disturbances in a Mach 6.8 flow along a sharp cone. The agreement between the fundamentally different PSE and DNS approaches is remarkable for the test cases presented.

Mars Pathfinder six-degree-of-freedom entry analysis

Abstract: The Mars Pathfinder mission provides the next opportunity for scientific exploration of the surface of Mars following a 7.6 km/s direct entry. In support of this effort, a six-degree-of-freedom trajectory analysis and aerodynamic characteristic assessment are performed to demonstrate vehicle flyability and to quantify the effect that each of numerous uncertainties has upon the nominal mission profile. The entry vehicle is shown to be aerodynamically stable over a large portion of its atmospheric flight. Two low angle-of-attack static instabilities (freestream velocities of about 6.5 and 3.5 km/s) and a low angle-of-attack dynamic instability (supersonic) are identified and shown to cause bounded increases in vehicle attitude. The effects of center-of-gravity placement, entry attitude, vehicle roll rate, aerodynamic misprediction, and atmospheric uncertainty on the vehicle attitude profile and parachute deployment conditions are quantified. A Monte Carlo analysis is performed to statistically assess the combined impact of multiple off-nominal conditions on the nominal flight characteristics. These results suggest that there is a 99.7% probability that the peak attitude throughout the entry will be less than 8.5 deg, the peak heating attitude will be below 6.2 deg, and the attitude at parachute deployment will be less than 3.9 deg.