Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “high order”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Scattering of impact wave by a crack in composite plate

The surface responses due to impact load on an infinite uniaxial graphite/epoxy plate with and without delamination cracks are investigated both in time and frequency domain by using a hybrid method combining the finite element discretization of the near-field with boundary integral representation of the field outside a contour enclosing completely the crack. This combined method leads to a set of linear unsymmetric complex matrix equations which are solved to obtain the response in the frequency domain by biconjugate gradient method. The time domain response is then obtained by using an FFT. In order to capture the time-domain characteristics accurately, high order finite elements have been used. Also, both the six node singular elements and eight node transition elements are used around the crack tips to model the crack-tip singularity. It is shown that from the numerical results for surface responses both depth and length of this crack can be identified.

Ju, T. H.↗
NASA NTRS1991

Design of an essentially non-oscillatory reconstruction procedure on finite-element type meshes

An essentially non-oscillatory reconstruction for functions defined on finite-element type meshes was designed. Two related problems are studied: the interpolation of possibly unsmooth multivariate functions on arbitrary meshes and the reconstruction of a function from its average in the control volumes surrounding the nodes of the mesh. Concerning the first problem, we have studied the behavior of the highest coefficients of the Lagrange interpolation function which may admit discontinuities of locally regular curves. This enables us to choose the best stencil for the interpolation. The choice of the smallest possible number of stencils is addressed. Concerning the reconstruction problem, because of the very nature of the mesh, the only method that may work is the so called reconstruction via deconvolution method. Unfortunately, it is well suited only for regular meshes as we show, but we also show how to overcome this difficulty. The global method has the expected order of accuracy but is conservative up to a high order quadrature formula only. Some numerical examples are given which demonstrate the efficiency of the method.

Abgrall, R.↗
NASA NTRS2005

Application of Hilbert-Huang Transform for Improved Defect Detection in Terahertz NDE of Shuttle Tiles

Terahertz NDE is being examined as a method to inspect the adhesive bond-line of Space Shuttle tiles for defects. Terahertz signals are generated and detected, using optical excitation of biased semiconductors with femtosecond laser pulses. Shuttle tile samples were manufactured with defects that included repair regions unbond regions, and other conditions that occur in Shuttle structures. These samples were inspected with a commercial terahertz NDE system that scanned a tile and generated a data set of RF signals. The signals were post processed to generate C-scan type images that are typically seen in ultrasonic NDE. To improve defect visualization the Hilbert-Huang Transform, a transform that decomposes a signal into oscillating components called intrinsic mode functions, was applied to test signals identified as being in and out of the defect regions and then on a complete data set. As expected with this transform, the results showed that the decomposed low-order modes correspond to signal noise while the high-order modes correspond to low frequency oscillations in the signal and mid-order modes correspond to local signal oscillations. The local oscillations compare well with various reflection interfaces and the defect locations in the original signal.

Anastasi, Robert F.↗
NASA NTRS1982

Self-tuning regulator design for adaptive control of aircraft wing/store flutter

The application of the self-tuning regulator concept to adaptively control aircraft wing/store flutter instability is described. A simple design based on a reduced-order aircraft model has been successfully tested on a high-order simulation of an advanced aircraft, and performance was found to be comparable to another design using on-line maximum likelihood identification of plant parameters. The main advantage of the self-tuning regulator is its simplicity, while the main disadvantage is the inadequacy of prior performance guarantees.

Johnson, T. L.↗
NASA NTRS1990

Robust Adaptive Control In Hilbert Space

Paper discusses generalization of scheme for adaptive control of finite-dimensional system to infinite-dimensional Hilbert space. Approach involves generalization of command-generator tracker (CGT) theory. Does not require reference model to be same order as that of plant, and knowledge of order of plant not needed. Suitable for application to high-order systems, main emphasis on adjustment of low-order feedback-gain matrix. Analysis particularly relevant to control of large, flexible structures.

Wen, John Ting-Yung↗
NASA NTRS1990

