Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximation methods”

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 1,621 records · Page 90

Polynomial approximation of Poincare maps for Hamiltonian system

Different methods are proposed and tested for transforming a non-linear differential system, and more particularly a Hamiltonian one, into a map without integrating the whole orbit as in the well-known Poincare return map technique. We construct piecewise polynomial maps by coarse-graining the phase-space surface of section into parallelograms and using either only values of the Poincare maps at the vertices or also the gradient information at the nearest neighbors to define a polynomial approximation within each cell. The numerical experiments are in good agreement with both the real symplectic and Poincare maps.

Froeschle, Claude↗

Calculations of molecular ionization energies using a self-consistent-charge Hartree-Fock-Slater method

A numerical-variational method for performing self-consistent molecular calculations in the Hartree-Fock-Slater (HFS) model is presented. Molecular wavefunctions are expanded in terms of basis sets constructed from numerical HFS solutions of selected one-center atomlike problems. Binding energies and wavefunctions for the molecules are generated using a discrete variational method for a given molecular potential. In the self-consistent-charge (SCC) approximation to the complete self-consistent-field (SCF) method, results of a Mulliken population analysis of the molecular eigenfunctions are used in each iteration to produce 'atomic' occupation numbers. The simplest SCC potential is then obtained from overlapping spherical atomlike charge distributions. Molecular ionization energies are calculated using the transition-state procedure; results are given for CO, H2O, H2S, AlCl, InCl, and the Ni5O surface complex. Agreement between experimental and theoretical ionization energies for the free-molecule valence levels is generally within 1 eV. The simple SCC procedure gives a reasonably good approximation to the molecular potential, as shown by comparison with experiment, and with complete SCF calculations for CO, H2O, and H2S.

Rosen, A.↗

A three-dimensional incompressible flow simulation method and its application to the Space Shuttle main engine. II Turbulent flow

An implicit finite difference code cast in general curvilinear coordinates is further developed for three-dimensional incompressible turbulent flows. The code is based on the method of pseudocompressibility and utilizes the Beam and Warming implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is further extended to include composite-grids to overcome the excessive computer memory required for solving turbulent flows in complex three-dimensional geometries. A simple turbulence model is proposed for internal flows. The code is being used for the Space Shuttle Main Engine (SSME) internal flow analyses.

Chang, J. L. C.↗

Numerical simulation of flow path in the oxidizer side hot gas manifold of the Space Shuttle main engine

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside the oxidizer side Hot Gas Manifold of the Space Shuttle Main Engine. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is developed to overcome the complexity of the geometry. In the present study, the laminar and turbulent flows in the oxidizer side Hot Gas Manifold have been computed. The study reveals that: (1) there exists large recirculation zones inside the bowl if no vanes are present; (2) strong secondary flows are observed in the transfer tube; and (3) properly shaped and positioned guide vanes are effective in eliminating flow separation.

Lin, S. J.↗

Theoretical Additional Span Loading Characteristics of Wings with Arbitrary Sweep, Aspect Ratio, and Taper Ratio

The Weissinger method for determining additional span loading has been used to find the lift-curve slope, spanwise center of pressure, aerodynamic center location, and span loading coefficients of untwisted and uncambered wings having a wide range of plan forms characterized by various combinations of sweep, aspect ratio, and taper ratio. The results are presented as variations of the aerodynamic characteristics with sweep angle for various values of aspect ratio and taper ratio. Methods are also included for determining induced drag and the approximate effects of compressibility. Despite the limitations of a lifting line method such as Weissinger's, the good agreement found between experimentally and theoretically determined characteristics warrants confidence in the method. In particular, it is believed that trends observed in results of the Weissinger method should be reliable. One of the most significant results showed that for each angle of sweep there is a taper ratio for which aspect ratio has little effect on the span loading and for which the loading is practically elliptical. This elliptic loading is approached at a taper ratio of 1.39 for 30 degree of sweepforward, 0.45 for zero degree of sweepback. (author)

Deyoung, John↗

On the dynamics of approximating schemes for dissipative nonlinear equations

