Search NASA⌕ Search

SEARCH · Search NASA

Results for “first 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 343 records · Page 19

A second-order theory of an artificial satellite under the influence of the oblateness of the earth

The general theory of the method of averaging is used to study the effect of the oblateness of the earth on the motion of an artificial satellite. The first-order harmonic J2 and the second-order harmonics J3 and J4 are included in the analysis. The first-order, periodic variations together with the second-order, averaged and secular variations are derived for all the six orbital elements. Alternative orbital elements and necessary relations are introduced to deal with circular orbits, orbits at critical inclination and equatorial orbits. Essential transformations and initialization procedures are also outlined.

Liu, J. J. F.↗

A natural conservative flux difference splitting for the hyperbolic systems of gasdynamics

Gasdynamic hyperbolic systems are treated by a novel conservative flux difference splitting upwind method which is applicable to explicit or implicit and iterative or direct schemes, for marching in time or space with the Euler or compressible Navier-Stokes equations. The method is able to capture sharp shocks correctly by maintaining global conservation. The embedded characteristics formulation is cast in the volumetric internal energy primitive variable, which is effective for real as well as perfect gases. Higher order upwind formulae are constructed from the distribution of pieces of the simple internal flux difference splitting, and the method reverts to first order at the appropriate points in order not to violate the domain of dependence by differencing across discontinuities. At insignificant computational overhead, switching is done by simple algebraic extensions of the truth functions that determine upwind direction for the first order scheme.

Lombard, C. K.↗

The dynamics and helium distribution in hydrogen-helium fluid planets

The simple case of a homogeneous planet without first-order phase transitions is considered and an investigation is conducted concerning a pure hydrogen planet in which a first-order phase transition takes place from fluid molecular hydrogen to fluid metallic hydrogen. Attention is also given to convection in the presence of a compositional gradient, the effects of helium insolubility in a cooling hydrogen-helium planet, a hydrogen-helium planet in its early evolution, and the case in which influence of phase transition occurs much later in the evolution of the planet.

Stevenson, D. J.↗

Lightweight Deployable Mirrors with Tensegrity Supports

The upper part of Figure 1 shows a small-scale prototype of a developmental class of lightweight, deployable structures that would support panels in precise alignments. In this case, the panel is hexagonal and supports disks that represent segments of a primary mirror of a large telescope. The lower part of Figure 1 shows a complete conceptual structure containing multiple hexagonal panels that hold mirror segments. The structures of this class are of the tensegrity type, which was invented five decades ago by artist Kenneth Snelson. A tensegrity structure consists of momentfree compression members (struts) and tension members (cables). The structures of this particular developmental class are intended primarily as means to erect large segmented primary mirrors of astronomical telescopes or large radio antennas in outer space. Other classes of tensegrity structures could also be designed for terrestrial use as towers, masts, and supports for general structural panels. An important product of the present development effort is the engineering practice of building a lightweight, deployable structure as an assembly of tensegrity modules like the one shown in Figure 2. This module comprises two octahedral tensegrity subunits that are mirror images of each other joined at their plane of mirror symmetry. In this case, the plane of mirror symmetry is both the upper plane of the lower subunit and the lower plane of the upper subunit, and is delineated by the midheight triangle in Figure 2. In the configuration assumed by the module to balance static forces under mild loading, the upper and lower planes of each sub-unit are rotated about 30 , relative to each other, about the long (vertical) axis of the structure. Larger structures can be assembled by joining multiple modules like this one at their sides or ends. When the module is compressed axially (vertically), the first-order effect is an increase in the rotation angle, but by virtue of the mirror arrangement, the net first-order rotation between the uppermost and lowermost planes is zero. The need to have zero net rotation between these planes under all loading conditions in a typical practical structure is what prompts the use of the mirror configuration. Force and moment loadings other than simple axial compression produce only second-order deformations through strains in the struts and cables.

Zeiders, Glenn W.↗

Least-squares finite element solution of 3D incompressible Navier-Stokes problems

