Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence 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 919 records · Page 51

Asymptotic boundary conditions for dissipative waves: General theory

An outstanding issue in the computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗

Asymptotic boundary conditions for dissipative waves - General theory

An outstanding issue in computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗

Improving the efficiency of aerodynamic shape optimization

The computational efficiency of an aerodynamic shape optimization procedure that is based on discrete sensitivity analysis is increased through the implementation of two improvements. The first improvement involves replacing a grid-point-based approach for surface representation with a Bezier-Bernstein polynomial parameterization of the surface. Explicit analytical expressions for the grid sensitivity terms are developed for both approaches. The second improvement proposes the use of Newton's method in lieu of an alternating direction implicit methodology to calculate the highly converged flow solutions that are required to compute the sensitivity coefficients. The modified design procedure is demonstrated by optimizing the shape of an internal-external nozzle configuration. Practically identical optimization results are obtained that are independent of the method used to represent the surface. A substantial factor of 8 decrease in computational time for the optimization process is achieved by implementing both of the design procedure improvements.

Burgreen, Greg W.↗

Design for Manufacturing of Structures with Automated Fiber Placement via Integration of Analysis and Process Planning

Under the Design for Manufacturing (DFM) task in the NASA HiCAM program, significant process has been made towards establishing a fully automated optimization process that spans the structural analysis, design, and manufacturing process planning for Automated Fiber Placement (AFP) structures. Previous efforts in this area established the data formats and mapping processes needed to exchange data between disciplines. The software tools used in the framework are HyperX (structural optimization), CAPP (process planning), and VCP (AFP path generation). The work has culminated with automation of data exchanges and an optimization process that drives the three software toward convergence. Within this framework, HyperX is used to generate mass-optimum composite designs which are iterated with VCP and CAPP to improve manufacturability of the designs. As the ply manufacturing process is adjusted by CAPP and VCP, HyperX continually re-assesses the structural integrity of the part and makes adjustments to the laminate as needed. For example, the position of fiber steering and tow overlaps/gaps can necessitate changes to the laminate design. Once the design converges, the end result is a stack of plies that meets both structural and manufacturing requirements. This paper presents the details of the framework and demonstrations the automated process on several parts with significant double-curvature.

HiCAM↗

Structural analysis of hollow blades: Torsional stress analysis of hollow fan blades for aircraft jet engines

A torsional stress analysis of hollow fans blades by the finite element method is presented. The fans are considered to be double circular arc blades, hollowed 30 percent, and twisted by a component of the centrifugal force by the rated revolution. The effects of blade hollowing on strength and rigidity are discussed. The effects of reinforcing webs, placed in the hollowed section in varying numbers and locations, on torsional rigidity and the convergence of stresses, are reported. A forecast of the 30 percent hollowing against torsional loadings is discussed.

Ogawa, A.↗

Meteorological Analysis of Icing Conditions Encountered in Low-Altitude Stratiform Clouds

Liquid-water content, droplet size, and temperature data measured during 22 flights in predominatly stratiform clouds through the 1948-49 and the 1949-50 winters are presented. Several icing encounters were of greater severity than those previously measured over the same geographical area, but were within the limits of similar measurements obtained over different terrain within the United States. An analysis of meteorological conditions existing during the 74 flights conducted for four winters indicated an inverse relation of liquid-water concentration to maximum horizontal extent of icing clouds. Data on the vertical extent of supercooled clouds are also presented. Icing conditions were most likely to occur in the southwest and northwest quadrants of a cyclone area, and least likely to occur in the southeast and northeast quadrants where convergent air flow and lifting over the associated warm frontal surface usually cause precipitation. Additional data indicated that, icing conditions were usually encountered in nonprecipitating clouds existing at subfreezing temperatures and were unlikely over areas where most weather observing stations reported the existence of precipitation. Measurements of liquid-water content obtained during 12 flights near the time and location of radiosonde observations were compared with theoretical values. The average liquid-water content of a cloud layer, as measured by the multicylinder technique, seldom exceeded two-thirds of that which could be released by adiabatic lifting. Local areas near the cloud tops equaled or occasionally exceeded the calculated maximum quantity of liquid water.

Kline, D. B.↗

Analysis of airfoil transitional separation bubbles

A previously developed local inviscid-viscous interaction technique for the analysis of airfoil transitional separation bubbles, ALESEP (Airfoil Leading Edge Separation) has been modified to utilize a more accurate windward finite difference procedure in the reversed flow region, and a natural transition/turbulence model has been incorporated for the prediction of transition within the separation bubble. Numerous calculations and experimental comparisons are presented to demonstrate the effects of the windward differencing scheme and the natural transition/turbulence model. Grid sensitivity and convergence capabilities of this inviscid-viscous interaction technique are briefly addressed. A major contribution of this report is that with the use of windward differencing, a second, counter-rotating eddy has been found to exist in the wall layer of the primary separation bubble.

Davis, R. L.↗