Applications of ultrashort laser pulses in science and technology; Proceedings of the Meeting, The Hague, Netherlands, Mar. 12, 13, 1990

Topics discussed are on the generation of high-intensity femtosecond lasers, the high-repetition and infrared femtosecond pulses, and physics of semiconductors and applications. Papers are presented on the femtosecond pulse generation at 193 nm; the generation of intense subpicosecond and femtosecond pulses; intense tunable subpicosecond and femtosecond pulses in the visible and infrared, generated by optical parametric oscillators; a high-efficiency high-energy optical amplifier for femtosecond pulses; and the generation of solitons, periodic pulsing, and nonlinearities in GaAs. Other papers are on ultrafast relaxation dynamics of photoexcited carriers in GaAs, high-order optical nonlinear susceptibilities of transparent glasses, subnanosecond risetime high-power pulse generation using photoconductive bulk GaAs devices, femtosecond studies of plasma formation in crystalline and amorphous silicon, and subpicosecond dynamics of hot carrier relaxation in InP and GaAs.

Antonetti, Andre↗
NASA NTRS2004

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel Benchmarks (NPB). In this paper, we present some interesting performance results of ow OpenMP parallel implementation on different architectures such as the SGI Origin2000, SGI Altix, and Cray MTA-2.

Feng, Huiyu↗
NASA NTRS2016

Phase-Retrieval Uncertainty Estimation and Algorithm Comparison for the JWST-ISIM Test Campaign

Phase retrieval, the process of determining the exitpupil wavefront of an optical instrument from image-plane intensity measurements, is the baseline methodology for characterizing the wavefront for the suite of science instruments (SIs) in the Integrated Science Instrument Module (ISIM) for the James Webb Space Telescope (JWST). JWST is a large, infrared space telescope with a 6.5-meter diameter primary mirror. JWST is currently NASA's flagship mission and will be the premier space observatory of the next decade. ISIM contains four optical benches with nine unique instruments, including redundancies. ISIM was characterized at the Goddard Space Flight Center (GSFC) in Greenbelt, MD in a series of cryogenic vacuum tests using a telescope simulator. During these tests, phase-retrieval algorithms were used to characterize the instruments. The objective of this paper is to describe the Monte-Carlo simulations that were used to establish uncertainties (i.e., error bars) for the wavefronts of the various instruments in ISIM. Multiple retrieval algorithms were used in the analysis of ISIM phase-retrieval focus-sweep data, including an iterativetransform algorithm and a nonlinear optimization algorithm. These algorithms emphasize the recovery of numerous optical parameters, including low-order wavefront composition described by Zernike polynomial terms and high-order wavefront described by a point-by-point map, location of instrument best focus, focal ratio, exit-pupil amplitude, the morphology of any extended object, and optical jitter. The secondary objective of this paper is to report on the relative accuracies of these algorithms for the ISIM instrument tests, and a comparison of their computational complexity and their performance on central and graphical processing unit clusters. From a phase-retrieval perspective, the ISIM test campaign includes a variety of source illumination bandwidths, various image-plane sampling criteria above and below the Nyquist- Shannon critical sampling value, various extended object sizes, and several other impactful effects.

Design Analysis↗
NASA NTRS1980

An airborne far-infrared spectrometer for astronomical observations

A Fabry-Perot spectrometer for airborne observations of astronomical spectral lines has been designed for use in the 50 to 200 micron region of the far-infrared. The spectrometer uses a fixed wavelength Fabry-Perot filter cooled to 4.2 K in conjunction with a high order scanning Fabry-Perot at room temperature. This arrangement provides high spectral resolution and also minimizes the thermal background radiation on the detector. The instrument has been successfully used to study all the presently detected fine structure lines between 50 and 120 microns.

Storey, J. W. V.↗
NASA NTRS1985

Numerical synthesis of tri-variate velocity realizations of turbulence

