Two-dimensional Spectral Analysis of SeaWiFS and Topex/Poseidon Data Supports Theory of Wave-Induced Variations of Passive Tracers
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Advanced formulations of boundary element method for periodic, transient transform domain and transient time domain solution of three-dimensional solids have been implemented using a family of isoparametric boundary elements. The necessary numerical integration techniques as well as the various solution algorithms are described. The developed analysis has been incorporated in a fully general purpose computer program BEST3D which can handle up to 10 subregions. A number of numerical examples are presented to demonstrate the accuracy of the dynamic analyses.
The role of turbulence in pollutant formation was investigated both experimentally and analytically with the objective of understanding the important mechanisms of a fairly basic level. Preliminary measurements of the cold-flow turbulence characteristics downstream of three perforated plate flame holders, of varying solidity ratio, were made. In the combustion experiment, measurements of the concentration distributions of CO, NO(x), HC, and CO2, downstream of the same flame holders, were made. The combustor configuration was tubular and constant-area, and the fuel was premixed, prevaporized propane. The experimental results were modeled by a one-dimensional analysis in which transport phenomena were neglected, and more successfully, by a quasi-one-dimensional analysis containing recirculation zones and turbulent transport effects.
Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic-plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intension of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels. In this regard, a series of fracture tests were conducted on curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. During the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).
The behavior of the eigenvalues of preconditioning matrices for the pseudospectral approximation to the derivative operator has been analyzed in one and two dimensions. The one-dimensional analysis resulted in real and positive eigenvalues for the selected tridiagonal matrices. In the two-dimensional analysis, the eigenvalues of the selected block-diagonal matrices behaved well, but the preconditioner is full and therefore not suitable for applications. The Richardson scheme has been applied in the unpreconditioned as well as the preconditioned version to find the solution of the model problem.
This study presents a numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed in the present study. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field of a subsonic turbine stage and the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid.
In this work, we present a data-driven strategy for effective construction of a surrogate model in high-dimensional parameter space for the ion energy-angle distribution (IEAD) output of hPIC simulations of plasma-surface interactions. The methodology is based on a bin-by-bin least-squares fitting of the IEAD in the parameter space. The fitting is performed in a transformed coordinate system to normalize the IEAD, and it employs sparse grids for sampling the parameter space to overcome sampling challenges in high dimensions. The surrogate model is significantly cheaper computationally than direct hPIC simulations yet maintains high fidelity to them, providing a fast emulator for hPIC simulations. Sensitivity analysis based on the surrogate model is utilized to characterize the dependence of the ion impact angle and energy moments on the physical parameters.
Analysis of a flexible beam with displaced end-located inertial masses is presented. The resulting three-dimensional mode shape is shown to consist of two one-plane bending modes and one torsional mode. These three components of the mode shapes are shown to be linear combinations of trigonometric and hyperbolic sine and cosine functions. Boundary conditions are derived to obtain nonlinear algebraic equations through kinematic coupling of the general solutions of the three governing partial differential equations. A method of solution which takes these boundary conditions into account is also presented. A computer program has been written to obtain unique solutions to the resulting nonlinear algebraic equations. This program, which calculates natural frequencies and three-dimensional mode shapes for any number of modes, is presented and discussed.
The sizing and efficiency of an aircraft is largely determined by the performance of its high-lift system. Subsonic civil transports most often use deployable multi-element airfoils to achieve the maximum-lift requirements for landing, as well as the high lift-to-drag ratios for take-off. However, these systems produce very complex flow fields which are not fully understood by the scientific community. In order to compete in today's market place, aircraft manufacturers will have to design better high-lift systems. Therefore, a more thorough understanding of the flows associated with these systems is desired. Flight and wind-tunnel experiments have been conducted on NASA Langley's B737-100 research aircraft to obtain detailed full-scale flow measurements on a multi-element high-lift system at various flight conditions. As part of this effort, computational aerodynamic tools are being used to provide preliminary flow-field information for instrumentation development, and to provide additional insight during the data analysis and interpretation process. The purpose of this paper is to demonstrate the ability and usefulness of a three-dimensional low-order potential flow solver, PMARC, by comparing computational results with data obtained from 1/8 scale wind-tunnel tests. Overall, correlation of experimental and computational data reveals that the panel method is able to predict reasonably well the pressures of the aircraft's multi-element wing at several spanwise stations. PMARC's versatility and usefulness is also demonstrated by accurately predicting inviscid three-dimensional flow features for several intricate geometrical regions.
The Department of Energy, the Stirling Technology Company (STC), and the NASA Glenn Research Center are developing Stirling convertors for Stirling radioisotope generators to provide electrical power for future NASA deep space missions. STC is developing the 55-We technology demonstration convertor (TDC) under contract to the Department of Energy. The Department of Energy recently named Lockheed Martin as the system integration contractor for the Stirling radioisotope generator development project. Lockheed Martin will develop the Stirling radioisotope generator engineering unit and has contract options to develop the qualification unit and the first flight unit. Glenn s role includes an in-house project to provide convertor, component, and materials testing and evaluation in support of the overall power system development. As a part of this work, Glenn has established an in-house Stirling research laboratory for testing, analyzing, and evaluating Stirling machines. STC has built four 55-We convertors for NASA, and these are being tested at Glenn. A cross-sectional view of the 55-We TDC is shown in the figure. Of critical importance to the successful development of the Stirling convertor for space power applications is the development of a lightweight and highly efficient linear alternator. In support, Glenn has been developing finite element analysis and finite element method tools for performing various linear alternator thermal and electromagnetic analyses and evaluating design configurations. A three-dimensional magnetostatic finite element model of STC's 55-We TDC linear alternator was developed to evaluate the demagnetization fields affecting the alternator magnets. Since the actual linear alternator hardware is symmetric to the quarter section about the axis of motion, only a quarter section of the alternator was modeled. The components modeled included the mover laminations, the neodymium-iron-boron magnets, the stator laminations, and the copper coils. The three-dimensional magnetostatic model was then coupled with a circuit simulator model of the alternator load and convertor controller. The coupled model was then used to generate alternator terminal voltage and current predictions. The predicted voltage and current waveforms agreed well with the experimental data, which tended to validate the accuracy of the coupled model. The model was then used to generate predictions of the demagnetization fields acting on the alternator magnets for the alternator under load. The preliminary model predictions indicate that the highest potential for demagnetization is along the inside surface of the uncovered magnets. The demagnetization field for the uncovered magnets when the mover is positioned at the end of a stroke is higher than it is when the mover is at the position of maximum induced voltage or maximum alternator current. Assuming normal load conditions, the model predicted that the onset of demagnetization is most likely to occur for magnet temperatures above 101 C.
A computer code has been developed to analyze the inviscid flow field in a supersonic combustion ramjet (scramjet) inlet. The code uses the three-dimensional Euler equations in full conservation form to describe the inlet flow. An algebraic numerical coordinate transformation is used to generate a set of boundary-fitted curvilinear coordinates. The governing equations are solved by a time-asymptotic, unsplit, two-step, finite-difference method. This method is highly efficient on the vector processing computers for which the current code is written. Detailed results are presented for two scramjet inlet configurations over a range of Mach numbers. The calculated results are compared with the available experimental and theoretical results.
Recent interest in environmental/thermal barrier coatings (EBC/TBCs) has prompted research to develop life-prediction methodologies for the coating systems of advanced high-temperature ceramic matrix composites (CMCs). Heat-transfer analysis of EBC/TBCs for CMCs is an essential part of the effort. It helps establish the resulting thermal profile through the thickness of the CMC that is protected by the EBC/TBC system. This report documents the results of a one-dimensional analysis of an advanced high-temperature CMC system protected with an EBC/TBC system. The one-dimensional analysis was used for tradeoff studies involving parametric variation of the conductivity; the thickness of the EBC/TBCs, bond coat, and CMC substrate; and the cooling requirements. The insight gained from the results will be used to configure a viable EBC/TBC system for CMC liners that meet the desired hot surface, cold surface, and substrate temperature requirements.
The importance of natural convection and other fluid motions in low-g space processing is now well recognized. Recent space experiments in the areas of natural convection and material processing, as well as results of theoretical studies, have yielded much needed information on fluid behavior in low-g environments. The state of knowledge of fluid motions in low-g environments is reviewed and the dimensional analysis approach used to assess the relative importances of various driving forces for fluid flow in four of the Skylab material processing experiments outlined. Results of dimensional analyses for the Skylab experiments, subsequently confirmed by actual space data, are presented. Finally, the limits of dimensional analysis in assessment studies are indicated.
In the course of previous composite structures test programs, the need for and the feasibility of developing analyses for scale-up effects has been demonstrated. The analysis techniques for scale-up effects fall into two categories. The first category pertains to developing analysis methods independently for a single, unique failure mode in composites, and using this compendium of analysis methods together with a global structural model to identify and predict the response and failure mode of full-scale built-up structures. The second category of scale-up effects pertains to similitude in structural validation testing. In this latter category, dimensional analysis is used to develop scale-up laws that enable extrapolation of sub-scale component test data to full-scale structures. This viewgraph presentation describes the approach taken and some developments accomplished in the first category of analysis for scale-up effects. Layup dependence of composite material properties severely limits the use of the dimensional analysis approach and these limitations are illustrated by examples.
To overcome the limitations of classical propeller theory, a computer program, NASPROP-E, was developed which solves for the flow field surrounding a multibladed propeller and axisymmetric nacelle combination using a finite difference method. The governing equations are the three dimensional unsteady Euler equations written in a cylindrical coordinate system. They are marched in time until a steady state solution is obtained. The Euler equations require no special treatment to model the blade work vorticity. The equations are solved using an implicit approximate factorization method. Numerical results are presented which have greatly increased the understanding of high speed propeller flow fields. Numerical results for swirl angle downstream of the propeller and propeller power coefficient are higher than experimental results. The radial variation of coefficient are higher than experimental results. The radial variation of swirl angle, however, is in reasonable agreement with the experimental results. The predicted variation of power coefficient with blade angle agrees very well with data.
We perform Bayesian model calibration of two-nucleon (NN) low-energy constants (LECs) appearing in an NN interaction based on pionless effective field theory (πEFT). The calibration is carried out for potentials constructed using naive dimensional analysis in NN relative momenta (p) up to next-to-leading order [NLO, O(p 2 )] and next-to-next-to-next-to-leading order [N3LO, O(p 4 )]. We consider two classes of pionless πEFT potential: one that acts in all partial waves and another that is dominated by s-wave physics. The two classes produce broadly similar results for calibrations to NN data up to E lab = 5 MeV. Our analysis accounts for the correlated uncertainties that arise from the truncation of the pionless πEFT. We simultaneously estimate both the πEFT LECs and the parameters that quantify the truncation error. This permits the first quantitative estimates of the pionless πEFT breakdown scale, Λ b : the 95% intervals are Λ b ∈[50.11,63.03] MeV at NLO and Λ b ∈[72.27,88.54] MeV at N3LO. Furthermore, invoking naive dimensional analysis for the NN potential, therefore, does not lead to consistent results across orders in pionless πEFT. This exemplifies the possible use of Bayesian tools to identify inconsistencies in a proposed EFT power counting.
This paper presents representative results from an aeroelastic code (TURBO-AE) based on an Euler/Navier-Stokes unsteady aerodynamic code (TURBO). Unsteady pressure, lift, and moment distributions are presented for a helical fan test configuration which is used to verify the code by comparison to two-dimensional linear potential (flat plate) theory. The results are for pitching and plunging motions over a range of phase angles, Good agreement with linear theory is seen for all phase angles except those near acoustic resonances. The agreement is better for pitching motions than for plunging motions. The reason for this difference is not understood at present. Numerical checks have been performed to ensure that solutions are independent of time step, converged to periodicity, and linearly dependent on amplitude of blade motion. The paper concludes with an evaluation of the current state of development of the TURBO-AE code and presents some plans for further development and validation of the TURBO-AE code.
Two currently available coupled inviscid/viscous multielement computational codes, including a relatively simple panel method and an Euler method, are used to analyze a high-lift system. The results are compared with two-dimensional wind-tunnel test results and then with the three-dimensional flight-test results obtained from the NASA Langley Transport Systems Research Vehicle five-element high-lift wing section. Comparisons were also made between the panel method, the Euler method, and flight data for two high-lift configurations, one representing a take-off configuration and the other an approach configuration. For the take-off configuration, both codes agreed reasonably well with experimental data, but both codes were found to overpredict the flap upper-surface pressures for the approach configuration.