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 181 records · Page 10

A numerical method for computing three dimensional viscous supersonic flow fields about slender bodies

A numerical method is presented for the calculation of steady, three-dimensional, viscous, compressible flow fields about slender bodies at angle of attack and at supersonic speeds. Approximations are introduced in modeling the flow in the longitudinal direction. Accordingly, the flow fields calculated with the program were computed with a model that permits viscous crossflow together with inviscid axial flow. An analysis of the errors introduced by such a treatment is presented. Numerical calculations were made and compared with experimental results for an ogive-cylinder and an airplane fuselage configuration. Generally, good agreement with experiment was obtained. However, boundary layer separation and body vortex positions differed from experimental locations on the ogive-cylinder, and the shock induced by the fuselage canopy was predicted at a slightly different location.

Walitt, L.↗

Turbulent jet mixing at high supersonic speeds.

The results of three related studies of turbulent jet mixing are presented and discussed. The major effort was the collection of mean-flow data at a much higher free stream Mach number than was previously available. The effects of initial swirl in the injectant were found to be minimal. Numerical predictions of the non-swirling experimental case obtained using an eddy viscosity approximation are found to be in reasonable agreement with the data. Lastly, an approximate analysis of the decay of the swirling velocity is developed. This analysis is based upon an eddy viscosity model for transport in the circumferential direction adapted from wing-tip vortex work.

Schetz, J. A.↗

The efficient solution of transonic wing flow fields

An evaluation of the transonic-wing-analysis computer code TWING is presented. TWING utilizes a fully implicit, approximate-factorization iteration scheme to solve the full-potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations have been analyzed, and comparisons of computed results have been made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative, full-potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations, including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.

Holst, T. L.↗

Time Dependence of Collision Probabilities During Satellite Conjunctions

The NASA Conjunction Assessment Risk Analysis (CARA) team has recently implemented updated software to calculate the probability of collision (P (sub c)) for Earth-orbiting satellites. The algorithm can employ complex dynamical models for orbital motion, and account for the effects of non-linear trajectories as well as both position and velocity uncertainties. This “3D P (sub c)” method entails computing a 3-dimensional numerical integral for each estimated probability. Our analysis indicates that the 3D method provides several new insights over the traditional “2D P (sub c)” method, even when approximating the orbital motion using the relatively simple Keplerian two-body dynamical model. First, the formulation provides the means to estimate variations in the time derivative of the collision probability, or the probability rate, R (sub c). For close-proximity satellites, such as those orbiting in formations or clusters, R (sub c) variations can show multiple peaks that repeat or blend with one another, providing insight into the ongoing temporal distribution of risk. For single, isolated conjunctions, R (sub c) analysis provides the means to identify and bound the times of peak collision risk. Additionally, analysis of multiple actual archived conjunctions demonstrates that the commonly used “2D P (sub c)” approximation can occasionally provide inaccurate estimates. These include cases in which the 2D method yields negligibly small probabilities (e.g., P (sub c)) is greater than 10 (sup -10)), but the 3D estimates are sufficiently large to prompt increased monitoring or collision mitigation (e.g., P (sub c) is greater than or equal to 10 (sup -5)). Finally, the archive analysis indicates that a relatively efficient calculation can be used to identify which conjunctions will have negligibly small probabilities. This small-P (sub c) screening test can significantly speed the overall risk analysis computation for large numbers of conjunctions.

Hall, Doyle T.↗

Trajectories for Flyby Sample Return at Saturn's Moons

Ballistic trajectories are computed which would enable a sample return mission to Titan or Enceladus without capturing, descending, or landing. The low-cost mission concept utilizes a free return trajectory that also involves a close flyby of the moon. This work extends the concept, and related trajectory analysis methodology, previously applied to a Europa mission. Specifically, a broad search algorithm is employed to systematically locate potentially feasible itineraries over an entire Saturn period. High-quality approximate solutions are then optimized to be continuous using high-fidelity dynamics. Techniques and software from the Europa analysis, were readily adapted and able to find numerous mission enabling trajectories. A direct mission to Titan is possible with flight time under 16 years and Earth-relative speeds below 11.0 km/sec. The VEEGA option is shown to substantially reduce launch C3, but flight times exceed 21 years. Unfortunately, an Enceladus mission requires a flight time of 25 years or more, and incurs fairly high relative speeds. Nevertheless, an optimized reference mission is computed.

Jones, Drew Ryan↗

Analysis of Numerical Simulation Database for Pressure Fluctuations Induced by High-Speed Turbulent Boundary Layers

