Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariant”

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 1,027 records · Page 57

The Physics of Turbulence in the Boundary Layer

The geometry of the velocity field in a numerically simulated incompressible turbulent boundary layer over a flat plate at Re theta=670 has been studied using the invariants of the velocity gradient tensor. These invariants are computed at every grid point in the flow and used to form the discriminant. Of primary interest are those regions in the flow where the discriminant is positive; regions where, according to the characteristic equation, the eigenvalues of the velocity gradient tensor are complex. An observer moving with a frame of reference which is attached to a fluid particle lying within such a region would see a local flow pattern of the type stable-focus-stretching or unstable-focus-compressing. When the flow is visualized this way, continuous, connected, large-scale structures are revealed that extend from the point just below the buffer layer out to the beginning of the wake region. These structures are aligned with the mean shear close to the wall and arch in the cross-stream direction away from the wall. In some cases the structures observed are very similar to to the hairpin eddy vision of boundary layer structure proposed by Theodorsen. That the structure of the flow is revealed more effectively by the discriminant rather than by the vorticity is important and adds support to recent observations of the discriminant in a channel flow simulation. Of particular importance is the fact that the procedure does not require the use of an arbitrary threshold in the discriminant. Further analysis using computer flow visualization shows a high degree of spatial correlation between regions of positive discriminant, extreme negative pressure fluctuations and large instantaneous values of Reynolds shear stress.

Kline, Stephen↗

Robust Stabilization of a Class of passive Nonlinear Systems

The problem of feedback stabilization is considered for a class of nonlinear, finite dimensional, time invariant passive systems that are affine in control. Using extensions of the Kalman-Yakubovch lemma, it is shown that such systems can be stabilized by a class of finite demensional, linear, time-invariant controllers which are strictly positive real in the weak or marginal sense. The stability holds regardless of model uncertainties, and is therefore, robust.

Joshi, Suresh M.↗

A Simple Theory of Capillary-Gravity Wave Turbulence

Employing a recently proposed 'multi-wave interaction' theory, inertial spectra of capillary gravity waves are derived. This case is characterized by a rather high degree of nonlinearity and a complicated dispersion law. The absence of scale invariance makes this and some other problems of wave turbulence (e.g., nonlinear inertia gravity waves) intractable by small-perturbation techniques, even in the weak-turbulence limit. The analytical solution obtained in the present work for an arbitrary degree of nonlinearity is shown to be in reasonable agreement with experimental data. The theory explains the dependence of the wave spectrum on wind input and describes the accelerated roll-off of the spectral density function in the narrow sub-range separating scale-invariant regimes of purely gravity and capillary waves, while the appropriate (long- and short-wave) limits yield power laws corresponding to the Zakharov-Filonenko and Phillips spectra.

Glazman, Roman E.↗

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

On the Asymptotic Regimes and the Strongly Stratified Limit of Rotating Boussinesq Equations

Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating Boussinesq equations are analyzed in the asymptotic limit, of strong stratification in the Burger number of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It is shown that in both regimes horizontally averaged buoyancy variable is an adiabatic invariant for the full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with vertical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes on inertial-gravity waves are evidenced. The 'split' of the energy transfer of the vortical and the wave components is established in the Craya-Herring cyclic basis. As the Burger number increases from zero to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. The energy spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. Intermediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure ) are obtained.

Babin, A.↗

A New View of the Dynamics of Reynolds Stress Generation in Turbulent Boundary Layers

The structure of a numerically simulated turbulent boundary layer over a flat plate at Re(theta) = 670 was studied using the invariants of the velocity gradient tensor (Q and R) and a related scalar quantity, the cubic discriminant (D = 27R(exp 2)/4 + Q(exp 3)). These invariants have previously been used to study the properties of the small-scale motions responsible for the dissipation of turbulent kinetic energy. In addition, these scalar quantities allow the local flow patterns to be unambiguously classified according to the terminology proposed by Chong et al. The use of the discriminant as a marker of coherent motions reveals complex, large-scale flow structures that are shown to be associated with the generation of Reynolds shear stress -u'v'(bar). These motions are characterized by high spatial gradients of the discriminant and are believed to be an important part of the mechanism that sustains turbulence in the near-wall region.

Cantwell, Brian J.↗

Self-Similar Compressible Free Vortices

