Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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.

At least 685 records · Page 38

The measurement of winds over the ocean from Skylab with application to measuring and forecasting typhoons and hurricanes

On Skylab, a combination microwave radar-radiometer (S193) made measurements in a tropical hurricane (AVA), a tropical storm, and various extratropical wind systems. The winds at each cell scanned by the instrument were determined by objective numerical analysis techniques. The measured radar backscatter is compared to the analyzed winds and shown to provide an accurate method for measuring winds from space. An operational version of the instrument on an orbiting satellite will be able to provide the kind of measurements in tropical cyclones available today only by expensive and dangerous aircraft reconnaissance. Additionally, the specifications of the wind field in the tropical boundary layer should contribute to improved accuracy of tropical cyclone forecasts made with numerical weather predictions models currently being applied to the tropical atmosphere.

Cardone, V. J.↗

Seasonal-scale Observational Data Analysis and Atmospheric Phenomenology for the Cold Land Processes Experiment

The Cold Land Processes Experiment (CLPX) experiment emphasized the development of a strong synergism between process-oriented understanding, land surface models and microwave remote sensing. Our work sought to investigate which topographically- generated atmospheric phenomena are most relevant to the CLPX MSA's for the purpose of evaluating their climatic importance to net local moisture fluxes and snow transport through the use of high-resolution data assimilation/atmospheric numerical modeling techniques. Our task was to create three long-term, scientific quality atmospheric datasets for quantitative analysis (for all CLPX researchers) and provide a summary of the meteorologically-relevant phenomena of the three MSAs (see Figure) over northern Colorado. Our efforts required the ingest of a variety of CLPX datasets and the execution an atmospheric and land surface data assimilation system based on the Navier-Stokes equations (the Local Analysis and Prediction System, LAPS, and an atmospheric numerical weather prediction model, as required) at topographically- relevant grid spacing (approx. 500 m). The resulting dataset will be analyzed by the CLPX community as a part of their larger research goals to determine the relative influence of various atmospheric phenomena on processes relevant to CLPX scientific goals.

Poulos, Gregory S.↗

Global Infrared Observations of Roughness Induced Transition on the Space Shuttle Orbiter

High resolution infrared observations made from a mobile ground based optical system captured the laminar-to-turbulent boundary layer transition process as it occurred during Space Shuttle Endeavour's return to earth following its final mission in 2011. The STS-134 imagery was part of a larger effort to demonstrate an emerging and reliable non-intrusive global thermal measurement capability and to complement a series of boundary layer transition flight experiments that were flown on the Shuttle. The STS-134 observations are believed to be the first time that the development and movement of a hypersonic boundary layer transition front has been witnessed in flight over the entire vehicle surface and in particular, at unprecedented spatial resolution. Additionally, benchmark surface temperature maps of the Orbiter lower surface collected over multiple flights and spanning a Mach range of 18 to 6 are now available and represent an opportunity for collaborative comparison with computational techniques focused on hypersonic transition and turbulence modeling. The synergy of the global temperature maps with the companion in-situ thermocouple measurements serve as an example of the effective leveraging of resources to achieve a common goal of advancing our understanding of the complex nature of high Mach number transition. It is shown that quantitative imaging can open the door to a multitude of national and international opportunities for partnership associated with flight-testing and subsequent validation of numerical simulation techniques. The quantitative imaging applications highlighted in this paper offer unique and complementary flight measurement alternatives and suggest collaborative instrumentation opportunities to advance the state of the art in transition prediction and maximize the return on investment in terms of developmental flight tests for future vehicle designs.

Horvath, Thomas J.↗

A hybrid perturbation-Galerkin technique for partial differential equations

A two-step hybrid perturbation-Galerkin technique for improving the usefulness of perturbation solutions to partial differential equations which contain a parameter is presented and discussed. In the first step of the method, the leading terms in the asymptotic expansion(s) of the solution about one or more values of the perturbation parameter are obtained using standard perturbation methods. In the second step, the perturbation functions obtained in the first step are used as trial functions in a Bubnov-Galerkin approximation. This semi-analytical, semi-numerical hybrid technique appears to overcome some of the drawbacks of the perturbation and Galerkin methods when they are applied by themselves, while combining some of the good features of each. The technique is illustrated first by a simple example. It is then applied to the problem of determining the flow of a slightly compressible fluid past a circular cylinder and to the problem of determining the shape of a free surface due to a sink above the surface. Solutions obtained by the hybrid method are compared with other approximate solutions, and its possible application to certain problems associated with domain decomposition is discussed.

