Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximate projection method”

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 55 records · Page 3

Numerical Simulations of an Unsteady Rocket Launch from the AH-64D Apache Longbow Helicopter

Rocket and missile firings from attack helicopters can cause main engine compressor stall. Studies of this phenomenon suggest that the main engine ingests either the plume from the rockets or the rocket blast waves. This creates surges at the inlet face, causing a loss of power in the main engine. The objective of this project is to set-up a computational fluid dynamics (CFD) simulation of the AH-64D Apache Longbow helicopter during a rocket launch, in order to qualitatively study the fluid dynamics of the problem. This project presents a progression of three unsteady Navier-Stokes solutions. The first unsteady solution involves only a rocket launch from its launch canister. The second solution is a launch from a canister mounted on the Apache's wing-pylon assembly. The last solution includes the Apache main engine and fuselage. The computations use a series of structured, overset grid systems, which allow for a rocket moving in a prescribed path. The method implements a Roe upwind scheme with LU-SGS (lower-upper factored symmetric Gauss-Seidel). A rotor pressure disk model approximates the helicopter rotor, while the rocket engine exit properties are applied as a prescribed boundary condition. Although the project is only at the half-way point, the first and second CFD simulations suggest the possibility of pressure wave interference. Sudden surges in pressure occur from two sources: at rocket start-up, and as the rocket leaves the canister. Wave patterns set-up by these sources appear to propagate to the location of the engine inlet. However the simplified geometry simulation with the main engine needs to be performed before coming to a conclusion.

Okamoto, Kevin↗

A Runge-Kutta discontinuous finite element method for high speed flows

A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.

Bey, Kim S.↗

Implementation of the 64-meter-diameter Antennas at the Deep Space Stations in Australia and Spain

The management and construction aspects of the Overseas 64-m Antenna Project in which two 64-m antennas were constructed at the Tidbinbilla Deep Space Communications Complex in Australia, and at the Madrid Deep Space Communications Complex in Spain are described. With the completion of these antennas the Deep Space Network is equipped with three 64-m antennas spaced around the world to maintain continuous coverage of spacecraft operations. These antennas provide approximately a 7-db gain over the capabilities of the existing 26-m antenna nets. The report outlines the project organization and management, resource utilization, fabrication, quality assurance, and construction methods by which the project was successfully completed. Major problems and their solutions are described as well as recommendations for future projects.

Bartos, K. P.↗

Climate Hazard Assessment for Stakeholder Adaptation Planning in New York City

This paper describes a time-sensitive approach to climate change projections, developed as part of New York City's climate change adaptation process, that has provided decision support to stakeholders from 40 agencies, regional planning associations, and private companies. The approach optimizes production of projections given constraints faced by decision makers as they incorporate climate change into long-term planning and policy. New York City stakeholders, who are well-versed in risk management, helped pre-select the climate variables most likely to impact urban infrastructure, and requested a projection range rather than a single 'most likely' outcome. The climate projections approach is transferable to other regions and consistent with broader efforts to provide climate services, including impact, vulnerability, and adaptation information. The approach uses 16 Global Climate Models (GCMs) and three emissions scenarios to calculate monthly change factors based on 30-year average future time slices relative to a 30- year model baseline. Projecting these model mean changes onto observed station data for New York City yields dramatic changes in the frequency of extreme events such as coastal flooding and dangerous heat events. Based on these methods, the current 1-in-10 year coastal flood is projected to occur more than once every 3 years by the end of the century, and heat events are projected to approximately triple in frequency. These frequency changes are of sufficient magnitude to merit consideration in long-term adaptation planning, even though the precise changes in extreme event frequency are highly uncertain

Horton, Radley M.↗

On High-Order Upwind Methods for Advection

