Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Nonlinear structural analysis on distributed-memory computers

A computational strategy is presented for the nonlinear static and postbuckling analyses of large complex structures on massively parallel computers. The strategy is designed for distributed-memory, message-passing parallel computer systems. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a nested dissection (or multilevel substructuring) ordering scheme; (3) parallel assembly of global matrices; and (4) a parallel sparse equation solver. The effectiveness of the strategy is assessed by applying it to thermo-mechanical postbuckling analyses of stiffened composite panels with cutouts, and nonlinear large-deflection analyses of HSCT models on Intel Paragon XP/S computers. The numerical studies presented demonstrate the advantages of nested dissection-based solvers over traditional skyline-based solvers on distributed memory machines.

Watson, Brian C.↗
NASA NTRS1995

Postbuckling and large-deflection nonlinear analyses on distributed-memory computers

A computational strategy is presented for postbuckling and nonlinear static analyses of large complex structures on distributed-memory parallel computers. The strategy is designed for message-passing parallel computer systems. The key elements of the proposed strategy are: (1) a multiple-parameter reduced basis technique; (2) a nested dissection (or multilevel substructuring) ordering scheme; (3) parallel assembly of global matrices; and (4) a parallel sparse equation solver. The effectiveness of the strategy is assessed by performing thermomechanical postbuckling analyses of stiffened composite panels with cutouts, and nonlinear large-deflection analyses of High Speed Civil Transport models on three distributed-memory computers. The numerical studies presented demonstrate the advantages of nested dissection-based solvers over traditional skyline-based solvers on distributed-memory machines.

Watson, Brian C.↗
NASA NTRS1994

Transient nucleate pool boiling in microgravity: Some initial results

Variable gravity provides an opportunity to test the understanding of phenomena which are considered to depend on buoyancy, such as nucleate pool boiling. The active fundamental research in nucleate boiling has sought to determine the mechanisms or physical processes responsible for its high effectiveness, manifested by the high heat flux levels possible with relatively low temperature differences. Earlier research on nucleate pool boiling at high gravity levels under steady conditions demonstrated quantitatively that the heat transfer is degraded as the buoyancy normal to the heater surfaced increases. Correspondingly, it was later shown, qualitatively for short periods of time only, that nucleate boiling heat transfer is enhanced as the buoyancy normal to the heater surface is reduced. It can be deduced that nucleate pool boiling can be sustained as a quasi-steady process provided that some means is available to remove the vapor generated from the immediate vicinity of the heater surface. One of the objectives of the research, the initial results of which are presented here, is to quantify the heat transfer associated with boiling in microgravity. Some quantitative results of nucleate pool boiling in high quality microgravity (a/g approximately 10(exp -5)) of 5s duration, obtained in an evacuated drop tower, are presented here. These experiments were conducted as precursors of longer term space experiments. A transient heating technique is used, in which the heater surface is a transparent gold film sputtered on a qua rtz substrate, simultaneously providing the mean surface temperature from resistance thermometry and viewing of the boiling process both from beneath and across the surface. The measurement of the transient mean heater surface temperature permits the computation, by numerical means, of the transient mean heat transfer coefficient. The preliminary data obtained demonstrates that a quasi-steady boiling process can occur in microgravity if the bulk liquid subcooling is sufficiently high and if the imposed heat flux is sufficiently low. This is attributed to suface tension effects at the liquid-vapor-solid junction causing rewetting to take place, sustaining the nucleate boiling. Otherwise, dryout at the heater surface will occur, as observed.

Merte, Herman, Jr.↗
NASA NTRS1995

An Analytic Approximation to Very High Specific Impulse and Specific Power Interplanetary Space Mission Analysis

A simple, analytic approximation is derived to calculate trip time and performance for propulsion systems of very high specific impulse (50,000 to 200,000 seconds) and very high specific power (10 to 1000 kW/kg) for human interplanetary space missions. The approach assumed field-free space, constant thrust/constant specific power, and near straight line (radial) trajectories between the planets. Closed form, one dimensional equations of motion for two-burn rendezvous and four-burn round trip missions are derived as a function of specific impulse, specific power, and propellant mass ratio. The equations are coupled to an optimizing parameter that maximizes performance and minimizes trip time. Data generated for hypothetical one-way and round trip human missions to Jupiter were found to be within 1% and 6% accuracy of integrated solutions respectively, verifying that for these systems, credible analysis does not require computationally intensive numerical techniques.

Williams, Craig Hamilton↗
NASA NTRS1995

Tensoral Present and Future