Direct numerical simulations (DNS) of Mach 6 turbulent boundary layer with nominal freestream Mach number of 6 and Reynolds number of Re(sub T) approximately 460 are conducted at two wall temperatures (Tw/Tr = 0.25, 0.76) to investigate the generated pressure fluctuations and their dependence on wall temperature. Simulations indicate that the influence of wall temperature on pressure fluctuations is largely limited to the near-wall region, with the characteristics of wall-pressure fluctuations showing a strong temperature dependence. Wall temperature has little influence on the propagation speed of the freestream pressure signal. The freestream radiation intensity compares well between wall-temperature cases when normalized by the local wall shear; the propagation speed of the freestream pressure signal and the orientation of the radiation wave front show little dependence on the wall temperature.

Duan, Lian↗

Thermal stresses in a spherical pressure vessel having temperature-dependent, transversely isotropic, elastic properties

Rayleigh-Ritz and modified Rayleigh-Ritz procedures are used to construct approximate solutions for the response of a thick-walled sphere to uniform pressure loads and an arbitrary radial temperature distribution. The thermoelastic properties of the sphere are assumed to be transversely isotropic and nonhomogeneous; variations in the elastic stiffness and thermal expansion coefficients are taken to be an arbitrary function of the radial coordinate and temperature. Numerical examples are presented which illustrate the effect of the temperature-dependence upon the thermal stress field. A comparison of the approximate solutions with a finite element analysis indicates that Ritz methods offer a simple, efficient, and relatively accurate approach to the problem.

Tauchert, T. R.↗

Analysis of solar radiation pressure induced coupled librations of a gravity stabilized axisymmetric satellite

This paper presents an analysis of solar radiation pressure induced coupled librations of gravity stabilized cylindrical spacecraft with a special reference to geostationary communication satellites. The Lagrangian approach is used to obtain the corresponding equations of motion. The solar induced torques are assumed to be free of librational angles and are represented by their Fourier expansion. The response and periodic solutions are obtained through linear and nonlinear analyses, using the method of harmonic balance in the latter case. The stability conditions are obtained using Routh-Hurwitz criteria. To establish the ranges of validity the analytic response is compared with the numerical solution. Finally, values of the system parameters are suggested to make the satellite behave as desired. Among these is a possible approach to subdue the solar induced roll resonance. It is felt that the approximate analysis presented here should significantly reduce the computational efforts involved in the design and stability analysis of the systems.

Shrivastava, S. K.↗

A new methodology for free wake analysis using curved vortex elements

A method using curved vortex elements was developed for helicopter rotor free wake calculations. The Basic Curve Vortex Element (BCVE) is derived from the approximate Biot-Savart integration for a parabolic arc filament. When used in conjunction with a scheme to fit the elements along a vortex filament contour, this method has a significant advantage in overall accuracy and efficiency when compared to the traditional straight-line element approach. A theoretical and numerical analysis shows that free wake flows involving close interactions between filaments should utilize curved vortex elements in order to guarantee a consistent level of accuracy. The curved element method was implemented into a forward flight free wake analysis, featuring an adaptive far wake model that utilizes free wake information to extend the vortex filaments beyond the free wake regions. The curved vortex element free wake, coupled with this far wake model, exhibited rapid convergence, even in regions where the free wake and far wake turns are interlaced. Sample calculations are presented for tip vortex motion at various advance ratios for single and multiple blade rotors. Cross-flow plots reveal that the overall downstream wake flow resembles a trailing vortex pair. A preliminary assessment shows that the rotor downwash field is insensitive to element size, even for relatively large curved elements.

Bliss, Donald B.↗

A rotor optimization using regression analysis

The design and development of helicopter rotors is subject to the many design variables and their interactions that effect rotor operation. Until recently, selection of rotor design variables to achieve specified rotor operational qualities has been a costly, time consuming, repetitive task. For the past several years, Kaman Aerospace Corporation has successfully applied multiple linear regression analysis, coupled with optimization and sensitivity procedures, in the analytical design of rotor systems. It is concluded that approximating equations can be developed rapidly for a multiplicity of objective and constraint functions and optimizations can be performed in a rapid and cost effective manner; the number and/or range of design variables can be increased by expanding the data base and developing approximating functions to reflect the expanded design space; the order of the approximating equations can be expanded easily to improve correlation between analyzer results and the approximating equations; gradients of the approximating equations can be calculated easily and these gradients are smooth functions reducing the risk of numerical problems in the optimization; the use of approximating functions allows the problem to be started easily and rapidly from various initial designs to enhance the probability of finding a global optimum; and the approximating equations are independent of the analysis or optimization codes used.

Giansante, N.↗

Zero potential vorticity envelopes for the zonal-mean velocity of the Venus/Titan atmospheres

