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 397 records · Page 22

Second-order solution for determining density and velocity fields of galaxies

In this Letter, I use second-order Lagrangian perturbation theory to calculate an analytical expression relating density to velocity in a gravitating system. This solution can be used to compare peculiar velocity field measurements with observations of large-scale structure. The predictions of both linear theory and second-order theory are compared with the results of N-body simulations. While linear theory systematically overestimates the velocity flows in high-density regions, the second-order corrections calculated herein remove this systematic error.

Gramann, Mirt↗

Reliable design of H-2 optimal reduced-order controllers via a homotopy algorithm

Due to control processor limitations, the design of reduced-order controllers is an active area of research. Suboptimal methods based on truncating the order of the corresponding linear-quadratic-Gaussian (LQG) compensator tend to fail if the requested controller dimension is sufficiently small and/or the requested controller authority is sufficiently high. Also, traditional parameter optimization approaches have only local convergence properties. This paper discusses a homotopy algorithm for optimal reduced-order control that has global convergence properties. The exposition is for discrete-time systems. The algorithm has been implemented in MATLAB and is applied to a benchmark problem.

Collins, Emmanuel G.↗

Solution of the linear shallow water equations by the fourth-order leapfrog scheme

Numerical schemes of the first and second order of approximation introduce numerical distortion when the wave propagation over a long distance is investigated. To alleviate this problem, the fourth-order leapfrog scheme is constructed. The standard leapfrog method is based on the truncated Taylor series expansion which depicts an error proportional to the second-order terms. In the proposed method the numerical solution is corrected for these terms. The space and time corrections work well in diminishing numerical dispersion and dissipation.

Kowalik, Z.↗

Regenerator performance in a Vuilleumier refrigerator compared with a third-order numerical model

A three-stage Vuilleumier refrigerator was used to measure the performance of various third stage regenerators. The refrigerator operates between 2.5 and 5.0 Hz and, depending on the material used in the third stage regenerator, achieves temperatures of 8 to 20 K at the cold end of the third stage. This paper presents a comparison of regenerator performance for four regenerator materials: 229 micron diameter spheres of Pb(+)5 pct Sb, 229 micron diameter spheres of brass, 216 micron irregularly-shaped GdRh powder, and a mixture of 229 microns and 762 microns diameter spheres of Pb(+5) pct Sb. The experimental results are compared with a first-order model that neglects the void volume within the regenerator and with a third-order model that considers the effect of pressure oscillations in the regenerator void volume. Experimental results indicate that regenerator losses are dominated by the pressure oscillation in the void volume rather than the mass flow through the temperature gradient in the regenerator. These results are consistent with the third-order numerical model. This model shows that the heat capacity of the gas in the void space as well as the heat capacity of the matrix influences the regenerator performance.

Bradley, P. E.↗

A Study of Second-Order Supersonic Flow Theory

Second-order solutions of supersonic-flow problems are sought by iteration, using the linearized solution as the first step. For plane and axially symmetric flows, particular solutions of the iteration equation are discovered which reduce the second-order problem to an equivalent linearized problem. Comparison of second-order solutions with exact and numerical results shows great improvement over linearized theory. For full three-dimensional flow, only a partial particular solution is found. The inclined cone is solved, and the possibility of treating more general problems is considered.

Van Dyke, Milton D↗

A wind-tunnel investigation of the effects of thrust-axis inclination on propeller first-order vibration

Data on the aerodynamic excitation of first-order vibration occurring in a representative three-blade propeller having its thrust axis inclined to the air stream at angles of 0 degrees, 4.55 degrees, and 9.8 degrees are included in this paper. For several representative conditions the aerodynamic excitation has been computed and compared with the measured values. Blade stresses also were measured to permit the evaluation of the blade stress resulting from a given blade aerodynamic excitation. It was concluded that the section aerodynamic exciting force of a pitched propeller may be computed accurately at low rotational speeds. As section velocities approach the speed of sound, the accuracy of computation of section aerodynamic exciting force is not always so satisfactory. First-order blade vibratory stresses were computed with satisfactory accuracy from untilted-propeller loading data. A stress prediction which assumes a linear relation between first-order vibratory stress and the product of pitch angle and dynamic pressure and which is based on stresses at low rotational speeds will be conservative when the outer portions of the blade are in the transonic and low supersonic speed range.