Developing and utilizing an Euler computational method for predicting the airframe/propulsion effects for an aft-mounted turboprop transport. Volume 1: Theory document

An Euler flow solver was developed for predicting the airframe/propulsion integration effects for an aft-mounted turboprop transport. This solver employs a highly efficient multigrid scheme, with a successive mesh-refinement procedure to accelerate the convergence of the solution. A new dissipation model was also implemented to render solutions that are grid insensitive. The propeller power effects are simulated by the actuator disk concept. An embedded flow solution method was developed for predicting the detailed flow characteristics in the local vicinity of an aft-mounted propfan engine in the presence of a flow field induced by a complete aircraft. Results from test case analysis are presented. A user's guide for execution of computer programs, including format of various input files, sample job decks, and sample input files, is provided in an accompanying volume.

Chen, H. C.↗

An improved fully stressed design algorithm for plate/shell structures

The stress-ratio algorithm associated with fully stressed design philosophy has been used as a convenient tool to achieve minimum weight design of strength-limited structures. The algorithm is effective and converges quickly for many cases. However, it presents extremely slow oscillatory iteration histories for plate-thickness design problems that involve transverse bending loads. Modification of the basic algorithm presented in this paper provides an effective remedy to this problem when both membrane and bending loads are present. The modified resizing algorithm requires numerical solutions of a fourth-order algebraic equation. No additional data, beyond the ordinary static analysis results, are required.

Miura, H.↗

Surface wind convergence as a short-term predictor of cloud-to-ground lightning at Kennedy Space Center

Cloud-to-ground lightning is a significant forecast problem at the Kennedy Space Center (KSC) in Florida. In this study, cloud-to-ground lightning is related in time and space to surface convergence for 244 days during the convective seasons of 1985 and 1986 over a 790 sq-km network at KSC. The method uses surface convergence, particularly the average over the area, to identify the potential for new, local thunderstorm growth, and it can be used to specify the likely time and location of lightning during the life cycle of the convection. A threshold of 75 x 10 to the -6th/s change in divergence is the main criterion used to define a convergence event, and a set of flashes less than 30 min apart defines a lightning event. Time intervals are found from the study to be approximately 1 h from beginning convergence to first flash, and another hour to the end of lightning. The influences of low-level winds and midlevel moisture in determining the location and intensity of convection are discussed. This is the first known dynamically-based forecast method for lightning prediction. The technique, currently in use at KSC, has been shown to be a systematic, quantitative tool for predicting lightning onset in situations where conventional analysis tools such as radar and satellite are limited.

Watson, Andrew I.↗

Improving the efficiency of aerodynamic shape optimization procedures

The computational efficiency of an aerodynamic shape optimization procedure which is based on discrete sensitivity analysis is increased through the implementation of two improvements. The first improvement involves replacing a grid point-based approach for surface representation with a Bezier-Bernstein polynomial parameterization of the surface. Explicit analytical expressions for the grid sensitivity terms are developed for both approaches. The second improvement proposes the use of Newton's method in lieu of an alternating direction implicit (ADI) methodology to calculate the highly converged flow solutions which are required to compute the sensitivity coefficients. The modified design procedure is demonstrated by optimizing the shape of an internal-external nozzle configuration. A substantial factor of 8 decrease in computational time for the optimization process was achieved by implementing both of the design improvements.

Burgreen, Greg W.↗

Maximum von Mises Stress in the Loading Environment of Mass Acceleration Curve

Method for calculating stress due to acceleration loading: 1) Part has been designed by FEA and hand calculation in one critical loading direction judged by the analyst; 2) Maximum stress can be due to loading in another direction; 3) Analysis procedure to be presented determines: a) The maximum Mises stress at any point; and b) The direction of maximum loading associated with the "stress". Concept of Mass Acceleration Curves (MAC): 1) Developed by JPL to perform preliminary structural sizing (i.e. Mariners, Voyager, Galileo, Pathfinder, MER,...MSL); 2) Acceleration of physical masses are bounded by a curve; 3) G-levels of vibro-acoustic and transient environments; 4) Convergent process before the couple loads cycle; and 5) Semi-empirical method to effectively bound the loads, not a simulation of the actual response.

von Mises stress↗

High Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

Taylor series integration is implemented in a spacecraft trajectory analysis code-the Spacecraft N-body Analysis Program (SNAP) - and compared with the code s existing eighth-order Runge-Kutta Fehlberg time integration scheme. Nine trajectory problems, including near Earth, lunar, Mars and Europa missions, are analyzed. Head-to-head comparison at five different error tolerances shows that, on average, Taylor series is faster than Runge-Kutta Fehlberg by a factor of 15.8. Results further show that Taylor series has superior convergence properties. Taylor series integration proves that it can provide rapid, highly accurate solutions to spacecraft trajectory problems.

Scott, James R.↗

An Analysis of the Origin and Propagation of the Multiple Coronal Mass Ejections of 2010 August 1