Geer, James F.↗

A new approach for synthesizing dual reflector antennas

It is pointed out that the offset geometry for dual reflector antennas is attractive for satellite communications where low side lobes are necessary to achieve good isolation between adjacent main lobes operating over the same frequency band. A description is presented of a new, numerically rigorous technique for synthesizing dual offset reflectors. The objective of the synthesis procedure is to transform a given primary feed pattern into a desired aperture distribution through two reflections in an offset dual reflector system. It is found that a given specification for the aperture distribution does not always lead to a smooth solution for the two surfaces. In this event, the so-called energy condition on the specified amplitude distribution is relaxed, and an approximate solution which satisfies the energy condition only within a certain numerical tolerance is derived.

Mittra, R.↗

Simulation and analysis of antennas radiating in a complex environment

A numerical procedure for computing the high-frequency radiation patterns of antennas mounted on curved surfaces is described. The procedure utilizes the uniform geometrical theory of diffraction to examine the antenna system's performance, which is dependent on antenna radiation patterns. Composite ellipsoid models of fuselage shapes are developed and the formation of geodesic paths on the models is studied; the shape of the fuselage affects the radiation patterns. The actual field radiated by the source and scattered by the structure is calculated using the ray field technique. The numerical solution is applied to the analysis of the antenna radiation patterns of a military aircraft, private aircraft, and the Space Shuttle orbiter. Good correlation between the calculated and measured radiation patterns is noted verifying the usefulness and accuracy of the numerical procedure.

Kim, J. J.↗

Toward Optimum Configuration of Circumferential Groove Casing Treatment in Transonic Compressor Rotors

The current paper first reviews experimental and numerical investigations to understand flow physics and to develop optimum configurations of circumferential grooves in compressor rotors. Circumferential grooves are used mainly to increase stall margin in axial compressors with small decrease in aerodynamic efficiency. Although circumferential groove casing treatment has been used widely, flow mechanisms of the circumferential grooves at near stall conditions are not well understood yet. Detailed time-dependent flow measurement inside tip gap in a high speed compressor is still a big challenge even though significant advance has been made in non-intrusive flow measurement technique. Therefore numerical approaches have been used to study relevant flow physics. However, optimum design of circumferential grooves to a given compressor with the computational tools is not practical yet. In the present paper, various investigations to study flow physics of circumferential groove casing treatment in axial compressor are reviewed first. Possible missing flow physics are identified and future research efforts for the optimum design are discussed.

Hah, Chunill↗

Trajectory optimization for an asymmetric launch vehicle

A numerical optimization technique is used to fully automate the trajectory design process for an symmetric configuration of the proposed Advanced Launch System (ALS). The objective of the ALS trajectory design process is the maximization of the vehicle mass when it reaches the desired orbit. The trajectories used were based on a simple shape that could be described by a small set of parameters. The use of a simple trajectory model can significantly reduce the computation time required for trajectory optimization. A predictive simulation was developed to determine the on-orbit mass given an initial vehicle state, wind information, and a set of trajectory parameters. This simulation utilizes an idealized control system to speed computation by increasing the integration time step. The conjugate gradient method is used for the numerical optimization of on-orbit mass. The method requires only the evaluation of the on-orbit mass function using the predictive simulation, and the gradient of the on-orbit mass function with respect to the trajectory parameters. The gradient is approximated with finite differencing. Prelaunch trajectory designs were carried out using the optimization procedure. The predictive simulation is used in flight to redesign the trajectory to account for trajectory deviations produced by off-nominal conditions, e.g., stronger than expected head winds.

Sullivan, Jeanne Marie↗

Cyclic Symmetry Finite Element Forced Response Analysis of a Distortion-Tolerant Fan with Boundary Layer Ingestion

