Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approximate computing”

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 505 records · Page 28

Computing region moments from boundary representations

The class of all possible formulas for computing arbitrary moments of a region from the region's boundary is derived. The selection of a particular formula depends on the choice of an independent parameter. Several choices of this parameter are explored for region boundaries approximated by polygons. The parameter choice that minimizes computation time for boundaries represented by chain code is derived. Algorithms are presented for computing arbitrary moments for a region from a polygonal approximation of its boundary and for computing low order moments from chain encoded boundaries.

Wilf, J. M.↗

Computational strategies for tire monitoring and analysis

Computational strategies are presented for the modeling and analysis of tires in contact with pavement. A procedure is introduced for simple and accurate determination of tire cross-sectional geometric characteristics from a digitally scanned image. Three new strategies for reducing the computational effort in the finite element solution of tire-pavement contact are also presented. These strategies take advantage of the observation that footprint loads do not usually stimulate a significant tire response away from the pavement contact region. The finite element strategies differ in their level of approximation and required amount of computer resources. The effectiveness of the strategies is demonstrated by numerical examples of frictionless and frictional contact of the space shuttle Orbiter nose-gear tire. Both an in-house research code and a commercial finite element code are used in the numerical studies.

Danielson, Kent T.↗

Integrated aerodynamic-structural design of a transport wing

The integrated aerodynamic-structural design of a subsonic transport wing for minimum weight subject to required range is formulated and solved. The problem requires large computational resources, and two methods are used to alleviate the computational burden. First, a modular sensitivity method that permits the usage of black-box disciplinary software packages, is used to reduce the cost of sensitivity derivatives. In particular, it is shown that derivatives of the aeroelastic response and divergence speed can be calculated without the costly computation of derivatives of aerodynamic influence coefficient and structural stiffness matrices. A sequential approximate optimization is used to further reduce computational cost. The optimization procedure is shown to require a relatively small number of analysis and sensitivity calculations.

Grossman, B.↗

Finite-Element Methods for Real-Time Simulation of Surgery

Two finite-element methods have been developed for mathematical modeling of the time-dependent behaviors of deformable objects and, more specifically, the mechanical responses of soft tissues and organs in contact with surgical tools. These methods may afford the computational efficiency needed to satisfy the requirement to obtain computational results in real time for simulating surgical procedures as described in Simulation System for Training in Laparoscopic Surgery (NPO-21192) on page 31 in this issue of NASA Tech Briefs. Simulation of the behavior of soft tissue in real time is a challenging problem because of the complexity of soft-tissue mechanics. The responses of soft tissues are characterized by nonlinearities and by spatial inhomogeneities and rate and time dependences of material properties. Finite-element methods seem promising for integrating these characteristics of tissues into computational models of organs, but they demand much central-processing-unit (CPU) time and memory, and the demand increases with the number of nodes and degrees of freedom in a given finite-element model. Hence, as finite-element models become more realistic, it becomes more difficult to compute solutions in real time. In both of the present methods, one uses approximate mathematical models trading some accuracy for computational efficiency and thereby increasing the feasibility of attaining real-time up36 NASA Tech Briefs, October 2003 date rates. The first of these methods is based on modal analysis. In this method, one reduces the number of differential equations by selecting only the most significant vibration modes of an object (typically, a suitable number of the lowest-frequency modes) for computing deformations of the object in response to applied forces.

Basdogan, Cagatay↗

Strategy for reflector pattern calculation: Let the computer do the work

Using high frequency approximations, the secondary pattern of a reflector antenna can be calculated by numerically evaluating a radiation integral I(u,v). In recent years, tremendous effort has been expended to reducing I(u,v) to Fourier integrals. These reduction schemes are invariably reflector geometry dependent. Hence, different analyses/computer software development must be carried out for different reflector shapes/boundaries. it is pointed out, that, as the computer power improves, these reduction schemes are no longer necessary. Comparable accuracy and computation time can be achieved by evaluating I(u,v) by a brute force FFT described in this note. Furthermore, there is virtually no restriction on the reflector geometry by using the brute force FFT.

