SEARCH · Search NASA
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.
An all metal energy dissipator construction is disclosed for dissipating kinetic energy force (F) by rolling balls which are forced by a tapered surface on an expandable sleeve to frictionally load a force rod. The balls are maintained in an initial position by a plate member which is biased by a spring member. A spring member returns the force rod to its initial position after a loading force is removed.
Stability of the closed-loop system is of vital importance when controller synthesis is being performed using a reduced order model of the open-loop plant. For passive systems such as flexible space structures, closed-loop stability is often ensured by low-authority dissipative controllers using collocated and compatible sensors and actuators. This paper considers the problem of eigensystem assignment for flexible structures using constant-gain output feedback, under the dissipativity constraints requiring that the gain matrices be nonnegative definite. A novel method is developed in which a sequential procedure is used to assign one complex-conjugate pair of closed-loop eigenvalues in each step. At each step, an iterative technique is employed for generating a nonnegative definite gain matrix that assigns a new pair while keeping the previously assigned pairs unchanged. This method can assign up to m closed-loop system eigenvalues, where m is the number of control inputs. A numerical example of damping enhancement for a large flexible structure is presented to demonstrate the approach.
The stability characteristics of dynamic dissipative compensators are investigated for multibody flexible space structures having nonlinear dynamics. The problem addressed is that of proving asymptotic stability of dynamic dissipative compensators. The stability proof uses the Liapunov approach and exploits the inherent passivity of such systems. For such systems these compensators are shown to be robust to parametric uncertainties and unmodeled dynamics. The results are applicable to a large class of structures such as flexible space structures with articulated flexible appendages.
Dissipative energy flow controllers are designed for interconnected modal subsystems. Active feedback controllers for vibration suppression are then viewed as either an additional subsystem or a dissipative coupling. These controllers, which are designed by the LQG positive real control approach, maximize energy flow from a specified modal subsystem.
Nonlocality of the triad interactions in the dissipation range of developed turbulence is investigated by numerical simulation and the quasi-normal theories. It is found that the energy transfer is dominated by nonlocal triad interactions over the wavenumber range 0.1 is less than k/k(sub d) is less than 4, where k(sub d) is the Kolmogorov wave number. The nonlocality of the interaction has a close relation with the power of an algebraic prefactor of the exponential decay of the energy spectrum in the far-dissipation range.
Since one can rarely write down the analytical solutions to nonlinear dissipative partial differential equations (PDE's), it is important to understand whether, and in what sense, the behavior of approximating schemes to these equations reflects the true dynamics of the original equations. Further, because standard error estimates between approximations of the true solutions coming from spectral methods - finite difference or finite element schemes, for example - and the exact solutions grow exponentially in time, this analysis provides little value in understanding the infinite time behavior of a given approximating scheme. The notion of the global attractor has been useful in quantifying the infinite time behavior of dissipative PDEs, such as the Navier-Stokes equations. Loosely speaking, the global attractor is all that remains of a sufficiently large bounded set in phase space mapped infinitely forward in time under the evolution of the PDE. Though the attractor has been shown to have some nice properties - it is compact, connected, and finite dimensional, for example - it is in general quite complicated. Nevertheless, the global attractor gives a way to understand how the infinite time behavior of approximating schemes such as the ones coming from a finite difference, finite element, or spectral method relates to that of the original PDE. Indeed, one can often show that such approximations also have a global attractor. We therefore only need to understand how the structure of the attractor for the PDE behaves under approximation. This is by no means a trivial task. Several interesting results have been obtained in this direction. However, we will not go into the details. We mention here that approximations generally lose information about the system no matter how accurate they are. There are examples that show certain parts of the attractor may be lost by arbitrary small perturbations of the original equations.
A variety of modifications to the modeled dissipation rate transport equation that have been proposed during the past two decades to account for rotational strains are examined. The models are subjected to two crucial test cases: the decay of isotropic turbulence in a rotating frame and homogeneous shear flow in a rotating frame. It is demonstrated that these modifications do not yield substantially improved predictions for these two test cases and in many instances give rise to unphysical behavior. An alternative proposal, based on the use of the tensor dissipation rate, is made for the development of improved models.
In the first part of this contribution we show that the master equation derived from the generalized version of the nonlinear Doebner-Goldin equation leads to the squeezing of one of the quadratures. In the second part we consider two familiar Hamiltonians, the Bateman- Caldirola-Kanai and the optical parametric oscillator; going back to their classical Lagrangian form we introduce a stochastic force and a dissipative factor. From this new Lagrangian we obtain a modified Hamiltonian that treats adequately the simultaneous amplification and dissipation phenomena, presenting squeezing, too.
For a dissipative system with Ohmic friction, we obtain a simple and exact solution for the wave function of the system plus the bath. It is described by the direct product in two independent Hilbert space. One of them is described by an effective Hamiltonian, the other represents the effect of the bath, i.e., the Brownian motion, thus clarifying the structure of the wave function of the system whose energy is dissipated by its interaction with the bath. No path integral technology is needed in this treatment. The derivation of the Weisskopf-Wigner line width theory follows easily.
A model is presented to compute the turbulent kinetic energy dissipation length scale l(sub epsilon) in a stably stratified shear flow. The expression for l(sub epsilon) is derived from solving the spectral balance equation for the turbulent kinetic energy. The buoyancy spectrum entering such equation is constructed using a Lagrangian timescale with modifications due to stratification. The final result for l(sub epsilon) is given in algebraic form as a function of the Froude number Fr and the flux Richardson number R(sub f), l(sub epsilon) = l(sub epsilon)(Fr, R(sub f). The model predicts that for R(sub f) less than R(sub fc), l(sub epsilon) decreases with stratification. An attractive feature of the present model is that it encompasses, as special cases, some seemingly different models for l(sub epsilon) that have been proposed in the past by Deardorff, Hunt et al., Weinstock, and Canuto and Minotti. An alternative form for the dissipation rate epsilon is also discussed that may be useful when one uses a prognostic equation for the heat flux. The present model is applicable to subgrid-scale models, which are needed in large eddy simulations (LES), as well as to ensemble average models. The model is applied to predict the variation of l(sub epsilon) with height z in the planetary boundary layer. The resulting l(sub epsilon) versus z profile reproduces very closely the nonmonotonic profile of l(sub epsilon) exhibited by many LES calculations, beginning with the one by Deardorff in 1974.
Studies of inclined rings inside galaxy potentials have mostly considered the influence of self-gravity and viscous dissipation separately. In this study, we construct models of highly inclined ('polar') rings in an external potential including both self-gravity and dissipation due to a drag force. We do not include pressure forces and thus ignore shock heating that dominates the evolution of gaseous rings inside strongly nonspherical potentials. We adopt an oblate spheroidal scale-free logarithmic potential with axis ratio q = 0.85 and an initial inclination of 80 deg for the self-gravitating rings. We find that stellar (dissipationless) rings suffer from mass loss during their evolution. Mass loss also drives a secular change of the mean inclination toward the poles of the potential. As much as half of the ring mass escapes in the process and forms an inner and an outer shell of precessing orbits. If the remaining mass is more than approximately 0.02 of the enclosed galaxy mass, rings remain bound and do not fall apart from differential precession. The rings precess at a constant rate for more than a precession period tau(sub p) finding the configuration predicted by Sparke in 1986 which warps at larger radii toward the poles of the potential. We model shear viscosity with a velocity-dependent drag force and find that nuclear inflow dominates over self-gravity if the characteristic viscous inflow time scale tau(sub vi) is shorter than approximately 25(tau(sub p)). Rings with (tau(sub vi))/(tau(sub p)) less than or approximately equal to 25 collapse toward the nucleus of the potential within one precession period independent of the amount of self-gravity. Our results imply that stars and gas in real polar rings exhibit markedly different dynamical evolutions.
The kinetic structure of slow magnetosonic shocks obtained from 1- and 2-D hybrid (fluid electrons, particle ions) simulations is discussed. The magnetic structure of the shock varies considerably with the Mach number. At low Mach numbers, the shock may be laminar with no associated waves, whereas, at high Mach numbers the shock may be unsteady with Alfvenic turbulence in the upstream and the downstream regions. The nature of ion dissipation also varies with the Mach number. At weak laminar shocks, ion dissipation consists of adiabatic heating and reflection of a small number of ions back into the upstream. At stronger shocks, particle scattering by Alfven/ion-cyclotron waves in the upstream contributes to the heating which continues well into the downstream region.
Two equation and higher order closures for compressible turbulence fail to capture the compressible wall layers' log scaling. Accounting for the distinction between Favre and Reynolds averaged variables in the compressible moment equations indicate that turbulent transport expressions obtained using the 'variable density approximation' are in error. The error is related to the enstrophy, a Reynolds averaged variable appearing in the equation for the Favre averaged k; recognizing this fact an expression for the transport of dissipation consistent with simple mixing length arguments is obtained. Within the (limited) context of a gradient transport hypothesis a rational form for the turbulent transport of the dissipation is found. Modestly better agreement with the well established compressible Van Driest log scaling is found in k - epsilon calculation.
Several algorithms for introducing artificial dissipation into a central difference approximation to the Euler and Navier Stokes equations are considered. The focus of the paper is on the convective upwind and split pressure (CUSP) scheme, which is designed to support single interior point discrete shock waves. This scheme is analyzed and compared in detail with scalar and matrix dissipation (MATD) schemes. Resolution capability is determined by solving subsonic, transonic, and hypersonic flow problems. A finite-volume discretization and a multistage time-stepping scheme with multigrid are used to compute solutions to the flow equations. Numerical results are also compared with either theoretical solutions or experimental data. For transonic airfoil flows the best accuracy on coarse meshes for aerodynamic coefficients is obtained with a simple MATD scheme.
The dissipation rate transport equation remains the most uncertain part of turbulence modeling. The difficulties arc increased when external agencies like rotation prevent straightforward dimensional analysis from determining the correct form of the modelled equation. In this work, the dissipation rate transport equation and subgrid scale models for rotating turbulence are derived from an analytical statistical theory of rotating turbulence. In the strong rotation limit, the theory predicts a turbulent steady state in which the inertial range energy spectrum scales as k(sup -2) and the turbulent time scale is the inverse rotation rate. This scaling has been derived previously by heuristic arguments.
Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.
Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.
Optimal controller design is addressed for a class of linear, time-invariant systems which are dissipative with respect to a quadratic power function. The system matrices are assumed to be affine functions of uncertain parameters confined to a convex polytopic region in the parameter space. For such systems, a method is developed for designing a controller which is dissipative with respect to a given power function, and is simultaneously optimal in the linear-quadratic-Gaussian (LQG) sense. The resulting controller provides robust stability as well as optimal performance. Three important special cases, namely, passive, norm-bounded, and sector-bounded controllers, which are also LQG-optimal, are presented. The results give new methods for robust controller design in the presence of parametric uncertainties.