Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 NTRS1986

A fourth order accurate finite difference scheme for the computation of elastic waves

A finite difference for elastic waves is introduced. The model is based on the first order system of equations for the velocities and stresses. The differencing is fourth order accurate on the spatial derivatives and second order accurate in time. The model is tested on a series of examples including the Lamb problem, scattering from plane interf aces and scattering from a fluid-elastic interface. The scheme is shown to be effective for these problems. The accuracy and stability is insensitive to the Poisson ratio. For the class of problems considered here it is found that the fourth order scheme requires for two-thirds to one-half the resolution of a typical second order scheme to give comparable accuracy.

Bayliss, A.↗
NASA NTRS1987

Order reduction in linear state estimation under performance constraints

The design and analysis of minimal-order state estimators for possibly time-varying linear systems, under constraints on the maximal allowable mean-square error, are considered. A global lower bound on the optimal error is derived, along with a lower bound on the minimal estimator order, needed for meeting the performance constraint. The ideal reduced-order estimator which satisfies the lower bound is derived, along with conditions for its realizability. When the ideal estimator is not realizable, its structure forms a suboptimal estimator, which maintains, in some sense, a local optimality property and is called the pseudoideal estimator. The mean-square error of the pseudoideal estimator defines upper bounds on the optimal error and on the estimator order needed for meeting the performance constraint. The lower and the upper bounds on the order define a reduced search set for the design problem. When the distance between the ideal and the pseudoideal estimators is sufficiently small in a certain numerical sense, the pseudoideal estimator may be considered optimal for practical purposes.

Baram, Yoram↗
NASA NTRS1989

On the state estimation of structures with second order observers

The use of Linear Quadratic Regulator (LQR) control synthesis techniques implies the availability of full state feedback. For vibration control of structures, usually only a limited number of states are measured from which an observer model reconstructs the full state. It is shown that using second order observers is a viable technique for reconstructing the unmeasured states of structures under mildly restrictive conditions. Moreover, the computational advantages of the second order observer as compared to a first order observer indicate that significantly larger observer models may be utilized. Numerical examples are used to demonstrate the performance of second order observers. The implications of second order observers in the development of Controls-Structures Interaction (CSI) technology is discussed.

Belvin, W. Keith↗
NASA NTRS1989

On the state estimation of structures with second order observers

The use of the linear quadratic regulator control synthesis techniques implies the availability of full state feedback. For vibration control of structures, usually only a limited number of states are measured from which an observer model reconstructs the full state. This paper shows that using second-order observers is a viable technique for reconstructing the unmeasured states of structures under mildly restrictive conditions. Moreover, the computational advantages of the second-order observer, as compared to a first-order observer, indicate that significantly larger observer models may be utilized. Numerical examples are used to demonstrate the performance of second-order observers. The implications of second-order observers in the development of the control/structures interaction technology is discussed.

Belvin, W. Keith↗
NASA NTRS1989

First and second order convex approximation strategies in structural optimization

In this paper, various methods based on convex approximation schemes are discussed that have demonstrated strong potential for efficient solution of structural optimization problems. First, the convex linearization method (Conlin) is briefly described, as well as one of its recent generalizations, the method of moving asymptotes (MMA). Both Conlin and MMA can be interpreted as first-order convex approximation methods that attempt to estimate the curvature of the problem functions on the basis of semiempirical rules. Attention is next directed toward methods that use diagonal second derivatives in order to provide a sound basis for building up high-quality explicit approximations of the behavior constraints. In particular, it is shown how second-order information can be effectively used without demanding a prohibitive computational cost. Various first-order and second-order approaches are compared by applying them to simple problems that have a closed form solution.

Fleury, C.↗
NASA NTRS1989

The role of order in distributed programs

The role of order in building distributed systems is discussed. It is the belief that a principle of event ordering underlies the wide range of operating systems mechanisms that were put forward for building robust distributed software. Stated concisely, this principle achieves correct distributed behavior by ordering classes of distributed events that conflict with one another. By focusing on order, simplified descriptions can be obtained and convincingly correct solutions to problems that might otherwise have looked extremely complex. Moreover, it is observed that there are a limited number of ways to obtain order, and that the choice made impacts greatly on performance.