Since one can rarely write down the analytical solutions to nonlinear dissipative partial differential equations (PDE's), it is important to understand whether, and in what sense, the behavior of approximating schemes to these equations reflects the true dynamics of the original equations. Further, because standard error estimates between approximations of the true solutions coming from spectral methods - finite difference or finite element schemes, for example - and the exact solutions grow exponentially in time, this analysis provides little value in understanding the infinite time behavior of a given approximating scheme. The notion of the global attractor has been useful in quantifying the infinite time behavior of dissipative PDEs, such as the Navier-Stokes equations. Loosely speaking, the global attractor is all that remains of a sufficiently large bounded set in phase space mapped infinitely forward in time under the evolution of the PDE. Though the attractor has been shown to have some nice properties - it is compact, connected, and finite dimensional, for example - it is in general quite complicated. Nevertheless, the global attractor gives a way to understand how the infinite time behavior of approximating schemes such as the ones coming from a finite difference, finite element, or spectral method relates to that of the original PDE. Indeed, one can often show that such approximations also have a global attractor. We therefore only need to understand how the structure of the attractor for the PDE behaves under approximation. This is by no means a trivial task. Several interesting results have been obtained in this direction. However, we will not go into the details. We mention here that approximations generally lose information about the system no matter how accurate they are. There are examples that show certain parts of the attractor may be lost by arbitrary small perturbations of the original equations.

Jones, Donald A.↗

Application of Dorodnitsyn's Integral Method to Nonequilibrium Flows over Pointed Bodies

Dorodnitsyn's integral method is used to obtain an approximate solution to the supersonic nonequilibrium flow over pointed bodies with attached shock waves. The partial differential equations governing the flow are converted to an approximate set of ordinary equations, which are solved by numerical integration starting at the body tip. Detailed analytical and numerical results for the first approximation are presented, considering the vibrational relaxation of a diatomic gas over a wedge or cone. It is shown that the first approximation yields: (1) The exact flow-variable gradients at the wedge tip (2) Expressions for the flow-variable gradients at the cone tip which are in agreement with extrapolations of characteristics calculations (3) A good approximate algebraic solution for frozen or equilibrium conical flow (4) An approximate expression for the nonequilibrium-flow stream function which affords a means of obtaining variations across the shock layer of the temperature and vibrational energy. Numerical results for both the wedge and cone compare favorably with identical cases computed by the method of characteristics.

NONEQUILIBRIUM FLOW↗

Receptivity and Transition of Supersonic Boundary Layers Over Swept Wings

The receptivity, stability, and transition of three-dimensional supersonic boundary layers over (1) a swept cylinder, (2) a swept wing with a sharp leading edge, and (3) a swept wing with a blunt leading edge are numerically investigated for a free-stream Mach number of 3. These computations are compared to an earlier experimental and computational study performed by Archambaud et al.1 The steady flow fields with and without roughness elements are obtained by solving the full Navier-Stokes equations. The N-factors computed in this study at the transition onset locations reported in Ref. 1 for flow over the swept cylinder are approximately 16.5 for traveling crossflow disturbances and 9 for stationary disturbances. The N-factors for the traveling crossflow are high based on our past experiences. However, they are comparatively smaller than those reported by Archambaud et al., who found N-factor values in the range of 20 to 25 for traveling disturbances and 13 to 20 for stationary disturbances. Similarly, the N-factors computed in this study for the traveling and stationary disturbances for the flow over the sharp wing are approximately 7 and 2.5, respectively, and for the flow over the blunt wing are 6.5 and 4.8, respectively. Using the envelope method, Archambaud et al. obtained values of approximately 8.0 and 4.0 for the sharp wing case and 16.0 and 12.0 for the blunt wing case.

Balakumar, Ponnampalam↗

Influence of Desert Dust Intrusions on Ground-based and Satellite Derived Ultraviolet Irradiance in Southeastern Spain

The desert dust aerosols strongly affect propagation of solar radiation through the atmosphere, reducing surface irradiance available for photochemistry and photosynthesis. This paper evaluates effects of desert dust on surface UV erythemal irradiance (UVER), as measured by a ground-based broadband UV radiometer and retrieved from the satellite Ozone Monitoring Instrument (OMI) at Granada (southern Spain) from January 2006 to December 2010. The dust effects are characterized by the transmittance ra tio of the measured UVER to the corresponding modeled clear sky value. The transmittance has an exponential dependency on aerosol optical depth (AOD), with minimum values of approximately 0.6 (attenuation of approximately 40%). The OMI UVER algorithm does not account for UV aerosol absorption, which results in overestimation of the ground-based UVER especially during dust episodes with a mean relative difference up to 40%. The application of aerosol absorption post-correction method reduces OMI bias up to approximately 13%. The results highlight great effect of desert dust on the surface UV irradiance in regions like southern Spain, where dust intrusions from Sahara region are very frequent.

Krotkov, Nickolay A.↗

Algorithm Plans Collision-Free Path for Robotic Manipulator

An algorithm has been developed to enable a computer aboard a robot to autonomously plan the path of the manipulator arm of the robot to avoid collisions between the arm and any obstacle, which could be another part of the robot or an external object in the vicinity of the robot. In simplified terms, the algorithm generates trial path segments and tests each segment for potential collisions in an iterative process that ends when a sequence of collision-free segments reaches from the starting point to the destination. The main advantage of this algorithm, relative to prior such algorithms, is computational efficiency: the algorithm is designed to make minimal demands upon the limited computational resources available aboard a robot. This path-planning algorithm utilizes a modified version of the collision-detection method described in "Improved Collision-Detection Method for Robotic Manipulator" (NPO-30356), NASA Tech Briefs, Vol. 27, No. 3 (June 2003), page 72. The method involves utilization of mathematical models of the robot constructed prior to operation and similar models of external objects constructed automatically from sensory data acquired during operation. This method incorporates a previously developed method, known in the art as the method of oriented bounding boxes (OBBs), in which an object is represented approximately, for computational purposes, by a box that encloses its outer boundary. Because many parts of a robotic manipulator are cylindrical, the OBB method has been extended in this method to enable the approximate representation of cylindrical parts by use of octagonal or other multiple-OBB assemblies denoted oriented bounding prisms (OBPs). A multiresolution OBB/OBP representation of the robot and its manipulator arm and a multiresolution OBB representation of external objects (including terrain) are constructed and used in a process in which collisions at successively finer resolutions are detected through computational detection of overlaps between the corresponding OBB and OBP models. For computational efficiency, the process is started at the coarsest resolution and stopped as soon as possible, preferably before reaching the finest resolution. At the coarsest resolution, there is a single OBB enclosing all relevant external objects and a single OBB enclosing the entire robot. At the next finer level of resolution, the coarsest-resolution OBB is divided into two OBBs, and so forth. If no collision is detected at the coarsest resolution, then there is no need for further computation to detect collisions. If a collision is detected at the coarsest resolution, then tests for collisions are performed at the next finer level of resolution. This process is continued to successively finer resolutions until either no more collisions are detected or the finest resolution is reached.

Backes, Paul↗

Curing Composite Materials Using Lower-Energy Electron Beams

In an improved method of fabricating composite-material structures by laying up prepreg tapes (tapes of fiber reinforcement impregnated by uncured matrix materials) and then curing them, one cures the layups by use of beams of electrons having kinetic energies in the range of 200 to 300 keV. In contrast, in a prior method, one used electron beams characterized by kinetic energies up to 20 MeV. The improved method was first suggested by an Italian group in 1993, but had not been demonstrated until recently. With respect to both the prior method and the present improved method, the impetus for the use of electron- beam curing is a desire to avoid the high costs of autoclaves large enough to effect thermal curing of large composite-material structures. Unfortunately, in the prior method, the advantages of electron-beam curing are offset by the need for special walls and ceilings on curing chambers to shield personnel from x rays generated by impacts of energetic electrons. These shields must be thick [typically 2 to 3 ft (about 0.6 to 0.9 m) if made of concrete] and are therefore expensive. They also make it difficult to bring large structures into and out of the curing chambers. Currently, all major companies that fabricate composite-material spacecraft and aircraft structures form their layups by use of automated tape placement (ATP) machines. In the present improved method, an electron-beam gun is attached to an ATP head and used to irradiate the tape as it is pressed onto the workpiece. The electron kinetic energy between 200 and 300 keV is sufficient for penetration of the ply being laid plus one or two of the plies underneath it. Provided that the electron-beam gun is properly positioned, it is possible to administer the required electron dose and, at the same time, to protect personnel with less shielding than is needed in the prior method. Adequate shielding can be provided by concrete walls 6 ft (approximately equal to 1.8 m) high and 16 in. (approximately equal to 41 cm) thick, without a ceiling. The success of the present method depends on the use of a cationic epoxy as the matrix material in the prepreg tape, heating the prepreg tape to a temperature of 50 C immediately prior to layup, and exposing the workpiece to an electron-beam dose of approximately 2 Mrad. Experiments have shown that structures fabricated by the present method have the same mechanical properties as those of nominally identical structures fabricated by the prior method with electron beams of 3 to 4 MeV.

Byrne, Catherine A.↗

The anisotropic radiative transfer problem in optically thick, strongly magnetized plasma - A comparison of results

The accuracy of four different methods for solving the coherent radiative transfer problem in a strongly magnetized plasma is compared. Specific attention is given to the case of a semi-infinite homogeneous plasma with a magnetic field perpendicular to the surface. The four approaches evaluated were: diffusion approximation schemes; the generalized H-function method of Silant'ev (1982) the Feautrier method of Nagel (1981) and the integral equation approach of Meszares and Bonazzola (1981). It is found that the coupled diffusion approximation schemes of Nagel (1981) and Kaminker et al. (1982) provided satisfactory results and were more efficient in many cases than direct numerical methods.

Pavlov, G. G.↗

Electron-Hydrogen Elastic Scattering

Scattering by single-electron systems is always of interest because the wave function of the target is known exactly. Various approximations have been employed to take into account distortion produced in the target. Among them are the method of polarized orbitals and the close coupling approximation. Recently, e-H and e-He+ S-wave scattering in the elastic region has been studied using the Feshbach projection operator formalism. In this approach, the usual Hartree-Fock and exchange potentials are augmented by an optical potential and the resulting phase shifts have rigorous lower bounds. Now this method is being applied to the e-H P-wave scattering in the elastic region. The number of terms in the Hylleraas-type wave function for the 1,3 P phase shifts is 84 and the resulting phase shifts (preliminary) are given. The results have been given up to five digits because to that accuracy they are rigorous lower bounds. They are in general agreement with the variational (VAR) results of Armstead, and those obtained from the intermediate energy R-matrix method (RM) of Scholz et al., and the finite element method (FEM) of Botero and Shertzer. The later two methods do not provide any bounds on phase shifts.

Bhatia, A. K.↗

The Problem of Torsion in Prismatic Members of Circular Segmental Cross Section

The problem is solved by approximation, by setting up a function complying with the differential equation of the stress function, and determining the coefficients appearing in it in such a way that the boundary condition is fulfilled as nearly as possible. For the semicircle, for which the solution is known, the method yields very accurate values; the approximated stress distribution is in good agreement with the accurately computed distribution. Stress and strain measurements indicate that the approximate solution is in sufficiently exact agreement with reality for segmental cross sections.

Weigand, A.↗

Spacecraft Modeling, Attitude Determination, and Control: Quaternion-Based Approach

Spacecraft attitude determination and control is an important part of a spacecraft to achieve its designed mission. As of today, many spacecrafts have been successfully launched, and most of them have performed well as they were designed. Many research papers have been published to address the attitude determination and control design problems. Several text books are available for students to learn the technology and for engineers to use as references. The most popular spacecraft models for attitude determination algorithms and control design methods are the Euler angle models and the quaternion models . The Euler angle models have been proved very efficient as the linearized models are controllable, and all standard linear control system design methods are directly applicable. The drawbacks related to the Euler angle methods are (a) the designs based on linearized models may not globally stabilize the original nonlinear spacecraft, i.e., the design may not work when the attitude of the spacecraft is far away from the point where the linearization is performed; (b) the models depend on the rotational sequences, this can be error prone if several teams work on the same project and they use different rotational sequences; (c) for any rotational sequence, there is a singular point where the model is not applicable; and (d) since most attitude determination methods use quaternion to represent the spacecraft attitude, there is a need to transform quaternion into Euler angles. On the other hand, for quaternion models, people have found controllers that can globally stabilize nonlinear spacecraft systems; the models do not depend on rotational sequences and they have no singular point; and the quaternion is provided by attitude determination system and ready to use. The main problem with the quaternion model based control system design is that the linearized quaternion model is not controllable. Therefore, most published design methods heavily rely on Lyapunov functions for the nonlinear spacecraft system. But there is no systematic way to obtain a desired Lyapunov functions. Moreover, the Lyapunov function based designs focus on the closed-loop system stability but pay little attention to the closed-loop system performance. In a series of papers, the author proposed some reduced quaternion models which lead to some controllable linearized spacecraft models. Therefore, all standard linear system theory can be directly applied to analyze and design the spacecraft control systems. We showed that, in some cases, the designed control system is not only optimal for the linearized system, but also globally stabilize the original nonlinear system . Clearly, the reduced quaternion models do not depend on rotational sequences. Due to the special structure of the linearized spacecraft model, some most important design methods, such as LQR design and robust pole assignment design are very simple, enjoy the analytical solutions for some problems, have direct connection to the performance measures, such as settling time, rising time , and percentage of overshoot . All these features are attractive for high quality control system designs. The idea mentioned above is then extended to more spacecraft control problems using specific actuators such as magnetic torque bars and control momentum gyroscopes. These types of actuators may not provide exactly desired torques. Most existing methods use different conversions to get approximate solutions, meaning that these actuators may generate a torque close to but not equal to the desired one. Using the reduced quaternion models that incorporate the actuators into the system model, the control inputs are not torques but the operational parameters. The main benefit of this idea is that the control actions are not approximate but accurate. As all actuators have their operational limit, design with input constraints are also considered in this book by using recently developed interior-point optimization techniques. This book grows up from my research on the spacecraft attitude determination and control design methods in more than a decade which is focused on using reduced quaternion models because of their merits stated above. The book provides all necessary background materials on orbital dynamics, rotations and quaternion, frequently used reference frames, transformations between reference frames, space environment and disturbance torques, ephemeris astronomical vector calculations and measurement instruments, spacecraft control actuators and their models, so that the readers will get a global picture and can apply all these information into the spacecraft system modeling, attitude determination, and spacecraft control system designs, which is the main purpose of this book. This book is different from existing books in that we focus on quaternion based spacecraft control system designs and we consider only attitude control system design related problems, from spacecraft modeling, to attitude determination and estimation, to control system design method selection, to control algorithm development, and to the simulation of the control system designs. Moreover, this book addresses different attitude control tasks in the spacecraft life cycle, including spacecraft maneuver, orbit raising, attitude control, and rendezvous. Finally, this book emphasizes the state space design methods rather than the classical frequency design methods.

Attitude determination and control↗

Minimizing Weight of Structural Designs

Approximation Concepts Code for Efficient Structural Synthesis (ACCESS3) combines methods to offer structural system designer powerful, efficient tool for synthesis of minimum-weight designs. Approximation concepts convert general structural synthesis problem into sequence of explicit problems of separable algebraic form; dual method exploits separable form to construct sequence of explicit dual functions.

Schmit, L. A., Jr.↗

Determination of the pure silicon monocarbide content of silicon carbide and products based on silicon carbide

Experience has shown that different methods of analysis of SiC products give different results. Methods identified as AFNOR, FEPA, and manufacturer P, currently used to detect SiC, free C, free Si, free Fe, and SiO2 are reviewed. The AFNOR method gives lower SiC content, attributed to destruction of SiC by grinding. Two products sent to independent labs for analysis by the AFNOR and FEPA methods showed somewhat different results, especially for SiC, SiO2, and Al2O3 content, whereas an X-ray analysis showed a SiC content approximately 10 points lower than by chemical methods.

Prost, L.↗

The large discretization step method for time-dependent partial differential equations

A new method for the acceleration of linear and nonlinear time dependent calculations is presented. It is based on the Large Discretization Step (LDS) approximation, defined in this work, which employs an extended system of low accuracy schemes to approximate a high accuracy discrete approximation to a time dependent differential operator. Error bounds on such approximations are derived. These approximations are efficiently implemented in the LDS methods for linear and nonlinear hyperbolic equations, presented here. In these algorithms the high and low accuracy schemes are interpreted as the same discretization of a time dependent operator on fine and coarse grids, respectively. Thus, a system of correction terms and corresponding equations are derived and solved on the coarse grid to yield the fine grid accuracy. These terms are initialized by visiting the fine grid once in many coarse grid time steps. The resulting methods are very general, simple to implement and may be used to accelerate many existing time marching schemes.

Haras, Zigo↗