Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PHASE SPACE”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

A direct Vlasov code to study the non-stationary current collection by a cylindrical Langmuir probe

The time-dependent current collection by a cylindrical Langmuir probe, whose bias is suddenly changed from zero to a positive or negative finite value, is studied with a novel direct Vlasov code. The numerical algorithm is based on finite-difference formulas to approximate spatial and velocity derivatives and the time integration is carried out with an explicit Runge-Kutta method, or in the case of probe radius small compared with the Debye length, by using the unconditionally stable backward Euler scheme. Both electrons and ions are treated kinetically by the code, which implements initial and boundary conditions that are consistent with the presence of the probe. Within the considered parameter range, the plasma sheath around the probe exhibited an overshoot and it later recovered a steady state. Phase space diagrams of the particle trajectories revealed the presence of a trapped population of particles. The dependence of this population as a function of the probe radius is presented as well as a comparison with the stationary theory. The performance of the code and a comparison with previously used particle-in-cell algorithms are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Lower hybrid wavepacket stochasticity revisited

Analysis is presented in support of the explanation in Ref. [1] for the observation of relativistic electrons during Lower Hybrid (LH) operation in EC pre-heated plasma at the WEGA stellarator [1,2]. LH power from the WEGA TE11 circular waveguide, 9 cm diameter, un-phased, 2.45 GHz antenna, is radiated into a B≅0.5 T, Ðœ„n{sub e}≅5×10{sup 17} 1/m{sup 3} plasma at T{sub e}≅10 eV bulk temperature with an EC generated 50 keV component [1]. The fast electrons cycle around flux or drift surfaces with few collisions, sufficient for randomizing phases but insufficient for slowing fast electrons down, and thus repeatedly interact with the rf field close to the antenna mouth, gaining energy in the process. Our antenna calculations reveal a standing electric field pattern at the antenna mouth, with which we formulate the electron dynamics via a relativistic Hamiltonian. A simple approximation of the equations of motion leads to a relativistic generalization of the area-preserving Fermi-Ulam (F-U) map [3], allowing phase-space global stochasticity analysis. At typical WEGA plasma and antenna conditions, the F-U map predicts an LH driven current of about 230 A, at about 225 W of dissipated power, in good agreement with the measurements and analysis reported in [1].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Canonical transformation for trapped/passing guiding-center orbits in axisymmetric tokamak geometry

The generating function for the canonical transformation from the parallel canonical coordinates (s,p{sub ||}) to the action-angle coordinates (ζ, J) for trapped/passing guiding-center orbits in axisymmetric tokamak geometry is presented. Drawing on the analogy between the phase-space portraits of the librating/rotating pendulum and the trapped/passing guiding-center orbits, the generating function is expressed in terms of the Jacobi zeta function, which can then readily be used to obtain an explicit expression for the bounce-center transformation for trapped/passing-particle guiding-center orbits in axisymmetric tokamak geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A study of the early-stage evolution of relativistic electron-ion shock using three-dimensional particle-in-cell simulations

We report the results of a 3D particle-in-cell simulation carried out to study the early-stage evolution of the shock formed when an unmagnetized relativistic jet interacts with an ambient electron-ion plasma. Full-shock structures associated with the interaction are observed in the ambient frame. When open boundaries are employed in the direction of the jet, the forward shock is seen as a hybrid structure consisting of an electrostatic shock combined with a double layer, while the reverse shock is seen as a double layer. The ambient ions show two distinct features across the forward shock: a population penetrating into the shocked region from the precursor region and an accelerated population escaping from the shocked region into the precursor region. This behavior is a signature of a combination of an electrostatic shock and a double layer. Jet electrons are seen to be electrostatically trapped between the forward and reverse shock structures showing a ring-like distribution in a phase-space plot, while ambient electrons are thermalized and become essentially isotropic in the shocked region. The magnetic energy density grows to a few percent of the jet kinetic energy density at both the forward and the reverse shock transition layers in a rather short time scale. We see little disturbance of the jet ions over this time scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Various sizes of sliding event bursts in the plastic flow of metallic glasses based on a spatiotemporal dynamic model

This paper presents a spatiotemporal dynamic model based on the interaction between multiple shear bands in the plastic flow of metallic glasses during compressive deformation. Various sizes of sliding events burst in the plastic deformation as the generation of different scales of shear branches occurred; microscopic creep events and delocalized sliding events were analyzed based on the established model. This paper discusses the spatially uniform solutions and traveling wave solution. The phase space of the spatially uniform system applied in this study reflected the chaotic state of the system at a lower strain rate. Moreover, numerical simulation showed that the microscopic creep events were manifested at a lower strain rate, whereas the delocalized sliding events were manifested at a higher strain rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

On the simulation of indistinguishable fermions in the many-body Wigner formalism