Birman, Kenneth P.↗
NASA NTRS1989

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Chun, H. M.↗
NASA NTRS1990

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗
NASA NTRS1990

Discrete integration of continuous Kalman filtering equations for time invariant second-order structural systems

A general form for the first-order representation of the continuous second-order linear structural-dynamics equations is introduced to derive a corresponding form of first-order continuous Kalman filtering equations. Time integration of the resulting equations is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete Kalman filtering equations involving only symmetric sparse N x N solution matrices.

Park, K. C.↗
NASA NTRS1991

Finite-volume application of high-order ENO schemes to two-dimensional boundary-value problems

Finite-volume applications of high-order accurate ENO schemes to two-dimensional boundary-value problems are studied. These schemes achieve high-order spatial accuracy, in smooth regions, by a piecewise polynomial approximation of the solution from cell averages. In addition, this spatial operation involves an adaptive stencil algorithm in order to avoid the oscillatory behavior that is associated with interpolation across steep gradients. High-order TVD Runge-Kutta methods are employed for time integration, thus making these schemes best suited for unsteady problems. Fifth- and sixth-order accurate applications are validated through a grid refinement study involving the solutions of scalar hyperbolic equations. A previously proposed extension for the Euler equations of gas dynamics is tested, including its application to solutions of boundary-value problems involving solid walls and curvilinear coordinates.

Casper, Jay↗
NASA NTRS1991

Ordered fast Fourier transforms on a massively parallel hypercube multiprocessor

