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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1989

New free-wake analysis of rotorcraft hover performance using influence coefficients

Free-wake analyses of helicopter rotor wakes in hover using time stepping have been shown to encounter instabilities which preclude convergence to valid free-vortex solutions for rotor-wake geometries. Previous work has demonstrated that these convergence difficulties can be overcome by implementing a new free-wake analysis method based on the use of influence coefficients. The present paper reviews this approach and documents its incorporation into a hover performance analysis called Evaluation of Hover Performance using Influence Coefficients (EHPIC). The technical principles underlying the EHPIC code are described with emphasis on steps taken to develop the single-filament wake models used in previous work into a multifilament wake valid for realistic hover performance predictions. The coupling of the wake model to a lifting surface loads analysis is described, and sample problems are solved that illustrate the robustness of the method. Performance calculations are also undertaken for hover to illustrate the utility of EHPIC in the analysis of rotorcraft performance.

Quackenbush, Todd R.↗

NASA NTRS1983

Improving the convergence rate of parabolic ADI methods

The rate of convergence to steady state of parabolic Alternating Direction Implicit (ADI) solvers is analyzed in terms of the L(2)-norms of the residuals. The analysis allows one to predict the number of iterations necessary for convergence as function of the Courant number, Lambda. A simple modification of existing ADI codes is devised. It improves the convergence rate substantially and is insensitive to the Courant number in a large range of Lambda.

Abarbanel, S. S.↗

NASA NTRS1985

The Structure and Dynamics of an Observed Moist Front

The structure and dynamics of the moist cold front of 25 and 26 April 1979, the third observing day of the Severe Environmental Storms and Mesoscale experiment (SESAME), are investigated through the use of a three-dimensional mesoscale numerical model. This work is one of the first studies in which model results are compared, in a one-to-one manner, with a detailed observational analysis, namely that of Ogura and Portis (1982) as taken from the SESAME observations. In addition, frontogenetical effects, both adiabatic and diabatic, are studied on a vertical cross-section through the front; similarities and differences with the adiabatic analysis of Ogura and Portis are discussed. Many similarities exist between the modeled and analyzed fields, although the analysis tends to have weaker horizontal gradients due to the coarseness of the observational network. Vorticity and convergence near the surface were found to have the same magnitude in both the model solution and the analysis, in contrast to idealized frontogenesis models which predict vorticity to be much larger than convergence.

Orlanski, I.↗

NASA NTRS1987

Convergence properties of modal costs for certain distributed parameter systems

A complete modal cost analysis is presented for the vibration of various kinds of simple continua with different boundary conditions. Explicit formulas for the norm of the response (called the 'cost') are derived for these distributed parameter systems. The convergence theorems developed are useful in the model reduction of the equivalent continuum models of large space structures as well as in the selection of finite-element code for control design.

Hu, A.↗

NASA NTRS1993

Helicopter trim analysis by shooting and finite element methods with optimally damped Newton iterations

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used, and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Achar, N. S.↗

NASA NTRS2005

Numerical Evaluation of P-Multigrid Method for the Solution of Discontinuous Galerkin Discretizations of Diffusive Equations

This paper describes numerical experiments with P-multigrid to corroborate analysis, validate the present implementation, and to examine issues that arise in the implementations of the various combinations of relaxation schemes, discretizations and P-multigrid methods. The two approaches to implement P-multigrid presented here are equivalent for most high-order discretization methods such as spectral element, SUPG, and discontinuous Galerkin applied to advection; however it is discovered that the approach that mimics the common geometric multigrid implementation is less robust, and frequently unstable when applied to discontinuous Galerkin discretizations of di usion. Gauss-Seidel relaxation converges 40% faster than block Jacobi, as predicted by analysis; however, the implementation of Gauss-Seidel is considerably more expensive that one would expect because gradients in most neighboring elements must be updated. A compromise quasi Gauss-Seidel relaxation method that evaluates the gradient in each element twice per iteration converges at rates similar to those predicted for true Gauss-Seidel.

Atkins, H. L.↗

NASA NTRS1993

Reliability enhancement of Navier-Stokes codes through convergence enhancement