In the fourth installment of the celebrated series of five papers entitled "Towards the ultimate conservative difference scheme", Van Leer (1977) introduced five schemes for advection, the first three are piecewise linear, and the last two, piecewise parabolic. Among the five, scheme I, which is the least accurate, extends with relative ease to systems of equations in multiple dimensions. As a result, it became the most popular and is widely known as the MUSCL scheme (monotone upstream-centered schemes for conservation laws). Schemes III and V have the same accuracy, are the most accurate, and are closely related to current high-order methods. Scheme III uses a piecewise linear approximation that is discontinuous across cells, and can be considered as a precursor of the discontinuous Galerkin methods. Scheme V employs a piecewise quadratic approximation that is, as opposed to the case of scheme III, continuous across cells. This method is the basis for the on-going "active flux scheme" developed by Roe and collaborators. Here, schemes III and V are shown to be equivalent in the sense that they yield identical (reconstructed) solutions, provided the initial condition for scheme III is defined from that of scheme V in a manner dependent on the CFL number. This equivalence is counter intuitive since it is generally believed that piecewise linear and piecewise parabolic methods cannot produce the same solutions due to their different degrees of approximation. The finding also shows a key connection between the approaches of discontinuous and continuous polynomial approximations. In addition to the discussed equivalence, a framework using both projection and interpolation that extends schemes III and V into a single family of high-order schemes is introduced. For these high-order extensions, it is demonstrated via Fourier analysis that schemes with the same number of degrees of freedom 𝐾 per cell, in spite of the different piecewise polynomial degrees, share the same sets of eigenvalues and thus, have the same stability and accuracy. Moreover, these schemes are accurate to order 2𝐾−1, which is higher than the expected order of 𝐾.

high-order methods↗

Analysis of defect structure in silicon. Silicon sheet growth development for the large area silicon sheet task of the Low-Cost Solar array Project

One hundred ninety-three silicon sheet samples, approximately 880 square centimeters, were analyzed for twin boundary density, dislocation pit density, and grain boundary length. One hundred fifteen of these samples were manufactured by a heat exchanger method, thirty-eight by edge defined film fed growth, twenty-three by the silicon on ceramics process, and ten by the dendritic web process. Seven solar cells were also step-etched to determine the internal defect distribution on these samples. Procedures were developed or the quantitative characterization of structural defects such as dislocation pits, precipitates, twin & grain boundaries using a QTM 720 quantitative image analyzing system interfaced with a PDP 11/03 mini computer. Characterization of the grain boundary length per unit area for polycrystalline samples was done by using the intercept method on an Olympus HBM Microscope.

Natesh, R.↗

Performance of FFT methods in local gravity field modelling

Fast Fourier transform (FFT) methods provide a fast and efficient means of processing large amounts of gravity or geoid data in local gravity field modelling. The FFT methods, however, has a number of theoretical and practical limitations, especially the use of flat-earth approximation, and the requirements for gridded data. In spite of this the method often yields excellent results in practice when compared to other more rigorous (and computationally expensive) methods, such as least-squares collocation. The good performance of the FFT methods illustrate that the theoretical approximations are offset by the capability of taking into account more data in larger areas, especially important for geoid predictions. For best results good data gridding algorithms are essential. In practice truncated collocation approaches may be used. For large areas at high latitudes the gridding must be done using suitable map projections such as UTM, to avoid trivial errors caused by the meridian convergence. The FFT methods are compared to ground truth data in New Mexico (xi, eta from delta g), Scandinavia (N from delta g, the geoid fits to 15 cm over 2000 km), and areas of the Atlantic (delta g from satellite altimetry using Wiener filtering). In all cases the FFT methods yields results comparable or superior to other methods.

Forsberg, Rene↗

Evaluating Light Rain Drop Size Estimates from Multiwavelength Micropulse Lidar Network Profiling