The coding of high-performance fluids simulations requires significant knowledge of both numerical and computational details. The magnitude and complexity of low-level details is often enough to discourage many users of turbulence data wishing to study more important, higher-level fluid dynamical questions. These same complexities are often a practical barrier to simulation experts who develop, verify and maintain the codes which generate this data. Future fluids codes, with high resolution and complex geometries, are likely to involve far more coding complexity. The research--the design and implementation of the Tensoral computer language--aims to greatly ease the coding of today's simulation and post-processing codes and at the same time provide a general computational tool for future simulations.

Dresselhaus, Eliot↗
NASA NTRS1996

Multigrid Methods for EHL Problems

In many bearings and contacts, forces are transmitted through thin continuous fluid films which separate two contacting elements. Objects in contact are normally subjected to friction and wear which can be reduced effectively by using lubricants. If the lubricant film is sufficiently thin to prevent the opposing solids from coming into contact and carries the entire load, then we have hydrodynamic lubrication, where the lubricant film is determined by the motion and geometry of the solids. However, for loaded contacts of low geometrical conformity, such as gears, rolling contact bearings and cams, this is not the case due to high pressures and this is referred to as Elasto-Hydrodynamic Lubrication (EHL) In EHL, elastic deformation of the contacting elements and the increase in fluid viscosity with pressure are very significant and cannot be ignored. Since the deformation results in changing the geometry of the lubricating film, which in turn determines the pressure distribution, an EHL mathematical model must simultaneously satisfy the complex elasticity (integral) and the Reynolds lubrication (differential) equations. The nonlinear and coupled nature of the two equations makes numerical calculations computationally intensive. This is especially true for highly loaded problems found in practice. One novel feature of these problems is that the solution may exhibit sharp pressure spikes in the outlet region. To this date both finite element and finite difference methods have been used to solve EHL problems with perhaps greater emphasis on the use of the finite difference approach. In both cases, a major computational difficulty is ensuring convergence of the nonlinear equations solver to a steady state solution. Two successful methods for achieving this are direct iteration and multigrid methods. Direct iteration methods (e.g Gauss Seidel) have long been used in conjunction with finite difference discretizations on regular meshes. Perhaps one of the best examples of the application of such methods is the recent Effective Influence Method of Dowson and Wang. Multigrid methods have also been used with great success by Venner and Venner and Lubrecht with a good summary being given by Venner. As both these finite difference discretization based approaches appear to provide an efficient way of solving EHL problems, it is important to understand their relative merits. This paper is a first attempt at providing such an understanding in the context of EHL point contact problem, (contact of two spheres), in which the contact zone is a point and an ellipse or circle for unloaded and loaded dry contacts respectively. Since the film thickness and the contact width are generally small compared to the local radius of curvature of the two surfaces, the reduced geometry of the surfaces in the contact area can be accurately approximated to the contact between a paraboloid and a flat surface. The layout of the remainder of this paper is as follows. In section 2 we introduce the form of the equations to be solved. The Effective Influence Newton Method is described in Section 3 while Section 4 describes the Multigrid method to be used. Sections 5 and 6 describe the test problems to be used in the comparison between the two methods and compare the performance of the two methods. Section 7 concludes the paper with an argument of the two methods and suggests some future research directions.

Nurgat, Elyas↗
NASA NTRS1998

New Parallel Algorithms for Structural Analysis and Design of Aerospace Structures

Subspace and Lanczos iterations have been developed, well documented, and widely accepted as efficient methods for obtaining p-lowest eigen-pair solutions of large-scale, practical engineering problems. The focus of this paper is to incorporate recent developments in vectorized sparse technologies in conjunction with Subspace and Lanczos iterative algorithms for computational enhancements. Numerical performance, in terms of accuracy and efficiency of the proposed sparse strategies for Subspace and Lanczos algorithm, is demonstrated by solving for the lowest frequencies and mode shapes of structural problems on the IBM-R6000/590 and SunSparc 20 workstations.

Nguyen, Duc T.↗
NASA NTRS1999

A New Compendium of Unsteady Aerodynamic Test Cases for CFD: Summary of AVT WG-003 Activities