Lam, P. T.↗

Strategy for reflector pattern calculation - Let the computer do the work

Using high frequency approximations, the secondary pattern of a reflector antenna can be calculated by numerically evaluating a radiation integral I(u,v). In recent years, tremendous effort has been expended to reducing I(u,v) to Fourier integrals. These reduction schemes are invariably reflector geometry dependent. Hence, different analyses/computer software development must be carried out for different reflector shapes/boundaries. It is pointed out, that, as the computer power improves, these reduction schemes are no longer necessary. Comparable accuracy and computation time can be achieved by evaluating I(u,v) by a brute force FFT described in this note. Furthermore, there is virtually no restriction on the reflector geometry by using the brute force FFT.

Lam, P. T.↗

Further Examination of a Simplified Model for Positronium-Helium Scattering

While carrying out investigations on Ps-He scattering we realized that it would be possible to improve the results of a previous work on zero-energy scattering of ortho-positronium by helium atoms. The previous work used a model to account for exchange and also attempted to include the effect of short-range Coulomb interactions in the close-coupling approximation. The 3 terms that were then included did not produce a well-converged result but served to give some justification to the model. Now we improve the calculation by using a simple variational wave function, and derive a much better value of the scattering length. The new result is compared with other computed values, and when an approximate correction due to the van der Waals potential is included the total is consistent with an earlier conjecture.

DiRienzi, J.↗

Newton's method applied to finite-difference approximations for the steady-state compressible Navier-Stokes equations

Finite-difference approximations for steady-state compressible Navier-Stokes equations, whose two spatial dimensions are written in generalized curvilinear coordinates and strong conservation-law form, are presently solved by means of Newton's method in order to obtain a lifting-airfoil flow field under subsonic and transonnic conditions. In addition to ascertaining the computational requirements of an initial guess ensuring convergence and the degree of computational efficiency obtainable via the approximate Newton method's freezing of the Jacobian matrices, attention is given to the need for auxiliary methods assessing the temporal stability of steady-state solutions. It is demonstrated that nonunique solutions of the finite-difference equations are obtainable by Newton's method in conjunction with a continuation method.

Bailey, Harry E.↗

Modeling composite and fluffy grains: The effects of porosity

Recent studies of interplanetary and interstellar dust provide evidence that cosmic dust grains are fluffy, composite objects, highlighting the need for models of the electromagnetic scattering by these grains. Effective medium theory (EMT) with Mie-type series solutions has been used to explore the effects of porosity which would be important in composite dust particles. While this indirect approach is both flexible and computationally efficient, it is not necessarily a good approximation. The need for EMT and its rather restrictive assumptions may be circumvented through a direct computation of the scattering properties via finite element methods, such as the discrete dipole approximation (DDA). Recently, the utility of the DDA method has been advanced significantly through improvements in theory, in numerical algorithms, and in computer hardware. Extensive calculations with the DDA method are used here to examine more directly the effects of porosity. A particular emphasis is placed upon developing a valid methodology. For both solid and porous targets we establish both numerical and physical convergence properties over the range of size parameter that is required for our study. DDA cross sections for grains with a range of porosity are compared to those computed by the EMT/series expansion technique to examine the applicability of several mixing rules, including two extensions of the Bruggeman rule. We show that for particles with Rayleigh vacuum inclusions, the extension proposed by Rouleau & Martin is quite successful. We also investigate the effects of larger, non-Rayleigh vacuum inclusions for various levels of porosity and find that they can be significant.

Wolff, Michael J.↗

Comparison of results obtained by solving the radiative transfer equation with an iterative method and a spherical harmonics method

