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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations

A space-time element method is presented for solving the compressible Euler and Navier-Stokes equations. The proposed formulation includes the variational equation, predictor multi-corrector algorithms and boundary conditions. The variational equation is based on the time-discontinuous Galerkin method, in which the physical entropy variables are employed. A least-squares operator and a discontinuity-capturing operator are added, resulting in a high-order accurate and unconditionally stable method. Implicit/explicit predictor multi-corrector algorithms, applicable to steady as well as unsteady problems, are presented; techniques are developed to enhance their efficiency. Implementation of boundary conditions is addressed; in particular, a technique is introduced to satisfy nonlinear essential boundary conditions, and a consistent method is presented to calculate boundary fluxes. Numerical results are presented to demonstrate the performance of the method.

Shakib, Farzin↗
NASA NTRS1992

A component modes projection and assembly model reduction methodology for articulated, multi-flexible body structures

A two-stage model reduction methodology, combining the classical Component Mode Synthesis (CMS) method and the newly developed Enhanced Projection and Assembly (EP&A) method is proposed in this research. The first stage of this methodology, called COmponent Mode Projection and Assembly model REduction (COMPARE) method, involves the generations of CMS mode sets, such as the MacNeal-Rubin mode sets. These mode sets are then used to reduce the order of each component model in the Rayleigh-Ritz sense. The resultant component models are then combined to generate reduced-order system models at various system configurations. A composite mode set which retains important system modes at all system configurations is then selected from these reduced-order system models. In the second stage, the EP&A model reduction method is employed to reduce further the order of the system model generated in the first stage. The effectiveness of the COMPARE methodology has been successfully demonstrated on a high-order, finite-element model of the cruise-configured Galileo spacecraft.

Lee, Allan Y.↗
NASA NTRS1993

Optical communications and a comparison of optical technologies for a high data rate return link from Mars

The important principles of direct- and heterodyne-detection optical free-space communications are reviewed. Signal-to-noise-ratio (SNR) and bit-error-rate (BER) expressions are derived for both the direct-detection and heterodyne-detection optical receivers. For the heterodyne system, performance degradation resulting from received-signal and local oscillator-beam misalignment and laser phase noise is analyzed. Determination of interfering background power from local and extended background sources is discussed. The BER performance of direct- and heterodyne-detection optical links in the presence of Rayleigh-distributed random pointing and tracking errors is described. Finally, several optical systems employing Nd:YAG, GaAs, and CO2 laser sources are evaluated and compared to assess their feasibility in providing high-data-rate (10- to 1000-Mbps) Mars-to-Earth communications. It is shown that the root mean square (rms) pointing and tracking accuracy is a critical factor in defining the system transmitting laser-power requirements and telescope size and that, for a given rms error, there is an optimum telescope aperture size that minimizes the required power. The results of the analysis conducted indicate that, barring the achievement of extremely small rms pointing and tracking errors (less than 0.2 microrad), the two most promising types of optical systems are those that use an Nd:YAG laser (lambda = 1.064 microns) and high-order pulse position modulator (PPM) and direct detection, and those that use a CO2 laser (lambda = 10.6 microns) and phase shifting keying homodyne modulation and coherent detection. For example, for a PPM order of M = 64 and an rms pointing accuracy of 0.4 microrad, an Nd:YAG system can be used to implement a 100-Mbps Mars link with a 40-cm transmitting telescope, a 20-W laser, and a 10-m receiving photon bucket. Under the same conditions, a CO2 system would require 3-m transmitting and receiving telescopes and a 32-W laser to implement such a link. Other types of optical systems, such as a semiconductor laser systems, are impractical in the presence of large rms pointing errors because of the high power requirements of the 100-Mbps Mars link, even when optimal-size telescopes are used.

Spence, Rodney L.↗
NASA NTRS1992

An exploratory investigation of the flight dynamics effects of rotor rpm variations and rotor state feedback in hover

This paper presents the results of an analytical study conducted to investigate airframe/engine interface dynamics, and the influence of rotor speed variations on the flight dynamics of the helicopter in hover, and to explore the potential benefits of using rotor states as additional feedback signals in the flight control system. The analytical investigation required the development of a parametric high-order helicopter hover model, which included heave/yaw body motion, the rotor speed degree of freedom, rotor blade motion in flapping and lead-lag, inflow dynamics, a drive train model with a flexible rotor shaft, and an engine/rpm governor. First, the model was used to gain insight into the engine/drive train/rotor system dynamics and to obtain an improved simple formula for easy estimation of the dominant first torsional mode, which is important in the dynamic integration of the engine and airframe system. Then, a linearized version of the model was used to investigate the effects of rotor speed variations and rotor state feedback on helicopter flight dynamics. Results show that, by including rotor speed variations, the effective vertical damping decreases significantly from that calculated with a constant speed assumption, thereby providing a better correlation with flight test data. Higher closed-loop bandwidths appear to be more readily achievable with rotor state feedback. The results also indicate that both aircraft and rotor flapping responses to gust disturbance are significantly attenuated when rotor state feedback is used.

