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 865 records · Page 48

Current Grid Generation Strategies and Future Requirements in Hypersonic Vehicle Design, Analysis and Testing

Recent advances in computational power enable computational fluid dynamic modeling of increasingly complex configurations. A review of grid generation methodologies implemented in support of the computational work performed for the X-38 and X-33 are presented. In strategizing topological constructs and blocking structures factors considered are the geometric configuration, optimal grid size, numerical algorithms, accuracy requirements, physics of the problem at hand, computational expense, and the available computer hardware. Also addressed are grid refinement strategies, the effects of wall spacing, and convergence. The significance of grid is demonstrated through a comparison of computational and experimental results of the aeroheating environment experienced by the X-38 vehicle. Special topics on grid generation strategies are also addressed to model control surface deflections, and material mapping.

Papadopoulos, Periklis↗

Shear-flexible finite-element models of laminated composite plates and shells

Several finite-element models are applied to the linear static, stability, and vibration analysis of laminated composite plates and shells. The study is based on linear shallow-shell theory, with the effects of shear deformation, anisotropic material behavior, and bending-extensional coupling included. Both stiffness (displacement) and mixed finite-element models are considered. Discussion is focused on the effects of shear deformation and anisotropic material behavior on the accuracy and convergence of different finite-element models. Numerical studies are presented which show the effects of increasing the order of the approximating polynomials, adding internal degrees of freedom, and using derivatives of generalized displacements as nodal parameters.

Noor, A. K.↗

Hybrid analytical technique for the nonlinear analysis of curved beams

The application of a two-step hybrid technique to the geometrically nonlinear analysis of curved beams is used to demonstrate the potential of hybrid analytical techniques in nonlinear structural mechanics. The hybrid technique is based on successive use of the perturbation method and a classical direct variational procedure. The functions associated with the various-order terms in the perturbation expansion of the fundamental unknowns, and their sensitivity derivatives with respect to material and geometric parameters of the beam, are first obtained by using the perturbation method. These functions are selected as coordinate functions (or modes) and the classical direct variational technique is then used to compute their amplitudes. The potential of the proposed hybrid technique for nonlinear analysis of structures is discussed. The effectiveness of the hybrid technique is demonstrated by means of numerical examples. The symbolic computation system Mathematica is used in the present study. The tasks performed on Mathematica include: (1) generation of algebraic expressions for the perturbation functions of the different response quantities and their sensitivity derivatives: and (2) determination of the radius of convergence of the perturbation series.

Noor, A. K.↗

Multigrid schemes for viscous hypersonic flows

Several multigrid schemes are considered for the numerical computation of viscous hypersonic flows. For each scheme, the basic solution algorithm employs upwind spatial discretization with explicit multistage time stepping. Two-level versions of the various multigrid algorithms are applied to the two-dimensional advection equation, and Fourier analysis is used to determine their damping properties. The capabilities of the multigrid methods are assessed by solving two different hypersonic flow problems. Some new multigrid schemes, based on semicoarsening strategies, are shown to be quite effective in relieving the stiffness caused by the high-aspect-ratio cells required to resolve high Reynolds number flows. These schemes exhibit good convergence rates for Reynolds numbers up to 200 x 10(exp 6).

Swanson, R. C.↗

Shocked clouds in the Cygnus Loop

This grant covers the analysis of ROSAT PSPC and HRI images of the Cygnus Loop, an elderly supernova remnant. The project, as proposed, includes not only the usual analysis of ROSAT data; the ROSAT data is being combined with optical and UV data, and new model calculations are being performed. The status is reported on optical imagery, echelle data, IUE data, ROSAT data, and the grain model. The major question being addressed is whether the blastwave-cloud interaction in the feature known as XA is basically a converging shock in a fairly large cloud or turbulent stripping of material from the edges of a smaller cloud.

Raymond, John C.↗

Evaluating the Large-Scale Hydrological Cycle Response within the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2) Ensemble

