Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Structural Damage Prediction and Analysis for Hypervelocity Impact: Properties of Largest Fragment Produced by Hypervelocity Impact of Aluminum Spheres with Thin Aluminum Sheets

Results of a series of hypervelocity impact tests are presented. In these tests, 1.275-g, 9.53-mm-diameter, 2017-T4 aluminum spheres were fired at normal incidence at eight thicknesses of 6061-T6 aluminum sheet. Bumper thickness to projectile diameter (t/D) ratio ranged from 0.026 to 0.424. Nominal impact velocity was 6.7 km/s. Results of five tests using 6.35, 9.53, and 12.70-mm-diameter aluminum spheres and other aluminum alloy bumpers are also given. A large chunky fragment of projectile was observed at the center of the debris clouds produced by the impacts. The equivalent diameter of this large fragment ranged from 5.5 mm for the lowest t/D ratio to a minimum of 0.6 mm for the case where maximum breakup of the projectile occurred (t/D approximately 0.2 to 0.3). When the t/D ratio was 0.42, numerous large flaky fragments were evenly distributed in the external bubble of bumper debris. Velocity of the large central fragments decreased continuously with increasing t/D ratio, ranging from about 99 percent to less than 80 percent of the impact velocity. The change in the velocity of small fragments spalling from the rear of the projectile was used to obtain a relationship showing a linear increase in the size of the central projectile fragment with decrease in the shock-induced stress in the projectile.

Piekutowski, Andrew J.↗
NASA NTRS1995

Pseudo-time methods for constrained optimization problems governed by PDE

In this paper we present a novel method for solving optimization problems governed by partial differential equations. Existing methods are gradient information in marching toward the minimum, where the constrained PDE is solved once (sometimes only approximately) per each optimization step. Such methods can be viewed as a marching techniques on the intersection of the state and costate hypersurfaces while improving the residuals of the design equations per each iteration. In contrast, the method presented here march on the design hypersurface and at each iteration improve the residuals of the state and costate equations. The new method is usually much less expensive per iteration step since, in most problems of practical interest, the design equation involves much less unknowns that that of either the state or costate equations. Convergence is shown using energy estimates for the evolution equations governing the iterative process. Numerical tests show that the new method allows the solution of the optimization problem in a cost of solving the analysis problems just a few times, independent of the number of design parameters. The method can be applied using single grid iterations as well as with multigrid solvers.

Taasan, Shlomo↗
NASA NTRS1997

Comparative study of energetic heliospheric electrons at wind and geotail

The current fleet of both heliospheric and magnetospheric spacecraft provides a unique opportunity of studying spatio-temporal plasma phenomena. Among the numerous topics that can be addressed by such a fleet, figures the 3D study of energetic (greater than 20 keV) solar electron) solar electron events. WIND and GEOTAIL are particulary interesting for such analysis: more than two years of data; up to 3/4 continuous days spent by GEOTAIL in the solar wind; and comparable experiments. During solar electron events, rapid electron flux changes can be associated with these structures. Sometimes they are seen at both spacecraft with a time difference corresponding to the convection time, but sometimes they are observed at one spacecraft but not at the other. Several preliminary conclusions can be drawn from such events: (1) there is direct evidence that, within distances lower than 1.5 x 10(exp 6) km at 1 A (approximately 7000 km at the sun), there are significant spatial variations of (a) the magnetic field line connections to the source region and probably of (b) the particle propagation features; (2) magnetic field lines inside and outside particle propagation structures can keep distinct access to the source region for more than one hour; (3) flux waves, of probable local origin, were observed; and (4) within one hour, spatial flux discrepancies, separated by magnetic structures, may remain or disappear.

Chaizy, P. A.↗
NASA NTRS2014

S-NPP ATMS Instrument Prelaunch and On-Orbit Performance Evaluation

