Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hamiltonian systems”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗
NASA NTRS1973

Statistical mechanics of negative temperature states

The dynamics of two-dimensional interacting line vortices is identical to that of the two-dimensional electrostatic guiding center plasma. Both are Hamiltonian systems and are therefore susceptible to statistical mechanical treatments. The predictions of the microcanonical ensemble are explored for this system. Interest focuses primarily on the regime of total positive interaction energy, which should be above the Onsager negative temperature threshold. Calculations of the probability distribution for a component by means of the central limit theorem are carried out in the manner of Khinchin. The probability distribution of a component reduced to the usual Gibbs distribution in the regime of positive temperatures, and is still explicitly calculable for negative temperatures. The negative temperature states are neither quiescent nor spatially uniform. Expressions for the temperature are explicitly provided in terms of the total particle energy and particle number. A BBGKY hierarchy can be derived for both temperature regimes. Numerical simulations involving solutions of the equations of motion of 4008 particles are presented.

Montgomery, D. C.↗
NASA NTRS1984

On the topological stability of magnetostatic equilibria

The topological stability of MHD equilibria is investigated by exploring the formal analogy, in the ideal MHD limit, between the topology of magnetic lines of force in coordinate space and the topology of integral surfaces of one- and two-dimensional Hamiltonian systems in phase space. It is demonstrated that in an astrophysical setting, symmetric magnetostatic equilibria satisfying the ideal MHD equations are exceptional. The principal result of the study is that previous infinitesimal perturbation theory calculations can be generalized to include finite-amplitude and symmetry-breaking effects. The effect of the ergodicity of perturbed symmetric equilibria on heat dispersal in magnetically dominated plasmas is discussed.

Tsinganos, K. C.↗
NASA NTRS1985

Meteorites may follow a chaotic route to earth

A numerical integration was employed to show that the 3/1 Kirkwood gap region between Jupiter and the sun can serve as the source for asteroidal debris that strikes the earth as chondrites. A previous study by Wetherill (1968) considered planetary close encounters, secular resonances, collisional destruction in space and subsequent atmospheric abrasion of the chondrites. The factors would have resulted in chaotic trajectories for asteroids in the asteroid belt, i.e., eccentricities large enough to make the objects earth-crossing. A sampling of trajectories which could reach large eccentricities is integrated in terms of deterministic solutions of Hamiltonian systems of equations which behave erratically and are dependent on the initial conditions. The examinations successfully predicted a trajectory with a peak eccentricity of 0.652, a time interval of 340,000 yr, and a semimajor axis between 2.4-2.53 AU that was earth-crossing.

Wisdom, J.↗
NASA NTRS1985

Pattern formation by particles settling in viscous flows

Results are reported of several simulations of the motions of small objects in two-dimensional cellular flows performed to study the evolution of the distributions of particle flows over time in processes such as sedimentation. Stream functions are defined in the form of a Hamiltonian system for the particle motions, with inertial and particle acceleration being neglected. Conditions are found for which the particles will become trapped in an oscillatory motion or the particle orbits attain an oscillary motion. Various simulation images are provided to illustrate the usefulness of the technique for modeling fractal swarms.

Smith, L. A.↗
NASA NTRS1990

Fractal boundaries in magnetotail particle dynamics

It has been recently established that particle dynamics in the magnetotail geometry can be described as a nonintegrable Hamiltonian system with well-defined entry and exit regions through which stochastic orbits can enter and exit the system after repeatedly crossing the equatorial plane. It is shown that the phase space regions occupied by orbits of different numbers of equatorial crossings or different exit modes are separated by fractal boundaries. The fractal boundaries in an entry region for stochastic orbits are examined and the capacity dimension is determined.

Chen, J.↗
NASA NTRS1991

A family of asymptotically stable control laws for flexible robots based on a passivity approach