The simulation of quantum systems consisting of interacting, indistinguishable fermions is an incredible mathematical problem which poses formidable numerical challenges. Many sophisticated methods addressing this problem are available which are based on the many-body Schrödinger formalism. Recently a Monte Carlo technique for the resolution of the many-body Wigner equation has been introduced and successfully applied to the simulation of distinguishable, spinless particles. This numerical approach presents several advantages over other methods. Indeed, it is based on an intuitive formalism in which quantum systems are described in terms of a quasi-distribution function, and highly scalable due to its Monte Carlo nature. In this work, we extend the many-body Wigner Monte Carlo method to the simulation of indistinguishable fermions. To this end, we first show how fermions are incorporated into the Wigner formalism. Then we demonstrate that the Pauli exclusion principle is intrinsic to the formalism. As a matter of fact, a numerical simulation of two strongly interacting fermions (electrons) is performed which clearly shows the appearance of a Fermi (or exchange–correlation) hole in the phase-space, a clear signature of the presence of the Pauli principle. To conclude, we simulate 4, 8 and 16 non-interacting fermions, isolated in a closed box, and show that, as the number of fermions increases, we gradually recover the Fermi–Dirac statistics, a clear proof of the reliability of our proposed method for the treatment of indistinguishable particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

A kinetic model of radiating electrons

A kinetic theory is developed to describe radiating electrons whose motion is governed by the Lorentz-Dirac equation. This gives rise to a generalized Vlasov equation coupled to an equation for the evolution of the physical submanifold of phase space. The pathological solutions of the 1-particle theory may be removed by expanding the latter equation in powers of {tau}:=q{sup 2}/6{pi}m. The radiation-induced change in entropy is explored and its physical origin is discussed. As a simple demonstration of the theory, the radiative damping rate of longitudinal plasma waves is calculated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The velocity distribution of pickup He{sup +} measured at 0.3 AU by MESSENGER

During its interplanetary trajectory in 2007-2009, the MErcury Surface, Space ENvrionment, GEochemistry, and Ranging (MESSENGER) spacecraft passed through the gravitational focusing cone for interstellar helium multiple times at a heliocentric distance R ≈ 0.3 AU. Observations of He{sup +} interstellar pickup ions made by the Fast Imaging Plasma Spectrometer sensor on MESSENGER during these transits provide a glimpse into the structure of newly formed inner heliospheric pickup-ion distributions. This close to the Sun, these ions are picked up in a nearly radial interplanetary magnetic field. Compared with the near-Earth environment, pickup ions observed near 0.3 AU will not have had sufficient time to be energized substantially. Such an environment results in a nearly pristine velocity distribution function that should depend only on pickup-ion injection velocities (related to the interstellar gas), pitch-angle scattering, and cooling processes. From measured energy-per-charge spectra obtained during multiple spacecraft observational geometries, we have deduced the phase-space density of He{sup +} as a function of magnetic pitch angle. Our measurements are most consistent with a distribution that decreases nearly monotonically with increasing pitch angle, rather than the more commonly modeled isotropic or hemispherically symmetric forms. These results imply that pitch-angle scattering of He{sup +} may not be instantaneous, as is often assumed, and instead may reflect the velocity distribution of initially injected particles. In a slow solar wind stream, we find a parallel-scattering mean free path of λ {sub ||} ∼ 0.1 AU and a He{sup +} production rate of ∼0.05 m{sup –3} s{sup –1} within 0.3 AU.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Analysis and prediction of aperiodic hydrodynamic oscillatory time series by feed-forward neural networks, fuzzy logic, and a local nonlinear predictor

Forecasting of aperiodic time series is a compelling challenge for science. In this work, we analyze aperiodic spectrophotometric data, proportional to the concentrations of two forms of a thermoreversible photochromic spiro-oxazine, that are generated when a cuvette containing a solution of the spiro-oxazine undergoes photoreaction and convection due to localized ultraviolet illumination. We construct the phase space for the system using Takens' theorem and we calculate the Lyapunov exponents and the correlation dimensions to ascertain the chaotic character of the time series. Finally, we predict the time series using three distinct methods: a feed-forward neural network, fuzzy logic, and a local nonlinear predictor. We compare the performances of these three methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Fock expansion of multimode pure Gaussian states

The Fock expansion of multimode pure Gaussian states is derived starting from their representation as displaced and squeezed multimode vacuum states. The approach is new and appears to be simpler and more general than previous ones starting from the phase-space representation given by the characteristic or Wigner function. Fock expansion is performed in terms of easily evaluable two-variable Hermite–Kampé de Fériet polynomials. A relatively simple and compact expression for the joint statistical distribution of the photon numbers in the different modes is obtained. In particular, this result enables one to give a simple characterization of separable and entangled states, as shown for two-mode and three-mode Gaussian states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Varying alpha and cosmic acceleration in Brans-Dicke-BSBM theory: stability analysis and observational tests

By integration of generalized BSBM and Brans-Dicke cosmological models, in this article, we investigate the theoretical framework of fine structure constant variation and current cosmic acceleration. We first develop a mathematical formalism to analyse the stability of the model. By employing observational data to constrain the model parameters, phase space is performed and the attractor solutions of the model are detected. We then examine the model against observational data such as observational Hubble parameter dataset and quasar absorption spectra. The results confirm current universe acceleration and also predicts fine structure constant variation. Furthermore, extrapolation of the best fitted model in high redshift (z > 15) illustrates a significantly larger variation of fine structure constant in earlier epoch of the universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Investigating the Relation between Galaxy Properties and the Gaussianity of the Velocity Distribution of Groups and Clusters