Accurate prediction of the blade vibration stress is required to determine overall durability of fan blade design under Boundary Layer Ingestion (BLI) distorted flow environments. Traditional single blade modeling technique is incapable of representing accurate modeling for the entire rotor blade system subject to complex dynamic loading behaviors and vibrations in distorted flow conditions. A particular objective of our work was to develop a high-fidelity full-rotor aeromechanics analysis capability for a system subjected to a distorted inlet flow by applying cyclic symmetry finite element modeling methodology. This reduction modeling method allows computationally very efficient analysis using a small periodic section of the full rotor blade system. Experimental testing by the use of the 8-foot by 6-foot Supersonic Wind Tunnel Test facility at NASA Glenn Research Center was also carried out for the system designated as the Boundary Layer Ingesting Inlet/Distortion-Tolerant Fan (BLI2DTF) technology development. The results obtained from the present numerical modeling technique were evaluated with those of the wind tunnel experimental test, toward establishing a computationally efficient aeromechanics analysis modeling tool facilitating for analyses of the full rotor blade systems subjected to a distorted inlet flow conditions. Fairly good correlations were achieved hence our computational modeling techniques were fully demonstrated. The analysis result showed that the safety margin requirement set in the BLI2DTF fan blade design provided a sufficient margin with respect to the operating speed range.

forced response analysis↗

Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties

A semi-implicit fully spectral collocation method for the simulation of three-dimensional mantle convection with depth-dependent thermo-dynamic and transport properties is presented. The variable property Navier-Stokes equation expressed in terms of the primitive variable velocity and pressure is solved with the mass continuity and temperature equations. The periodic horizontal boundary conditions allow a Fourier expansion for the two horizontal directions. The stress-free, impermeable isothermal boundary conditions along with the depth dependent coefficients are handled with a Chebyshev expansion in the vertical direction. In the limit of an infinite Prandtl number appropriate to mantle convection, the inertial terms in the momentum equation are unimportant. In this case an explicit solution of a Poisson equation for pressure can be avoided; instead a fourth-order equation for vertical velocity can be solved. Simultaneous imposition of both impermeable and continuity boundary conditions during the vertical velocity evaluation is discussed. The pressure distributions on the top and bottom bounding planes were determined by means of an influence matrix technique. The numerical method employed here avoids time-splitting errors and enforces velocity boundary conditions and continuity over the entire domain, including the boundaries, to machine accuracy. Strongly time-dependent three-dimensional solutions up to a surface Rayleigh number of 1 x 10(exp 7) have been obtained. Strong upwellings, pulsating chaotically, are formed by the collective merging of cylindrical plumes.

Balachandar, S.↗

Flow and Turbulence Modeling and Computation of Shock Buffet Onset for Conventional and Supercritical Airfoils

Flow and turbulence models applied to the problem of shock buffet onset are studied. The accuracy of the interactive boundary layer and the thin-layer Navier-Stokes equations solved with recent upwind techniques using similar transport field equation turbulence models is assessed for standard steady test cases, including conditions having significant shock separation. The two methods are found to compare well in the shock buffet onset region of a supercritical airfoil that involves strong trailing-edge separation. A computational analysis using the interactive-boundary layer has revealed a Reynolds scaling effect in the shock buffet onset of the supercritical airfoil, which compares well with experiment. The methods are next applied to a conventional airfoil. Steady shock-separated computations of the conventional airfoil with the two methods compare well with experiment. Although the interactive boundary layer computations in the shock buffet region compare well with experiment for the conventional airfoil, the thin-layer Navier-Stokes computations do not. These findings are discussed in connection with possible mechanisms important in the onset of shock buffet and the constraints imposed by current numerical modeling techniques.

Bartels, Robert E.↗

A computerized symbolic integration technique for development of triangular and quadrilateral composite shallow-shell finite elements

Computerized symbolic integration was used in conjunction with group-theoretic techniques to obtain analytic expressions for the stiffness, geometric stiffness, consistent mass, and consistent load matrices of composite shallow shell structural elements. The elements are shear flexible and have variable curvature. A stiffness (displacement) formulation was used with the fundamental unknowns consisting of both the displacement and rotation components of the reference surface of the shell. The triangular elements have six and ten nodes; the quadrilateral elements have four and eight nodes and can have internal degrees of freedom associated with displacement modes which vanish along the edges of the element (bubble modes). The stiffness, geometric stiffness, consistent mass, and consistent load coefficients are expressed as linear combinations of integrals (over the element domain) whose integrands are products of shape functions and their derivatives. The evaluation of the elemental matrices is divided into two separate problems - determination of the coefficients in the linear combination and evaluation of the integrals. The integrals are performed symbolically by using the symbolic-and-algebraic-manipulation language MACSYMA. The efficiency of using symbolic integration in the element development is demonstrated by comparing the number of floating-point arithmetic operations required in this approach with those required by a commonly used numerical quadrature technique.

Anderson, C. M.↗