A general family of asymptotically stabilizing control laws is introduced for a class of nonlinear Hamiltonian systems. The inherent passivity property of this class of systems and the Passivity Theorem are used to show the closed-loop input/output stability which is then related to the internal state space stability through the stabilizability and detectability condition. Applications of these results include fully actuated robots, flexible joint robots, and robots with link flexibility.

Lanari, Leonardo↗
NASA NTRS1993

Nonlinear modes of the tensor Dirac equation and CPT violation

Recently, it has been shown that Dirac's bispinor equation can be expressed, in an equivalent tensor form, as a constrained Yang-Mills equation in the limit of an infinitely large coupling constant. It was also shown that the free tensor Dirac equation is a completely integrable Hamiltonian system with Lie algebra type Poisson brackets, from which Fermi quantization can be derived directly without using bispinors. The Yang-Mills equation for a finite coupling constant is investigated. It is shown that the nonlinear Yang-Mills equation has exact plane wave solutions in one-to-one correspondence with the plane wave solutions of Dirac's bispinor equation. The theory of nonlinear dispersive waves is applied to establish the existence of wave packets. The CPT violation of these nonlinear wave packets, which could lead to new observable effects consistent with current experimental bounds, is investigated.

Reifler, Frank J.↗
NASA NTRS1994

Modal decomposition of Hamiltonian variational equations

Over any finite arc of trajectory, the variational equations of a Hamiltonian system can be separated into 'normal' modes. This transformation is canonical, and the Lyapunov exponents over the trajectory arc occur as positive/negative pairs for conjugate modes, while the modal vectors remain unit vectors. This decomposition effectively solves the variational equations for any canonical, linear-dependent system. As an example, we study the Voyager I trajectory. In an interplanetary flyby, some of the modal variables increase by very large multiplicative factors, but this means that their conjugate modal variables decrease by those same very large multiplicative vectors. Maneuver strategies for this case are explored, and the minimum delta upsilon maneuver is found.

Wiesel, William E.↗
NASA NTRS2012

Tuning the Sensitivity of an Optical Cavity with Slow and Fast Light

We have measured mode pushing by the dispersion of a rubidium vapor in a Fabry-Perot cavity and have shown that the scale factor and sensitivity of a passive cavity can be strongly enhanced by the presence of such an anomalous dispersion medium. The enhancement is the result of the atom-cavity coupling, which provides a positive feedback to the cavity response. The cavity sensitivity can also be controlled and tuned through a pole by a second, optical pumping, beam applied transverse to the cavity. Alternatively, the sensitivity can be controlled by the introduction of a second counter-propagating input beam that interferes with the first beam, coherently increasing the cavity absorptance. We show that the pole in the sensitivity occurs when the sum of the effective group index and an additional cavity delay factor that accounts for mode reshaping goes to zero, and is an example of an exceptional point, commonly associated with coupled non-Hermitian Hamiltonian systems. Additionally we show that a normal dispersion feature can decrease the cavity scale factor and can be generated through velocity selective optical pumping

Smith, David D.↗
NASA NTRS1992

Non-Hamiltonian perturbations of integrable systems and resonance trapping

This paper studies general, non-Hamiltonian perturbations of integrable systems with two degrees of freedom and derives conditions for temporary and permanent resonance trapping. The analysis involves a noncanonical transformation of variables near the resonant manifold and averaging with respect to the fast phase to investigate oscillatory behavior on the intermediate timescale. The resulting reduced system is Hamiltonian to leading order and permits, after averaging on the intermediate, or libration, timescale, a canonical transformation to action-angle variables in the oscillation zone. The final system so obtained reveals the possible existence of two- and three-dimensional invariant tori in the vicinity of the resonant manifold. An explicit divergence condition for general perturbations to be dissipative on the slow timescale follows from the analysis. An application of this approach to the problem of resonant trapping and escape is outlined for the restricted problem of three bodies subject to dissipative perturbations with a radial symmetry.

Ghil, M.↗
NASA NTRSDate not supplied

Hamiltonian Optimal Control of Distributed Lagrangian Systems