The mid-Pliocene (∼3 Ma) is one of the most recent warm periods with high CO2 concentrations in the atmosphere and resulting high temperatures, and it is often cited as an analog for near-term future climate change. Here, we apply a moisture budget analysis to investigate the response of the large-scale hydrological cycle at low latitudes within a 13-model ensemble from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The results show that increased atmospheric moisture content within the mid-Pliocene ensemble (due to the thermodynamic effect) results in wetter conditions over the deep tropics, i.e., the Pacific intertropical convergence zone (ITCZ) and the Maritime Continent, and drier conditions over the subtropics. Note that the dynamic effect plays a more important role than the thermodynamic effect in regional precipitation minus evaporation (PmE) changes (i.e., northward ITCZ shift and wetter northern Indian Ocean). The thermodynamic effect is offset to some extent by a dynamic effect involving a northward shift of the Hadley circulation that dries the deep tropics and moistens the subtropics in the Northern Hemisphere (i.e., the subtropical Pacific). From the perspective of Earth's energy budget, the enhanced southward cross-equatorial atmospheric transport (0.22 PW), induced by the hemispheric asymmetries of the atmospheric energy, favors an approximately 1∘ northward shift of the ITCZ. The shift of the ITCZ reorganizes atmospheric circulation, favoring a northward shift of the Hadley circulation. In addition, the Walker circulation consistently shifts westward within PlioMIP2 models, leading to wetter conditions over the northern Indian Ocean. The PlioMIP2 ensemble highlights that an imbalance of interhemispheric atmospheric energy during the mid-Pliocene could have led to changes in the dynamic effect, offsetting the thermodynamic effect and, hence, altering mid-Pliocene hydroclimate.

mid-Pliocene↗

Integration of Structural Analysis and Manufacturing Process Planning for Global Optimization with Automated Fiber Placement

Design of mass-efficient composite structures intended for Automated Fiber Placement (AFP) requires close interaction between structural analysis and manufacturing process planning. Tools exist for each of these disciplines, but software interplay has been insufficient for rapid and efficient design iteration. Within the NASA Advanced Composites Consortium (ACC), the Design for Manufacturing (DFM) task has made significant progress towards linking these disciplines and respective software – HyperX (design), CAPP (process planning), and VCP (tool path generation). The initial focus in previous work was on data exchange between disciplines. The ability to both export and consume composite design and manufacturing data to and from each tool. This paper focuses on the effort to automate and streamline the connection between the tools listed above, with the goal of being able to automatically generate a composite AFP design that is mass-efficient and manufacturable. The optimization method being pursued is a bi-level approach, where each tool performs optimization within its discipline. The optimization in HyperX is focused on mass and laminate strength, while CAPP is focused on maximizing manufacturability. VCP is used to generate fiber paths for each design iteration. These sub-processes are wrapped with a global level optimization, driven by HyperX, used to converge the design. This paper describes the current state of this effort, which is a completed HyperX-VCP iteration loop and initial work on the HyperX-CAPP iteration loop. Additionally, example results are shown for a wind blade structure with double curvature.

Automated Fiber Placement↗

Theory and implementation of high-order adaptive hp methods for analysis of incompressible viscous flows

An account is given of 'smart' algorithms for CFD which change in structure and performance with during flow calculations to accommodate changing properties of the solution. Such algorithms prominently include adaptive FEM methods, which are designed to adjust mesh parameters for the control of numerical error; attention is presently given to those which change the mesh size h and the local spectral order p in order to achieve high accuracies with minimal numbers of degrees of freedom. These 'hp methods' produce exponentially convergent approximations through which flow features are resolved by automatically distributing element sizes and spectral orders. This leads to calculation of the local (elementwise) error in the approximation.

Oden, J. T.↗

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.↗

Processes of Equatorial Thermal Structure: An Analysis of Galileo Temperature Profile with 3-D Model

The Jupiter Thermosphere General Circulation Model (JTGCM) calculates the global dynamical structure of Jupiter's thermosphere self-consistently with its global thermal structure and composition. The main heat source that drives the thermospheric flow is high-latitude Joule heating. A secondary source of heating is the auroral process of particle precipitation. Global simulations of Jovian thermospheric dynamics indicate strong neutral outflows from the auroral ovals with velocities up to approximately 2 kilometers per second and subsequent convergence and downwelling at the Jovian equator. Such circulation is shown to be an important process for transporting significant amounts of auroral energy to equatorial latitudes and for regulating the global heat budget in a manner consistent with the high thermospheric temperatures observed by the Galileo probe. Adiabatic compression of the neutral atmosphere resulting from downward motion is an important source of equatorial heating (less than 0.06 microbar). The adiabatic heating continues to dominate between 0.06 and 0.2 microbar, but with an addition of comparable heating due to horizontal advection induced by the meridional flow. Thermal conduction plays an important role in transporting heat down to lower altitudes (greater than 0.2microbar) where it is balanced by the cooling associated with the wind transport processes. Interestingly, we find that radiative cooling caused by H3(+), CH4, and C2H2 emissions does not play a significant role in interpreting the Galileo temperature profile.

