Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Based Procedures”

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,189 records · Page 66

Time variant analysis of large scale constrained rotorcraft systems dynamics - An exploitation of IBM-3090 vector-processor's pipe-lining feature

A generalized algorithmic procedure is presented for handling the constraints in transmissions, which are treated as a multibody system of interconnected rigid/flexible bodies. The type of constraints are classified based on the interconnection of the bodies, assuming one or more points of contact to exist between them. The method is explained through flow charts and configuration/interaction tables. A significant increase in speed of execution is achieved by vectorizing the developed code in computationally intensive areas. The study of an example consisting of two meshing disks rotating at high angular velocity is carried out. The dynamic behavior of the constraint forces associated with the generalized coordinates of the system are plotted by selecting various modes. Applications are intended for the study of dynamic and subsequent prediction of constraint forces at the gear teeth contacting points in helicopter transmissions with the aim of improving performance dependability.

Amirouche, F. M. L.↗

A PDE Sensitivity Equation Method for Optimal Aerodynamic Design

The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Borggaard, Jeff↗

On the Use of CAD-Native Predicates and Geometry in Surface Meshing

Several paradigms for accessing CAD geometry during surface meshing for CFD are discussed. File translation, inconsistent geometry engines and non-native point construction are all identified as sources of non-robustness. The paper argues in favor of accessing CAD parts and assemblies in their native format, without translation, and for the use of CAD-native predicates and constructors in surface mesh generation. The discussion also emphasizes the importance of examining the computational requirements for exact evaluation of triangulation predicates during surface meshing. The native approach is demonstrated through an algorithm for the generation of closed manifold surface triangulations from CAD geometry. CAD parts and assemblies are used in their native format, and a part's native geometry engine is accessed through a modeler-independent application programming interface (API). In seeking a robust and fully automated procedure, the algorithm is based on a new physical space manifold triangulation technique specially developed to avoid robustness issues associated with poorly conditioned mappings. In addition, this approach avoids the usual ambiguities associated with floating-point predicate evaluation on constructed coordinate geometry in a mapped space. The technique is incremental, so that each new site improves the triangulation by some well defined quality measure. The algorithm terminates after achieving a prespecified measure of mesh quality and produces a triangulation such that no angle is less than a given angle bound, a or greater than pi - 2alpha. This result also sets bounds on the maximum vertex degree, triangle aspect-ratio and maximum stretching rate for the triangulation. In addition to the output triangulations for a variety of CAD parts, the discussion presents related theoretical results which assert the existence of such an angle bound, and demonstrate that maximum bounds of between 25 deg and 30 deg may be achieved in practice.

Aftosmis, M. J.↗

Practical gust load alleviation and flutter suppression control laws based on a LQG methodology

A modified linear quadratic Gaussian (LQG) synthesis procedure has been used to design low-order robust multiloop controllers for a flexible airplane. The introduction of properly constructed fictitious Gauss-Markov processes in the control loops allowed meeting classical frequency-domain stability criteria using the direct synthesis procedures of modern time-domain control theory. Model reduction was used to simplify the control laws to the point where they could be easily implemented on onboard flight computers. These control laws provided excellent gust load and flutter mode control with good stability margins and compared very favorably to other control laws synthesized by the classical root-locus technique.

Gangsaas, D.↗

Parabolized reduced Navier-Stokes computational techniques

A review is presented of methods in which composite or reduced Navier-Stokes (RNS) equations are treated with a pressure-gradient-based flux-vector splitting. The methods are similar to large Re asymptotic formulations, and streamwise diffusion terms are ignored in favor of an explicit deferred corrector based on higher-order diffusion terms. The methods can be used for 2D and 3D supersonic flows in which the effects of real gas are incorporated. Several examples of the procedure are given, and subsonic and supersonic flows are handled well with relaxation procedures that incorporate multigrid acceleration. Effective solutions are described for problems ranging from incompressible flows and supersonic flows to sharp shocks and reverse-flow capturing.

