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 127 records · Page 7

Multidimensional difference schemes with fourth-order accuracy

An explicit finite-difference algorithm is presented for the solution of quasilinear divergence free multidimensional hyperbolic systems. The method consists of four steps per time level. The resulting scheme is fourth-order accurate in both space and time, though the intermediate steps are only first-order accurate. The family of schemes introduced is dissipative, and hence, suitable for both smooth flows and flows containing shocks. This method is compared, in several numerical examples, with both second-order schemes and others that are fourth order in space, but second order in time.

Turkel, E.↗

An algorithm to compute the sequency ordered Walsh transform

A fast sequency-ordered Walsh transform algorithm is presented; this sequency-ordered fast transform is complementary to the sequency-ordered fast Walsh transform introduced by Manz (1972) and eliminating gray code reordering through a modification of the basic fast Hadamard transform structure. The new algorithm retains the advantages of its complement (it is in place and is its own inverse), while differing in having a decimation-in time structure, accepting data in normal order, and returning the coefficients in bit-reversed sequency order. Applications include estimation of Walsh power spectra for a random process, sequency filtering and computing logical autocorrelations, and selective bit reversing.

Larsen, H.↗

A fourth-order box method for solving the boundary layer equations

A fourth order box method for calculating high accuracy numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations is presented. The method is the natural extension of the second order Keller Box scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary layer equations. Numerical results for high accuracy test cases show the method to be significantly faster than other higher order and second order methods.

Wornom, S. F.↗

A critical study of higher-order numerical methods for solving the boundary-layer equations

A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗

Principles of equivalence, Eotvos experiments, and gravitational red-shift experiments - The free fall of electromagnetic systems to post-post-Coulombian order

Free fall in an external static gravitational field of a composite test body consisting of electromagnetically interacting charged particles is analyzed using the TH-epsilon-mu formalism developed by Lightman and Lee (1973). The Lightman and Lee calculation that was carried out to 'post-Coulombian' order is extended to 'post-post-Coulombian' order, and violations of the weak equivalence principle are found which result from the magnetostatic self-energy of the test body. The current generation of Eotvos experiments is shown to test the cited formalism at first order in the external gravitational potential to an accuracy of 4 by 10 to the -10th power for electrostatic energies and 6 by 10 to the -6th power for magnetostatic energies. It is concluded that although the formalism can be tested to an accuracy of 0.04 for electrostatic energies at second order in the gravitational potential, a significant second-order test for magnetostatic energies would require Eotvos experiments accurate to better than one part in 10 to the 15th power.

Haugan, M. P.↗

A second-order closure analysis of turbulent diffusion flames

A complete second-order closure computer program for the investigation of compressible, turbulent, reacting shear layers was developed. The equations for the means and the second order correlations were derived from the time-averaged Navier-Stokes equations and contain third order and higher order correlations, which have to be modeled in terms of the lower-order correlations to close the system of equations. In addition to fluid mechanical turbulence models and parameters used in previous studies of a variety of incompressible and compressible shear flows, a number of additional scalar correlations were modeled for chemically reacting flows, and a typical eddy model developed for the joint probability density function for all the scalars. The program which is capable of handling multi-species, multistep chemical reactions, was used to calculate nonreacting and reacting flows in a hydrogen-air diffusion flame.

Varma, A. K.↗

Application of higher-order numerical methods to the boundary-layer equations

A fourth-order method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations for both attached and separated flows. The efficiency of the present method is compared with other higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, the three-point spline methods, and a modified finite-element method. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗

Intracrystalline cation order in a lunar crustal troctolite

Lunar sample 76535 appears to be one of the most slowly cooled bits of silicate material yet studied. It provides, therefore, a unique opportunity for the study of ordering processes in the minerals present. A better understanding of these processes may permit better characterization of the thermal history of this and similar rocks. The cation ordering in the olivine is consistent with terrestrial olivines favoring the interpretation that ordering in olivines increases with increasing temperature. In low bronzite, the deviations from the common orthopyroxene space group appear to be caused by cation order on the basis of four M sites instead of two. The degree of cation order in each of these minerals is consistent with the rock having been excavated from its depth of formation by tectonic or impact processes without being reheated above 300 C.

Smyth, J. R.↗

Critical study of higher order numerical methods for solving the boundary-layer equations

A fourth order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method, which is the natural extension of the second order box scheme to fourth order, was demonstrated with application to the incompressible, laminar and turbulent, boundary layer equations. The efficiency of the present method is compared with two point and three point higher order methods, namely, the Keller box scheme with Richardson extrapolation, the method of deferred corrections, a three point spline method, and a modified finite element method. For equivalent accuracy, numerical results show the present method to be more efficient than higher order methods for both laminar and turbulent flows.

Wornom, S. F.↗

A study of the applicability of nucleation theory to quasi-thermodynamic transitions of second and higher Ehrenfest-order

The applicability of classical nucleation theory to second (and higher) order thermodynamic transitions in the Ehrenfest sense has been investigated and expressions have been derived upon which the qualitative and quantitative success of the basic approach must ultimately depend. The expressions describe the effect of temperature undercooling, hydrostatic pressure, and tensile stress upon the critical parameters, the critical nucleus size, and critical free energy barrier, for nucleation in a thermodynamic transition of any general order. These expressions are then specialized for the case of first and second order transitions. The expressions for the case of undercooling are then used in conjunction with literature data to estimate values for the critical quantities in a system undergoing a pseudo-second order transition (the glass transition in polystyrene). Methods of estimating the interfacial energy gamma in systems undergoing a first and second order transition are also discussed.

Barker, R. E., Jr.↗

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.↗

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↗

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↗

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↗

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.↗

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.↗

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.↗

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↗