Lie group methods are used to find both exact and numerical similarity solutions for compressible perturbations to all incompressible, two-dimensional, axisymmetric vortex reference flow. The reference flow vorticity satisfies an eigenvalue problem for which the solutions are a set of two-dimensional, self-similar, incompressible vortices. These solutions are augmented by deriving a conserved quantity for each eigenvalue, and identifying a Lie group which leaves the reference flow equations invariant. The partial differential equations governing the compressible perturbations to these reference flows are also invariant under the action of the same group. The similarity variables found with this group are used to determine the decay rates of the velocities and thermodynamic variables in the self-similar flows, and to reduce the governing partial differential equations to a set of ordinary differential equations. The ODE's are solved analytically and numerically for a Taylor vortex reference flow, and numerically for an Oseen vortex reference flow. The solutions are used to examine the dependencies of the temperature, density, entropy, dissipation and radial velocity on the Prandtl number. Also, experimental data on compressible free vortex flow are compared to the analytical results, the evolution of vortices from initial states which are not self-similar is discussed, and the energy transfer in a slightly-compressible vortex is considered.

vonEllenrieder, Karl↗

Feasibility of Decentralized Linear-Quadratic-Gaussian Control of Autonomous Distributed Spacecraft

A distributed satellite formation, modeled as an arbitrary number of fully connected nodes in a network, could be controlled using a decentralized controller framework that distributes operations in parallel over the network. For such problems, a solution that minimizes data transmission requirements, in the context of linear-quadratic-Gaussian (LQG) control theory, was given by Speyer. This approach is advantageous because it is non-hierarchical, detected failures gracefully degrade system performance, fewer local computations are required than for a centralized controller, and it is optimal with respect to the standard LQG cost function. Disadvantages of the approach are the need for a fully connected communications network, the total operations performed over all the nodes are greater than for a centralized controller, and the approach is formulated for linear time-invariant systems. To investigate the feasibility of the decentralized approach to satellite formation flying, a simple centralized LQG design for a spacecraft orbit control problem is adapted to the decentralized framework. The simple design uses a fixed reference trajectory (an equatorial, Keplerian, circular orbit), and by appropriate choice of coordinates and measurements is formulated as a linear time-invariant system.

Carpenter, J. Russell↗

Asymptotic Stability of Interconnected Passive Non-Linear Systems

This paper addresses the problem of stabilization of a class of internally passive non-linear time-invariant dynamic systems. A class of non-linear marginally strictly passive (MSP) systems is defined, which is less restrictive than input-strictly passive systems. It is shown that the interconnection of a non-linear passive system and a non-linear MSP system is globally asymptotically stable. The result generalizes and weakens the conditions of the passivity theorem, which requires one of the systems to be input-strictly passive. In the case of linear time-invariant systems, it is shown that the MSP property is equivalent to the marginally strictly positive real (MSPR) property, which is much simpler to check.

Isidori, A.↗

Feasibility of Decentralized Linear-Quadratic-Gaussian Control of Autonomous Distributed Spacecraft

A distributed satellite formation, modeled as an arbitrary number of fully connected nodes in a network, could be controlled using a decentralized controller framework that distributes operations in parallel over the network. For such problems, a solution that minimizes data transmission requirements, in the context of linear-quadratic-Gaussian (LQG) control theory, was given by Speyer. This approach is advantageous because it is non-hierarchical, detected failures gracefully degrade system performance, fewer local computations are required than for a centralized controller, and it is optimal with respect to the standard LQG cost function. Disadvantages of the approach are the need for a fully connected communications network, the total operations performed over all the nodes are greater than for a centralized controller, and the approach is formulated for linear time-invariant systems. To investigate the feasibility of the decentralized approach to satellite formation flying, a simple centralized LQG design for a spacecraft orbit control problem is adapted to the decentralized framework. The simple design uses a fixed reference trajectory (an equatorial, Keplerian, circular orbit), and by appropriate choice of coordinates and measurements is formulated as a linear time-invariant system.

Carpenter, J. Russell↗

Sliding Mode Thermal Control System for Space Station Furnace Facility

The decoupled control of the nonlinear, multiinput-multioutput, and highly coupled space station furnace facility (SSFF) thermal control system is addressed. Sliding mode control theory, a subset of variable-structure control theory, is employed to increase the performance, robustness, and reliability of the SSFF's currently designed control system. This paper presents the nonlinear thermal control system description and develops the sliding mode controllers that cause the interconnected subsystems to operate in their local sliding modes, resulting in control system invariance to plant uncertainties and external and interaction disturbances. The desired decoupled flow-rate tracking is achieved by optimization of the local linear sliding mode equations. The controllers are implemented digitally and extensive simulation results are presented to show the flow-rate tracking robustness and invariance to plant uncertainties, nonlinearities, external disturbances, and variations of the system pressure supplied to the controlled subsystems.