This paper investigates multiwavelength retrievals of median equivolumetric drop diameter D(sub 0) suitable for drizzle and light rain, through collocated 355-/527-nm Micropulse Lidar Network (MPLNET) observations collected during precipitation occurring 9 May 2012 at the Goddard Space Flight Center (GSFC) project site. By applying a previously developed retrieval technique for infrared bands, the method exploits the differential backscatter by liquid water at 355 and 527 nm for water drops larger than approximately 50 micrometers. In the absence of molecular and aerosol scattering and neglecting any transmission losses, the ratio of the backscattering profiles at the two wavelengths (355 and 527 nm), measured from light rain below the cloud melting layer, can be described as a color ratio, which is directly related to D(sub 0). The uncertainty associated with this method is related to the unknown shape of the drop size spectrum and to the measurement error. Molecular and aerosol scattering contributions and relative transmission losses due to the various atmospheric constituents should be evaluated to derive D(sub 0) from the observed color ratio profiles. This process is responsible for increasing the uncertainty in the retrieval. Multiple scattering, especially for UV lidar, is another source of error, but it exhibits lower overall uncertainty with respect to other identified error sources. It is found that the total error upper limit on D(sub 0) approaches 50%. The impact of this retrieval for long-term MPLNET monitoring and its global data archive is discussed.

light-rain drop size estimates↗

Design of robust state feedback controllers via eigenvalue/eigenvector assignment

A new method of selecting a multivariable state feedback controller is presented. The resulting controller simultaneously realizes arbitrary closed-loop eigenvalues, approximates specified modal behavior and achieves low eigensystem sensitivity with respect to plant parameter variations. The method characterizes a vector space slightly larger than the null space of the closed-loop system matrix and projects a desired modal matrix onto this space. Sensitivity of eigenvalues and eigenvectors is then minimized local to the desired modal matrix using a gradient search technique. A tutorial example to illustrate the design procedure is given.

Mielke, R. R.↗

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗

Collisional redistribution of radiation. I - The density matrix

As a first step in deriving accurate, calculable expressions for collisional redistribution of radiation, valid outside the impact regime a method for calculating an accurate binary-collision density operator for an atom in the presence of a driving field and perturbers is presented. Projection-operator techniques to establish first the precise validity of the Markoff approximation for radiative relaxation are used. Using the same techniques, it is shown how the collisional-relaxation problem may be analyzed in a practical manner outside the Markoff approximation. The effect of correlations between radiative and collisional events is included in a consistent way and the physical implications of these correlation effects are demonstrated for a simple example, the two-level atom with nondegenerate levels.

Burnett, K.↗

Positronium formation in e+ plus H- collisions

Cross sections for positronium formation by capture from the negative hydrogen ion are given. Orthogonalization corrections to the Coulomb (First) Born Approximation (CBA) differential and total cross sections are calculated using approximate H- wave functions of both Lowdin and Chandrasekhar. Various methods of orthogonalizing the unbound projectile to the possible bound states are considered. It is found that treating the atomic nuclei as if they were isotopic spin projections of a single type of nucleon gives cross sections that are an improvement over the CBA.

Straton, Jack C.↗

Krylov subspace methods - Theory, algorithms, and applications

Projection methods based on Krylov subspaces for solving various types of scientific problems are reviewed. The main idea of this class of methods when applied to a linear system Ax = b, is to generate in some manner an approximate solution to the original problem from the so-called Krylov subspace span. Thus, the original problem of size N is approximated by one of dimension m, typically much smaller than N. Krylov subspace methods have been very successful in solving linear systems and eigenvalue problems and are now becoming popular for solving nonlinear equations. The main ideas in Krylov subspace methods are shown and their use in solving linear systems, eigenvalue problems, parabolic partial differential equations, Liapunov matrix equations, and nonlinear system of equations are discussed.

Sad, Youcef↗

Photographic Volume Estimation of CPAS Main Parachutes

Capsule Parachute Assembly System (CPAS) flight tests regularly stage a helicopter to observe inflation of 116 ft D o ringsail Main parachutes. These side views can be used to generate 3-D models of inflating canopies to estimate enclosed volume. Assuming a surface of revolution is inadequate because reefed canopies in a cluster are elongated due to mutual aerodynamic interference. A method was developed to combine the side views with upward looking HD video to account for non-circular cross sections. Approximating the cross sections as elliptical greatly improves accuracy. But since that correction requires manually tracing projected outlines, the actual irregular shapes can be used to generate high fidelity models. Compensation is also made for apparent tilt angle. Validation was accomplished by comparing perimeter and projected area with known line lengths and/or high quality photogrammetry.