An approach for synthesizing trivariate turbulence velocity field spatial realizations is presented. Some of the spatial frequency characteristics of the random velocity field are described by the von Karman spectrum. The simulation algorithm is based on an efficient autoregressive-moving average (ARMA) scheme involving coefficient square matrices of order three. The determination of the efficient low order ARMA algorithm is preceded by the determination of a suitable high order autoregressive (AR) simulation algorithm. The numerical results are presented in a dimensionless form. Thus, they are applicable for any scale of turbulence.

Spanos, P.-T. D.↗
NASA NTRS2000

An Application of the Quadrature-Free Discontinuous Galerkin Method

The process of generating a block-structured mesh with the smoothness required for high-accuracy schemes is still a time-consuming process often measured in weeks or months. Unstructured grids about complex geometries are more easily generated, and for this reason, methods using unstructured grids have gained favor for aerodynamic analyses. The discontinuous Galerkin (DG) method is a compact finite-element projection method that provides a practical framework for the development of a high-order method using unstructured grids. Higher-order accuracy is obtained by representing the solution as a high-degree polynomial whose time evolution is governed by a local Galerkin projection. The traditional implementation of the discontinuous Galerkin uses quadrature for the evaluation of the integral projections and is prohibitively expensive. Atkins and Shu introduced the quadrature-free formulation in which the integrals are evaluated a-priori and exactly for a similarity element. The approach has been demonstrated to possess the accuracy required for acoustics even in cases where the grid is not smooth. Other issues such as boundary conditions and the treatment of non-linear fluxes have also been studied in earlier work This paper describes the application of the quadrature-free discontinuous Galerkin method to a two-dimensional shear layer problem. First, a brief description of the method is given. Next, the problem is described and the solution is presented. Finally, the resources required to perform the calculations are given.

Lockard, David P.↗
NASA NTRS1988

Robust control with structured perturbations

Two important problems in the area of control systems design and analysis are discussed. The first is the robust stability using characteristic polynomial, which is treated first in characteristic polynomial coefficient space with respect to perturbations in the coefficients of the characteristic polynomial, and then for a control system containing perturbed parameters in the transfer function description of the plant. In coefficient space, a simple expression is first given for the l(sup 2) stability margin for both monic and non-monic cases. Following this, a method is extended to reveal much larger stability region. This result has been extended to the parameter space so that one can determine the stability margin, in terms of ranges of parameter variations, of the closed loop system when the nominal stabilizing controller is given. The stability margin can be enlarged by a choice of better stabilizing controller. The second problem describes the lower order stabilization problem, the motivation of the problem is as follows. Even though the wide range of stabilizing controller design methodologies is available in both the state space and transfer function domains, all of these methods produce unnecessarily high order controllers. In practice, the stabilization is only one of many requirements to be satisfied. Therefore, if the order of a stabilizing controller is excessively high, one can normally expect to have a even higher order controller on the completion of design such as inclusion of dynamic response requirements, etc. Therefore, it is reasonable to have a lowest possible order stabilizing controller first and then adjust the controller to meet additional requirements. The algorithm for designing a lower order stabilizing controller is given. The algorithm does not necessarily produce the minimum order controller; however, the algorithm is theoretically logical and some simulation results show that the algorithm works in general.

Keel, Leehyun↗
NASA NTRS1977

Determination of pole sensitivities by Danilevskii's method

In control theory, a synonymous term for pole sensitivity is eigenvalue sensitivity. Existing methods of calculating eigenvalues are cumbersome, and cannot be trusted for systems roughly greater than tenth order. The method proposed in the present paper is applicable to high-order system. (It has been routinely used to generate eigenvalue sensitivities for systems up to 26th order, using a UNIVAC 1106.) Danilevskii's method is shown to be suitable for performing the necessary evaluations. The result is a rational function that can be used to evaluate the sensitivities for all distinct poles.

Nail, J. B.↗
NASA NTRS1992

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.

Pruett, C. D.↗
NASA NTRS2006

Track-Level-Compensation Look-Up Table Improves Antenna Pointing Precision