The first of a new generation of microwave sounders was launched aboard the Suomi-National Polar-Orbiting Partnership satellite in October 2011. The Advanced Technology Microwave Sounder (ATMS) combines the capabilities and channel sets of three predecessor sounders into a single package to provide information on the atmospheric vertical temperature and moisture profiles that are the most critical observations needed for numerical weather forecast models. Enhancements include size/mass/power approximately one third of the previous total, three new sounding channels, the first space-based, Nyquist-sampled cross-track microwave temperature soundings for improved fusion with infrared soundings, plus improved temperature control and reliability. This paper describes the ATMS characteristics versus its predecessor, the advanced microwave sounding unit (AMSU), and presents the first comprehensive evaluation of key prelaunch and on-orbit performance parameters. Two-year on-orbit performance shows that the ATMS has maintained very stable radiometric sensitivity, in agreement with prelaunch data, meeting requirements for all channels (with margins of ~40% for channels 1-15), and improvements over AMSU-A when processed for equivalent spatial resolution. The radiometric accuracy, determined by analysis from ground test measurements, and using on-orbit instrument temperatures, also shows large margins relative to requirements (specified as <1.0K for channels 1, 2, and 16-22 and <0.75 K for channels 3-15). A thorough evaluation of the performance of ATMS is especially important for this first proto-flight model unit of what will eventually be a series of ATMS sensors providing operational sounding capability for the U.S. and its international partners well into the next decade.

temperature↗
NASA NTRS1987

A higher-order theory for geometrically nonlinear analysis of composite laminates

A third-order shear deformation theory of laminated composite plates and shells is developed, the Navier solutions are derived, and its finite element models are developed. The theory allows parabolic description of the transverse shear stresses, and therefore the shear correction factors of the usual shear deformation theory are not required in the present theory. The theory also accounts for the von Karman nonlinear strains. Closed-form solutions of the theory for rectangular cross-ply and angle-ply plates and cross-ply shells are developed. The finite element model is based on independent approximations of the displacements and bending moments (i.e., mixed finite element model), and therefore, only C sup o -approximation is required. The finite element model is used to analyze cross-ply and angle-ply laminated plates and shells for bending and natural vibration. Many of the numerical results presented here should serve as references for future investigations. Three major conclusions resulted from the research: First, for thick laminates, shear deformation theories predict deflections, stresses and vibration frequencies significantly different from those predicted by classical theories. Second, even for thin laminates, shear deformation effects are significant in dynamic and geometrically nonlinear analyses. Third, the present third-order theory is more accurate compared to the classical and firt-order theories in predicting static and dynamic response of laminated plates and shells made of high-modulus composite materials.

Reddy, J. N.↗
NASA NTRS1982

Nonlinear zero-sum differential game analysis by singular perturbation methods

A class of nonlinear, zero-sum differential games, exhibiting time-scale separation properties, can be analyzed by singular-perturbation techniques. The merits of such an analysis, leading to an approximate game solution, as well as the 'well-posedness' of the formulation, are discussed. This approach is shown to be attractive for investigating pursuit-evasion problems; the original multidimensional differential game is decomposed to a 'simple pursuit' (free-stream) game and two independent (boundary-layer) optimal-control problems. Using multiple time-scale boundary-layer models results in a pair of uniformly valid zero-order composite feedback strategies. The dependence of suboptimal strategies on relative geometry and own-state measurements is demonstrated by a three dimensional, constant-speed example. For game analysis with realistic vehicle dynamics, the technique of forced singular perturbations and a variable modeling approach is proposed. Accuracy of the analysis is evaluated by comparison with the numerical solution of a time-optimal, variable-speed 'game of two cars' in the horizontal plane.

Sinar, J.↗
NASA NTRS1986

Adaptive computational methods for SSME internal flow analysis

Adaptive finite element methods for the analysis of classes of problems in compressible and incompressible flow of interest in SSME (space shuttle main engine) analysis and design are described. The general objective of the adaptive methods is to improve and to quantify the quality of numerical solutions to the governing partial differential equations of fluid dynamics in two-dimensional cases. There are several different families of adaptive schemes that can be used to improve the quality of solutions in complex flow simulations. Among these are: (1) r-methods (node-redistribution or moving mesh methods) in which a fixed number of nodal points is allowed to migrate to points in the mesh where high error is detected; (2) h-methods, in which the mesh size h is automatically refined to reduce local error; and (3) p-methods, in which the local degree p of the finite element approximation is increased to reduce local error. Two of the three basic techniques have been studied in this project: an r-method for steady Euler equations in two dimensions and a p-method for transient, laminar, viscous incompressible flow. Numerical results are presented. A brief introduction to residual methods of a-posterior error estimation is also given and some pertinent conclusions of the study are listed.

