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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Lunar Soil Erosion Physics for Landing Rockets on the Moon

To develop a lunar outpost, we must understand the blowing of soil during launch and landing of the new Altair Lander. For example, the Apollo 12 Lunar Module landed approximately 165 meters from the deactivated Surveyor Ill spacecraft, scouring its surfaces and creating numerous tiny pits. Based on simulations and video analysis from the Apollo missions, blowing lunar soil particles have velocities up to 2000 m/s at low ejection angles relative to the horizon, reach an apogee higher than the orbiting Command and Service Module, and travel nearly the circumference of the Moon [1-3]. The low ejection angle and high velocity are concerns for the lunar outpost.

Clegg, Ryan N.↗
NASA NTRS1990

A numerical and analytical study of nonlinear bifurcations associated with the morphological stability of two-dimensional single crystals

The nonlinear stability of a two-dimensional single crystal of pure material in an undercooled melt is studied both analytically and numerically. The quasi-steady state approximation is used for the thermal fields, and the effects of different solid and liquid thermal conductivities, isotropic interfacial growth kinetics, and isotropic surface tension are included. The bifurcation analysis is performed by calculating the instantaneous value of the fundamental component of the local normal growth speed for an interface perturbed by a single Fourier shape component. Numerically, the fundamental component of the interfacial growth speed is found by Fourier analysis of the solution to an integrodifferential equation obeyed at the interface. Analytically, an expansion technique is used to derive a solvability condition defining each of these bifurcation points. The analytical and numerical results are in very close agreement. Almost all of the bifurcations are subcritical, and the results are presented by giving values of the Landau coefficient as a function of the different dimensionless parameters used in the model.

Brush, L. N.↗
NASA NTRS1982

Numerical analysis of confined turbulent flow

The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.

Lin, A.↗
NASA NTRS1984

Applications of finite element and wave envelope element approximations to turbofan engine noise radiation including flight effects

The problem of acoustic radiation from turbofan engine inlets in flow has not lent itself fully to analysis by numerical means because of the large domains and high frequencies involved. The current work has extended the use of finite elements and wave envelope elements, elements which simulate decay and wavelike behaviour in their interpolation functions, from the no-flow case in which they have been proven, to cases incorporating mean flow. By employing an irrotational mean flow assumption, the acoustics problem has been posed in an axisymmetric formulation in terms of acoustic velocity potential, thus minimizing computer solution storage requirements. The results obtained from the numerical procedures agree well with known analytical solutions, static experimental jet engines inflow data, and also with flight test results.

Parrett, A. V.↗
NASA NTRS1981

On the utility of finite element theory for computational fluid dynamics

An implicit finite element numerical solution algorithm is derived for the compressible Navier-Stokes equations expressed in generalized coordinates. The theoretical basis utilizes a Galerkin-Weighted Residuals formulation, and extremization of approximation error within the context of a multipole expansion. A von Neumann analysis for a simplified form indicates the algorithm fourth- to sixth-order phase accurate, with third-order dissipation for the elementary linear element construction. Performance is improved for the algorithm constructed using quadratic interpolation. Numerical experiments for shocked duct flows are employed to optimize the several algorithm parameters. Additional numerical solutions validate algorithm accuracy and utility for aerodynamics applications.

Baker, A. J.↗
NASA NTRS1988

The quasi-rigid rotation of coronal magnetic fields

Spherical harmonic analysis and numerical simulations are used to study the rotational properties of the coronal magnetic field under the assumption that it can be approximated by a current-free extension of the photospheric field. It is found that the rotation rate in the outer corona is determined, principally, by coronal filtering, the global averages of the photospheric rotation rate, and ongoing source eruptions. The present model is able to account for observationally inferred rotational properties. It is suggested that the coronal rotation rate accelerates gradually due to the equatorward migration of sunspots, and that the 27-day equatorial period is approached toward sunspot minimum as the decaying photospheric flux becomes localized near the equator.

Wang, Y.-M.↗
NASA NTRS2014

Mission Design and Analysis for Suborbital Intercept and Fragmentation of an Asteroid with Very Short Warning Time

