Search NASA⌕ Search

SEARCH · Search NASA

Results for “ITERATIVE SOLUTION”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Fast methods to numerically integrate the Reynolds equation for gas fluid films

The alternating direction implicit (ADI) method is adopted, modified, and applied to the Reynolds equation for thin, gas fluid films. An efficient code is developed to predict both the steady-state and dynamic performance of an aerodynamic journal bearing. An alternative approach is shown for hybrid journal gas bearings by using Liebmann's iterative solution (LIS) for elliptic partial differential equations. The results are compared with known design criteria from experimental data. The developed methods show good accuracy and very short computer running time in comparison with methods based on an inverting of a matrix. The computer codes need a small amount of memory and can be run on either personal computers or on mainframe systems.

Dimofte, Florin↗
NASA NTRS1993

Analyzing Transonic Flow Over A Wing And Fuselage

TAWFIVE computer program calculates transonic flow over transport-type wing and fuselage. Combines multigrid acceleration technique in iterative solution of potential equation with use of integral-form boundary-layer equations to provide software tool computationally efficient, yet accurate enough for use in design. Simplifies solution process by breaking problem into loosely coupled set of modified equations. Inviscid-flow method, using standard inviscid-flow equations, valid in "outer" region away from wing, whereas boundary-layer equations valid in thin region near solid surface of wing. Written in FORTRAN 77.

Melson, N. Duane↗
NASA NTRS1992

Structural optimization of thin shells using finite element method

The objective of the present work was the structural optimization of thin shell structures that are subjected to stress and displacement constraints. In order to accomplish this, the structural optimization computer program DESAP1 was modified and improved. In the static analysis part of the DESAP1 computer program the torsional spring elements, which are used to analyze thin, shallow shell structures, were eliminated by modifying the membrane stiffness matrix of the triangular elements in the local coordinate system and adding a fictitious rotational stiffness matrix. This simplified the DESAP1 program input, improved the accuracy of the analysis, and saved computation time. In the optimization part of the DESAP1 program the stress ratio formula, which redesigns the thickness of each finite element of the structure, was solved by an analytical method. This scheme replaced the iterative solution that was previously used in the DESAP1 program, thus increasing the accuracy and speed of the design. The modified program was used to design a thin, cylindrical shell structure with optimum weight, and the results are reported in this paper.

Gotsis, Pascal K.↗
NASA NTRS1992

A microwave scattering model for layered vegetation

A microwave scattering model was developed for layered vegetation based on an iterative solution of the radiative transfer equation up to the second order to account for multiple scattering within the canopy and between the ground and the canopy. The model is designed to operate over a wide frequency range for both deciduous and coniferous forest and to account for the branch size distribution, leaf orientation distribution, and branch orientation distribution for each size. The canopy is modeled as a two-layered medium above a rough interface. The upper layer is the crown containing leaves, stems, and branches. The lower layer is the trunk region modeled as randomly positioned cylinders with a preferred orientation distribution above an irregular soil surface. Comparisons of this model with measurements from deciduous and coniferous forests show good agreements at several frequencies for both like and cross polarizations. Major features of the model needed to realize the agreement include allowance for: (1) branch size distribution, (2) second-order effects, and (3) tree component models valid over a wide range of frequencies.

Karam, Mostafa A.↗
NASA NTRS1993

High Temperature Composite Analyzer (HITCAN) theoretical manual, version 1.0

This manual outlines some of the theoretical aspects embedded in the computer code, HITCAN (High Temperature Composite ANalyzer). HITCAN is a general purpose computer program for predicting nonlinear global structural and local stress-strain response of arbitrarily oriented, multilayered high temperature metal matrix composite structures. This code combines composite mechanics and laminate theory with an internal data base of the constituents (matrix, fiber and interface) material properties. The thermal and mechanical properties of the constituents are considered to be nonlinearly dependent on several parameters including temperature, stress and stress rate. The computational procedure for the analysis of the composite structure uses the finite element method. HITCAN consists of three modules: COBSTRAN, METCAN and MHOST. COBSTRAN generates the geometry (pre-processor) and defines the layup of the different plies. METCAN computes the material behavior of the composite and of the constituents. Finally, MHOST is a finite element program based on the mixed iterative solution technique. MHOST has a library for 2D and 3D isoparametric elements. HITCAN is written in FORTRAN 77 computer language and has been configured and executed on the NASA Lewis Research Center CRAY XMP and YMP computers.