This lecture presents a Hamiltonian control method and a distributed optimal control method for distributed Lagrangian systems. The distributed optimal control theory is formulated using a semi-group abstraction resulting in an integro-differential Riccati equation.

Distributed Optimal Control↗
NASA NTRS1980

Hamiltonian indices and rational spectral densities

Several (global) topological properties of various spaces of linear systems, particularly symmetric, lossless, and Hamiltonian systems, and multivariable spectral densities of fixed McMillan degree are announced. The study is motivated by a result asserting that on a connected but not simply connected manifold, it is not possible to find a vector field having a sink as its only critical point. In the scalar case, this is illustrated by showing that only on the space of McMillan degree = /Cauchy index/ = n, scalar transfer functions can one define a globally convergent vector field. This result holds both in discrete-time and for the nonautonomous case. With these motivations in mind, theorems of Bochner and Fogarty are used in showing that spaces of transfer functions defined by symmetry conditions are, in fact, smooth algebraic manifolds.

Byrnes, C. I.↗
NASA NTRS2006

Using Quasiparticle Poisoning To Detect Photons

According to a proposal, a phenomenon associated with excitation of quasiparticles in certain superconducting quantum devices would be exploited as a means of detecting photons with exquisite sensitivity. The phenomenon could also be exploited to perform medium-resolution spectroscopy. The proposal was inspired by the observation that Coulomb blockade devices upon which some quantum logic gates are based are extremely sensitive to quasiparticles excited above the superconducting gaps in their leads. The presence of quasiparticles in the leads can be easily detected via the charge states. If quasiparticles could be generated in the leads by absorption of photons, then the devices could be used as very sensitive detectors of electromagnetic radiation over the spectral range from x-rays to submillimeter waves. The devices in question are single-Cooper-pair boxes (SCBs), which are mesoscopic superconducting devices developed for quantum computing. An SCB consists of a small superconducting island connected to a reservoir via a small tunnel junction and connected to a voltage source through a gate capacitor. An SCB is an artificial two-level quantum system, the Hamiltonian of which can be controlled by the gate voltage. One measures the expected value of the charge of the eigenvectors of this quantum system by use of a radio-frequency single-electron transistor. A plot of this expected value of charge as a function of gate voltage resembles a staircase that, in the ideal case, consists of steps of height 2 e (where e is the charge of one electron). Experiments have shown that depending on the parameters of the device, quasiparticles in the form of "broken" Cooper pairs present in the reservoir can tunnel to the island, giving rise to steps of 1 e. This effect is sometimes called "poisoning." Simulations have shown that an extremely small average number of quasiparticles can generate a 1-e periodic signal. In a device according to the proposal, this poisoning would be turned to advantage. Depending on the wavelength, an antenna or other component would be used to couple radiation into the reservoir, wherein the absorption of photons would break Cooper pairs, thereby creating quasiparticles that, in turn, would tunnel to the island, creating a 1-e signal. On the basis of conservative estimates of device parameters derived from experimental data and computational simulations that fit the data, it has been estimated that the noise equivalent power of a device according to the proposal could be as low as 6 10(exp -22) W/Hz(exp 1/2). It has also been estimated that the spectroscopic resolution (photon energy divided by increment of photon energy) of such a device in visible light would exceed 100.

Echternach, Pierre↗
NASA NTRS1985

Vortex studies relating to boundary layer turbulence and noise

The present study considers the two-dimensional case of an array of N rectilinear, like-sign vortices above an infinite flat boundary. The method of images can be employed with this configuration to reduce the problem to that of 2N vortices in free space, constrained by 2N symmetry relations. This system is Hamiltonian and therefore certain invariants of the motion are known. Futher, from the Hamiltonian constant, the equations of motion are readily derived and may be integrated numerically to determine the vortex trajectories. This knowledge of the time-dependent vortex motion then allows the resulting noise radiation to be computed by standard aeroacoustic techniques.

Adelman, J. L.↗