On 2010 August 1, the northern solar hemisphere underwent significant activity that involved a complex set of active regions near central meridian with, nearby, two large prominences and other more distant active regions. This activity culminated in the eruption of four major coronal mass ejections (CMEs), effects of which were detected at Earth and other solar system bodies. Recognizing the unprecedented wealth of data from the wide range of spacecraft that were available-providing the potential for us to explore methods for CME identification and tracking, and to assess issues regarding onset and planetary impact-we present a comprehensive analysis of this sequence of CMEs.We show that, for three of the four major CMEs, onset is associated with prominence eruption, while the remaining CME appears to be closely associated with a flare. Using instrumentation on board the Solar Terrestrial Relations Observatory spacecraft, three of the CMEs could be tracked out to elongations beyond 50◦; their directions and speeds have been determined by various methods, not least to assess their potential for Earth impact. The analysis techniques that can be applied to the other CME, the first to erupt, are more limited since that CME was obscured by the subsequent, much faster event before it had propagated far from the Sun; we discuss the speculation that these two CMEs interact. The consistency of the results, derived from the wide variety of methods applied to such an extraordinarily complete data set, has allowed us to converge on robust interpretations of the CME onsets and their arrivals at 1 AU.

Harrison, R. A.↗

2D axisymmetric analysis of SRM ignition transient

To analyze ignition transient of Space Shuttle solid rocket motor, a transient two-dimensional numerical model based on turbulent compressible Navier-Stokes equations in a generalized coordinate system was developed. One-dimensional numerical models (Peretz et al., 1973; Han, 1992; Pardue and Han, 1992) with empirical correlations data obtained from steady turbulent boundary layer flows agrees reasonably well with test rocket data by adjusting a few parameters. However, a 1D model can not provide a physical insight into the complex multidimensional thermal fields and flowfields in the chamber and the converging-diverging rocket nozzle. As an interim step, a 2D model was developed and compared with test data. A modified version of SIMPLE algorithm was used for the numerical model, and the standard k-epsilon model with a wall function was used for turbulence closure. Transient flowfields and thermal fields in the combustion chamber and the attached nozzle were obtained for a selected rocket geometry and propellant. Transient behaviors of the flow and thermal fields were analyzed, and were found to be in good agreement with physical expectations.

Bai, S. D.↗

Buckling of simply supported plate assemblies subject to shear loading

An analysis is presented of the shear buckling load of an arbitrary simply-supported prismatic plate assembly or stiffened panel configuration. The method utilizes stiffness matrices developed by Wittrick (1968). For the restraint of simple support along edges perpendicular to the longitudinal direction, various wavelengths of response are combined using an energy method that involves only a boundary integral. Numerical results for plate and panel configurations show the accuracy and convergence of the method.

Anderson, M. S.↗

On the Hilbert-Huang Transform Theoretical Developments

One of the main heritage tools used in scientific and engineering data spectrum analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as linearity, of being stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectrum analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposition data, the HHT allows spectrum analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a near orthogonal adaptive basis, a basis that is derived from the data. The IMFs can be further analyzed for spectrum interpretation by the classical Hilbert Transform. A new engineering spectrum analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications post additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs near orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the developments of new HHT processing options, such as real-time and 2-D processing using Field Programmable Array (FPGA) computational resources, enhanced HHT synthesis, and broaden the scope of HHT applications for signal processing.

Kizhner, Semion↗

On Certain Theoretical Developments Underlying the Hilbert-Huang Transform

One of the main traditional tools used in scientific and engineering data spectral analysis is the Fourier Integral Transform and its high performance digital equivalent - the Fast Fourier Transform (FFT). Both carry strong a-priori assumptions about the source data, such as being linear and stationary, and of satisfying the Dirichlet conditions. A recent development at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC), known as the Hilbert-Huang Transform (HHT), proposes a novel approach to the solution for the nonlinear class of spectral analysis problems. Using a-posteriori data processing based on the Empirical Mode Decomposition (EMD) sifting process (algorithm), followed by the normalized Hilbert Transform of the decomposed data, the HHT allows spectral analysis of nonlinear and nonstationary data. The EMD sifting process results in a non-constrained decomposition of a source real-value data vector into a finite set of Intrinsic Mode Functions (IMF). These functions form a nearly orthogonal derived from the data (adaptive) basis. The IMFs can be further analyzed for spectrum content by using the classical Hilbert Transform. A new engineering spectral analysis tool using HHT has been developed at NASA GSFC, the HHT Data Processing System (HHT-DPS). As the HHT-DPS has been successfully used and commercialized, new applications pose additional questions about the theoretical basis behind the HHT and EMD algorithms. Why is the fastest changing component of a composite signal being sifted out first in the EMD sifting process? Why does the EMD sifting process seemingly converge and why does it converge rapidly? Does an IMF have a distinctive structure? Why are the IMFs nearly orthogonal? We address these questions and develop the initial theoretical background for the HHT. This will contribute to the development of new HHT processing options, such as real-time and 2-D processing using Field Programmable Gate Array (FPGA) computational resources,

Kizhner, Semion↗