Lackney, J. J.↗
NASA NTRS1992

A finite element study of incompressible flows past oscillating cylinders and aerofoils

We present our numerical results for certain unsteady flows past oscillating cylinders and aerofoils. The computations are based on the stabilized space-time finite element formulation. The implicit equation systems resulting from the space-time finite element discretizations are solved using iterative solution techniques. One of the problems studied is flow past a cylinder which is forced to oscillate in the horizontal direction. In this case we observe a change from an unsymmetric mode of vortex shedding to a symmetric one. An extensive study was carried out for the case in which a cylinder is mounted on lightly damped springs and allowed to oscillate in the vertical direction. In this case the motion of the cylinder needs to be determined as part of the solution, and under certain conditions this motion changes the vortex-shedding pattern of the flow field significantly. This nonlinear fluid-structure interaction exhibits certain interesting behavior such as 'lock-in' and 'hysteresis', which are in good agreement with the laboratory experiments carried out by other researchers in the past. Preliminary results for flow past a pitching aerofoil are also presented.

Mittal, S.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗
NASA NTRS1995

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 2: Algorithmic developments and implementation

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Li, Wei↗
NASA NTRS1994

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 1: Formulation

The present work is part of a research program for the numerical simulation of electromagnetic (EM) fields within conventional Ginzburg-Landau (GL) superconductors. The final goal of this research is to formulate, develop and validate finite element (FE) models that can accurately capture electromagnetic thermal and material phase changes in a superconductor. The formulations presented here are for a time-independent Ginzburg-Landau superconductor and are derived from a potential-based variational principle. We develop an appropriate variational formulation of time-independent supercontivity for the general three-dimensional case and specialize it to the one-dimensional case. Also developed are expressions for the material-dependent parameters alpha and beta of GL theory and their dependence upon the temperature T. The one-dimensional formulation is then discretized for finite element purposes and the first variation of these equations is obtained. The resultant Euler equations contain nonlinear terms in the primary variables. To solve these equations, an incremental-iterative solution method is used. Expressions for the internal force vector, external force vector, loading vector and tangent stiffness matrix are therefore developed for use with the solution procedure.

Schuler, James J.↗
NASA NTRS1996

Correction of WindScat Scatterometric Measurements by Combining with AMSR Radiometric Data

The Seawinds scatterometer on the advanced Earth observing satellite-2 (ADEOS-2) will determine surface wind vectors by measuring the radar cross section. Multiple measurements will be made at different points in a wind-vector cell. When dense clouds and rain are present, the signal will be attenuated, thereby giving erroneous results for the wind. This report describes algorithms to use with the advanced mechanically scanned radiometer (AMSR) scanning radiometer on ADEOS-2 to correct for the attenuation. One can determine attenuation from a radiometer measurement based on the excess brightness temperature measured. This is the difference between the total measured brightness temperature and the contribution from surface emission. A major problem that the algorithm must address is determining the surface contribution. Two basic approaches were developed for this, one using the scattering coefficient measured along with the brightness temperature, and the other using the brightness temperature alone. For both methods, best results will occur if the wind from the preceding wind-vector cell can be used as an input to the algorithm. In the method based on the scattering coefficient, we need the wind direction from the preceding cell. In the method using brightness temperature alone, we need the wind speed from the preceding cell. If neither is available, the algorithm can work, but the corrections will be less accurate. Both correction methods require iterative solutions. Simulations show that the algorithms make significant improvements in the measured scattering coefficient and thus is the retrieved wind vector. For stratiform rains, the errors without correction can be quite large, so the correction makes a major improvement. For systems of separated convective cells, the initial error is smaller and the correction, although about the same percentage, has a smaller effect.