Although significant progress has been made in the finite element solution of incompressible viscous flow problems. Development of more efficient methods is still needed before large-scale computation of 3D problems becomes feasible. This paper presents such a development. The most popular finite element method for the solution of incompressible Navier-Stokes equations is the classic Galerkin mixed method based on the velocity-pressure formulation. The mixed method requires the use of different elements to interpolate the velocity and the pressure in order to satisfy the Ladyzhenskaya-Babuska-Brezzi (LBB) condition for the existence of the solution. On the other hand, due to the lack of symmetry and positive definiteness of the linear equations arising from the mixed method, iterative methods for the solution of linear systems have been hard to come by. Therefore, direct Gaussian elimination has been considered the only viable method for solving the systems. But, for three-dimensional problems, the computer resources required by a direct method become prohibitively large. In order to overcome these difficulties, a least-squares finite element method (LSFEM) has been developed. This method is based on the first-order velocity-pressure-vorticity formulation. In this paper the LSFEM is extended for the solution of three-dimensional incompressible Navier-Stokes equations written in the following first-order quasi-linear velocity-pressure-vorticity formulation.

Jiang, Bo-Nan↗

Guidance law development for aeroassisted transfer vehicles using matched asymptotic expansions

This report addresses and clarifies a number of issues related to the Matched Asymptotic Expansion (MAE) analysis of skip trajectories, or any class of problems that give rise to inner layers that are not associated directly with satisfying boundary conditions. The procedure for matching inner and outer solutions, and using the composite solution to satisfy boundary conditions is developed and rigorously followed to obtain a set of algebraic equations for the problem of inclination change with minimum energy loss. A detailed evaluation of the zeroth order guidance algorithm for aeroassisted orbit transfer is performed. It is shown that by exploiting the structure of the MAE solution procedure, the original problem, which requires the solution of a set of 20 implicit algebraic equations, can be reduced to a problem of 6 implicit equations in 6 unknowns. A solution that is near optimal, requires a minimum of computation, and thus can be implemented in real time and on-board the vehicle, has been obtained. Guidance law implementation entails treating the current state as a new initial state and repetitively solving the zeroth order MAE problem to obtain the feedback controls. Finally, a general procedure is developed for constructing a MAE solution up to first order, of the Hamilton-Jacobi-Bellman equation based on the method of characteristics. The development is valid for a class of perturbation problems whose solution exhibits two-time-scale behavior. A regular expansion for problems of this type is shown to be inappropriate since it is not valid over a narrow range of the independent variable. That is, it is not uniformly valid. Of particular interest here is the manner in which matching and boundary conditions are enforced when the expansion is carried out to first order. Two cases are distinguished-one where the left boundary condition coincides with, or lies to the right of, the singular region, and another one where the left boundary condition lies to the left of the singular region. A simple example is used to illustrate the procedure where the obtained solution is uniformly valid to O(Epsilon(exp 2)). The potential application of this procedure to aeroassisted plane change is also described and partially evaluated.

Calise, Anthony J.↗

The rainout parameterization in a photochemical model

Photochemical models of tropospheric trace species typically simulate the effects of wet removal of soluble species from the atmosphere through a first-order rainout parameterization, which is independent of the species solubility and constant in time. In this work a modified first-order rainout parameterization is proposed that is a function of the species' solubility and takes into account the intermittent nature of atmospheric wet removal. In the case of a highly soluble species, such as HNO3, the inclusion of the effects of intermittent rainout in the parameterization leads to averaged rainout lifetimes considerably higher than ones used in most photochemical models. For illustrative purposes this parameterization is applied to a one-dimensional steady state model with parameterized eddy diffusion, although in principle it can be used in more sophisticated multidimensional models. The model results are compared to analagous calculations obtained with a 3-D tracer transport model.

Giorgi, F.↗

Cosmic rays at fluid discontinuities

Cosmic-ray transport near discontinuities in the background fluid velocity is considered. Matching conditions for the cosmic-ray distribution are derived for both shear and compressive (shock) discontinuities, keeping terms to second order in the ratio of fluid speed to energetic-particle speed. Acceleration is found at shear discontinuities, which is not present in the first-order theory, and a modification of the matching condition at shocks. If there is no particle source concentrated at a shock, the new condition reduces to that obtained from first-order theory. Monte Carlo simulations show good agreement with the theory.