Majeed, T.↗

Characterization and Detailed Analysis of Regression Behavior for HTPB Solid Fuels Containing High Aluminum Loadings

NASA Marshall Space Flight Center's Materials and Processes Department, with support from the Propulsion Systems Department, has renewed the development and maintenance of a hybrid test bed for exposing ablative thermal protection materials to an environment similar to that seen in solid rocket motors (SRM). The Solid Fuel Torch (SFT), operated during the Space Shuttle program, utilized gaseous oxygen for oxidizer and an aluminized hydroxyl-terminated polybutadiene (HTPB) fuel grain to expose a converging section of phenolic material to a 400 psi, 2-phase flow combustion environment. The configuration allows for up to a 2 foot long, 5 inch diameter fuel grain cartridge. Wanting to now test rubber insulation materials with a turn-back feature to mimic the geometry of an aft dome being impinged by alumina particles, the throat area has now been increased by several times to afford flow similarity. Combined with the desire to maintain a higher operating pressure, the oxidizer flow rate is being increased by a factor of 10. Out of these changes has arisen the need to characterize the fuel/oxidizer combination in a higher mass flux condition than has been previously tested at MSFC, and at which the literature has little to no reporting as well. For (especially) metalized fuels, hybrid references have pointed out possible dependence of fuel regression rate on a number of variables: mass flux, G - oxidizer only (G0), or - total mass flux (Gtot), Length, L, Pressure, P, and Diameter, D.

Kibbey, Timothy P.↗

Nozzle flow of laser-heated radiating hydrogen with application to a laser-heated rocket

This paper presents a model for the steady heating of flowing hydrogen by a CW 10.6 micron laser, to consider the feasibility of a laser-heated rocket. The hydrogen flow and the laser beam are parallel, and move into a converging-diverging nozzle. The absorption of laser energy is initiated by a laser-supported combustion wave. The hydrogen is in chemical equilibrium, absorbs laser energy by inverse Bremsstrahlung, and loses energy by radiation. The hydrogen flow was calculated from the rear of the LSC wave to the throat. Estimates of convective heat losses were made using a hydrogen boundary layer analysis. Specific impulse, obtained by expanding isentropically from the throat to 1 atm or a vacuum, varies from 1400 to 3000 s. Radiation losses are 5 to 20%, though the energy fluxes to the walls are quite high. Convective loss estimates are high enough to indicate that coupling to the hot gas flow is required for a 10 kW engine, but not for a 5 MW engine.

Kemp, N. H.↗

Stress-intensity factors for two symmetric corner cracks

This paper presents stress-intensity factors, calculated by a three-dimensional finite-element analysis, for shallow or deep quarter-elliptical corner cracks at the edge of a hole in a finite-thickness plate. The plate was subjected to remote uniform tension, remote bending, or simulated pin loading in the hole; a wide range of configuration parameters was investigated. To verify the accuracy of the three-dimensional finite-element models employed, convergence was studied by varying the numbers of degrees of freedom. The stress-intensity factor variations along the crack front are compared with solutions from the literature.

Raju, I. S.↗

Efficient computation of the stability of three-dimensional compressible boundary layers

Methods for the computer analysis of the stability of three-dimensional compressible boundary layers are discussed and the user-oriented Compressible Stability Analysis (COSAL) computer code is described. The COSAL code uses a matrix finite-difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.

Malik, M. R.↗

Unsteady transonic algorithm improvements for realistic aircraft applications

Improvements to a time-accurate approximate factorization (AF) algorithm have been implemented for steady and unsteady transonic analysis of realistic aircraft configurations. The algorithm improvements include: an Engquist-Osher (E-O) type-dependent switch to more accurately and efficiently treat regions of supersonic flow, extension of the E-O switch for second-order spatial accuracy in these regions, nonreflecting far field boundary conditions for more accurate unsteady applications, and several modifications which accelerate convergence to steady-state. Calculations are presented for several configurations to evaluate the algorithm modifications. The modifications have significantly improved the stability of the AF algorithm and hence the reliability of the CAP-TSD code in general.

Batina, John T.↗

The effect of horizontal resolution on systematic errors of the GLA forecast model