Fluxes and intensities of light scattered by a model atmosphere are computed by a spherical harmonics approximation and by an iterative method of solving the radiative transfer equation and are compared. The large differences in the net fluxes and intensities reported by Dave and Armstrong (1974) for the two methods are reduced here by making a few changes in the iterative routine. Decreasing the polar angle increment from 2 to 1 deg in the iterative method of computing the source function does not improve the results as suggested by Dave and Armstrong.

Bahethi, O. P.↗

Recent CFD Simulations of Shuttle Orbiter Contingency Abort Aerodynamics

Modern Computational Fluid Dynamics (CFD) techniques were used to compute aerodynamic forces and moments of the Space Shuttle Orbiter in specific portions of contingency abort trajectory space. The trajectory space covers a Mach number range of 3.5-15, an angle-of-attack range of 20-60 degrees, an altitude range of 100-190 kft, and several different settings of the control surfaces (elevons, body flap, and speed brake). While approximately 40 cases have been computed, only a sampling of the results is presented here. The computed results, in general, are in good agreement with the Orbiter Operational Aerodynamic Data Book (OADB) data (i.e., within the uncertainty bands) for almost all the cases. However, in a limited number of high angle-of-attack cases (at Mach 15), there are significant differences between the computed results, especially the vehicle pitching moment, and the OADB data. A preliminary analysis of the data from the CFD simulations at Mach 15 shows that these differences can be attributed to real-gas/Mach number effects.

Papadopoulos, Periklis↗

Eigenvalue and eigenvector sensitivity and approximate analysis for repeated eigenvalue problems

A set of computationally efficient equations for eigenvalue and eigenvector sensitivity analysis are derived, and a method for eigenvalue and eigenvector approximate analysis in the presence of repeated eigenvalues is presented. The method developed for approximate analysis involves a reparamaterization of the multivariable structural eigenvalue problem in terms of a single positive-valued parameter. The resulting equations yield first-order approximations of changes in both the eigenvalues and eigenvectors associated with the repeated eigenvalue problem. Examples are given to demonstrate the application of such equations for sensitivity and approximate analysis.

Hou, Gene J. W.↗

Simultaneous quaternion estimation (QUEST) and bias determination

Tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function are presented. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and Earth-pointing reference attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. They converge to the same estimates at the same rate for inertially-fixed attitude, but the new algorithm converges more slowly than the differential correction for Earth-pointing attitude. The slower convergence of the new method for non-zero attitude rates is believed to be due to the use of an inadequate approximation for a partial derivative matrix. The new method requires about twice the computational effort of the differential correction. Improving the approximation for the partial derivative matrix in the new method is expected to improve its convergence at the cost of increased computational effort.

Markley, F. Landis↗

Computational structures technology and UVA Center for CST

Rapid advances in computer hardware have had a profound effect on various engineering and mechanics disciplines, including the materials, structures, and dynamics disciplines. A new technology, computational structures technology (CST), has recently emerged as an insightful blend between material modeling, structural and dynamic analysis and synthesis on the one hand, and other disciplines such as computer science, numerical analysis, and approximation theory, on the other hand. CST is an outgrowth of finite element methods developed over the last three decades. The focus of this presentation is on some aspects of CST which can impact future airframes and propulsion systems, as well as on the newly established University of Virginia (UVA) Center for CST. The background and goals for CST are described along with the motivations for developing CST, and a brief discussion is made on computational material modeling. We look at the future in terms of technical needs, computing environment, and research directions. The newly established UVA Center for CST is described. One of the research projects of the Center is described, and a brief summary of the presentation is given.

Noor, Ahmed K.↗

Feasibility of Actively Cooled Silicon Nitride Airfoil for Turbine Applications Demonstrated