Reduction of total computing time required by an iterative algorithm for solving Navier-Stokes equations is an important aspect of making the existing and future analysis codes more cost effective. Several attempts have been made to accelerate the convergence of an explicit Runge-Kutta time-stepping algorithm. These acceleration methods are based on local time stepping, implicit residual smoothing, enthalpy damping, and multigrid techniques. Also, an extrapolation procedure based on the power method and the Minimal Residual Method (MRM) were applied to the Jameson's multigrid algorithm. The MRM uses same values of optimal weights for the corrections to every equation in a system and has not been shown to accelerate the scheme without multigriding. Our Distributed Minimal Residual (DMR) method based on our General Nonlinear Minimal Residual (GNLMR) method allows each component of the solution vector in a system of equations to have its own convergence speed. The DMR method was found capable of reducing the computation time by 10-75 percent depending on the test case and grid used. Recently, we have developed and tested a new method termed Sensitivity Based DMR or SBMR method that is easier to implement in different codes and is even more robust and computationally efficient than our DMR method.

Choi, K.-Y.↗

NASA NTRS1986

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations, the steady Euler equations, using Newton's linearization procedure is presented. A theorem indicating quadratic convergence for the case of differential equations is demonstrated. A condition for the domain of quadratic convergence Omega(2) is obtained which indicates that whether an approximation lies in Omega(2) depends on the rate of change and the smoothness of the flow vectors, and hence is problem-dependent. The choice of spatial differencing, of particular importance for the present method, is discussed. The treatment of boundary conditions is addressed, and the system of equations resulting from the foregoing analysis is summarized and solution strategies are discussed. The convergence of calculated solutions is demonstrated by comparing them with exact solutions to one and two-dimensional problems.

Liou, M.-S.↗

NASA NTRS2017

A Dose of Reality: Radiation Analysis for Realistic Human Spacecraft

INTRODUCTION As with most computational analyses, a tradeoff exists between problem complexity, resource availability and response accuracy when modeling radiation transport from the source to a detector. The largest amount of analyst time for setting up an analysis is often spent ensuring that any simplifications made have minimal impact on the results. The vehicle shield geometry of interest is typically simplified from the original CAD design in order to reduce computation time, but this simplification requires the analyst to "re-draw" the geometry with a limited set of volumes in order to accommodate a specific radiation transport software package. The resulting low-fidelity geometry model cannot be shared with or compared to other radiation transport software packages, and the process can be error prone with increased model complexity. The work presented here demonstrates the use of the DAGMC (Direct Accelerated Geometry for Monte Carlo) Toolkit from the University of Wisconsin, to model the impacts of several space radiation sources on a CAD drawing of the US Lab module. METHODS The DAGMC toolkit workflow begins with the export of an existing CAD geometry from the native CAD to the ACIS format. The ACIS format file is then cleaned using SpaceClaim to remove small holes and component overlaps. Metadata is then assigned to the cleaned geometry file using CUBIT/Trelis from csimsoft (Registered Trademark). The DAGMC plugin script removes duplicate shared surfaces, facets the geometry to a specified tolerance, and ensures that the faceted geometry is water tight. This step also writes the material and scoring information to a standard input file format that the analyst can alter as desired prior to running the radiation transport program. The scoring results can be transformed, via python script, into a 3D format that is viewable in a standard graphics program. RESULTS The CAD model of the US Lab module of the International Space Station, inclusive of all the racks and components, was simplified to remove holes and volume overlaps. Problematic features within the drawing were also removed or repaired to prevent runtime issues. The cleaned drawing was then run through the DAGMC workflow to prepare for analysis. Pilot tests modeling transport of 1GeV proton and 800MeV/A oxygen sources show that reasonable results are converged upon in an acceptable amount of overall computation time from drawing preparation to data analysis. The FLUKA radiation transport code will next be used to model both a GCR and a trapped radiation source. These results will then be compared with measurements that have been made by the radiation instrumentation deployed inside the US Lab module. DISCUSSION Early analyses have indicated that the DAGMC workflow is a promising toolkit for running vehicle geometries of interest to NASA through multiple radiation transport codes. In addition, recent work has shown that a realistic human phantom, provided via a subcontract with the University of Florida, can be placed inside any vehicle geometry for a combinatorial analysis. This added functionality gives the user the ability to score various parameters at the organ level, and the results can then be used as input for cancer risk models.

Barzilla, J. E.↗

NASA NTRS1974