The present evaluation of alternative, massively parallel hypercube processor-applicable designs for ordered radix-2 decimation-in-frequency FFT algorithms gives attention to the reduction of computation time-dominating communication. A combination of the order and computational phases of the FFT is accordingly employed, in conjunction with sequence-to-processor maps which reduce communication. Two orderings, 'standard' and 'cyclic', in which the order of the transform is the same as that of the input sequence, can be implemented with ease on the Connection Machine (where orderings are determined by geometries and priorities. A parallel method for trigonometric coefficient computation is presented which does not employ trigonometric functions or interprocessor communication.

Tong, Charles↗
NASA NTRS1993

Conservative high-order-accurate finite-difference methods for curvilinear grids

Two fourth-order-accurate finite-difference methods for numerically solving hyperbolic systems of conservation equations on smooth curvilinear grids are presented. The first method uses the differential form of the conservation equations; the second method uses the integral form of the conservation equations. Modifications to these schemes, which are required near boundaries to maintain overall high-order accuracy, are discussed. An analysis that demonstrates the stability of the modified schemes is also provided. Modifications to one of the schemes to make it total variation diminishing (TVD) are also discussed. Results that demonstrate the high-order accuracy of both schemes are included in the paper. In particular, a Ringleb-flow computation demonstrates the high-order accuracy and the stability of the boundary and near-boundary procedures. A second computation of supersonic flow over a cylinder demonstrates the shock-capturing capability of the TVD methodology. An important contribution of this paper is the dear demonstration that higher order accuracy leads to increased computational efficiency.

Rai, Man M.↗
NASA NTRS1993

A consistent model reduction of measured modal parameters for reduced-order active control

The problem of synthesizing reduced-order linear models of vibrating structures for the design of fixed-order dynamic feedback control is investigated. The present technique builds on a recently developed procedure for constructing an objective set of mass and stiffness matrices from measured modal parameters that are akin to the Craig-Bampton synthesized ones obtained from finite element models. The constructed mass and stiffness matrices are determined directly from the identification of experimental data, however, rather than through correlation or reconciliation of a finite element model. A model truncation criterion is then applied to the identified minimum-order mass and stiffness model to satisfy certain observability/controllability requirements for the reduced model. Numerical examples illustrate the effectiveness of the proposed technique for synthesizing reduced-order controllers from system realizations of experimental data. The dynamic performance of the resulting closed-loop models is assessed using the known full-order structural dynamics and compared with existing model reduction techniques.

Alvin, K. F.↗
NASA NTRS1993

Dynamic localization and second-order subgrid-scale models in large eddy simulations of channel flow

The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Cabot, William H.↗
NASA NTRS1994

Technical report series on global modeling and data assimilation. Volume 2: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere

High order horizontal diffusion of the form K Delta(exp 2m) is widely used in spectral models as a means of preventing energy accumulation at the shortest resolved scales. In the spectral context, an implicit formation of such diffusion is trivial to implement. The present note describes an efficient method of implementing implicit high order diffusion in global finite difference models. The method expresses the high order diffusion equation as a sequence of equations involving Delta(exp 2). The solution is obtained by combining fast Fourier transforms in longitude with a finite difference solver for the second order ordinary differential equation in latitude. The implicit diffusion routine is suitable for use in any finite difference global model that uses a regular latitude/longitude grid. The absence of a restriction on the timestep makes it particularly suitable for use in semi-Lagrangian models. The scale selectivity of the high order diffusion gives it an advantage over the uncentering method that has been used to control computational noise in two-time-level semi-Lagrangian models.

Max J. Suarez↗
NASA NTRS1994

Second-order reconstruction of the inflationary potential

To first order in the deviation from scale invariance the inflationary potential and its first two derivatives can be expressed in terms of the spectral indices of the scalar and tensor perturbations, n and n(sub T), and their contributions to the variance of the quadrupole CBR temperature anisotropy, S and T. In addition, there is a 'consistency relation' between these quantities: n(sub T) = (-1/ 7)(T/S). We derive the second-order expressions for the inflationary potential and its first two derivatives and the first-order expression for its third derivative, in terms, of n, n(sub T), S, T, and dn/d ln gamma. We also obtain the second-order consistency relation, n(sub T) = (-1/7)(T/S)(1 + 0.11(T/S) + 0.15(n-1)). As an example we consider the exponential potential, the only known case where exact analytic solutions for the perturbation spectra exist. We reconstruct the potential via Taylor expansion (with coefficients calculated at both first and second order), and introduce the Pade approximate as a greatly improved alternative.

Liddle, Andrew R.↗
NASA NTRS1995

A high performance totally ordered multicast protocol

This paper presents the Reliable Multicast Protocol (RMP). RMP provides a totally ordered, reliable, atomic multicast service on top of an unreliable multicast datagram service such as IP Multicasting. RMP is fully and symmetrically distributed so that no site bears un undue portion of the communication load. RMP provides a wide range of guarantees, from unreliable delivery to totally ordered delivery, to K-resilient, majority resilient, and totally resilient atomic delivery. These QoS guarantees are selectable on a per packet basis. RMP provides many communication options, including virtual synchrony, a publisher/subscriber model of message delivery, an implicit naming service, mutually exclusive handlers for messages, and mutually exclusive locks. It has commonly been held that a large performance penalty must be paid in order to implement total ordering -- RMP discounts this. On SparcStation 10's on a 1250 KB/sec Ethernet, RMP provides totally ordered packet delivery to one destination at 842 KB/sec throughput and with 3.1 ms packet latency. The performance stays roughly constant independent of the number of destinations. For two or more destinations on a LAN, RMP provides higher throughput than any protocol that does not use multicast or broadcast.

Montgomery, Todd↗
NASA NTRS1994

Higher order Godunov schemes for isothermal hydrodynamics

In this paper we construct higher order Godunov schemes for isothermal flow. Isothermal hydrodynamics serves as a good representation for several systems of astrophysical interest. The schemes designed here have second-order accuracy in space and time and some are third-order accurate for advection. Moreover, several ingredients of these schemes are essential components of even higher order. The methods designed here have excellent ability to represent smooth flow yet capture shocks with high resolution. Several test problems are presented. The algorithms presented here are compared with other algorithms having a comparable formal order of accuracy.

Balsara, Dinshaw S.↗