Small near-Earth objects (NEOs) approximately 50-150 m in size are far more numerous (hundreds of thousands to millions yet to be discovered) than larger NEOs. Small NEOs, which are mostly asteroids rather than comets, are very faint in the night sky due to their small sizes, and are, therefore, difficult to discover far in advance of Earth impact. Furthermore, even small NEOs are capable of creating explosions with energies on the order of tens or hundreds of megatons (Mt). We are, therefore, motivated to prepare to respond effectively to short warning time, small NEO impact scenarios. In this paper we explore the lower bound on actionable warning time by investigating the performance of notional upgraded Intercontinental Ballistic Missiles (ICBMs) to carry Nuclear Explosive Device (NED) payloads to intercept and disrupt a hypothetical incoming NEO at high altitudes (generally at least 2500 km above Earth). We conduct this investigation by developing optimal NEO intercept trajectories for a range of cases and comparing their performances. Our results show that suborbital NEO intercepts using Minuteman III or SM-3 IIA launch vehicles could achieve NEO intercept a few minutes prior to when the NEO would strike Earth. We also find that more powerful versions of the launch vehicles (e.g., total delta V of approximately 9.5-11 km/s) could intercept incoming NEOs several hours prior to when the NEO would strike Earth, if launched at least several days prior to the time of intercept. Finally, we discuss a number of limiting factors and practicalities that affect whether the notional systems we describe could become feasible.

near- Earth object↗
NASA NTRS2006

Consistency of Post-Newtonian Waveforms with Numerical Relativity

General relativity predicts the gravitational radiation signatures of mergers of compact binaries, such as coalescing binary black hole systems. Derivations of waveform predictions for such systems are required for optimal scientific analysis of observational gravitational wave data, and have so far been achieved primarily with the aid of the post-Newtonian (PN) approximation. The quality of this treatment is unclear, however, for the important late inspiral portion. We derive late-inspiral waveforms via a complementary approach, direct numerical simulation of Einstein's equations, which has recently matured sufficiently for such applications. We compare waveform phasing from simulations covering the last approximately 14 cycles of gravitational radiation from an equal-mass binary system of nonspinning black holes with the corresponding 3PN and 3.5PN orbital phasing. We find agreement consistent with internal error estimates based on either approach at the level of one radian over approximately 10 cycles. The result suggests that PN waveforms for this system are effective roughly until the system reaches its last stable orbit just prior to the final merger/

Baker, John G.↗
NASA NTRS2007

Consistency of Post-Newtonian Waveforms with Numerical Relativity

General relativity predicts the gravitational radiation signatures of mergers of compact binaries,such as coalescing binary black hole systems. Derivations of waveform predictions for such systems are required for optimal scientific analysis of observational gravitational wave data, and have so far been achieved primarily with the aid of the post-Newtonian (PN) approximation. The quaIity of this treatment is unclear, however, for the important late inspiral portion. We derive late-inspiral wave forms via a complementary approach, direct numerical simulation of Einstein's equations, which has recently matured sufficiently for such applications. We compare waveform phasing from simulations covering the last approximately 14 cycles of gravitational radiation from an equal-mass binary system of nonspinning black holes with corresponding 3PN and 3.5PN waveforms. We find phasing agreement consistent with internal error estimates based in either approach, at the level of one radian over approximately 10 cycles. The result suggests that PN waveforms for this system are effective roughly until the system reaches its last stable orbit just prior to the final merger.

Baker, John G.↗
NASA NTRS2003

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗
NASA NTRS2004

Residual Strength Characterization of a Curved Integrally-Stiffened Panel

Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic-plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intension of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels. In this regard, a series of fracture tests were conducted on curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. During the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).

Seshadri, B. R.↗
NASA NTRS1992

Numerical analysis of turbine blade tip treatments

Three-dimensional solutions of the Navier-Stokes equations for a turbine blade with a turning angle of 180 degrees have been computed, including blade tip treatments involving cavities. The geometry approximates a preliminary design for the GGOT (Generic Gas Oxidizer Turbine). The data presented here will be compared with experimental data to be obtained from a linear cascade using original GGOT blades. Results have been computed for a blade with 1 percent clearance, based on chord, and three different cavity sizes. All tests were conducted at a Reynolds number of 4 x 10 exp 7. The grid contains 39,440 points with 10 spanwise planes in the tip clearance region of 5.008E-04 m. Streamline plots and velocity vectors together with velocity divergence plots reveal the general flow behavior in the clearance region. Blade tip temperature calculations suggest placement of a cavity close to the upstream side of the blade tip for reduction of overall blade tip temperature. The solutions do not account for the relative motion between the endwall and the turbine blade. The solutions obtained are generally consistent with previous work done in this area,

Gopalaswamy, Nath S.↗
NASA NTRS1995

Analysis of discretization errors in LES