We investigate the dependence of stellar population properties of galaxies on group dynamical stage for a subsample of the Yang catalog. We classify groups according to their galaxy velocity distribution into Gaussian (G) and Non-Gaussian (NG). Using two totally independent approaches, we have shown that our measurement of Gaussianity is robust and reliable. Our sample covers Yang’s groups in the redshift range 0.03 ≤ z ≤ 0.1, with mass ≥ 10{sup 14}M{sub ⊙}. The new method, called Hellinger Distance, to determine whether a group has a velocity distribution Gaussian or NG is very effective in distinguishing between the two families. NG groups present halo masses higher than the G ones, confirming previous findings. Examining the skewness and kurtosis of the velocity distribution of G and NG groups, we find that faint galaxies in NG groups are mainly infalling, for the first time, into the groups. We show that considering only faint galaxies in the outskirts; those in NG groups are older and more metal-rich than those in G groups. Also, examining the Projected Phase Space of cluster galaxies, we see that bright and faint galactic systems in G groups are in dynamical equilibrium—which does not seem to be the case in NG groups. These findings suggest that NG systems have a higher infall rate, assembling more galaxies that have experienced preprocessing before entering the group.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

E × B shear pattern formation by radial propagation of heat flux waves

A novel theory to describe the formation of E×B flow patterns by radially propagating heat flux waves is presented. A model for heat avalanche dynamics is extended to include a finite delay time between the instantaneous heat flux and the mean flux, based on an analogy between heat avalanche dynamics and traffic flow dynamics. The response time introduced here is an analogue of the drivers' response time in traffic dynamics. The microscopic foundation for the time delay is the time for mixing of the phase space density. The inclusion of the finite response time changes the model equation for avalanche dynamics from Burgers equation to a nonlinear telegraph equation. Based on the telegraph equation, the formation of heat flux jams is predicted. The growth rate and typical interval of jams are calculated. The connection of the jam interval to the typical step size of the E×B staircase is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Investigation of effect of solenoid magnet on emittances of ion beam from laser ablation plasma

A magnetic field can increase an ion current of a laser ablation plasma and is expected to control the change of the plasma ion current. However, the magnetic field can also make some fluctuations of the plasma and the effect on the beam emittance and the emission surface is not clear. To investigate the effect of a magnetic field, we extracted the ion beams under three conditions where without magnetic field, with magnetic field, and without magnetic field with higher laser energy to measure the beam distribution in phase space. Then we compared the relations between the plasma ion current density into the extraction gap and the Twiss parameters with each condition. We observed the effect of the magnetic field on the emission surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Experimental investigation of quasiperiodic-chaotic-quasiperiodic-chaotic transition in a direct current magnetron sputtering plasma

Floating potential fluctuations from a direct current magnetron sputtering plasma have been analysed using time series analysis techniques like phase space plots, power spectra, frequency bifurcation plot, etc. The system exhibits quasiperiodic-chaotic-quasiperiodic-chaotic transitions as the discharge voltage was increased. The transitions of the fluctuations, quantified using the largest Lyapunov exponent, have been corroborated by Hurst exponent and the Shannon entropy. The Shannon entropy is high for quasiperiodic and low for chaotic oscillations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Dirac open-quantum-system dynamics: Formulations and simulations

In this work, we present an open-system interaction formalism for the Dirac equation. Overcoming a complexity bottleneck of alternative formulations, our framework enables efficient numerical simulations (utilizing a typical desktop) of relativistic dynamics within the von Neumann density matrix and Wigner phase-space descriptions. Employing these instruments, we gain important insights into the effect of quantum dephasing for relativistic systems in many branches of physics. In particular, the conditions for robustness of Majorana spinors against dephasing are established. Using the Klein paradox and tunneling as examples, we show that quantum dephasing does not suppress negative energy particle generation. Hence, the Klein dynamics is also robust to dephasing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Simultaneous multiparticle emissions in hot nuclei evaporation process

This work presents a new mechanism for the evaporation with simultaneous particles emission mechanism in the evaporation chain as new channels opened to high excitation energy regime of the compound nucleus. The probability of multiple simultaneous emissions is determined based on phase space approach. A Monte Carlo simulation is employed to compute the final average yield of emitted particles after the decay chain. The neutron, proton, alpha and fission yields are obtained and compared to the conventional calculation with sequential simple particles emission and the relevance of the different channels in competition is also analyzed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

Direct numerical simulation of electrokinetic instability and transition to chaotic motion

A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3 -step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic “room disturbances” and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the experimentally observed absence of regular (near-sinusoidal) oscillations for the overlimiting regimes. There is a qualitative agreement of the experimental and the theoretical values of the threshold of instability, the dominant size of the observed coherent structures, and the experimental and theoretical volt–current characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