Jackson Mark E.↗

Possible Statistics of Two Coupled Random Fields: Application to Passive Scalar

We use the relativity postulate of scale invariance to derive the similarity transformations between two coupled scale-invariant random elds at different scales. We nd the equations leading to the scaling exponents. This formulation is applied to the case of passive scalars advected i) by a random Gaussian velocity field; and ii) by a turbulent velocity field. In the Gaussian case, we show that the passive scalar increments follow a log-Levy distribution generalizing Kraichnan's solution and, in an appropriate limit, a log-normal distribution. In the turbulent case, we show that when the velocity increments follow a log-Poisson statistics, the passive scalar increments follow a statistics close to log-Poisson. This result explains the experimental observations of Ruiz et al. about the temperature increments.

Dubrulle, B.↗

Some Thermodynamic Considerations on the Physical and Quantum Nature of Space and Time

It is suggested that the Planck h = m(sub k)c Lambda(sub k) and the Boltzmann k = m(sub k)c nu(sub k)Constants have stochastic foundation. It is further suggested that a body of fluid at equilibrium is composed of a spectrum of molecular clusters (energy levels) the size of which are governed by the Maxwell-Boltzmann distribution function. Brownian motions are attributed to equilibrium between suspensions and molecular clusters. Atomic (molecular) transition between different size atomic- (molecular-) clusters (energy levels) is shown to result in emission/absorption of energy in accordance with Bohr's theory of atomic spectra. Physical space is identified as a tachyonic fluid that is Dirac's stochastic ether or de Broglie's hidden thermostat. Compressibility of physical space, in accordance with Planck's compressible ether, is shown to result in the Lorentz-Fitzgerald contraction, thus providing a causal explanation of relativistic effect in accordance with the perceptions of Poincare and Lorentz. The invariant Schrodinger equation is derived from the invariant Bernoulli equation for incompressible potential flow. Following Heisenberg a temporal uncertainty relation is introduced as Delta(nu(sub Beta)) Delta(Rho(sub Beta)) > = k.

Sohrab, Siavash H.↗

Theory of the Shape of the Heliospheric Termination Shock: Axial Symmetry

The heliospheric termination shock must exhibit asymmetry in its shape, due in part to internal latitude variations in the solar wind, and in part to the special directions defined by the external interstellar flow and/or the galactic magnetic field. This asymmetry shows up naturally in numerical simulations of the interaction between the heliosphere and local interstellar medium. However, to date only one analytical treatment of the theory has appeared, which discussed the modification of the shock shape due to solar wind latitude variations in the presence of spherically symmetric outer boundary conditions. In the present report, we discuss an extension and generalization of this theory to an axially symmetric gasdynamic system, in which departure from spherical symmetry may be due either to internal solar wind variations or to the directional properties of the external interstellar medium. It is shown that for steady flow the post-shock region is characterized by an infinite set of quantities that are conserved along streamlines; among these invariants are the stagnation pressure and a quantity closely related to vorticity. Moreover, for a given latitude profile of the (supersonic) solar wind, the geometry of the termination shock uniquely determines the valued of these invariants at the points where the streamlines emerge from the shock. A self-consistent solution of the problem thus requires a match between the shock geometry and distant boundary conditions. The theory is applied to give estimates of termination-shock asymmetry for various internal and external conditions.

Barnes, Aaron↗

Discrete Abstractions of Hybrid Systems: Verification of Safety and Application to User-Interface Design