Analysis of complex elastic structures by a Rayleigh-Ritz component modes method using Lagrange multipliers

The free vibrations of elastic structures of arbitrary complexity were analyzed in terms of their component modes. The method was based upon the use of the normal unconstrained modes of the components in a Rayleigh-Ritz analysis. The continuity conditions were enforced by means of Lagrange Multipliers. Examples of the structures considered are: (1) beams with nonuniform properties; (2) airplane structures with high or low aspect ratio lifting surface components; (3) the oblique wing airplane; and (4) plate structures. The method was also applied to the analysis of modal damping of linear elastic structures. Convergence of the method versus the number of modes per component and/or the number of components is discussed and compared to more conventional approaches, ad-hoc methods, and experimental results.

Klein, L. R.↗

NASA NTRS1984

Aerospace engineering design by systematic decomposition and multilevel optimization

A method for systematic analysis and optimization of large engineering systems, by decomposition of a large task into a set of smaller subtasks that is solved concurrently is described. The subtasks may be arranged in hierarchical levels. Analyses are carried out in each subtask using inputs received from other subtasks, and are followed by optimizations carried out from the bottom up. Each optimization at the lower levels is augmented by analysis of its sensitivity to the inputs received from other subtasks to account for the couplings among the subtasks in a formal manner. The analysis and optimization operations alternate iteratively until they converge to a system design whose performance is maximized with all constraints satisfied. The method, which is still under development, is tentatively validated by test cases in structural applications and an aircraft configuration optimization.

Sobieszczanski-Sobieski, J.↗

NASA NTRS1984

Automated mesoscale winds derived from GOES multispectral imagery

An automated technique for extracting mesoscale winds from sequences of GOES VISSR image pairs was developed, tested and configured for quasi-real time/research applications on a computing system which gives mesoscale wind estimates at the highest spatial/temporal resolution possible from the VISSR imagery down to a wind vector separation of 10 km. Preprocessing of imagery using IR resampling, VIS edge preserving filtering, and reduced VIS resolution averaging improved height assignments and vector extraction for 10, 15, and 30 min imagery. An objective quality control system provides much greater than 99% accuracy in eliminating questionable wind estimates. Automated winds generally have better spatial coverage and density, and have random error estimates half as large as the manual winds. Dynamical analysis of cloud wind divergence revealed temporally consistent convergence centers on the meso beta scale that are highly correlated with on going and future developing convective storms. The entire system of computer codes was successfully vectorized for execution on an array processor resulting in job turnaround in less than one hour.

Wilson, G. S.↗

NASA NTRS1984

Aerospace engineering design by systematic decomposition and multilevel optimization

This paper describes a method for systematic analysis and optimization of large engineering systems, e.g., aircraft, by decomposition of a large task into a set of smaller, self-contained subtasks that can be solved concurrently. The subtasks may be arranged in many hierarchical levels with the assembled system at the top level. Analyses are carried out in each subtask using inputs received from other subtasks, and are followed by optimizations carried out from the bottom up. Each optimization at the lower levels is augmented by analysis of its sensitivity to the inputs received from other subtasks to account for the couplings among the subtasks in a formal manner. The analysis and optimization operations alternate iteratively until they converge to a system design whose performance is maximized with all constraints satisfied. The method, which is still under development, is tentatively validated by test cases in structural applications and an aircraft configuration optimization. It is pointed out that the method is intended to be compatible with the typical engineering organization and the modern technology of distributed computing.

Sobieszczanski-Sobieski, J.↗

NASA NTRS1992

Magnetic flipping - Reconnection in three dimensions without null points

In three dimensions, magnetic reconnection may take place in a sheared magnetic field at any singular field line, where the nearby field has X-type topology in planes perpendicular to the field line and where an electric field is present parallel to the field line. In the ideal region around the singular line there will, in general, be singularities in the plasma flow and electric field, both at the singular line and at 'magnetic flipping layers', which are remnants of local magnetic separatrices. In the absence of a three-dimensional magnetic neutral point or null point, reconnection of field lines can still occur by a process of magnetic flipping, in which the plasma crosses the flipping layers but the field lines rapidly flip along them by magnetic diffusion. Depending on the boundary conditions, there may be two or four flipping layers which converge on the singular line. A boundary layer analysis of a flipping layer is given, in which the magnetic field parallel to the layer decreases as one crosses it while the plasma pressure (or magnetic pressure associated with the field along the singular line) increases. The width of the flipping layer decreases with distance from the singular line.

