Search NASA⌕ Search

SEARCH · Search NASA

Results for “equation solving”

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 163 records · Page 9

Ablation and radiation coupled viscous hypersonic shock layers, volume 2

A computer program can be used to determine the radiative flux and flux divergence through a nonisothermal planar slab of gas is described. The program considers species typical of air at high temperatures and nylon or carbon phenolic ablation products. The equations solved are for a variable optical depth line and continuum gas model using species cross sections and line widths as the basic data.

Engel, C. D.↗

Elliptical cutouts in cylindrical shells

The stress concentrations due to an elliptical cutout in a shallow cylindrical shell are determined by means of a system of singular integral equations, solved numerically. Three loading conditions are considered: tension, and internal pressure. Results confirm those obtained by previous authors by other methods and extend the range of parameters over which results are available.

Tsai, C.-J.↗

Normal and shear impact of layered composite with a crack - Dynamic stress intensification

The dynamic response of a layered composite under normal and shear impact is analyzed by assuming that the composite contains an initial flaw in the matrix material. One of the objectives was to develop an analytical method for determining dynamic stress solutions which should lead to a numerical method which utilizes Fourier transform for the space variable and Laplace transform for the time variable. The time-dependent angle loading is separated into two parts: a symmetric and a skew-symmetric with reference to the crack plane. By superposition, the transient boundary conditions consist of applying normal and shear tractions to a crack embedded in a layered composite; one phase of the composite could represent the fiber while the other could be the matrix. Mathematically, these conditions reduce the problem to a system of dual integral equations solved in the transform plane for the transform of the dynamic stress-intensity factor.

Sih, G. C.↗

Numerical simulation of three-dimensional boattail afterbody flow fields

The thin shear layer approximations of the three-dimensional, compressible Navier-Stokes equations are solved for subsonic, transonic, and supersonic flow over axisymmetric boattail bodies at moderate angles of attack. The plume is modeled by a solid body configuration identical to those used in experimental tests. An implicit algorithm of second-order accuracy is used to solve the equations on the ILLIAC IV computer. The turbulence is expressed by an algebraic model applicable to three-dimensional flow fields with moderate separation. The computed results compare favorably with three different sets of experimental data reported by Reubush, Shrewsbury, and Benek, respectively

Deiwert, G. S.↗

The optimization of convergence for Chebyshev polynomial methods in an unbounded domain

Grosch and Orszag (1977) have performed a numerical analysis of the problem of solving differential equations in a semiinfinite or infinite domain using Chebyshev polynomials. The principal limitation of the conducted study was that it was entirely empirical. Various differential equations were solved in different ways and the numbers were compared. The present investigation has the objective to extend the studies conducted by Grosch and Orszag by deriving asymptotic approximations to the Chebyshev coefficients of simple model functions. This approach makes it possible to conduct more systematic comparisons of different methods, extend the range of comparisons, and, perhaps most important, give simple analytic formulas for choosing the optimum domain size or mapping parameter L for various situations.

Boyd, J. P.↗

Three-Dimensional, Subsonic, Turbulent Juncture Region Flow

CMC3DPNS predicts three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Versatile program for comprehensive flow-field analysis applies finite-element methodology to nonlinear field problems. Order-ofmagnitude analysis of subsonic, three-dimensional, steady time-averaged NavierStokes equations for semibounded aerodynamic juncture geometries yields para bolic Navier-Stokes equations. CMC3DPNS solves these equations to predict juncture-region flow.

Baker, A. J.↗

Thundercloud current determination from measurements at the earth's surface

For the analysis of the actual dominant mechanism or mechanisms of electrification in a given storm cell, it is of great importance to make observations which can be used to determine the configuration and the magnitudes of the currents and impedances of the generator systems in active thunderclouds. Krider and Musser (1982) have measured Maxwell currents below Florida thunderclouds and shown that their patterns are consistent with simple generator geometries. It is, therefore, of great interest to study the extent to which a set of Maxwell current measurements can be inverted to obtain the generator currents and locations. Nisbet (1983) has developed an electrical model for active thunderclouds which solves equations for current densities and electric fields. The present investigation is concerned with the extent to which such a model can be used to derive data about the cloud generator magnitudes and locations, taking into account also the conductivity profile from ground-based measurements or estimates of the Maxwell current density.

Nisbet, J. S.↗

Predicting False Lock in Phase-Locked Loops

Theoretical paper discusses advances in mathematical analysis of phase-locked loops. Presents new results in prediction of false locking. Interest to users of phase-lock circuits and to researchers seeking ways to detect or avoid false lock. Equations solved numerically by Newton-Raphson technique, with Jacobian computed by finite-difference scheme. Algorithm produces results limited only by precision of computer on which executed.

Reed, Bill↗

Calculating Transport Of Moisture Through Honeycomb Panels

Diffusion equations solved for steady-state and transient conditions. HUMID is computer program for predicting moisture gradients and internal bursting pressure in composite honeycomb-sandwich panels. Designed to monitor ebb and flow of moisture in Space Shuttle Orbiter payload bay doors, but applicable for general use. Written in FORTRAN.

Martin, R. T.↗

Wave phenomena in a high Reynolds number compressible boundary layer