Human interaction with a complex control system involves the user, the automation’s discrete mode logic, and the underlying continuous dynamics of the physical system. The user-interface of such systems always displays a reduced set of information about the entire system. Designing interfaces such that all the pertinent information is available and assuring that this information is correct is important for any user-interface, but especially so for safety-critical systems such as automotive systems and autopilots. Here we describe a methodology for the analysis of hybrid control systems that incorporate user interaction, with the goal of assuring that the information provided to the user is correct. That is, the user-interface must contain all information necessary to safely complete a desired procedure or task. We begin with a hybrid system model which incorporates discrete mode logic as well as nonlinear continuous dynamics. Using a hybrid computational tool for reachability, we find the largest region of the state-space in which we can guarantee the state of the system can always remain – this is the safe region of operation. By implementing a controller for safety which arises from this computation, we mathematically guarantee that this safe region is invariant, meaning that the system will always remain within the safe region if the determined controller is used on the boundary of the safe region. Verification within a hybrid framework allows us to account for the continuous dynamics underlying the discrete representations displayed to the user. Using the computed invariant regions as discrete states, we can abstract a discrete event system from this hybrid system with safety restrictions. This abstraction can be used to determine what information must be provided on the display. Furthermore, in cases in which an interface already exists, the abstraction provides the necessary input into existing interface verification methods. We provide two examples: a car traveling through a yellow light at an intersection and an aircraft autopilot in an automatic landing/go-around maneuver. The examples demonstrate the applicability of this methodology to hybrid systems that have operational constraints we can pose in terms of safety. This methodology differs from existing work in hybrid system verification in that we directly account for the user’s interactions with the system.

Meeko Oishi↗

Radio Frequency Trap for Containment of Plasmas in Antimatter Propulsion Systems Using Rotating Wall Electric Fields

A containment apparatus for containing a cloud of charged particles comprises a cylindrical vacuum chamber having a longitudinal axis. Within the vacuum chamber is a containment region. A magnetic field is aligned with the longitudinal axis of the vacuum chamber. The magnetic field is time invariant and uniform in strength over the containment region. An electric field is also aligned with the longitudinal axis of the vacuum chamber and the magnetic field. The electric field is time invariant, and forms a potential well over the containment region. One or more means are disposed around the cloud of particles for inducing a rotating electric field internal to the vacuum chamber. The rotating electric field imparts energy to the charged particles within the containment region and compress the cloud of particles. The means disposed around the outer surface of the vacuum chamber for inducing a rotating electric field are four or more segments forming a segmented ring, the segments conforming to the outer surface of the vacuum chamber. Each of the segments is energized by a separate alternating voltage. The sum of the voltages imposed on each segment establishes the rotating field. When four segments form a ring, the rotating field is obtained by a signal generator applying a sinusoidal signal phase delayed by 90,180 and 270 degrees in sequence to the four segments.

Sims, William Herbert, III↗

Inclusive dielectron cross sections in p + p and p + d interactions at beam energies from 1.04 to 4.88 GeV

Measurements of dielectron production in p + p and p + d collisions with beamkinetic energies from 1.04 to 4.88 GeV are presented. The differential cross section is presented as a function of invariant pair mass, transverse momentum, and rapidity. The shapes of the mass spectra and their evolution with beam energy provide information about the relative importance of the various dielectron production mechanisms in this energy regime. The p + d to p + p ratio of the dielectron yield is also presented as a function of invariant pair mass, transverse momentum, and rapidity. The shapes of the transverse momentum and rapidity spectra from the p + d and p + p systems are found to be similar to one another for each of the beam energies studied. The beam energy dependence of the integrated cross sections is also presented.

NASA Discipline Radiation Health↗

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Number

Present simulation of turbulent flows involving shock wave/boundary layer interaction invariably overestimates heat flux by almost a factor of two. One possible reason for such a performance is a result of the fact that the turbulence models employed make use of Morkovin's hypothesis. This hypothesis is valid for non-hypersonic Mach numbers and moderate rates of heat transfer. At hypersonic Mach numbers, high rates of heat transfer exist in regions where shock wave/boundary layer interactions are important. As a result, one should not expect traditional turbulence models to yield accurate results. The goal of this investigation is to explore the role of a variable Prandtl number formulation in predicting heat flux in flows dominated by strong shock wave/boundary layer interactions. The intended applications involve external flows in the absence of combustion such as those encountered in supersonic inlets. This can be achieved by adding equations for the temperature variance and its dissipation rate. Such equations can be derived from the exact Navier-Stokes equations. Traditionally, modeled equations are based on the low speed energy equation where the pressure gradient term and the term responsible for energy dissipation are ignored. It is clear that such assumptions are not valid for hypersonic flows. The approach used here is based on the procedure used in deriving the k-zeta model, in which the exact equations that governed k, the variance of velocity, and zeta, the variance of vorticity, were derived and modeled. For the variable turbulent Prandtl number, the exact equations that govern the temperature variance and its dissipation rate are derived and modeled term by term. The resulting set of equations are free of damping and wall functions and are coordinate-system independent. Moreover, modeled correlations are tensorially consistent and invariant under Galilean transformation. The final set of equations will be given in the paper.

Xiao, X.↗