Oden, J. T.↗
NASA NTRS1990

Stability analysis of the accretion line

The stability of the flow along the accretion line in the Bondi-Hoyle-Lyttleton-type accretion flow is studied. The WKB approximation is used to analyze short wavelength linear perturbations of the axisymmetric steady state flow. The viscosity terms in the gasdynamical equations are included in the analysis. The radial modes are studied both in two- and three-dimensional accretion flow. Since a displaced accretion line cannot be defined in three-dimensional flow, the tangential modes are only analyzed in the two-dimensional flow. The analysis of the perturbed gasdynamical equations are made simpler as a result of the radial and tangential modes being decoupled in the linear regime. As was shown before, the three-dimensional accretion flow is unstable against radial modes. The two-dimensional accretion flow is found to be unstable against tangential modes, as well as against radial modes. Using typical numerical parameters which are in common use, the wavelength which gives the maximum growth rate of a small perturbation is found to be about 0.5 r(a), where r(a) is the accretion radius.

Soker, Noam↗
NASA NTRS1959

Numerical Solution of the Flow of a Perfect Gas Over A Circular Cylinder at Infinite Mach Number

A solution for the two-dimensional flow of an inviscid perfect gas over a circular cylinder at infinite Mach number is obtained by numerical methods of analysis. Nonisentropic conditions of curved shock waves and vorticity are included in the solution. The analysis is divided into two distinct regions, the subsonic region which is analyzed by the relaxation method of Southwell and the supersonic region which was treated by the method of characteristics. Both these methods of analysis are inapplicable on the sonic line which is therefore considered separately. The shapes of the sonic line and the shock wave are obtained by iteration techniques. The striking result of the solution is the strong curvature of the sonic line and of the other lines of constant Mach number. Because of this the influence of the supersonic flow on the sonic line is negligible. On comparison with Newtonian flow methods, it is found that the approximate methods show a larger variation of surface pressure than is given by the present solution.

Hamaker, Frank M.↗
NASA NTRSDate not supplied

Towards Real-time Global Assessment of Post-fire Debris Flow Hazard

As the risk of wildfires increases worldwide, burned steeplands are vulnerable to the secondary hazard of widespread sediment mobilization through debris flows. Following an initial burn, sediment and soil previously restrained by vegetation are no longer consolidated, allowing for easy mobilization into channels and along steep hillslopes through runoff. Sufficiently powerful rainfall incorporates entrained material into turbulent flows and serves as the primary trigger for debris flow initiation. There is thus an ongoing need to establish the relationship between rainfall and debris flow initiation based on a variety of spatiotemporal preconditions. Previous work establishes regional and local thresholds to constrain the effect of rainfall in recently burned areas, but no empirical or numerical solution has worldwide application. Building from regionally-based efforts in the U.S., this work considers how remote sensing data can be applied to better approximate the post-fire debris flow hazards worldwide using freely available global datasets and software. Our work assesses the utility of remote sensing resources for analyzing burn characteristics, topography, rainfall intensity/duration, and, thus, debris flow initiation. Early results show that global observations are sufficient to delineate background rainfall rates from storms likely to cause debris flows across a variety of burn severity and topographic conditions. However, the dearth of publicly-available post-fire debris flow inventories globally limit the ability to test how the model framework performs within different climatologic and morphologic areas. This work will present preliminary analysis over the Western United States and demonstrate the feasibility of a global, near-real time model to provide situational awareness of potential hazards within recently burned areas worldwide. Future work will also consider how global or regional precipitation forecasts may increase the lead time for improved early warning of these hazards.

Elijah Orland↗
NASA NTRS1987

Down-to-earth supercomputing at NASA Ames

The requirements for computers for aerospace vehicle design problems, in particular fluid dynamics, are examined. The levels of approximation to Navier-Stokes equations are: (1) linearized inviscid form, (2) nonlinear inviscid equation, (3) Reynolds-averaged Navier-Stokes equation, (4) large eddy simulation with models for subgrid scale turbulence, and (5) full Navier-Stokes equations to calculate laminar/turbulent transition and turbulence dissipation. The Cray 2 with 256 megaword memory and 4.1 nsec clock and the Cray X-MP/48 with a 128 megaword memory and a 9.5 nsec clock of the Ames facility are described, and an example revealing their applicability is presented. The use of zonal methods to simulate aerodynamic flow fields numerically is discussed. Consideration is given to the development of a program for simulation of aerodynamic flow fields, the compilers, symbolic debuggers, the scheduling of multiple CPUs, and the analysis and storage of the data.