This article presents the improvement of the beam-waveguide antenna pointing accuracy due to the implementation of the track-level-compensation look-up table. It presents the development of the table, from the measurements of the inclinometer tilts to the processing of the measurement data and the determination of the threeaxis alidade rotations. The table consists of three axis rotations of the alidade as a function of the azimuth position. The article also presents the equations to determine the elevation and cross-elevation errors of the antenna as a function of the alidade rotations and the antenna azimuth and elevation positions. The table performance was verified using radio beam pointing data. The pointing error decreased from 4.5 mdeg to 1.4 mdeg in elevation and from 14.5 mdeg to 3.1 mdeg in cross-elevation. I. Introduction The Deep Space Station 25 (DSS 25) antenna shown in Fig. 1 is one of NASA s Deep Space Network beam-waveguide (BWG) antennas. At 34 GHz (Ka-band) operation, it is necessary to be able to track with a pointing accuracy of 2-mdeg root-mean-square (rms). Repeatable pointing errors of several millidegrees of magnitude have been observed during the BWG antenna calibration measurements. The systematic errors of order 4 and lower are eliminated using the antenna pointing model. However, repeatable pointing errors of higher order are out of reach of the model. The most prominent high-order systematic errors are the ones caused by the uneven azimuth track. The track is shown in Fig. 2. Manufacturing and installation tolerances, as well as gaps between the segments of the track, are the sources of the pointing errors that reach over 14-mdeg peak-to-peak magnitude, as reported in [1,2]. This article presents a continuation of the investigations and measurements of the pointing errors caused by the azimuth-track-level unevenness that were presented in [1] and [2], and it presents the implementation results. Track-level-compensation (TLC) look-up tables were created for the DSS 25, DSS 26, DSS 34, and DSS 55 antennas. To date, the most complete and detailed results were obtained for the DSS 25 and DSS 55 antennas. In this article, for brevity of presentation, we present the DSS 25 antenna results only. 1 Communications Ground Systems Section. The research described in this publication was carried out by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.

Gawronski, W.↗
NASA NTRS1983

Vibration characteristics of hexagonal radial rib and hoop platforms

Experiment and analysis have been used to characterize the modes of vibration of planar radial rib and hoop hexagonal platforms. Finite element analysis correlated very well with experimental results. The sensitivity of mode shapes and frequencies to cable stiffness and initial tension is presented. Threshold values have been identified, above which changes in cable stiffness do not affect the first few platform vibration modes. Primary vibration modes of the radial rib platform involve beam bending. Vibration modes of the hoop platform exhibit both beam bending and frame bending and torsion. Results indicate for low order polygonal structures, the radial rib concept produced a higher fundamental frequency. For high order polygonal structures, the hoop concept has the potential to achieve a higher fundamental frequency than the radial rib concept.

Belvin, W. K.↗
NASA NTRS1993

A homotopy algorithm for digital optimal projection control GASD-HADOC

The linear-quadratic-gaussian (LQG) compensator was developed to facilitate the design of control laws for multi-input, multi-output (MIMO) systems. The compensator is computed by solving two algebraic equations for which standard closed-loop solutions exist. Unfortunately, the minimal dimension of an LQG compensator is almost always equal to the dimension of the plant and can thus often violate practical implementation constraints on controller order. This deficiency is especially highlighted when considering control-design for high-order systems such as flexible space structures. This deficiency motivated the development of techniques that enable the design of optimal controllers whose dimension is less than that of the design plant. A homotopy approach based on the optimal projection equations that characterize the necessary conditions for optimal reduced-order control. Homotopy algorithms have global convergence properties and hence do not require that the initializing reduced-order controller be close to the optimal reduced-order controller to guarantee convergence. However, the homotopy algorithm previously developed for solving the optimal projection equations has sublinear convergence properties and the convergence slows at higher authority levels and may fail. A new homotopy algorithm for synthesizing optimal reduced-order controllers for discrete-time systems is described. Unlike the previous homotopy approach, the new algorithm is a gradient-based, parameter optimization formulation and was implemented in MATLAB. The results reported may offer the foundation for a reliable approach to optimal, reduced-order controller design.

Collins, Emmanuel G., Jr.↗