Jokipii, J. R.↗

On the Local Convergence of Pattern Search

We examine the local convergence properties of pattern search methods, complementing the previously established global convergence properties for this class of algorithms. We show that the step-length control parameter which appears in the definition of pattern search algorithms provides a reliable asymptotic measure of first-order stationarity. This gives an analytical justification for a traditional stopping criterion for pattern search methods. Using this measure of first-order stationarity, we analyze the behavior of pattern search in the neighborhood of an isolated local minimizer. We show that a recognizable subsequence converges r-linearly to the minimizer.

Dolan, Elizabeth D.↗

Extensions of guiding center motion to higher order

In a static magnetic field, some well-known guiding-center equations maintain their form when extended to next order in gyroradius. In these cases, it is only necessary to include the next order term in the magnetic moment series. The differential equation for guiding-center motion which describes both the parallel and perpendicular velocities correctly through first order in gyroradius is given. The question of how to define the guiding center position through second order arises and is discussed, and second order drifts are derived for one usual definition. The toroidal canonical angular momentum, P-phi, of the guiding center in an axisymmetric field is shown to be conserved using the guiding center velocity correct through first order. When second-order motion is included, P-phi is no longer a constant. The above extensions of guiding-center theory help to resolve the different tokamak orbits obtained either by using the guiding-center equations of motion or by using conservation of P-phi.

Northrop, T. G.↗

Costas loop tracking performance in the presence of antenna switching transients

The present paper examines the impact of phase, amplitude and Doppler transients on the tracking performance of a Costas loop, both in the presence and absence of noise; the loop's lock maintenance capability is also examined. Specific considerations involve analyzing second order loop transient behavior in the absence of noise and the presence of the three transients noted above; here, the operational differential equations are formulated and numerically solved. In the context of noise presence, a first order loop is treated in conjunction with amplitude and phase transients; toward this end the time dependent version of the Fokker-Planck partial differential equation is developed and numerically solved. Selected performance curves describe second order loop behavior in the presence of transients, whose values are of varying severity. First order loop behavior in the presence of noise is characterized by the illustration of suitably defined time dependent 'truncated probability density functions' and associated loss of lock probabilities.

Weinberg, A.↗

Linear feedback rate bounds for regressive channels

Bounds for the linear feedback capacity of m-th order Gaussian autoregressive channels are derived. The upper bound is tighter than that found by Tiernan and Schalwijk (1974) for the feedback capacity of a first-order autoregressive Gaussian channel with not necessarily linear processing. The separation between the upper and lower bounds is small, and it is conjectured that the lower bound converges to the feedback capacity of the first-order channel as the number of signals tends to infinity.

Butman, S. A.↗

Reliability Sensitivity Analysis and Design Optimization of Composite Structures Based on Response Surface Methodology

This report discusses the development and application of two alternative strategies in the form of global and sequential local response surface (RS) techniques for the solution of reliability-based optimization (RBO) problems. The problem of a thin-walled composite circular cylinder under axial buckling instability is used as a demonstrative example. In this case, the global technique uses a single second-order RS model to estimate the axial buckling load over the entire feasible design space (FDS) whereas the local technique uses multiple first-order RS models with each applied to a small subregion of FDS. Alternative methods for the calculation of unknown coefficients in each RS model are explored prior to the solution of the optimization problem. The example RBO problem is formulated as a function of 23 uncorrelated random variables that include material properties, thickness and orientation angle of each ply, cylinder diameter and length, as well as the applied load. The mean values of the 8 ply thicknesses are treated as independent design variables. While the coefficients of variation of all random variables are held fixed, the standard deviations of ply thicknesses can vary during the optimization process as a result of changes in the design variables. The structural reliability analysis is based on the first-order reliability method with reliability index treated as the design constraint. In addition to the probabilistic sensitivity analysis of reliability index, the results of the RBO problem are presented for different combinations of cylinder length and diameter and laminate ply patterns. The two strategies are found to produce similar results in terms of accuracy with the sequential local RS technique having a considerably better computational efficiency.

Rais-Rohani, Masoud↗

A Comparison of Experimental and Theoretical Results for Labyrinth Gas Seals