Systematic prediction errors of the Goddard Laboratory for Atmospheres (GLA) forecast system are reduced when the higher-resolution (2 x 2.5 deg) model version is used. Based on a budget analysis of the 200-mb eddy streamfunction, the improvement of stationary eddy forecasting is seen to be caused by the following mechanism: by increasing the horizontal spatial resolution of the forecast model, atmospheric diabatic heating over the three tropical continents is changed in a way that intensifies the planetary-scale divergent circulations associated with the three pairs of divergent-convergent centers over these continents. The intensified divergent circulation results in an enhancement of vorticity sources in the Northern Hemisphere. The additional vorticity is advected eastward by a stationary wave train along 30 deg N, thereby reducing systematic errors in the lower-resolution (4 x 5 deg) GLA model.

Chen, Tsing-Chang↗

Improved inhomogeneous finite elements for fabric reinforced composite mechanics analysis

There is a need to do routine stress/failure analysis of fabric reinforced composite microstructures to provide additional confidence in critical applications and guide materials development. Conventional methods of 3-D stress analysis are time consuming to set up, run and interpret. A need exists for simpler methods of modeling these structures and analyzing the models. The principal difficulty is the discrete element mesh generation problem. Inhomogeneous finite elements are worth investigating for application to these problems because they eliminate the mesh generation problem. However, there are penalties associated with these elements. Their convergence rates can be slow compared to homogeneous elements. Also, there is no accepted method for obtaining detailed stresses in the constituent materials of each element. This paper shows that the convergence rate can be significantly improved by a simple device which substitutes homogeneous elements for the inhomogeneous ones. The device is shown to work well in simple one and two dimensional problems. However, demonstration of the application to more complex two and three dimensional problems remains to be done. Work is also progressing toward more realistic fabric microstructural geometries.

Foye, R. L.↗

The Climate Signal in Regional Moisture Fluxes: A Comparison of Three Global Data Assimilation Products

This study assesses the quality of estimates of climate variability in moisture flux and convergence from three assimilated data sets: two are reanalysis products generated at the Goddard Data Assimilation Office (DAO) and the National Centers for Environmental Prediction/National Centers for Atmospheric Research (NCEPJNCAR), and the third consists of the operational analyses generated at the European Center for Medium Range Forecasts (ECMWF). The regions under study (the United States Great Plains, the Indian monsoon region, and Argentina east of the Andes) are characterized by frequent low level jets (LLJs) and other interannual low level wind variations tied to the large-scale flow. While the emphasis is on the reanalysis products, the comparison with the operational product is provided to help assess the improvements gained from a fixed analysis system. All three analyses capture the main moisture flux anomalies associated with selected extreme climate (drought and flood) events during the period 1985-93. The correspondence is strongest over the Great Plains and weakest over the Indian monsoon region reflecting differences in the observational coverage. For the reanalysis products, the uncertainties in the lower tropospheric winds is by far the dominant source of the discrepancies in the moisture flux anomalies in the middle latitude regions. Only in the Indian Monsoon region, where interannual variability in the low level winds is comparatively small, does the moisture bias play a substantial role. In contrast, the comparisons with the operational product show differences in moisture which are comparable torhe differences in the wind in all three regions. Compared with the fluxes, the anomalous moisture convergences show substantially larger differences among the three products. The best agreement occurs over the Great Plains region where all three products show vertically-integrated moisture convergence during the floods and divergence during the drought with differences in magnitude of about 25%. The reanalysis products, in particular, show good agreement in depicting the different roles of the mean flow and transients during the flood and drought periods. Differences between the three products in the other two regions exceed 100% reflecting differences in the low level jets and the large scale circulation patterns. The operational product tends to have locally larger amplitude convergence fields which average out in area-mean budgets: this appears to be at least in part due to errors in the surface pressure fields and aliasing from the higher resolution of the original ECMWF fields. On average, the reanalysis products show higher coherence with each other than with the operational product in the estimates of interannual variability. This result is less clear in the Indian monsoon region where differences in the input observations appears to be an important factor. The agreement in the anomalous convergence patterns is, however, still rather poor even over relatively data dense regions such as the United States Great Plains. These differences are attributed to deficiencies in the assimilating GCM's representations of the planetary boundary layer and orography, and a global observing system incapable of resolving the highly confined low level winds associated with the climate anomalies.

Min, Wei↗