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 235 records · Page 13

Analysis and design of numerical schemes for gas dynamics 1: Artificial diffusion, upwind biasing, limiters and their effect on accuracy and multigrid convergence

The theory of non-oscillatory scalar schemes is developed in this paper in terms of the local extremum diminishing (LED) principle that maxima should not increase and minima should not decrease. This principle can be used for multi-dimensional problems on both structured and unstructured meshes, while it is equivalent to the total variation diminishing (TVD) principle for one-dimensional problems. A new formulation of symmetric limited positive (SLIP) schemes is presented, which can be generalized to produce schemes with arbitrary high order of accuracy in regions where the solution contains no extrema, and which can also be implemented on multi-dimensional unstructured meshes. Systems of equations lead to waves traveling with distinct speeds and possibly in opposite directions. Alternative treatments using characteristic splitting and scalar diffusive fluxes are examined, together with modification of the scalar diffusion through the addition of pressure differences to the momentum equations to produce full upwinding in supersonic flow. This convective upwind and split pressure (CUSP) scheme exhibits very rapid convergence in multigrid calculations of transonic flow, and provides excellent shock resolution at very high Mach numbers.

Jameson, Antony↗

Parallel Computational Environment for Substructure Optimization

Design optimization of large structural systems can be attempted through a substructure strategy when convergence difficulties are encountered. When this strategy is used, the large structure is divided into several smaller substructures and a subproblem is defined for each substructure. The solution of the large optimization problem can be obtained iteratively through repeated solutions of the modest subproblems. Substructure strategies, in sequential as well as in parallel computational modes on a Cray YMP multiprocessor computer, have been incorporated in the optimization test bed CometBoards. CometBoards is an acronym for Comparative Evaluation Test Bed of Optimization and Analysis Routines for Design of Structures. Three issues, intensive computation, convergence of the iterative process, and analytically superior optimum, were addressed in the implementation of substructure optimization into CometBoards. Coupling between subproblems as well as local and global constraint grouping are essential for convergence of the iterative process. The substructure strategy can produce an analytically superior optimum different from what can be obtained by regular optimization. For the problems solved, substructure optimization in a parallel computational mode made effective use of all assigned processors.

Gendy, Atef S.↗

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster↗

Nonlinear panel flutter - A finite-element approach

A finite-element approach has been developed for computing nonlinear flutter characteristics of rectangular isotropic panels with stream-alined side edges, based on aerodynamic forces from supersonic two-dimensional quasi-steady aerodynamic theory. Stress distributions and panel oscillation frequencies were determined from the analysis. The finite-element formulation, solution procedure, and convergence characteristics are presented. Comparisons are made with linear flutter and large-amplitude vibration results and demonstrate that good accuracy is obtained. Non-linear flutter results are presented for effects of aerodynamic damping, length-width ratio, initial in-plane forces and boundary-support conditions. Comparisons with experimental results are also presented.

Mei, C.↗

Analytical study of STOL Aircraft in ground effect. Part 1: Nonplanar, nonlinear wing/jet lifting surface method

A nonlinear, nonplanar three dimensional jet flap analysis, applicable to the ground effect problem, is presented. Lifting surface methodology is developed for a wing with arbitrary planform operating in an inviscid and incompressible fluid. The classical, infintely thin jet flap model is employed to simulate power induced effects. An iterative solution procedure is applied within the analysis to successively approximate the jet shape until a converged solution is obtained which closely satisfies jet and wing boundary conditions. Solution characteristics of the method are discussed and example results are presented for unpowered, basic powered and complex powered configurations. Comparisons between predictions of the present method and experimental measurements indicate that the improvement of the jet with the ground plane is important in the analyses of powered lift systems operating in ground proximity. Further development of the method is suggested in the areas of improved solution convergence, more realistic modeling of jet impingement and calculation efficiency enhancements.

Shollenberger, C. A.↗

Numerical simulation of steady transonic flow about airfoils

A computer code has been developed that couples a fast transonic full-potential AF2 solver with both an efficient integral boundary-layer method and a viscous wedge approximation of the shock/boundary-layer interaction. The efficiency of the coupled analysis methods and the method of coupling has resulted in a uniquely efficient analysis tool. The airfoil geometry is modified by the displacement thickness before the shock and the displacement thickness plus the viscous wedge thickness after the shock by considering the viscous effects as an equivalent transpiration boundary condition. The flow about conventional and supercritical airfoils under moderately strong shock situations has been calculated. Comparisons with experimental data indicate that this viscous correction method has improved the accuracy of the full-potential analysis. Furthermore, the computer time required to obtain a converged solution has been reduced.