Rubin, Stanley G.↗

F-16XL and F-18 High Speed Acoustic Flight Test Databases

This report presents the recorded acoustic data and the computed narrow-band and 1/3-octave band spectra produced by F-18 and F-16XL aircraft in subsonic flight over an acoustic array. Both broadband-shock noise and turbulent mixing noise are observed in the spectra. Radar and c-band tracking systems provided the aircraft position which enabled directivity and smear angles from the aircraft to each microphone to be computed. These angles are based on source emission time and thus give some idea about the directivity of the radiated sound field due to jet noise. A follow-on static test was also conducted where acoustic and engine data were obtained. The acoustic data described in the report has application to community noise analysis, noise source characterization and validation of prediction models. A detailed description of the signal processing procedures is provided. Follow-on static tests of each aircraft were also conducted for which engine data and far-field acoustic data are presented.

Kelly, J. J.↗

Safety Aboard the International Space Station

As with any task that NASA takes on, safety is of utmost importaqce. There are pages of safety codes and procedures that must be followed before any idea can be brought to life. Unfortunately, the International Space Station s (ISS) safety regulations and procedures are based on lg standards rather than on Og. To aide in making this space age home away from home a less hazardous environment, I worked on several projects revolving around the dangers of flammable items in microgravity. The first task I was assigned was to track flames. This involves turning eight millimeter video recordings, of tests run in the five second drop tower, into avi format on the computer. The footage is then compressed and altered so that the flame can be seen more clearly. Using another program called Spotlight, line profiles were used to collect data describing the luminescence of the flame at different points. These raw data are saved as text files and run trough a macro so that a Matlab program can analyze it. By fitting the data to a curve and determining the areas of brightest luminescence, the behavior of the flame can be recorded numerically. After entering the data into a database, researchers can come back later and easily get information on flames resulting from different gas and liquid mixtures in microgravity. I also worked on phase two of the FATE project, which deals with safety aboard the ISS. This phase involves igniting projected droplets and determining how they react with secondary materials. Such simulations represent, on a small scale, the spread of onboard fires due to the effervescence of burning primary materials. I set up existing hardware to operate these experiments and ran tests with it, photographing the results. I also made CAD drawings of the apparatus and the area available on the (SF)2 rig for it to fit into. The experiment will later be performed on the KC-135, and the results gathered will be used to reanalyze current safety standards for the ISS, including the distance of required separation for flammable materials. Additional information is included in the original extended abstract.

Mintz, Shauna M.↗

Expansion of Tabulated Scattering Matrices in Generalized Spherical Functions

An efficient way to solve the vector radiative transfer equation for plane-parallel turbid media is to Fourier-decompose it in azimuth. This methodology is typically based on the analytical computation of the Fourier components of the phase matrix and is predicated on the knowledge of the coefficients appearing in the expansion of the normalized scattering matrix in generalized spherical functions. Quite often the expansion coefficients have to be determined from tabulated values of the scattering matrix obtained from measurements or calculated by solving the Maxwell equations. In such cases one needs an efficient and accurate computer procedure converting a tabulated scattering matrix into the corresponding set of expansion coefficients. This short communication summarizes the theoretical basis of this procedure and serves as the user guide to a simple public-domain FORTRAN program.

scattering matrix↗

TPSAS-NF1676L-13122-DND

An approach for assessing the delamination propagation and growth capabilities in commercial finite element codes was developed and demonstrated for the Virtual Crack Closure Technique (VCCT) implementations in ABAQUS?. The Double Cantilever Beam (DCB) specimen was chosen as an example. First, benchmark results to assess delamination propagation capabilities under static loading were created using models simulating specimens with different delamination lengths. Second, starting from an initially straight front, the delamination was allowed to propagate based on the algorithms implemented in the commercial finite element software. The load-displacement relationship obtained from the propagation analysis results and the benchmark results were compared. Good agreements could be achieved by selecting the appropriate input parameters, which were determined in an iterative procedure. A benchmark case to assess the fatigue growth prediction capabilities was also created. Input control parameters were varied to study the effect on the computed delamination growth. The benchmarks enabled the selection of the appropriate input parameters that yielded good agreements between the results obtained from the automated analysis and the benchmark results. Overall, the results are encouraging but further assessment for mixed-mode delamination is required.