Priest, E. R.↗

NASA NTRS1995

Comparison between global latent heat flux computed from multisensor (SSM/I and AVHRR) and from in situ data

Comprehensive Ocean-Atmosphere Data Set (COADS) and satellite-derived parameters are input to a similarity theory-based model and treated in completely equivalent ways to compute global latent heat flux (LHF). In order to compute LHF exclusively from satellite measurements, an empirical relationship (Q-W relationship) is used to compute the air mixing ratio from Special Sensor Microwave/Imager (SSM/I) precipitable water W and a new one is derived to compute the air temperature also from retrieved W(T-W relationship). First analyses indicate that in situ and satellite LHF computations compare within 40%, but systematic errors increase the differences up to 100% in some regions. By investigating more closely the origin of the discrepancies, the spatial sampling of ship reports has been found to be an important source of error in the observed differences. When the number of in situ data records increases (more than 20 per month), the agreement is about 50 W/sq m rms (40 W/sq m rms for multiyear averages). Limitations of both empirical relationships and W retrieval errors strongly affect the LHF computation. Systematic LHF overestimation occurs in strong subsidence regions and LHF underestimation occurs within surface convergence zones and over oceanic upwelling areas. The analysis of time series of the different parameters in these regions confirms that systematic LHF discrepancies are negatively correlated with the differences between COADS and satellite-derived values of the air mixing ratio and air temperature. To reduce the systematic differences in satellite-derived LHF, a preliminary ship-satellite blending procedure has been developed for the air mixing ratio and air temperature.

Jourdan, Didier↗

NASA NTRS1992

Small Business Innovations (Helicopters)

The amount of engine power required for a helicopter to hover is an important, but difficult, consideration in helicopter design. The EHPIC program model produces converged, freely distorted wake geometries that generate accurate analysis of wake-induced downwash, allowing good predictions of rotor thrust and power requirements. Continuum Dynamics, Inc., the Small Business Innovation Research (SBIR) company that developed EHPIC, also produces RotorCRAFT, a program for analysis of aerodynamic loading of helicopter blades in forward flight. Both helicopter codes have been licensed to commercial manufacturers.

Source record↗

NASA NTRS2008

Bounded Linear Stability Margin Analysis of Nonlinear Hybrid Adaptive Control

This paper presents a bounded linear stability analysis for a hybrid adaptive control that blends both direct and indirect adaptive control. Stability and convergence of nonlinear adaptive control are analyzed using an approximate linear equivalent system. A stability margin analysis shows that a large adaptive gain can lead to a reduced phase margin. This method can enable metrics-driven adaptive control whereby the adaptive gain is adjusted to meet stability margin requirements.

Nguyen, Nhan T.↗

NASA NTRS2011

Convergence of Defect-Correction and Multigrid Iterations for Inviscid Flows

Convergence of multigrid and defect-correction iterations is comprehensively studied within different incompressible and compressible inviscid regimes on high-density grids. Good smoothing properties of the defect-correction relaxation have been shown using both a modified Fourier analysis and a more general idealized-coarse-grid analysis. Single-grid defect correction alone has some slowly converging iterations on grids of medium density. The convergence is especially slow for near-sonic flows and for very low compressible Mach numbers. Additionally, the fast asymptotic convergence seen on medium density grids deteriorates on high-density grids. Certain downstream-boundary modes are very slowly damped on high-density grids. Multigrid scheme accelerates convergence of the slow defect-correction iterations to the extent determined by the coarse-grid correction. The two-level asymptotic convergence rates are stable and significantly below one in most of the regions but slow convergence is noted for near-sonic and very low-Mach compressible flows. Multigrid solver has been applied to the NACA 0012 airfoil and to different flow regimes, such as near-tangency and stagnation. Certain convergence difficulties have been encountered within stagnation regions. Nonetheless, for the airfoil flow, with a sharp trailing-edge, residuals were fast converging for a subcritical flow on a sequence of grids. For supercritical flow, residuals converged slower on some intermediate grids than on the finest grid or the two coarsest grids.

Diskin, Boris↗