Lee, S. C.↗

Automated mesoscale winds derived from GOES multispectral imagery

An automated technique for extracting mesoscale winds from sequences of GOES visible infrared spin scan radiatiometer (VISSR) image pairs has been developed, tested extensively, and configured for quasi-real time research applications on the Atmospheric Sciences Division's research computing sytem. The entire system of computer codes was successfully vectorized for execution on an array processor resulting in job turnaround in less than 1 hour. An objective quality control system provides much greater than 99 percent accuracy in eliminating questionable wind estimates. Dynamical analysis of cloud wind divergence has revealed temporally consistent convergence centers on the meso-beta scale that are highly correlated with ongoing and future developing convective storms.

Wilson, G. S.↗

The Cycle of Cosmic Dust

The presence of dust in the interstellar medium of galaxies dramatically affects their spectral appearance, and the interpretation of their underlying physical properties. Consequently, the diagnostic of galaxy spectra depends crucially on our understanding of the nature and properties of these dust particles. Models of interstellar dust particles must be able to reproduce a basic set of observational constraints, including the general interstellar extinction and infrared emission observed in the diffuse interstellar medium (ISM). Recent analysis of the solar spectrum have resulted in a convergence between the solar and B star abundances. This development, and the steadiness in the solar abundance determination of the primary refractory elements Mg, Si, and Fe, strongly suggest that any viable dust model must also obey interstellar abundances constraints. I will present fifteen dust models that differ in composition and size distribution and that simultaneously satisfy the local extinction, infrared emission, and abundances constraints have been shown to exist. This multitude of viable dust models provides us with an increased flexibility in understanding dust evolution and the many variations in dust properties in different phases of the ISM and stellar environments.

Dwek, Eliahu↗

ACT Payload Shroud Structural Concept Analysis and Optimization

Aerospace structural applications demand a weight efficient design to perform in a cost effective manner. This is particularly true for launch vehicle structures, where weight is the dominant design driver. The design process typically requires many iterations to ensure that a satisfactory minimum weight has been obtained. Although metallic structures can be weight efficient, composite structures can provide additional weight savings due to their lower density and additional design flexibility. This work presents structural analysis and weight optimization of a composite payload shroud for NASA s Ares V heavy lift vehicle. Two concepts, which were previously determined to be efficient for such a structure are evaluated: a hat stiffened/corrugated panel and a fiber reinforced foam sandwich panel. A composite structural optimization code, HyperSizer, is used to optimize the panel geometry, composite material ply orientations, and sandwich core material. HyperSizer enables an efficient evaluation of thousands of potential designs versus multiple strength and stability-based failure criteria across multiple load cases. HyperSizer sizing process uses a global finite element model to obtain element forces, which are statistically processed to arrive at panel-level design-to loads. These loads are then used to analyze each candidate panel design. A near optimum design is selected as the one with the lowest weight that also provides all positive margins of safety. The stiffness of each newly sized panel or beam component is taken into account in the subsequent finite element analysis. Iteration of analysis/optimization is performed to ensure a converged design. Sizing results for the hat stiffened panel concept and the fiber reinforced foam sandwich concept are presented.

Zalewski, Bart B.↗

A Critical Study of Agglomerated Multigrid Methods for Diffusion

Agglomerated multigrid techniques used in unstructured-grid methods are studied critically for a model problem representative of laminar diffusion in the incompressible limit. The studied target-grid discretizations and discretizations used on agglomerated grids are typical of current node-centered formulations. Agglomerated multigrid convergence rates are presented using a range of two- and three-dimensional randomly perturbed unstructured grids for simple geometries with isotropic and stretched grids. Two agglomeration techniques are used within an overall topology-preserving agglomeration framework. The results show that multigrid with an inconsistent coarse-grid scheme using only the edge terms (also referred to in the literature as a thin-layer formulation) provides considerable speedup over single-grid methods but its convergence deteriorates on finer grids. Multigrid with a Galerkin coarse-grid discretization using piecewise-constant prolongation and a heuristic correction factor is slower and also grid-dependent. In contrast, grid-independent convergence rates are demonstrated for multigrid with consistent coarse-grid discretizations. Convergence rates of multigrid cycles are verified with quantitative analysis methods in which parts of the two-grid cycle are replaced by their idealized counterparts.