The behavior of spatially unstable waves in a high Reynolds number compressible laminar boundary layer is investigated by solution of the laminar two-dimensional compressible Navier-Stokes equations (solved to fourth-order accuracy) over a flat plate with a fluctuating disturbance generated at the inflow. A significant nonlinear distortion is produced, in qualitative agreement with experimental data. It is shown that increasing compressibility can significantly stabilize the flow over a flat plate, and that the mechanism of phase cancellation is a viable mechanism for the control of growing disturbances.

Bayliss, A.↗

Piezoviscosity In Lubrication Of Nonconformal Contacts

Developments in theory of lubrication. Analysis of piezoviscous-rigid regime of lubrication of two ellipsoidal contacts. Begins with Reynolds equation for point contact. Equation nondimensionalized using Roelands empirical formula and Dowson and Higginson formula. Equation solved numerically. Solutions obtained for full spectrum of conditions to find effects of dimensionless load, speed, parameters of lubricated and lubricating materials, and angle between direction of rolling and direction of entrainment of lubricant.

Jeng, Yeau-Ren↗

Comparison of two- and three-dimensional Navier-Stokes solutions with NASA experimental data for CAST-10 airfoil

The two-dimensional (2-D) and three-dimensional Navier-Stokes equations are solved for flow over a NAE CAST-10 airfoil model. Recently developed finite-volume codes that apply a multistage time stepping scheme in conjunction with steady state acceleration techniques are used to solve the equations. Two-dimensional results are shown for flow conditions uncorrected and corrected for wind tunnel wall interference effects. Predicted surface pressures from 3-D simulations are compared with those from 2-D calculations. The focus of the 3-D computations is the influence of the sidewall boundary layers. Topological features of the 3-D flow fields are indicated. Lift and drag results are compared with experimental measurements.

Swanson, R. Charles↗

Transient finite element computations on the transputer system

The aim was to study the solution of transient finite element problems on the Transputer system of parallel processors. The central difference time integration rule was used so that no equation solving was necessary. Also investigated was subcycling time integration which uses different time steps in different subdomains of the finite element mesh. A one-dimensional bar problem was analyzed using the parallel time integration algorithm. This involves subdividing the bar into subproblems which are assigned to different processors. Results show that the significant speed-up can be obtained through parallel processing. Also subcycling can give an additional speed-up in certain classes of problems. A two-dimensional problem was also examined to evaluate the effect of the communication to computation ratio on solution time.

Smolinski, Patrick J.↗

Knowledge representation to support reasoning based on multiple models

Model Based Reasoning is a powerful tool used to design and analyze systems, which are often composed of numerous interactive, interrelated subsystems. Models of the subsystems are written independently and may be used together while they are still under development. Thus the models are not static. They evolve as information becomes obsolete, as improved artifact descriptions are developed, and as system capabilities change. Researchers are using three methods to support knowledge/data base growth, to track the model evolution, and to handle knowledge from diverse domains. First, the representation methodology is based on having pools, or types, of knowledge from which each model is constructed. In addition information is explicit. This includes the interactions between components, the description of the artifact structure, and the constraints and limitations of the models. The third principle we have followed is the separation of the data and knowledge from the inferencing and equation solving mechanisms. This methodology is used in two distinct knowledge-based systems: one for the design of space systems and another for the synthesis of VLSI circuits. It has facilitated the growth and evolution of our models, made accountability of results explicit, and provided credibility for the user community. These capabilities have been implemented and are being used in actual design projects.

Gillam, April↗

Efficiency of group implicit concurrent algorithms for transient finite element analysis

The performance of group implicit algorithms is assessed on actual concurrent computers. It is shown that, as the number of subdomains is increased, performance enhancements are derived from two sources: the increased parallelism in the computations; and a reduction in equation solving effort. Moreover, these two performance enhancements are synergistic, in the sense that the corresponding speed-ups are multiplied, rather than merely added. Simulations on a 32-node hypercube are presented for which the interprocessor communications efficiencies obtained are consistently in excess of 90 percent.

Ortiz, M.↗

A multidimensional finite element method for CFD

A finite element method is used to solve the equations of motion for 2- and 3-D fluid flow. The time-dependent equations are solved explicitly using quadrilateral (2-D) and hexahedral (3-D) elements, mass lumping, and reduced integration. A Petrov-Galerkin technique is applied to the advection terms. The method requires a minimum of computational storage, executes quickly, and is scalable for execution on computer systems ranging from PCs to supercomputers.

Pepper, Darrell W.↗

Numerical simulation of fluid flow around a scramaccelerator projectile

Numerical simulations of the fluid motion and temperature distribution around a 'scramaccelerator' projectile are obtained for Mach numbers in the 5-10 range. A finite element method is used to solve the equations of motion for inviscid and viscous two-dimensional or axisymmetric compressible flow. The time-dependent equations are solved explicitly, using bilinear isoparametric quadrilateral elements, mass lumping, and a shock-capturing Petrov-Galerkin formulation. Computed results indicate that maintaining on-design performance for controlling and stabilizing oblique detonation waves is critically dependent on projectile shape and Mach number.

Pepper, Darrell W.↗

Unified Formulation Of Dynamics Of Robot Arms

Report presents analysis of variety of algorithms developed to solve equations of motion of multiple-rigid-link robot arms. Techniques of spatial-operator algebra used to derive algorithms concisely and study relationships among them.

Jain, Abhinandan↗