The basic equations are derived for a two control volume model for compressible flow in a labyrinth seal. The flow is assumed to be completely turbulent and isoenergetic. The wall friction factors are determined using the Blasius formula. Jet flow theory is used for the calculation of the recirculation velocity in the cavity. Linearized zeroth and first order perturbation equations are developed for small motion about a centered position by an expansion in the eccentricity ratio. The zeroth order pressure distribution is found by satisfying the leakage equation. The circumferential velocity distribution is determined by satisfying the momentum equations. The first order equations are solved by a separation of variable solution. Integration of the resultant pressure distribution along and around the seal defines the reaction force developed by the seal and the corresponding dynamic coefficients. The results of this analysis are compared to experimental test results.

Scharrer, Joseph Kirk↗

Some Aspects of Essentially Nonoscillatory (ENO) Formulations for the Euler Equations, Part 3

An essentially nonoscillatory (ENO) formulation is described for hyperbolic systems of conservation laws. ENO approaches are based on smart interpolation to avoid spurious numerical oscillations. ENO schemes are a superset of Total Variation Diminishing (TVD) schemes. In the recent past, TVD formulations were used to construct shock capturing finite difference methods. At extremum points of the solution, TVD schemes automatically reduce to being first-order accurate discretizations locally, while away from extrema they can be constructed to be of higher order accuracy. The new framework helps construct essentially non-oscillatory finite difference methods without recourse to local reductions of accuracy to first order. Thus arbitrarily high orders of accuracy can be obtained. The basic general ideas of the new approach can be specialized in several ways and one specific implementation is described based on: (1) the integral form of the conservation laws; (2) reconstruction based on the primitive functions; (3) extension to multiple dimensions in a tensor product fashion; and (4) Runge-Kutta time integration. The resulting method is fourth-order accurate in time and space and is applicable to uniform Cartesian grids. The construction of such schemes for scalar equations and systems in one and two space dimensions is described along with several examples which illustrate interesting aspects of the new approach.

Chakravarthy, Sukumar R.↗

Effect of resonance-oblateness coupling on a satellite orbit

Second-order effects of the coupling between geopotential resonance and oblateness on a satellite orbit are calculated. These effects arise from the interaction of resonance with the secular changes of the orbit's node, perigee, and mean anomaly. They have the same period and phase as first-order resonance perturbations. But their amplitudes are proportional to the square of the period and dominate the first-order effects as the orbit becomes commensurate. A striking example of this coupling is seen in the 18-day resonance variation of the node of the orbit of the first earth resources technology satellite. Analysis of this 1 arc sec (31 m) variation yields a strong fourteenth-order constraint to the geopotential field for odd degree terms. This constraint (lumped coefficient) is poorly predicted by most current models.

Wagner, C. A.↗

Effect of resonance-oblateness coupling on a satellite orbit

Second order effects of the coupling between geopotential resonance and oblateness on a satellite orbit are calculated. Results show that: (1) these effects arise from the interaction of resonance with the secular changes of the orbit's node, perigee, and mean anomaly; (2) they have the same period and phase as first order resonance perturbations; and (3) their amplitudes are proportional to the square of the period and dominate the first order effects as the orbit becomes commensurate. A striking example of this coupling is seen in the 18 day resonance variation of the node of the orbit of the first earth resources technology satellite. Analysis of this one arc second (31m) variation yielded a strong 14th order constraint to the geopotential for odd degree terms. This constraint is poorly predicted by current models.

Wagner, C. A.↗

VLBI determination of synchronous equatorial satellite orbits.

A first-order theory of satellite orbit determination using time-difference data is developed for the case of a synchronous, equatorial spacecraft. Also, a first order error analysis is presented that permits one to estimate how errors in assumed or measured quantities affect the values of the derived orbital parameters. The theory is applied to observational data obtained during a series of VLBI (Very-Long-Baseline Interferometry) experiments conducted between stations at Rosman, North Carolina, and Mojave, California, using the ATS-3 satellite as a radio source. Orbit parameters derived from the VLBI data show good agreement with those determined by the ATS project office using more conventional tracking techniques.

Ross, S.↗