Ronald Krueger↗

An Improved Ionospheric Correction Method for Wide-Area Augmentation Systems

A powerful approach for generating ionospheric corrections in wide area differential GPS applications has been developed, that can be applied to the Federal Aviation Administration's Wide Area Augmentation System (WAAS). This approach has been used to support ionospheric calibration for NASA's Deep Space Network and will be supporting real-time operations for global ionospheric specification and, possibly, single-frequency satellite altimeter calibration. It is a real-time, grid-based technique relying on a computationally efficient Kalman-type filter to produce accurate, smoothly varying ionospheric correction maps over the coverage area. Formal error maps are also computed, providing vertical delay errors over the WAAS grid, which are useful in integrity monitoring. This solution to the ionospheric correction problem is self-calibrating, since GPS transmitter inter-frequency biases are obtained as a by-product of the mapping procedure. This contrasts with other techniques in which bias values must be provided from some additional source.

Wide-Area↗

Analysis of high-aspect-ratio jet-flap wings of arbitrary geometry

An analytical technique to compute the performance of an arbitrary jet-flapped wing is developed. The solution technique is based on the method of Maskell and Spence in which the well-known lifting-line approach is coupled with an auxiliary equation providing the extra function needed in jet-flap theory. The present method is generalized to handle straight, uncambered wings of arbitrary planform, twist, and blowing (including unsymmetrical cases). An analytical procedure is developed for continuous variations in the above geometric data with special functions to exactly treat discontinuities in any of the geometric and blowing data. A rational theory for the effect of finite wing thickness is introduced as well as simplified concepts of effective aspect ratio for rapid estimation of performance.

Lissaman, P. B. S.↗

Analytical Derivation of the Vegetation Optical Depth from the Microwave Polarization Difference Index

A numerical solution for the canopy optical depth in an existing microwave-based land surface parameter retrieval model is presented. The optical depth is derived from the microwave polarization difference index and the dielectric constant of the soil. The original procedure used an approximation in the form of a logarithmic decay function to define this relationship, and was derived through a series of lengthy polynomials. These polynomials had to be recalculated when the scattering albedo or antenna incidence angle changes. The new procedure is computationally more efficient and accurate.

Meesters, Antoon G. C. A.↗

F-ANG+: A 3-D Augmented-Stencil Face-Averaged Nodal-Gradient Cell-Centered Finite-Volume Method for Hypersonic Flows

We describe the extension of a 2-D simplified face-averaged nodal-gradient (F-ANG) method to 3-D and demonstrate that the 3-D simplified F-ANG method is accomplished by augmenting the nodecentered gradient least squares stencil. This augmented stencil F-ANG method is shown to result in advection and diffusion schemes that are stable for hexahedral, prismatic, pyramidal and tetrahedral cells without having to resort to cell-averaged nodal gradients. In addition, we describe the modifications to the augmented stencil required to support the use of wall function boundary conditions. Finally we describe a consistent, face-stencil based multi-dimensional limiter procedure (MLP), and show it to be fully consistent and compatible with the linearity-preserving unstructured- MUSCL (LP-U-MUSCL) scheme for all values of kappa. These methods and schema are implemented in the cell-centered finite-volume code VULCAN-CFD, which is then used to investigate whether the robustness improvements demonstrated in 2-D carry over to 3-D by computing hypersonic flows using mixed-element grids as well as highly adapted tetrahedral grids.

Weighted Least-Squares↗