With the continuous progress in hardware and numerical schemes, Computational Unsteady Aerodynamics (CUA), that is, the application of Computational Fluid Dynamics (CFD) to unsteady flowfields, is slowly finding its way as a useful and reliable tool (turbulence and transition modeling permitting) in the aircraft, helicopter, engine and missile design and development process. Before a specific code may be used with confidence it is essential to validate its capability to describe the physics of the flow correctly, or at least to the level of approximation required, for which purpose a comparison with accurate experimental data is needed. Unsteady wind tunnel testing is difficult and expensive; two factors which dramatically limit the number of organizations with the capability and/or resources to perform it. Thus, unsteady experimental data is scarce, often classified and scattered in diverse documents. Additionally, access to the reports does not necessarily assure access to the data itself. The collaborative effort described in this paper was conceived with the aim of collecting into a single easily accessible document as much quality data as possible. The idea is not new. In the early 80's NATO's AGARD (Advisory Group for Aerospace Research & Development) Structures and Material Panel (SMP) produced AGARD Report No. 702 "Compendium of Unsteady Aerodynamic Measurements", which has found and continues to find extensive use within the CUA Community. In 1995 AGARD's Fluid Dynamics Panel (FDP) decided to update and expand the former database with new geometries and physical phenomena, and launched Working Group WG-22 on "Validation Data for Computational Unsteady Aerodynamic Codes". Shortly afterwards AGARD was reorganized as the RTO (Research and Technology Organization) and the WG was renamed as AVT (Applied Vehicle Technolology) WG-003. Contributions were received from AEDC, BAe, DLR, DERA, Glasgow University, IAR, NAL, NASA, NLR, and ONERA. The final publication with the results of the exercise is expected in the second part of 1999. The aim of the present paper is to announce and present the new database to the Aeroelasticity community. It is also intended to identify, together with one of the groups of end users it targets, deficiencies in the compendium that should be addressed by means of new wind tunnel tests or by obtaining access to additionally existing data.

Ruiz-Calavera, Luis P.↗
NASA NTRS1999

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗
NASA NTRS1999

A Global Ocean Tide Model From TOPEX/POSEIDON Altimetry: GOT99.2

Goddard Ocean Tide model GOT99.2 is a new solution for the amplitudes and phases of the global oceanic tides, based on over six years of sea-surface height measurements by the TOPEX/POSEIDON satellite altimeter. Comparison with deep-ocean tide-gauge measurements show that this new tidal solution is an improvement over previous global models, with accuracies for the main semidiurnal lunar constituent M2 now below 1.5 cm (deep water only). The new solution benefits from use of prior hydrodynamic models, several in shallow and inland seas as well as the global finite-element model FES94.1. This report describes some of the data processing details involved in handling the altimetry, and it provides a comprehensive set of global cotidal charts of the resulting solutions. Various derived tidal charts are also provided, including tidal loading deformation charts, tidal gravimetric charts, and tidal current velocity (or transport) charts. Finally, low-degree spherical harmonic coefficients are computed by numerical quadrature and are tabulated for the major short-period tides; these are useful for a variety of geodetic and geophysical purposes, especially in combination with similar estimates from satellite laser ranging.

Ray, Richard D.↗
NASA NTRS2001

Global MHD Models of the Solar Corona

Global magnetohydrodynamic (MHD) models of the solar corona are computationally intensive, numerically complex simulations that have produced important new results over the past few years. After a brief overview of how these models usually work, I will address three topics: (1) How these models are now routinely used to predict the morphology of the corona and analyze Earth and space-based remote observations of the Sun; (2) The direct application of these models to the analysis of physical processes in the corona and chromosphere and to the interpretation of in situ solar wind observations; and (3) The use of results from global models to validate the approximations used to make detailed studies of physical processes in the corona that are not otherwise possible using the global models themselves.

Suess, S. T.↗
NASA NTRS2001

D3: A Collaborative Infrastructure for Aerospace Design

DARWIN is a NASA developed, Internet-based system for enabling aerospace researchers to securely and remotely access and collaborate on the analysis of aerospace vehicle design data, primarily the results of wind-tunnel testing and numeric (e.g., computational fluid dynamics) model executions. DARWIN captures, stores and indexes data, manages derived knowledge (such as visualizations across multiple data sets) and provides an environment for designers to collaborate in the analysis of the results of testing. DARWIN is an interesting application because it supports high volumes of data, integrates multiple modalities of data display (e.g. images and data visualizations), and provides non-trivial access control mechanisms. DARWIN enables collaboration by allowing not only sharing visualizations of data, but also commentary about and view of data.

Walton, Joan↗
NASA NTRS2002

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗
NASA NTRS1986

Digital Imaging