Nishikawa, Hiroaki↗

Shape Servoing of Deformable Objects using Adaptive Deformation Model Estimation

In this paper, we propose an adaptive shape servoing method to deform a soft object into a desired 3-D shape. The high dimensional representation and the unknown deformation properties of the soft object pose a challenge to actively manipulate its shape. To address this issue, we develop a method to compute the deformation Jacobian matrix in real-time. The Jacobian is estimated using a set of basis functions and its corresponding parameters to capture the dynamics of the system and relate the applied input motion to changes in the soft object's shape. An integral concurrent learning (ICL) based adaptive update law is derived using Lyapunov analysis to estimate the deformation parameters and prove its convergence. A physics-based simulation is used to validate the proposed method and controller by performing manipulation tasks with different desired configurations. The performance is compared with a standard gradient update law to demonstrate the accuracy and robustness of our approach.

Vrithik Raj Guthikonda↗

Investigation of the Potential for Lidar & Polarimeter Combined Aerosol Retrievals Using Various Instrument Combinations

As a scientific community, we are motivated to achieve better observational capabilities for aerosols, to reduce uncertainties about aerosol’s role in climate. Encouraged by the National Academies Decadal Survey, many are looking to the synergy of lidar and polarimeter as a means of achieving improved measurements. Polarimetric observations of total-column radiances have excellent sensitivity to aerosol size distribution and absorption properties. Active lidar measurements have fewer wavelengths and view angles and therefore more limited sensitivity to aerosol microphysical properties such as complex refractive index. Yet the advantage of vertically resolved profiles is game-changing, especially given that models contain significant uncertainty about the vertical distribution of aerosol and that polarimeter and other passive retrievals must make assumptions about vertical distribution which can impact performance. To optimize the synergy, a joint polarimeter + lidar mission requires a joint retrieval. However, the exact instruments that will be launched on a future satellite are unknown, and the capabilities of polarimeters and lidars cover a large range. In this study, we perform information content analysis for a range of nominal instrument combinations. The methodology is ideal for assessing the impact of measurements and a priori information (assumptions and constraints) with respect to the desired retrieval products, and for comparing measurement systems. It has advantages that make it useful in advance of or in parallel with retrieval development, such as producing quantitatively useful comparisons even for retrieval scenarios that are significantly underdetermined and therefore impossible to retrieve without significant assumptions, and it is not impacted by retrieval-specific sources of error such as incomplete or incorrect convergence which can complicate the interpretation of a retrieval performance analysis. Here, we will derive the expected uncertainty in retrieved aerosol concentration, effective radius, fine mode fraction and single scattering albedo, for a combination of a polarimeter with various grades of lidar, from an elastic backscatter lidar, similar to the pathfinder CALIOP on the CALIPSO satellite, to a multi-wavelength high spectral resolution lidar.

Sharon P Burton↗

Pattern Search Algorithms for Bound Constrained Minimization

We present a convergence theory for pattern search methods for solving bound constrained nonlinear programs. The analysis relies on the abstract structure of pattern search methods and an understanding of how the pattern interacts with the bound constraints. This analysis makes it possible to develop pattern search methods for bound constrained problems while only slightly restricting the flexibility present in pattern search methods for unconstrained problems. We prove global convergence despite the fact that pattern search methods do not have explicit information concerning the gradient and its projection onto the feasible region and consequently are unable to enforce explicitly a notion of sufficient feasible decrease.

Lewis, Robert Michael↗

Gas Seal Pad With Herringbone-Grooved Rotor-Stiffness and Load Capacity

The principle of herringbone-grooved journal bearings has been applied to the case of a seal disc running under a finger seal pad. The inward pumping action of herringbone grooves on the disc generates load capacity and stiffness to maintain a fluid film and prevent contact of the pad and disc. This mechanism does not depend on a converging film under the pad, such as analyzed in previous works. Analysis shows that significant stiffness and load capacity can be supplied by herringbone grooves. In order for the grooves to be effective, the seal pressure drop must be taken outside of the grooved portion of the rotor, but this may be acceptable in order to gain freedom from maintaining a precise film convergence.

Flemming, David P.↗

Analysis of a swirl stabilized arcjet thruster

Viscous swirling flow through transonic converging-diverging nozzles (with and without centerbodies) is examined numerically. An implicit technique is used to solve the axisymmetric Thin Layer Naiver-Stokes equations. The effect of swirl and viscosity on low Reynolds number flows is shown by examining integral parameters and flowfield details.