Stevens, Kenneth G., Jr.↗
NASA NTRS1990

On the partitioning strategy based on symmetry transformations

A computational procedure is presented for the analysis of unsymmetric structures. The procedure is based on a modified version of the symmetry-transformation partitioning strategy, in which the response of the structure is approximated by a linear combination of symmetric/antisymmetric response vectors, each obtained by using only a fraction of the degrees of freedom of the original FEM model. The three key elements of the procedure are: (1) mixed (or primitive-variable) formulation with independent shape functions for the different fields; (2) restructuring of the governing discrete equations of the structure into uncoupled sets in the symmetric and antisymmetric response vectors; and (3) a stable and efficient iterative process for generating the response of the structure. The effectiveness of the proposed procedure and its advantages over classical substructuring are demonstrated by means of numerical examples.

Noor, Ahmed K.↗
NASA NTRS1999

Cloud Statistics for NASA Climate Change Studies

The Principal Investigator participated in two field experiments and developed a global data set on cirrus cloud frequency and optical depth to aid the development of numerical models of climate. Four papers were published under this grant. The accomplishments are summarized: (1) In SUCCESS (SUbsonic aircraft: Contrail & Cloud Effects Special Study) the Principal Investigator aided weather forecasters in the start of the field program. A paper also was published on the clouds studied in SUCCESS and the use of the satellite stereographic technique to distinguish cloud forms and heights of clouds. (2) In SHEBA (Surface Heat Budget in the Arctic) FIRE/ACE (Arctic Cloud Experiment) the Principal Investigator provided daily weather and cloud forecasts for four research aircraft crews, NASA's ER-2, UCAR's C-130, University of Washington's Convert 580, and the Canadian Atmospheric Environment Service's Convert 580. Approximately 105 forecasts were written. The Principal Investigator also made daily weather summaries with calculations of air trajectories for 54 flight days in the experiment. The trajectories show where the air sampled during the flights came from and will be used in future publications to discuss the origin and history of the air and clouds sampled by the aircraft. A paper discussing how well the FIRE/ACE data represent normal climatic conditions in the arctic is being prepared. (3) The Principal Investigator's web page became the source of information for weather forecasting by the scientists on the SHEBA ship. (4) Global Cirrus frequency and optical depth is a continuing analysis of global cloud cover and frequency distribution are being made from the NOAA polar orbiting weather satellites. This analysis is sensitive to cirrus clouds because of the radiative channels used. During this grant three papers were published which describe cloud frequencies, their optical properties and compare the Wisconsin FM Cloud Analysis to other global cloud data such as the International Satellite Cloud Climatology Program (ISCCP) and the Stratospheric Aerosol and Gas Experiment (SAGE). A summary of eight years of HIRS data will be published in late 1998. Important information from this study are: 1) cirrus clouds cover most of the earth, 2) they are found about 40% of the time globally, 3) in the tropics cirrus cloud frequencies are even higher, from 80-100%, 4) there is slight evidence that cirnis cloud cover is increasing in the northern hemisphere at about 0.5% per year, and 5) cirrus clouds have an average infrared transmittance of about 40% of the terrestrial radiation. (5) Global Cloud Frequency Statistics published on the Principal Investigator's web page have been used in the planning of the future CRYSTAL experiment and have been used for refinements of a global numerical model operated at the Colorado State University.

Wylie, Donald P.↗
NASA NTRS1986

Application of viscous-inviscid interaction methods to transonic turbulent flows