Digital Imaging is the computer processed numerical representation of physical images. Enhancement of images results in easier interpretation. Quantitative digital image analysis by Perceptive Scientific Instruments, locates objects within an image and measures them to extract quantitative information. Applications are CAT scanners, radiography, microscopy in medicine as well as various industrial and manufacturing uses. The PSICOM 327 performs all digital image analysis functions. It is based on Jet Propulsion Laboratory technology, is accurate and cost efficient.

Source record↗
NASA NTRS2003

A Discrete Analysis of Non-reflecting Boundary Conditions for Discontinuous Galerkin Method

We present a discrete analysis of non-reflecting boundary conditions for the discontinuous Galerkin method. The boundary conditions considered in this paper include the recently proposed Perfectly Matched Layer absorbing boundary condition for the linearized Euler equation and two non-reflecting boundary conditions based on the characteristic decomposition of the flux on the boundary. The analyses for the three boundary conditions are carried out in a unifled way. In each case, eigensolutions of the discrete system are obtained and applied to compute the numerical reflection coefficients of a specified out-going wave. The dependencies of the reflections at the boundary on the out-going wave angle and frequency as well as the mesh sizes arc? studied. Comparisons with direct numerical simulation results are also presented.

Hu, Fang Q.↗
NASA NTRS2003

Prompt and Precise Prototyping

For Sanders Design International, Inc., of Wilton, New Hampshire, every passing second between the concept and realization of a product is essential to succeed in the rapid prototyping industry where amongst heavy competition, faster time-to-market means more business. To separate itself from its rivals, Sanders Design aligned with NASA's Marshall Space Flight Center to develop what it considers to be the most accurate rapid prototyping machine for fabrication of extremely precise tooling prototypes. The company's Rapid ToolMaker System has revolutionized production of high quality, small-to-medium sized prototype patterns and tooling molds with an exactness that surpasses that of computer numerically-controlled (CNC) machining devices. Created with funding and support from Marshall under a Small Business Innovation Research (SBIR) contract, the Rapid ToolMaker is a dual-use technology with applications in both commercial and military aerospace fields. The advanced technology provides cost savings in the design and manufacturing of automotive, electronic, and medical parts, as well as in other areas of consumer interest, such as jewelry and toys. For aerospace applications, the Rapid ToolMaker enables fabrication of high-quality turbine and compressor blades for jet engines on unmanned air vehicles, aircraft, and missiles.

Source record↗
NASA NTRS2004

Injection Induced Mixing in Flows Separating From Smooth Surfaces

An analytic model for predicting the effect of unsteady local surface injection on the flow separating from a streamlined body at angle of attack is proposed. The model uses the premise that separation control results from enhanced mixing along the shear layer that develops between the main stream and the fluid in the underlying recirculation zone. High-Reynolds-number asymptotic methods are used to connect the unsteady surface injection to an instability wave propagating on the separating shear layer and then to the large-scale coherent structures that produce the increased mixing. The results is a tool that can guide the choice of fluid-actuator parameters to maximize flow-control effectiveness and may also facilitate computer-based numerical experiments.

Adamczyk, John J.↗
NASA NTRS1998

Efficient Jitter Analysis for Spacecraft

Typically in space missions, the science instruments require a specific degree of pointing accuracy as well as dynamical quietness. This dynamical quietness, which is needed to allow the instruments to make measurements (remote sensing applications) or perform other functions, is usually characterized in terms of jitter and stability specifications. In order to insure that the spacecraft meets the requirements of its instruments, several jitter analyses are performed throughout the design phase of the spacecraft and beyond as the models of the spacecraft, its components, and disturbances mature. Each such analysis involves the simulation of the spacecraft and instrument dynamical response to all known disturbance scenarios, followed by the computation of jitter values for each instrument based on the specified jitter time windows. The direct approach for computing jitter values by sweeping maxima and minima throughout the time history may be costly in the computational sense as the size and number of the time histories involved could be quite large. Keeping in mind that typical spacecraft simulation time histories may easily involve hundreds of thousands or millions of points, it is imperative that a jitter analysis algorithm be developed which is more efficient than the direct approach. This paper presents a vectorized algorithm for efficient computation of spacecraft jitter values. The algorithm identifies the extreme points in the time history, which are the points that may dominate the jitter values depending on the location of the jitter window along the time history. The span of influence of each extremum is then computed by the algorithm and used in an efficient and vectorized fashion to obtain the jitter values. The algorithm deals with the multiple jitter windows 2 sequentially, first computing jitter values for the smallest time window, then looping over the remaining time windows until all jitter values are computed. A numerical example is carried out to demonstrate the efficiency and feasibility of the proposed jitter analysis technique.

Maghami, Peiman G.↗