The diagnostic analysis of numerical simulations of the Venus/Titan wind regime reveals an overlooked constraint upon the latitudinal structure of their zonal-mean angular momentum. The numerical experiments, as well as the limited planetary observations, are approximately consistent with the hypothesis that within the latitudes bounded by the wind maxima the total Ertel potential vorticity associated with the zonal-mean motion is approximately well mixed with respect to the neutral equatorial value for a stable circulation. The implied latitudinal profile of angular momentum is of the form M equal to or less than M(sub e)(cos lambda)(exp 2/Ri), where lambda is the latitude and Ri the local Richardson number, generally intermediate between the two extremes of uniform angular momentum (Ri approaches infinity) and uniform angular velocity (Ri = 1). The full range of angular momentum profile variation appears to be realized within the observed meridional - vertical structure of the Venus atmosphere, at least crudely approaching the implied relationship between stratification and zonal velocity there. While not itself indicative of a particular eddy mechanism or specific to atmospheric superrotation, the zero potential vorticity (ZPV) constraint represents a limiting bound for the eddy - mean flow adjustment of a neutrally stable baroclinic circulation and may be usefully applied to the diagnostic analysis of future remote sounding and in situ measurements from planetary spacecraft.

Allison, Michael↗

Model-size reduction technique for the analysis of symmetric anisotropic structures