Ray, Eric S.↗

New approaches to virtual environment surgery

This research focused on two main problems: 1) low cost, high fidelity stereoscopic imaging of complex tissues and organs; and 2) virtual cutting of tissue. A further objective was to develop these images and virtual tissue cutting methods for use in a telemedicine project that would connect remote sites using the Next Generation Internet. For goal one we used a CT scan of a human heart, a desktop PC with an OpenGL graphics accelerator card, and LCD stereoscopic glasses. Use of multiresolution meshes ranging from approximately 1,000,000 to 20,000 polygons speeded interactive rendering rates enormously while retaining general topography of the dataset. For goal two, we used a CT scan of an infant skull with premature closure of the right coronal suture, a Silicon Graphics Onyx workstation, a Fakespace Immersive WorkBench and CrystalEyes LCD glasses. The high fidelity mesh of the skull was reduced from one million to 50,000 polygons. The cut path was automatically calculated as the shortest distance along the mesh between a small number of hand selected vertices. The region outlined by the cut path was then separated from the skull and translated/rotated to assume a new position. The results indicate that widespread high fidelity imaging in virtual environment is possible using ordinary PC capabilities if appropriate mesh reduction methods are employed. The software cutting tool is applicable to heart and other organs for surgery planning, for training surgeons in a virtual environment, and for telemedicine purposes.

NASA Discipline Neuroscience↗

The construction and use of divergence free vector expansions for incompressible fluid flow calculations

For incompressible fluids the law of mass conservation reduces to a constraint on the velocity vector, namely that it be divergence free. This constraint has long been a source of great difficulty to the numericist seeking to discretize the Navier-Stokes and Euler equations. A spectral method is discussed which overcomes this difficulty. Its efficacy is demonstrated on some simple problems. The velocity is approximated by a finite sum of divergence free vectors, each of which satisfies the same boundary conditions as the velocity. Projecting the governing equation onto the space of inviscid vector fields eliminates the pressure term and produces a set of ordinary differential equations that must be solved for the coefficents in the velocity. The pressure can then be recovered if it is needed.

Mhuiris, N. M. G.↗

Proposed methods for the reconstruction of the Pioneer/Venus probe entry trajectories

In connection with the Pioneer Venus Project, an orbiter and an array of atmospheric probes are sent to Venus in separate launches during 1978. The four probes, which are transported to the planet on an instrumented Bus are released approximately 20 days prior to encounters. There are three identical small probes and a large probe. Each is designed to make in-situ atmospheric measurements. One of the principal scientific objectives of the Pioneer Venus (P/V) multiprobe mission is the determination of the structure and composition of the Venusian atmosphere down to the planet surface. To enhance the science return it is of paramount importance to be able to accurately determine the spacecraft entry trajectories. In addition to the in-situ spacecraft measurements, earth-based Doppler data are available. The process felt to be most applicable for the P/V multiprobe entry trajectory reconstruction is considered, with emphasis placed on the highly dynamic upper entry segment of the trajectory.

Findlay, J. T.↗

Cruise noise of the 2/9th scale model of the Large-scale Advanced Propfan (LAP) propeller, SR-7A

Noise data on the Large-scale Advanced Propfan (LAP) propeller model SR-7A were taken in the NASA Lewis Research Center 8 x 6 foot Wind Tunnel. The maximum blade passing tone noise first rises with increasing helical tip Mach number to a peak level, then remains the same or decreases from its peak level when going to higher helical tip Mach numbers. This trend was observed for operation at both constant advance ratio and approximately equal thrust. This noise reduction or, leveling out at high helical tip Mach numbers, points to the use of higher propeller tip speeds as a possible method to limit airplane cabin noise while maintaining high flight speed and efficiency. Projections of the tunnel model data are made to the full scale LAP propeller mounted on the test bed aircraft and compared with predictions. The prediction method is found to be somewhat conservative in that it slightly overpredicts the projected model data at the peak.

Dittmar, James H.↗