Song, S.↗
NASA NTRS1996

Preconditioned Mixed Spectral Element Methods for Elasticity and Stokes Problems

Preconditioned iterative methods for the indefinite systems obtained by discretizing the linear elasticity and Stokes problems with mixed spectral elements in three dimensions are introduced and analyzed. The resulting stiffness matrices have the structure of saddle point problems with a penalty term, which is associated with the Poisson ratio for elasticity problems or with stabilization techniques for Stokes problems. The main results of this paper show that the convergence rate of the resulting algorithms is independent of the penalty parameter, the number of spectral elements Nu and mildly dependent on the spectral degree eta via the inf-sup constant. The preconditioners proposed for the whole indefinite system are block-diagonal and block-triangular. Numerical experiments presented in the final section show that these algorithms are a practical and efficient strategy for the iterative solution of the indefinite problems arising from mixed spectral element discretizations of elliptic systems.

Pavarino, Luca F.↗
NASA NTRS1996

Application of an Upwind High Resolution Finite-Differencing Scheme and Multigrid Method in Steady-State Incompressible Flow Simulations

The analysis and design of a submarine propulsor requires the ability to predict the characteristics of both laminar and turbulent flows to a higher degree of accuracy. This report presents results of certain benchmark computations based on an upwind, high-resolution, finite-differencing Navier-Stokes solver. The purpose of the computations is to evaluate the ability, the accuracy and the performance of the solver in the simulation of detailed features of viscous flows. Features of interest include flow separation and reattachment, surface pressure and skin friction distributions. Those features are particularly relevant to the propulsor analysis. Test cases with a wide range of Reynolds numbers are selected; therefore, the effects of the convective and the diffusive terms of the solver can be evaluated separately. Test cases include flows over bluff bodies, such as circular cylinders and spheres, at various low Reynolds numbers, flows over a flat plate with and without turbulence effects, and turbulent flows over axisymmetric bodies with and without propulsor effects. Finally, to enhance the iterative solution procedure, a full approximation scheme V-cycle multigrid method is implemented. Preliminary results indicate that the method significantly reduces the computational effort.

NAVIER STOKES EQUATIONS↗
NASA NTRS2000

Simulation of Aluminum Micro-mirrors for Space Applications at Cryogenic Temperatures

Closed form and finite element models are developed to predict the device response of aluminum electrostatic torsion micro-mirrors fabricated on silicon substrate for space applications at operating temperatures of 30K. Initially, closed form expressions for electrostatic pressure arid mechanical restoring torque are used to predict the pull-in and release voltages at room temperature. Subsequently, a detailed mechanical finite element model is developed to predict stresses and vertical beam deflection induced by the electrostatic and thermal loads. An incremental and iterative solution method is used in conjunction with the nonlinear finite element model and closed form electrostatic equations to solve. the coupled electro-thermo-mechanical problem. The simulation results are compared with experimental measurements at room temperature of fabricated micro-mirror devices.

Kuhn, J. L.↗
NASA NTRS2000

Aerothermodynamic Analysis of Stardust Sample Return Capsule with Coupled Radiation and Ablation

An aerothermodynamic analysis of the forebody aeroshell of the Stardust Sample Return Capsule is carried out by using the axisymmetric viscous shock-layer equations with and without fully coupled radiation and ablation. Formulation of the viscous shock-layer equations with shoulder radius as the length scale and implementation of the Vigneron pressure condition allow resolution of the flowfield over the shoulder. With a predominantly supersonic outflow over the shoulder, a globally iterated solution or viscous shock-layer equations can be obtained. The stagnation-point results are obtained along a specified trajectory, whereas detailed calculations along the body are provided at the peak-heating point. The equilibrium calculations with ablation injection are the focus of the present study because of the lack of a general chemical nonequilibrium analysis that accounts for both surface and flowfield effect. The equilibrium calculations also provide a simple way to conserve surface (and flowfield) elemental composition for the current small ablation injection rates, where the surface elemental composition is a mixture of freestream and ablator elements. Therefore, the coupled laminar and turbulent flow solutions with radiation and ablation are obtained by using the equilibrium flow chemistry, whereas a nonequilibrium chemistry model is used for solutions without ablation and turbulence. Various computed results are compared with those obtained by the other researchers.