Nickel-base superalloys currently limit gas turbine engine performance. Active cooling has extended the temperature range of service of nickel-base superalloys in current gas turbine engines, but the margin for further improvement appears modest. Therefore, significant advancements in materials technology are needed to raise turbine inlet temperatures above 2400 F to increase engine specific thrust and operating efficiency. Because of their low density and high-temperature strength and thermal conductivity, in situ toughened silicon nitride ceramics have received a great deal of attention for cooled structures. However, the high processing costs and low impact resistance of silicon nitride ceramics have proven to be major obstacles for widespread applications. Advanced rapid prototyping technology in combination with conventional gel casting and sintering can reduce high processing costs and may offer an affordable manufacturing approach. Researchers at the NASA Glenn Research Center, in cooperation with a local university and an aerospace company, are developing actively cooled and functionally graded ceramic structures. The objective of this program is to develop cost-effective manufacturing technology and experimental and analytical capabilities for environmentally stable, aerodynamically efficient, foreign-object-damage-resistant, in situ toughened silicon nitride turbine nozzle vanes, and to test these vanes under simulated engine conditions. Starting with computer aided design (CAD) files of an airfoil and a flat plate with internal cooling passages, the permanent and removable mold components for gel casting ceramic slips were made by stereolithography and Sanders machines, respectively. The gel-cast part was dried and sintered to final shape. Several in situ toughened silicon nitride generic airfoils with internal cooling passages have been fabricated. The uncoated and thermal barrier coated airfoils and flat plates were burner rig tested for 30 min without and with air cooling. Without cooling, the surface temperature of the flat plate reached approximately 2350 F. Starting with computer aided design (CAD) files of an airfoil and a flat plate with internal cooling passages, the permanent and removable mold components for gel casting ceramic slips were made by stereolithography and Sanders machines, respectively. The gel-cast part was dried and sintered to final shape. Several in situ toughened silicon nitride generic airfoils with internal cooling passages have been fabricated. The uncoated and thermal barrier coated airfoils and flat plates were burner rig tested for 30 min without and with air cooling. Without cooling, the surface temperature of the flat plate reached approximately 2350 F. With cooling, the surface temperature decreased to approximately 1910 F--a drop of approximately 440 F. This preliminary study demonstrates that a near-net-shape silicon nitride airfoil can be fabricated and that silicon nitride can sustain severe thermal shock and the thermal gradients induced by cooling and, thus, is a viable candidate for cooled components.

Bhatt, Ramakrishna T.↗

RKH space approximations for the feedback operator in a linear hereditary control system

Computational implementation of feedback control laws for linear hereditary systems requires the approximation of infinite dimensional feedback operators with finite dimensional operators. The dense subspaces of K-polygonal functions in reproducing kernel Hilbert spaces, RKH spaces, suggest finite dimensional approximations of the matrix representations of the control operators. A convergence theorem is developed for the approximations and the numerical implementation of the approximations is discussed.

Reneke, J. A.↗

Investigations on the amount of downwash behind rectangular and elliptical wings

The downwash behind a wing may be computed from the system of vortices created by the wing and which can be substituted for the latter. If the actual rolling-up process of the free vortex sheet is disregarded, the establishment of boundaries for the downwash factors is contingent upon the premise that no rolling up takes place or else that the developed tip vortices start direct on the effective line. The actual downwash factors for the rectangular wing are rather of the approximate magnitude of the values computable with the developed tip vortices starting direct on the wing.

Muttray, H↗

A Boundary Condition Relaxation Algorithm for Strongly Coupled, Ablating Flows Including Shape Change

Implementations of a model for equilibrium, steady-state ablation boundary conditions are tested for the purpose of providing strong coupling with a hypersonic flow solver. The objective is to remove correction factors or film cooling approximations that are usually applied in coupled implementations of the flow solver and the ablation response. Three test cases are considered - the IRV-2, the Galileo probe, and a notional slender, blunted cone launched at 10 km/s from the Earth's surface. A successive substitution is employed and the order of succession is varied as a function of surface temperature to obtain converged solutions. The implementation is tested on a specified trajectory for the IRV-2 to compute shape change under the approximation of steady-state ablation. Issues associated with stability of the shape change algorithm caused by explicit time step limits are also discussed.

Gnoffo, Peter A.↗