Chen, Robert T. N.↗
NASA NTRS1993

Numerical Simulation of a High Mach Number Jet Flow

The recent efforts to develop accurate numerical schemes for transition and turbulent flows are motivated, among other factors, by the need for accurate prediction of flow noise. The success of developing high speed civil transport plane (HSCT) is contingent upon our understanding and suppression of the jet exhaust noise. The radiated sound can be directly obtained by solving the full (time-dependent) compressible Navier-Stokes equations. However, this requires computational storage that is beyond currently available machines. This difficulty can be overcome by limiting the solution domain to the near field where the jet is nonlinear and then use acoustic analogy (e.g., Lighthill) to relate the far-field noise to the near-field sources. The later requires obtaining the time-dependent flow field. The other difficulty in aeroacoustics computations is that at high Reynolds numbers the turbulent flow has a large range of scales. Direct numerical simulations (DNS) cannot obtain all the scales of motion at high Reynolds number of technological interest. However, it is believed that the large scale structure is more efficient than the small-scale structure in radiating noise. Thus, one can model the small scales and calculate the acoustically active scales. The large scale structure in the noise-producing initial region of the jet can be viewed as a wavelike nature, the net radiated sound is the net cancellation after integration over space. As such, aeroacoustics computations are highly sensitive to errors in computing the sound sources. It is therefore essential to use a high-order numerical scheme to predict the flow field. The present paper presents the first step in a ongoing effort to predict jet noise. The emphasis here is in accurate prediction of the unsteady flow field. We solve the full time-dependent Navier-Stokes equations by a high order finite difference method. Time accurate spatial simulations of both plane and axisymmetric jet are presented. Jet Mach numbers of 1.5 and 2.1 are considered. Reynolds number in the simulations was about a million. Our numerical model is based on the 2-4 scheme by Gottlieb & Turkel. Bayliss et al. applied the 2-4 scheme in boundary layer computations. This scheme was also used by Ragab and Sheen to study the nonlinear development of supersonic instability waves in a mixing layer. In this study, we present two dimensional direct simulation results for both plane and axisymmetric jets. These results are compared with linear theory predictions. These computations were made for near nozzle exit region and velocity in spanwise/azimuthal direction was assumed to be zero.

Hayder, M. Ehtesham↗
NASA NTRS1991

Controls for space structures

Assembly and operation of large space structures (LSS) in orbit will require robot-assisted docking and berthing of partially-assembled structures. These operations require new solutions to the problems of controls. This is true because of large transient and persistent disturbances, controller-structure interaction with unmodeled modes, poorly known structure parameters, slow actuator/sensor dynamical behavior, and excitation of nonlinear structure vibrations during control and assembly. For on-orbit assembly, controllers must start with finite element models of LSS and adapt on line to the best operating points, without compromising stability. This is not easy to do, since there are often unmodeled dynamic interactions between the controller and the structure. The indirect adaptive controllers are based on parameter estimation. Due to the large number of modes in LSS, this approach leads to very high-order control schemes with consequent poor stability and performance. In contrast, direct model reference adaptive controllers operate to force the LSS to track the desirable behavior of a chosen model. These schemes produce simple control algorithms which are easy to implement on line. One problem with their use for LSS has been that the model must be the same dimension as the LSS - i.e., quite large. A control theory based on the command generator tracker (CGT) ideas of Sobel, Mabins, Kaufman and Wen, Balas to obtain very low-order models based on adaptive algorithms was developed. Closed-loop stability for both finite element models and distributed parameter models of LSS was proved. In addition, successful numerical simulations on several LSS databases were obtained. An adaptive controller based on our theory was also implemented on a flexible robotic manipulator at Martin Marietta Astronautics. Computation schemes for controller-structure interaction with unmodeled modes, the residual mode filters or RMF, were developed. The RMF theory was modified to compensate slow actuator/sensor dynamics. These new ideas are being applied to LSS simulations to demonstrate the ease with which one can incorporate slow actuator/sensor effects into our design. It was also shown that residual mode filter compensation can be modified for small nonlinearities to produce exponentially stable closed-loop control. A theory for disturbance accommodating controllers based on reduced order models of structures was developed, and stability results for these controllers in closed-loop with large-scale finite element models of structures were obtained.