Lunar Flight Study Series: Volume 7. Earth-Moon Transit Studies Based on Ephemeris Data and Using Best Available Computer Program: Principles for Reducing Earth-Moon Trajectory Analysis to Fundamentals - 2

This report presents an empirical investigation of earth-to-moon trajectories. The primary objective is to uncover relationships that result in the reduction of mission analysis or trajectory analysis problems to fundamentals. It is considered imperative that the results be accurate to the extent that any relationships that are uncovered are not brought about by simplifying or unrealistic assumptions. Consequently, the model used is as accurate as is presently available for use and the computations are performed under strict error controls. Procedures are reasonable from the engineering viewpoint. Some very helpful principles are uncovered. These may be briefly stated as follows: 1) Trajectories of constant flight time to the moon (in the 66 hour to 90 hour range), and arriving within a few minutes of the same time pass over a common point on the moon regardless of the arrival altitude and inclination, This common point of passage for such a family of trajectories is arbitrarily referred to as a VERTEX point; 2) These VERTEX points lie within a bounded region on the lunar surface as arrival time varies throughout the years, for a given flight time or flight time spread; 3) Increasing flight time results primarily in a longitude shift of the VERTEX point for a given arrival time; 4) Launch azimuth should be treated as a strong parameter regarding its influence on the vertex. Launch azimuth variations up to +/- 20 degrees about 90 degrees can shift the vertex as much as 3 degrees in latitude even when one is free to pick launch time, coast time, and S-IVB burn time appropriately with launch azimuth; and 5) The locus of periselenum, as arrival inclination takes on all possible values, is nearly circular about the vertex point for a given arrival altitude. A publication is now in preparation that applies these principles to much advantage in the solution of various trajectory analysis problems.

LUNAR FLIGHT↗

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗

The Mount Wilson solar magnetograph - Scanning and data system

The paper describes a computer-operated image-scanning and data-collection system for the magnetograph at the Mt. Wilson 150-foot Tower telescope. The system is based on a minicomputer with a 32K word core memory and a generalized interface unit for controlling image motion, a keyboard, and an associated television screen. Operation of the solar image guider and the data-collection assembly is outlined along with the observation and data-reduction procedures. Advantages of the system include the ability to move the image in almost any conceivable fashion, a wide choice of integration times, and increased accuracy in magnetic and Doppler calibrations as well as in setting of the magnetic zero level.

Howard, R.↗

Taxonomy of asteroids

A taxonomic system for asteroids is discussed which is based on seven directly observable parameters from polarimetry, spectrophotometry, radiometry, and UBV photometry. The classification scheme is entirely empirical and independent of specific mineralogical interpretations. Five broad classes (designated C, S, M, E, and R), as well as an 'unclassifiable' designation, are defined on the basis of observational data for 523 asteroids. Computer-generated type classifications and derived diameters are given for the 523 asteroids, and the application of the classification procedure is illustrated. Of the 523 asteroids classified, 190 are identified as C objects, 141 as S type, 13 as type M, three as type E, three as type R, 55 as unclassifiable, and 118 as ambiguous. The present taxonomic system is compared with several other asteroid classification systems.

Bowell, E.↗

Aerodynamic optimization by simultaneously updating flow variables and design parameters

The application of conventional optimization schemes to aerodynamic design problems leads to inner-outer iterative procedures that are very costly. An alternative approach is presented based on the idea of updating the flow variable iterative solutions and the design parameter iterative solutions simultaneously. Two schemes based on this idea are applied to problems of correcting wind tunnel wall interference and optimizing advanced propeller designs. The first of these schemes is applicable to a limited class of two-design-parameter problems with an equality constraint. It requires the computation of a single flow solution. The second scheme is suitable for application to general aerodynamic problems. It requires the computation of several flow solutions in parallel. In both schemes, the design parameters are updated as the iterative flow solutions evolve. Computations are performed to test the schemes' efficiency, accuracy, and sensitivity to variations in the computational parameters.

Rizk, M. H.↗