Keith, Theo G., Jr.↗

Error Estimates of the Ares I Computed Turbulent Ascent Longitudinal Aerodynamic Analysis

Numerical predictions of the longitudinal aerodynamic characteristics for the Ares I class of vehicles, along with the associated error estimate derived from an iterative convergence grid refinement, are presented. Computational results are based on an unstructured grid, Reynolds-averaged Navier-Stokes analysis. The validity of the approach to compute the associated error estimates, derived from a base grid to an extrapolated infinite-size grid, was first demonstrated on a sub-scaled wind tunnel model at representative ascent flow conditions for which the experimental data existed. Such analysis at the transonic flow conditions revealed a maximum deviation of about 23% between the computed longitudinal aerodynamic coefficients with the base grid and the measured data across the entire roll angles. This maximum deviation from the wind tunnel data was associated with the computed normal force coefficient at the transonic flow condition and was reduced to approximately 16% based on the infinite-size grid. However, all the computed aerodynamic coefficients with the base grid at the supersonic flow conditions showed a maximum deviation of only about 8% with that level being improved to approximately 5% for the infinite-size grid. The results and the error estimates based on the established procedure are also presented for the flight flow conditions.

Abdol-Hamid, Khaled S.↗

Wind-Stress Simulations and Equatorial Dynamics in an AGCM: Basic Results from a 1979-1999 Forced SST Experiment - Part 1

This is the first of a two part study examining the connection of the equatorial momentum budget in an AGCM (Atmospheric General Circulation Model), with simulated equatorial surface wind stresses over the Pacific. The AGCM used in this study forms part of a newly developed coupled forecasting system used at NASA's Seasonal- to-Interannual Prediction Project. Here we describe the model and present results from a 20-year (1979-1999) AMIP-type experiment forced with observed SSTs (Sea Surface Temperatures). Model results are compared them with available observational data sets. The climatological pattern of extra-tropical planetary waves as well as their ENSO-related variability is found to agree quite well with re-analysis estimates. The model's surface wind stress is examined in detail, and reveals a reasonable overall simulation of seasonal interannual variability, as well as seasonal mean distributions. However, an excessive annual oscillation in wind stress over the equatorial central Pacific is found. We examine the model's divergent circulation over the tropical Pacific and compare it with estimates based on re-analysis data. These comparisons are generally good, but reveal excessive upper-level convergence in the central Pacific. In Part II of this study a direct examination of individual terms in the AGCM's momentum budget is presented. We relate the results of this analysis to the model's simulation of surface wind stress.

Bacmeister, Julio T.↗

The Influence of Volcanic Eruptions on the Climate of Tropical South America During the Last Millennium in an Isotope-Enabled General Circulation Model

Currently, little is known on how volcanic eruptions impact large-scale climate phenomena such as South American paleo-intertropical Convergence Zone (ITCZ) position and summer monsoon behavior. In this paper, an analysis of observations and model simulations is employed to assess the influence of large volcanic eruptions on the climate of tropical South America. This problem is first considered for historically recent volcanic episodes for which more observations are available but where fewer events exist and the confounding effects of El Niño-Southern Oscillation (ENSO) lead to inconclusive interpretation of the impact of volcanic eruptions at the continental scale. Therefore, we also examine a greater number of reconstructed volcanic events for the period 850 CE to present that are incorporated into the NASA GISS ModelE2-R simulation of the last millennium. An advantage of this model is its ability to explicitly track water isotopologues throughout the hydrologic cycle and simulating the isotopic imprint following a large eruption. This effectively removes a degree of uncertainty associated with error-prone conversion of isotopic signals into climate variables, and allows for a direct comparison between GISS simulations and paleoclimate proxy records. Our analysis reveals that both precipitation and oxygen isotope variability respond with a distinct seasonal and spatial structure across tropical South America following an eruption. During austral winter, the heavy oxygen isotope in precipitation is enriched, likely due to reduced moisture convergence in the ITCZ domain and reduced rainfall over northern South America. During austral summer, however, more negative values of the precipitation isotopic composition are simulated over Amazonia, despite reductions in rainfall, suggesting that the isotopic response is not a simple function of the "amount effect". During the South American monsoon season, the amplitude of the temperature response to volcanic forcing is larger than the rather weak and spatially less coherent precipitation signal, complicating the isotopic response to changes in the hydrologic cycle.

Colose, Christopher M.↗