Gray, W H↗

Second-order design sensitivities for linear elastic problems

The subject of this paper is the efficient and accurate determination of second-order design sensitivities in elastic bodies. The approach being carried out here is the direct differentiation of the governing derivative boundary element method (DBEM) formulation of the problem. Second-order sensitivities of boundary stresses are obtained in an elegant manner. A numerical implementation of the method is carried out with isoparametric quadratic boundary elements and numerical results are presented for several sample problems. Considerable savings in computational effort for an optimization procedure is possible through the use of efficiently determined accurate second-order design sensitivities.

Zhang, Qing↗

High order hybrid numerical simulations of two dimensional detonation waves

In order to study multi-dimensional unstable detonation waves, a high order numerical scheme suitable for calculating the detailed transverse wave structures of multidimensional detonation waves was developed. The numerical algorithm uses a multi-domain approach so different numerical techniques can be applied for different components of detonation waves. The detonation waves are assumed to undergo an irreversible, unimolecular reaction A yields B. Several cases of unstable two dimensional detonation waves are simulated and detailed transverse wave interactions are documented. The numerical results show the importance of resolving the detonation front without excessive numerical viscosity in order to obtain the correct cellular patterns.

Cai, Wei↗

Failure modes of reduced-order orbit determination filters and their remedies

Ways in which failure can occur in reduced-order, orbit determination filter, error covariance calculations are discussed. In the context of this article, reduced-order filters denote nonoptimal filters which include fixed levels of uncertainty in some parameters of the measurement models or in the spacecraft dynamical model which are not explicitly estimated in the filter equations. Failure is defined as an increase in the orbit determination covariance with the addition of data or as an unreasonable growth in the covariance with time, i.e., nonasymptotic behavior of the covariance. Some simple, known cases of failure are discussed along with their traditional remedies. In addition, more modern remedies are discussed which are currently under development at the Jet Propulsion Laboratory. The article first describes the known problems of reduced-order filters when they are employed for orbit determination, and their traditional remedies. Then, having defined these, the relevancy and desirability of the more modern remedies are made apparent.

Scheeres, D. J.↗

Cation ordering in orthopyroxenes and cooling rates of meteorites: Low temperature cooling rates of Estherville, Bondoc and Shaw

The cooling rates of meteorites provide important constraints on the size of their parent bodies, and their accretionary and evolutionary histories. However, the cooling rates obtained so far from the commonly used metallographic, radiometric and fission-track methods have been sometimes quite controversial, such as in the case of the mesosiderites and the meteorite Shaw. We have undertaken a systematic study of the cooling rates of meteorites using a different approach, which involves single crystal x-ray determination of Fe(2+)-Mg ordering in orthopyroxenes (OP(x)) in meteorites, subject to bulk compositional constraints, and numerical simulation of the evolution of the ordering state as a function of cooling rate, within the framework of the thermodynamic and kinetic principles governing cation ordering. We report the results obtained for OP(x) crystals from Shaw and two mesosiderites, Estherville and Bondoc.

Ganguly, J.↗

GMM-1: A 50 degree and order gravitational field model for Mars

Knowledge of the gravitational field, in combination with surface topography, provides one of the principal means of inferring the internal structure of a planetary body. The highest resolution gravitational field for Mars published thus far was derived from Doppler tracking data from the Mariner 9 and Viking 1 and 2 spacecraft and is complete to degree and order 18 corresponding to a half wavelength resolution of approximately 600 km. This field, which is characterized by a spatial resolution that is slightly better than that of the highest resolution (16x16) topographic model, has been utilized extensively in analyses of the state of stress and isostatic compensation of the Martian lithosphere. However, the resolution and quality of current gravity and topographic fields are such that the origin and evolution of even the major physiographic features on Mars, such as the hemispheric dichotomy and Tharsis rise, are not well understood. We have re-analyzed the Viking and Mariner data sets and have derived a new gravitational field, which we designated GMM-1 (Goddard Mars Model-1). This model is complete to spherical harmonic degree and order 50 with a corresponding (half wavelength) spatial resolution of 200-300 km where the data permit. In contrast to previous models, GMM-1 was solved to as high degree and order as necessary to nearly exhaust the attenuated gravitational signal contained in the tracking data.

Smith, D. E.↗

Fourth-order 2N-storage Runge-Kutta schemes