Gupta, Roop N.↗
NASA NTRS2002

Fast Bound Methods for Large Scale Simulation with Application for Engineering Optimization

In this work, we have focused on fast bound methods for large scale simulation with application for engineering optimization. The emphasis is on the development of techniques that provide both very fast turnaround and a certificate of Fidelity; these attributes ensure that the results are indeed relevant to - and trustworthy within - the engineering context. The bound methodology which underlies this work has many different instantiations: finite element approximation; iterative solution techniques; and reduced-basis (parameter) approximation. In this grant we have, in fact, treated all three, but most of our effort has been concentrated on the first and third. We describe these below briefly - but with a pointer to an Appendix which describes, in some detail, the current "state of the art."

Patera, Anthony T.↗
NASA NTRS2004

Magnetic Circuit Model of PM Motor-Generator to Predict Radial Forces

A magnetic circuit model is developed for a PM motor for flywheel applications. A sample motor is designed and modeled. Motor configuration and selection of materials is discussed, and the choice of winding configuration is described. A magnetic circuit model is described, which includes the stator back iron, rotor yoke, permanent magnets, air gaps and the stator teeth. Iterative solution of this model yields flux linkages, back EMF, torque, power, and radial force at the rotor caused by eccentricity. Calculated radial forces are then used to determine motor negative stiffness.

McLallin, Kerry↗
NASA NTRS2004

Cirrus Clouds Optical, Microphysical and Radiative Properties Observed During Crystal-Face Experiment: I. A Radar-Lidar Retrieval System

A method of retrieving cloud microphysical properties using combined observations from both cloud radar and lidar is introduced. This retrieval makes use of an improvement to the traditional optimal estimation retrieval method, whereby a series of corrections are applied to the state vector during the search for an iterative solution. This allows faster convergence to a solution and is less processor intensive. The method is first applied to a synthetic cloud t o demonstrate its validity, and it is shown that the retrieval reliably reproduces vertical profiles of ice water content. The retrieval method is then applied to radar and lidar observations from the CRYSTAL-FACE experiment, and vertical profiles of ice crystal diameter, number concentration, and ice water content are retrieved for a cirrus cloud layers observed one day of that experiment. The validity of the relationship between visible extinction coefficient and radar reflectivity was examined. While synthetic tests showed such a functional relationship, the measured data only partially supported such a conclusion. This is due to errors in the forward model (as explained above) as well as errors in the data sets, including possible mismatch between lidar and radar profiles or errors in the optical depth. Empirical relationships between number concentrations and mean particle diameter were also examined. The results indicate that a distinct and robust relationship exists between these retrieved quantities and it is argued that such a relationship is more than an artifact of the retrieval process offering insight into the nature of the microphysical processes taking place in cirrus.

Mitrescu, C.↗
NASA NTRS1989

Analysis of an algae-based CELSS. Part 1: model development

A steady state chemical model and computer program have been developed for a life support system and applied to trade-off studies. The model is based on human demand for food and oxygen determined from crew metabolic needs. The model includes modules for water recycle, waste treatment, CO2 removal and treatment, and food production. The computer program calculates rates of use and material balance for food. O2, the recycle of human waste and trash, H2O, N2, and food production supply. A simple non-iterative solution for the model has been developed using the steady state rate equations for the chemical reactions. The model and program have been used in system sizing and subsystem trade-off studies of a partially closed life support system.

Models, Chemical↗