All numerical simulations of turbulence (DNS or LES) involve some discretization errors. The integrity of such simulations therefore depend on our ability to quantify and control such errors. In the classical literature on analysis of errors in partial differential equations, one typically studies simple linear equations (such as the wave equation or Laplace's equation). The qualitative insight gained from studying such simple situations is then used to design numerical methods for more complex problems such as the Navier-Stokes equations. Though such an approach may seem reasonable as a first approximation, it should be recognized that strongly nonlinear problems, such as turbulence, have a feature that is absent in linear problems. This feature is the simultaneous presence of a continuum of space and time scales. Thus, in an analysis of errors in the one dimensional wave equation, one may, without loss of generality, rescale the equations so that the dependent variable is always of order unity. This is not possible in the turbulence problem since the amplitudes of the Fourier modes of the velocity field have a continuous distribution. The objective of the present research is to provide some quantitative measures of numerical errors in such situations. Though the focus of this work is LES, the methods introduced here can be just as easily applied to DNS. Errors due to discretization of the time-variable are neglected for the purpose of this analysis.

Ghosal, Sandip↗
NASA NTRS1979

A numerical comparison with an exact solution for the transient response of a cylinder immersed in a fluid

The transient response of an elastic cylindrical shell immersed in an acoustic media that is engulfed by a plane wave is determined numerically. The method applies to the USA-STAGS code which utilizes the finite element method for the structural analysis and the doubly asymptotic approximation for the fluid-structure interaction. The calculations are compared to an exact analysis for two separate loading cases: a plane step wave and an exponentially decaying plane wave.

Giltrud, M. E.↗
NASA NTRS1990

Resolution of the 1D regularized Burgers equation using a spatial wavelet approximation

The Burgers equation with a small viscosity term, initial and periodic boundary conditions is resolved using a spatial approximation constructed from an orthonormal basis of wavelets. The algorithm is directly derived from the notions of multiresolution analysis and tree algorithms. Before the numerical algorithm is described these notions are first recalled. The method uses extensively the localization properties of the wavelets in the physical and Fourier spaces. Moreover, the authors take advantage of the fact that the involved linear operators have constant coefficients. Finally, the algorithm can be considered as a time marching version of the tree algorithm. The most important point is that an adaptive version of the algorithm exists: it allows one to reduce in a significant way the number of degrees of freedom required for a good computation of the solution. Numerical results and description of the different elements of the algorithm are provided in combination with different mathematical comments on the method and some comparison with more classical numerical algorithms.

Liandrat, J.↗
NASA NTRS1993

Autonomous vehicle motion control, approximate maps, and fuzzy logic

Progress on research on the control of actions of autonomous mobile agents using fuzzy logic is presented. The innovations described encompass theoretical and applied developments. At the theoretical level, results of research leading to the combined utilization of conventional artificial planning techniques with fuzzy logic approaches for the control of local motion and perception actions are presented. Also formulations of dynamic programming approaches to optimal control in the context of the analysis of approximate models of the real world are examined. Also a new approach to goal conflict resolution that does not require specification of numerical values representing relative goal importance is reviewed. Applied developments include the introduction of the notion of approximate map. A fuzzy relational database structure for the representation of vague and imprecise information about the robot's environment is proposed. Also the central notions of control point and control structure are discussed.

Ruspini, Enrique H.↗
NASA NTRSDate not supplied

Variable Mixing Nozzle Design with Slotted Vortex Generators for Jet Noise Reduction

A new variable geometry turbofan nozzle concept is presented with the dual goals of airport noise reduction and high propulsive efficiency at cruise, while employing only a single moving part. The nozzle utilizes a number of curved vanes to function as vortex generators (VGs) during takeoff and initial climb. Streamwise vortices provide an increase in mixing between the jet and external stream, with a resulting decrease in high frequency noise generation and a noise source distribution which is more amenable to airframe-based shielding. Unlike conventional VGs, the vanes are positioned without transverse incidence, and a pressure difference across the surface of each vane is created via an adjacent slot. During high altitude flight, when nozzle efficiency concerns outweigh any preference for jet noise reduction, all slots may be closed through rotation of a slotted ring inside the nozzle. RANS-based numerical analysis is performed to calculate propulsive metrics, determine noise source characteristics and understand various design parameter sensitivities. A new formulation is developed to approximately correct computed noise intensities for any throttle adjustments required to maintain takeoff thrust, and a rough estimate of shielding effectiveness is proposed by means of integration over a spatially distributed noise source. A total of 19 different nozzle geometries are considered in the present study. Results indicate significant noise benefits at takeoff, including a 2 kHz fly-over noise reduction of roughly 2-3 dB, while allowing for less than 15% of the net thrust reduction at supersonic cruise calculated for a fixed-penetration chevron nozzle.

Jet Noise↗