Balas, Mark↗
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 NTRS1992

Interference patterns in the Spacelab 2 plasma wave data - Oblique electrostatic waves generated by the electron beam

During the Spacelab 2 mission the University of Iowa's Plasma Diagnostics Package (PDP) explored the plasma environment around the shuttle. Wideband spectrograms of plasma waves were obtained from the PDP at frequencies of 0-30 kHz and at distances up to 400 m from the shuttle. Strong low-frequency (below 10 kHz) electric field noise was observed in the wideband data during two periods in which an electron beam was ejected from the shuttle. This noise shows clear evidence of interference patterns caused by the finite (3.89 m) antenna length. The low-frequency noise was the most dominant type of noise produced by the ejected electron beam. Analysis of antenna interference patterns generated by these waves permits a determination of the wavelength, the direction of propagation, and the location of the source region. The observed waves have a linear dispersion relation very similar to that of ion acoustic waves. The waves are believed to be oblique ion acoustic or high-order ion cyclotron waves generated by a current of ambient electrons returning to the shuttle in response to the ejected electron beam.

Feng, Wei↗
NASA NTRS1993

Gravitational tides in the outer planets. I - Implications of classical tidal theory. II - Interior calculations and estimation of the tidal dissipation factor

Classical tidal theory is applied to the atmospheres of the outer planets. The tidal geopotential due to satellites of the outer planets is discussed, and the solution of Laplace's tidal equation for Hough modes appropriate to tides on the outer planets is examined. The vertical structure of tidal modes is described, noting that only relatively high-order meridional mode numbers can propagate vertically with growing amplitude. Expected magnitudes for tides in the visible atmosphere of Jupiter are discussed. The classical theory is extended to planetary interiors taking the effects of spherically and self-gravity into account. The thermodynamic structure of Jupiter is described and the WKB theory of the vertical structure equation is presented. The regions for which inertial, gravity, and acoustic oscillations are possible are delineated. The case of a planet with a neutral interior is treated, discussing the various atmospheric boundary conditions and showing that the tidal response is small.

Ioannou, Petros J.↗
NASA NTRS1993

Effects of spatial order of accuracy on the computation of vortical flowfields

The effect of the order-of-accuracy, used for the spatial discretization, on the resolution of the leading edge vortices over sharp-edged delta wings is investigated. The flowfield is computed using a viscous/inviscid zonal approach. The viscous flow in the vicinity of the wing is computed using the conservative formulation of the compressible, thin-layer Navier-Stokes equations. The leeward-side vortical flowfield and the other flow regions away from the surface are computed as inviscid. The time integration is performed with both an explicit fourth-order Runge-Kutta scheme and an implicit, factorized, iterative scheme. High-order-accurate inviscid fluxes are computed using both a conservative and a non-conservative (primitive variable) formulation. The nonlinear, inviscid terms of the primitive variable form of the governing equations are evaluated with a finite-difference numerical scheme based on the sign of the eigenvalues. High-order, upwind-biased, finite difference formulas are used to evaluate the derivatives of the nonlinear convective terms. Computed results are compared with available experimental data, and comparisons of the flowfield in the vicinity of the vortex cores are presented.

Ekaterinaris, J. A.↗
NASA NTRS1988

Maximum likelihood identification for large space structures