A two-step computational procedure is presented for reducing the size of the analysis model for an anisotropic symmetric structure to that of the corresponding orthotropic structure. The key elements of the procedure are: (1) decomposition of the stiffness matrix into the sum of an orthotropic and nonorthotropic (anisotropic) parts; and (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The global approximation vectors are selected to be the solution corresponding to zero nonorthotropic matrix and its various-order derivatives with respect to an anisotropic tracing parameter (identifying the nonorthotropic material coefficients). The size of the analysis model used in generating the global approximation vectors is identical to that of the corresponding orthotropic structure. The effectiveness of the proposed technique is demonstrated by means of numerical examples and its potential for solving other quasi-symmetric problems is discussed.

Noor, A. K.↗

Improved 3-D turbomachinery CFD algorithm

The building blocks of a computer algorithm developed for the time-accurate flow analysis of rotating machines are described. The flow model is a finite volume method utilizing a high resolution approximate Riemann solver for interface flux definitions. This block LU implicit numerical scheme possesses apparent unconditional stability. Multi-block composite gridding is used to orderly partition the field into a specified arrangement. Block interfaces, including dynamic interfaces, are treated such as to mimic interior block communication. Special attention is given to the reduction of in-core memory requirements by placing the burden on secondary storage media. Broad applicability is implied, although the results presented are restricted to that of an even blade count configuration. Several other configurations are presently under investigation, the results of which will appear in subsequent publications.

Janus, J. Mark↗

The computation of induced drag with nonplanar and deformed wakes

The classical calculation of inviscid drag, based on far field flow properties, is reexamined with particular attention to the nonlinear effects of wake roll-up. Based on a detailed look at nonlinear, inviscid flow theory, it is concluded that many of the classical, linear results are more general than might have been expected. Departures from the linear theory are identified and design implications are discussed. Results include the following: Wake deformation has little effect on the induced drag of a single element wing, but introduces first order corrections to the induced drag of a multi-element lifting system. Far field Trefftz-plane analysis may be used to estimate the induced drag of lifting systems, even when wake roll-up is considered, but numerical difficulties arise. The implications of several other approximations made in lifting line theory are evaluated by comparison with more refined analyses.

Kroo, Ilan↗

The computation of induced drag with nonplanar and deformed wakes

The classical calculation of inviscid drag, based on far field flow properties, is re-examined with particular attention to the nonlinear effects of wake roll-up. Based on a detailed look at nonlinear, inviscid flow theory, the paper concludes that many of the classical, linear results are more general than might have been expected. Departures from the linear theory are identified and design implications are discussed. Results include the following: wake deformation has little effect on the induced drag of a single element wing, but introduces first order corrections to the induced drag of a multi-element lifting system. Far field Trefftz-plane analysis may be used to estimate the induced drag of lifting systems, even when wake roll-up is considered, but numerical difficulties arise. The implications of several other approximations made in lifting line theory are evaluated by comparison with more refined analyses.

Kroo, Ilan↗

Analysis and control of hourglass instabilities in underintegrated linear and nonlinear elasticity

Methods are described to identify and correct a bad finite element approximation of the governing operator obtained when under-integration is used in numerical code for several model problems: the Poisson problem, the linear elasticity problem, and for problems in the nonlinear theory of elasticity. For each of these problems, the reason for the occurrence of instabilities is given, a way to control or eliminate them is presented, and theorems of existence, uniqueness, and convergence for the given methods are established. Finally, numerical results are included which illustrate the theory.

Jacquotte, Olivier P.↗

Investigation of MMOD Impact on STS-115 Shuttle Payload Bay Door Radiator

The Orbiter radiator system consists of eight individual 4.6 m x 3.2 m panels located with four on each payload bay door. Forward panels #1 and #2 are 2.3 cm thick while the aft panels #3 and #4 have a smaller overall thickness of 1.3 cm. The honeycomb radiator panels consist of 0.028 cm thick Aluminum 2024-T81 facesheets and Al5056-H39 cores. The face-sheets are topped with 0.005 in. (0.127 mm) silver-Teflon tape. The radiators are located on the inside of the shuttle payload bay doors, which are closed during ascent and reentry, limiting damage to the on-orbit portion of the mission. Post-flight inspections at the Kennedy Space Center (KSC) following the STS-115 mission revealed a large micrometeoroid/orbital debris (MMOD) impact near the hinge line on the #4 starboard payload bay door radiator panel. The features of this impact make it the largest ever recorded on an orbiter payload bay door radiator. The general location of the damage site and the adjacent radiator panels can be seen in Figure 2. Initial measurements of the defect indicated that the hole in the facesheet was 0.108 in. (2.74 mm) in diameter. Figure 3 shows an image of the front side damage. Subsequent observations revealed exit damage on the rear facesheet. Impact damage features on the rear facesheet included a 0.03 in. diameter hole (0.76 mm), a approx.0.05 in. tall bulge (approx.1.3 mm), and a larger approx.0.2 in. tall bulge (approx.5.1 mm) that exhibited a crack over 0.27 in. (6.8 mm) long. A large approx.1 in. (25 mm) diameter region of the honeycomb core was also damaged. Refer to Figure 4 for an image of the backside damage to the panel. No damage was found on thermal blankets or payload bay door structure under the radiator panel. Figure 5 shows the front facesheet with the thermal tape removed. Ultrasound examination indicated a maximum facesheet debond extent of approximately 1 in. (25 mm) from the entry hole. X-ray examinations revealed damage to an estimated 31 honeycomb cells with an extent of 0.85 in. x 1.1 in. (21.6 x 27.9 mm). Pieces of the radiator at and surrounding the impact site were recovered during the repair procedures at KSC. They included the thermal tape, front facesheet, honeycomb core, and rear facesheet. These articles were examined at JSC using a scanning electron microscope (SEM) with an energy dispersive x-ray spectrometer (EDS). Figure 6 shows SEM images of the entry hole in the facesheet. The asymmetric height of the lip may be attributed to projectile shape and impact angle. Numerous instances of a glass-fiber organic matrix composite were observed in the facesheet tape sample. The fibers were approximately 10 micrometers in diameter and variable lengths. EDS analysis indicated a composition of Mg, Ca, Al, Si, and O. Figures 7 and 8 present images of the fiber bundles, which were believed to be circuit board material based on similarity in fiber diameter, orientation, consistency, and composition. A test program was initiated in an attempt to simulate the observed damage to the radiator facesheet and honeycomb. Twelve test shots were performed using projectiles cut from a 1.6 mm thick fiberglass circuit board substrate panel. Results from test HITF07017, shown in figures 9 and 10, correlates with the observed impact features reasonably well. The test was performed at 4.14 km/sec with an impact angle of 45 degrees using a cylindrical projectile with a diameter and length of 1.25 mm. The fiberglass circuit board material had a density of 1.65 g/cu cm, giving a projectile mass of 2.53 mg. An analysis was performed using the Bumper code to estimate the probability of impact to the shuttle from a 1.25 mm diameter particle. Table 1 shows a 1.6% chance (impact odds = 1 in 62) of a 1.25 mm or larger MMOD impact on the radiators of the vehicle during a typical ISS mission. There is a 0.4% chance (impact odds = 1 in 260) that a 1.25 mm or larger MMOD particle would impact the RCC wing leading edge and nose cap during a typical miion. Figure 11 illustrates the vulnerable areas of the wing leading edge reinforced carbon-carbon (RCC), an area of the vehicle that is very sensitive to impact damage. The highlighted red, orange, yellow, and light green areas would be expected to experience critical damage if impacted by an OD particle such as the one that hit the RH4 radiator panel on STS-115.

Hyde, J.↗

A New Higher-Order Composite Theory for Analysis and Design of High Speed Tilt-Rotor Blades

A higher-order theory is developed to model composite box beams with arbitrary wall thicknesses. The theory, based on a refined displacement field, represents a three-dimensional model which approximates the elasticity solution. Therefore, the cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are automatically included in the formulation. The model accurately captures the transverse shear stresses through the thickness of each wall while satisfying all stress-free boundary conditions. Several numerical results are presented to validate the present theory. The developed theory is then used to model the load carrying member of a tilt-rotor blade which has thick-walled sections. The composite structural analysis is coupled with an aerodynamic analysis to compute the aeroelastic stability of the blade. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt-rotor aircraft. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem and a hybrid approximate analysis is used to reduce the computational effort. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt-rotor blade.

McCarthy, Thomas Robert↗