Two different viscous-inviscid interaction schemes were developed for the analysis of steady, turbulent, transonic, separated flows over axisymmetric bodies. The viscous and inviscid solutions are coupled through the displacement concept using a transpiration velocity approach. In the semi-inverse interaction scheme, the viscous and inviscid equations are solved in an explicitly separate manner and the displacement thickness distribution is iteratively updated by a simple coupling algorithm. In the simultaneous interaction method, local solutions of viscous and inviscid equations are treated simultaneously, and the displacement thickness is treated as an unknown and is obtained as a part of the solution through a global iteration procedure. The inviscid flow region is described by a direct finite-difference solution of a velocity potential equation in conservative form. The potential equation is solved on a numerically generated mesh by an approximate factorization (AF2) scheme in the semi-inverse interaction method and by a successive line overrelaxation (SLOR) scheme in the simultaneous interaction method. The boundary-layer equations are used for the viscous flow region. The continuity and momentum equations are solved inversely in a coupled manner using a fully implicit finite-difference scheme.

Lee, D.↗
NASA NTRS1995

An instability in planetary rings due to ballistic transport

Ballistic transport in planetary rings is the net radial transport of mass and angular momentum due to exchanges of meteoroid impact ejecta between neighboring ring regions. The detailed linear stability analysis in this paper demonstrates that ballistic transport causes wavelike disturbances to grow and propagate in an otherwise uniform ring. The growth is strongest for intermediate values of the normal ring optical depth tau = 0.1 to 1.0 and goes to zero as tau approaches 0 and tau approaches infinity. For nominal values of various physical parameters, the minimum e-folding time is approximately 10(exp 5) years for tau approximately 0.4. The direction of propagation is opposite to the sense of any asymmetry that may exist in the ejecta direction distribution (inward for prograde ejecta and outward for retrograde ejecta). The additional effect of viscous transport tends to damp wavelike perturbations strongly at short wavelengths and at high values of tau. The quantitative agreement between this analytic work and numerical simulations reported elsewhere is generally quite good. As applied to Saturn's rings, the results in this paper strengthen the earlier conclusion from numerical calculations that the 100-km structure in the inner B Ring is caused by ballistic transport. However, it is also clear that ballistic transport cannot produce the complex structure seen in the outer two-thirds of the B Ring where tau greater than or approximately 1.5. Wavelike structures in the C Ring might also be attributed to ballistic transport; but this requires further study.

Durisen, Richard H.↗
NASA NTRS1971

Transient-forced convection film boiling on an isothermal flat plate

An approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model is developed. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady state condition are obtained. The influence of interfacial shear is seen to be negligible while the Prandtl Number and the ratio (C sub p delta T divided by h sub fg times Pr) have major influence. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5 percent. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100 percent higher than those at steady state. The average rate of vapor convected away is 10 to 15 percent lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger. Further, the total heat transferred during transient increases and the evaporation decreases for increasing values of C sub p delta T divided by h sub fg times Pr

Nagendra, H. R.↗
NASA NTRS1972

Smoothing of functions of range and range rate measurements from earth orbiting satellites

It is shown that for satellites in circular Earth orbits with altitudes of 500 kilometers to 1500 kilometers, and for satellites in elliptical orbits with an approximate 4000 kilometer height of perigee, a high degree least squares polynomial (e.g., a 9th or 10th degree in some cases) is required to smooth both range and range rate data for purposes of input to orbit determination programs. In order to circumvent this problem, functions of range and range rate are smoothed with lower degree least squares polynomials (e.g., 3rd and 4th degree) and it is shown that under the above geometric constraints the standard deviation of fit can be reduced to levels commensurate with typical S-band tracking system resolution which is 1 to 2 meters in range and 0.005 meters/second in range rate for a 1 per second data rate. Also shown are the effects of Gaussian random noise, biases, and periodic noise. This analysis includes numerous examples applied to the 44 point data smoothing interval currently used in much of the operational preprocessing at the Goddard Space Flight Center.

Grenchik, T. J.↗
NASA NTRS1987

Model-size reduction for the analysis of symmetric structures with asymmetric boundary conditions

A simple computational procedure is presented for reducing the size of the analysis model for a symmetric structure with asymmetric boundary conditions to that of the corresponding structure with symmetric boundary conditions. The procedure is based on approximating the asymmetric response of the structure by a linear combination of symmetric and antisymmetric global approximation vectors (or modes). The key elements of the procedure are (1) restructuring the governing finite-element equations to delineate the contributions to the symmetric and antisymmetric components of the asymmetric response, (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The effectiveness of the computational procedure is demonstrated by means of numerical examples of linear static problems of shells, and its potential for solving nonlinear problems is discussed.

Noor, Ahmed K.↗