This paper examines the use of on-orbit identification based on Maximum Likelihood Estimation (MLE) to provide these high-order, high-accuracy control design models for large space structures (LSS's). First, it outlines a general MLE identification algorithm, together with a covariance-analysis procedure to assess algorithm performance in terms of systematic and stochastic errors. Next, it examines various simplifications appropriate for the LSS identification application. Simplified analytical performance results are presented, as are numerical results to support these analyses. Finally, a graphical interpretation of these results is given.

Barrett, Michael F.↗
NASA NTRS1993

Summation by parts, projections, and stability

We have derived stability results for high-order finite difference approximations of mixed hyperbolic-parabolic initial-boundary value problems (IBVP). The results are obtained using summation by parts and a new way of representing general linear boundary conditions as an orthogonal projection. By slightly rearranging the analytic equations, we can prove strict stability for hyperbolic-parabolic IBVP. Furthermore, we generalize our technique so as to yield strict stability on curvilinear non-smooth domains in two space dimensions. Finally, we show how to incorporate inhomogeneous boundary data while retaining strict stability. Using the same procedure one can prove strict stability in higher dimensions as well.

Olsson, Pelle↗
NASA NTRS1993

A component modes projection and assembly model reduction methodology for articulated, multi-flexible body structures

A two-stage model reduction methodology, combining the classical Component Mode Synthesis (CMS) method and the newly developed Enhanced Projection and Assembly (EP&A) method, is proposed in this research. The first stage of this methodology, called the COmponent Modes Projection and Assembly model REduction (COMPARE) method, involves the generation of CMS mode sets, such as the MacNeal-Rubin mode sets. These mode sets are then used to reduce the order of each component model in the Rayleigh-Ritz sense. The resultant component models are then combined to generate reduced-order system models at various system configurations. A composite mode set which retains important system modes at all system configurations is then selected from these reduced-order system models. In the second stage, the EP&A model reduction method is employed to reduce further the order of the system model generated in the first stage. The effectiveness of the COMPARE methodology has been successfully demonstrated on a high-order, finite-element model of the cruise-configured Galileo spacecraft.

Lee, Allan Y.↗
NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗
NASA NTRS1993

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Lufkin, Eric A.↗
NASA NTRS1992

Numerical simulation of entry flow over blunt swept-wing planes

Viscous, heat-conducting flow with chemical and vibrational relaxation processes of the constituent gases surrounding winged spacecraft is considered in the continuum regime. The Navier-Stokes equations are appended by additional vibrational energy and species rate equations and supplemented by the equations of state and the phenomenological laws based on mixture rules or collisional cross sections. Numerical convective flux can be obtained from several forms of one-dimensional Riemann solver, with or without entropy correction. High-order accuracy is obtained from two types of reconstructive interpolation. A number of explicit and implicit numerical schemes have been implemented as a means to yield converged solutions. Both shock-fitting, finite-difference and shock-capturing, finite-volume techniques have been tested for configurations such as a sphere, double ellipsoid, blunt-edge delta wing, a European Hermes vehicle, and the U.S. Shuttle Orbiter. The shock-fitting code provides excellent results only for simple configurations, whereas the shock-capturing code leads to overall satisfying solutions for complex geometries.

Li, C. P.↗
NASA NTRS1993

Component modes damping assignment methodology for articulated, multiflexible body structures

To simulate the dynamical motion of articulated, multiflexible body structures, one can use multibody simulation packages such as DISCOS. To this end, one must supply appropriate reduced-order models for all of the flexible components involved. The component modes projection and assembly model reduction (COMPARE) methodology is one way to construct these reduced-order component models, which when reassembled capture important system input-to-output mapping of the full-order model at multiple system configurations of interest. In conjunction, we must also supply component damping matrices which when reassembled generate a system damping matrix that has certain desirable properties. The problem of determining the damping factors of components' modes to achieve a given system damping matrix is addressed here. To this end, we must establish from first principles a matrix-algebraic relation between the system's modal damping matrix and the components' modal damping matrices. An unconstrained/constrained optimization problem can then be formulated to determine the component modes' damping factors that best satisfy that matrix-algebraic relation. The effectiveness of the developed methodology, called ModeDamp, has been successfully demonstrated on a high-order, finite element model of the Galileo spacecraft.

Lee, Allan Y.↗
NASA NTRS1992

An exploratory investigation of the flight dynamics effects of rotor rpm variations and rotor state feedback in hover

This paper presents the results of an analytical study conducted to investigate airframe/engine interface dynamics, and the influence of rotor speed variations on the flight dynamics of the helicopter in hover, and to explore the potential benefits of using rotor states as additional feedback signals in the flight control system. The analytical investigation required the development of a parametric high-order helicopter hover model, which included heave/yaw body motion, the rotor speed degree of freedom, rotor blade motion in flapping and lead-lag, inflow dynamics, a drive train model with a flexible rotor shaft, and an engine/rpm governor. First, the model was used to gain insight into the engine/drive train/rotor system dynamics and to obtain an improved simple formula for easy estimation of the dominant first torsional mode, which is important in the dynamic integration of the engine and airframe system. Then, a linearized version of the model was used to investigate the effects of rotor speed variations and rotor state feedback on helicopter flight dynamics. Results show that, by including rotor speed variations, the effective vertical damping decreases significantly from that calculated with a constant speed assumption, thereby providing a better correlation with flight test data. Higher closed-loop bandwidths appear to be more readily achievable with rotor state feedback. The results also indicate that both aircraft and rotor flapping responses to gust disturbance are significantly attenuated when rotor state feedback is used.

Chen, Robert T. N.↗