A family of five-stage fourth-order Runge-Kutta schemes is derived; these schemes required only two storage locations. A particular scheme is identified that has desirable efficiency characteristics for hyperbolic and parabolic initial (boundary) value problems. This scheme is competitive with the classical fourth-order method (high-storage) and is considerably more efficient and accurate than existing third-order low-storage schemes.

Carpenter, Mark H.↗

Third-order 2N-storage Runge-Kutta schemes with error control

A family of four-stage third-order explicit Runge-Kutta schemes is derived that requires only two storage locations and has desirable stability characteristics. Error control is achieved by embedding a second-order scheme within the four-stage procedure. Certain schemes are identified that are as efficient and accurate as conventional embedded schemes of comparable order and require fewer storage locations.

Carpenter, Mark H.↗

Pseudospectral collocation methods for fourth order differential equations

Collocation schemes are presented for solving linear fourth order differential equations in one and two dimensions. The variational formulation of the model fourth order problem is discretized by approximating the integrals by a Gaussian quadrature rule generalized to include the values of the derivative of the integrand at the boundary points. Collocation schemes are derived which are equivalent to this discrete variational problem. An efficient preconditioner based on a low-order finite difference approximation to the same differential operator is presented. The corresponding multidomain problem is also considered and interface conditions are derived. Pseudospectral approximations which are C1 continuous at the interfaces are used in each subdomain to approximate the solution. The approximations are also shown to be C3 continuous at the interfaces asymptotically. A complete analysis of the collocation scheme for the multidomain problem is provided. The extension of the method to the biharmonic equation in two dimensions is discussed and results are presented for a problem defined in a nonrectangular domain.

Malek, Alaeddin↗

Development of homotopy algorithms for fixed-order mixed H2/H(infinity) controller synthesis

A major difficulty associated with H-infinity and mu-synthesis methods is the order of the resulting compensator. Whereas model and/or controller reduction techniques are sometimes applied, performance and robustness properties are not preserved. By directly constraining compensator order during the optimization process, these properties are better preserved, albeit at the expense of computational complexity. This paper presents a novel homotopy algorithm to synthesize fixed-order mixed H2/H-infinity compensators. Numerical results are presented for a four-disk flexible structure to evaluate the efficiency of the algorithm.

Whorton, M.↗

Probing the antiferromagnetic long-range order with Glauber spin states

It is well known that the ground state of low-dimensional antiferromagnets deviates from Neel states due to strong quantum fluctuations. Even in the presence of long-range order, those fluctuations produce a substantial reduction of the magnetic moment from its saturation value. Numerical simulations in anisotropic antiferromagnetic chains suggest that quantum fluctuations over Neel order appear in the form of localized reversal of pairs of neighboring spins. In this paper, we propose a coherent state representation for the ground state to describe the above situation. In the one-dimensional case, our wave function corresponds to a two-mode Glauber state, when the Neel state is used as a reference, while the boson fields are associated to coherent flip of spin pairs. The coherence manifests itself through the antiferromagnetic long-range order that survives the action of quantum fluctuations. The present representation is different from the standard zero-point spin wave state, and is asymptotically exact in the limit of strong anisotropy. The fermionic version of the theory, obtained through the Jordan-Wigner transformation, is also investigated.

Cabrera, Guillermo G.↗

Extensional flow convecting a reactant undergoing a first order homogeneous reaction and diffusional mass transfer from a sphere at low to intermediate Peclet and Damkohler numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Shah, N. Y.↗

High-order essentially non-oscillatory methods for computational aeroacoustics

The desire to obtain acoustic information from the numerical solution of a nonlinear system of equations is a demanding proposition for a computational algorithm. High-order accuracy is required for the propagation of high-frequency, low-amplitude waves. In addition, it is desirable to highly resolve discontinuities that can develop in the solutions of the Euler or Navier-Stokes equations. The class of essentially non-oscillatory (ENO) shock-capturing schemes has been designed to have both of these properties. The dual capacity of ENO schemes for high-order accuracy and non-oscillatory shock-capturing is achieved through the use of adaptive stenciling, which makes these schemes highly nonlinear. These schemes are briefly described and referenced herein. A fourth-order algorithm is then applied to the solution of an acoustic wave in a quasi-one-dimensional converging-diverging nozzle.

Casper, Jay↗