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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Theoretical and numerical modelling of chaotic electrostatic ion cyclotron (EIC) oscillations by Jerk equation

In the last few years, third order explicit autonomous differential equations, known as jerk equations, have generated great interest as they show features of regular and chaotic motion. In this paper, we have modelled chaotic electrostatic ion cyclotron oscillations using a third order nonlinear ordinary differential equation (ODE) and investigated its nonlinear dynamical properties. The nonlinear ODE has been derived for a plasma system from a two fluid model in the presence of a source term, under the influence of an external magnetic field, which is perpendicular to the direction of the wave vector. It is seen that the equation does not require an external forcing term to obtain chaotic behaviour. The stability of the solutions of the equation has been investigated analytically as well as numerically, and the bifurcation diagram obtained shows a number of interesting phenomena for various regimes of parameters. The coexisting attractors as well as their corresponding basins are shown and the phase space portraits at different conditions are obtained numerically and shown here. The results obtained here are in agreement with preliminary experiments conducted for a similar configuration of a plasma system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Why is Boris algorithm so good?

Due to its excellent long term accuracy, the Boris algorithm is the de facto standard for advancing a charged particle. Despite its popularity, up to now there has been no convincing explanation why the Boris algorithm has this advantageous feature. In this paper, we provide an answer to this question. We show that the Boris algorithm conserves phase space volume, even though it is not symplectic. The global bound on energy error typically associated with symplectic algorithms still holds for the Boris algorithm, making it an effective algorithm for the multi-scale dynamics of plasmas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Adiabatic electron response and solitary wave generation by trapped particle nonlinearity in a hydrogen plasma

The finite amplitude ion acoustic waves that trap electrons modify the structure of the evolving nonlinear soliton solutions. In the numerical simulations, self-consistently generated solitary waves are studied that emerge as a result of a current driven microinstability growing the ion acoustic mode in a collisionless Vlasov plasma. The growth saturates as a result of nonlinear effects governed by a combination of nonlinearities originating from the hydrodynamic model and kinetic particle trapping effects. The resulting solitary waves also coexist with a finite current and an electron plasma wave capable of perturbing the trapping potential. The results of multiscale simulation are analyzed and characterized following the kinetic prescription of undamped trapped particle mode in the form of phase space vortex solutions that are generalized form of Sagdeev's solitons and obey the solutions of a modified Korteweg-de Vries equation, accounting for a stronger nonlinearity originating from the electron trapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Observing electron spin resonance between 0.1 and 67 GHz at temperatures between 50 mK and 300 K using broadband metallic coplanar waveguides

We describe a fully broadband approach for electron spin resonance (ESR) experiments, where it is possible to tune not only the magnetic field but also the frequency continuously over wide ranges. Here, a metallic coplanar transmission line acts as compact and versatile microwave probe that can easily be implemented in different cryogenic setups. We perform ESR measurements at frequencies between 0.1 and 67 GHz and at temperatures between 50 mK and room temperature. Three different types of samples (Cr{sup 3+} ions in ruby, organic radicals of the nitronyl-nitroxide family, and the doped semiconductor Si:P) represent different possible fields of application for the technique. We demonstrate that an extremely large phase space in temperature, magnetic field, and frequency for ESR measurements, substantially exceeding the range of conventional ESR setups, is accessible with metallic coplanar lines.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Ion beam emittance from an ECRIS

Simulation of ion beam extraction from an Electron Cyclotron Resonance Ion Source (ECRIS) is a fully 3 dimensional problem, even if the extraction geometry has cylindrical symmetry. Because of the strong magnetic flux density, not only the electrons are magnetized but also the Larmor radius of ions is much smaller than the geometrical dimension of the plasma chamber (Ø 64 × 179 mm). If we assume that the influence of collisions is small on the path of particles, we can do particle tracking through the plasma if the initial coordinates of particles are known. We generated starting coordinates of plasma ions by simulation of the plasma electrons, accelerated stochastically by the 14.5 GHz radio frequency power fed to the plasma. With that we were able to investigate the influence of different electron energies on the extracted beam. Using these assumptions, we can reproduce the experimental results obtained 10 years ago, where we monitored the beam profile with the help of viewing targets. Additionally, methods have been developed to investigate arbitrary 2D cuts of the 6D phase space. To this date, we are able to discuss full 4D information. Currently, we extend our analysis tool towards 5D and 6D, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Quasi-classical description of molecular dynamics based on Egorov's theorem

Egorov's theorem on the classical propagation of quantum observables is related to prominent quasi-classical descriptions of quantum molecular dynamics as the linearized semiclassical initial value representation, the Wigner phase space method, or the statistical quasiclassical method. The error estimates show that different accuracies are achievable for the computation of expectation values and position densities. Numerical experiments for a Morse model of diatomic iodine and confined Henon–Heiles systems in various dimensions illustrate the theoretical results.

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

Testing the origin of cosmological magnetic fields through the large-scale structure consistency relations

We study the symmetries of the post-recombination cosmological magnetohydrodynamical equations which describe the evolution of dark matter, baryons and magnetic fields in a self-consistent way. This is done both at the level of fluid equations and of Vlasov-Poisson-Maxwell equations in phase space. We discuss some consistency relations for the soft limit of the (n+1)-correlator functions involving magnetic fields and matter overdensities. In particular, we stress that any violation of such consistency relations at equal-time would point towards an inflationary origin of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Rigorous Two-Step Activation for Fusion Systems with PyNE - Paper 127

In fusion systems, high energy neutrons activate system components. After shutdown, the biological dose rate from photons emitted from these components must be estimated as a function of position and time for maintenance planning. High resolution results can be obtained using the mesh-based Rigorous Two-Step (R2S) method. A neutron transport simulation provides mesh-based neutron fluxes. Using these fluxes and an appropriate irradiation scenario, a nuclear inventory analysis code is used to determine the photon source density within each mesh volume element at any decay time after shutdown. The photon source density mesh is used as a source for photon transport. Photon fluxes are tallied at positions of interest and flux-to-dose rate conversion factors are used to estimate the biological dose rate. Previous work by the authors involved the creation of R2S-ACT, an open source R2S workflow, using the Direct Accelerated Geometry Monte Carlo (DAGMC) version of MCNP5 (DAG-MCNP5) for both transport steps and the ALARA nuclear inventory analysis code. Various use cases may require different combinations of transport and activation codes, and nearly all of the capabilities of R2S-ACT have broader utility. For these reasons, the individual components of R2S-ACT were reimplemented within PyNE. PyNE is a collaborative, trans-institutional, open source project consisting of Python, Cython, and C++ code pertinent to nuclear simulations and analysis. This work involved additions to the PyNE mcnp, mesh, material, alara, dagmc, and sampling modules. A lightweight r2s module was then added to connect the code from these modules. The r2s and sampling modules, together with DAG-MCNP5 and ALARA form a complete R2S workflow for obtaining shutdown dose rate results. Mesh operations, geometry discretization, physics code input/output parsing, and mesh-based source sampling capabilities have been added to PyNE generically, eliminating the need to duplicate this effort across multiple applications. The PyNE r2s module uses these capabilities to couple MCNP and ALARA to form an R2S workflow, which has been validated with the FNG shutdown dose rate benchmark. Future work will involve benchmarking the workflow using tetrahedral mesh and also non-analog sampling options. added so that they can be used interchangeably. This will allow for comparisons between results obtained by different physics codes with minimal human effort. Work will also be done to create and implement an importance sampling method for biasing neutron transport towards regions of phase space of high importance to photon dose rate tallies. This will improve the efficiency of obtaining shutdown dose rate results in large, highly-attenuating geometries. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

Recent results from Super-Kamiokande

The Super-Kamiokande experiment has collected more than 11 live-years of atmospheric neutrino data. Atmospheric neutrinos cover a wide phase space in both energy and distance travelled, the parameters relevant for studying neutrino oscillations. We present here recent measurements of the three-flavor neutrino oscillation parameters using this atmospheric neutrino data, as well as new limits on mixing with a fourth sterile neutrino state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Semiclassical quantization of nonadiabatic systems with hopping periodic orbits

We present a semiclassical quantization condition, i.e., quantum–classical correspondence, for steady states of nonadiabatic systems consisting of fast and slow degrees of freedom (DOFs) by extending Gutzwiller’s trace formula to a nonadiabatic form. The quantum–classical correspondence indicates that a set of primitive hopping periodic orbits, which are invariant under time evolution in the phase space of the slow DOF, should be quantized. The semiclassical quantization is then applied to a simple nonadiabatic model and accurately reproduces exact quantum energy levels. In addition to the semiclassical quantization condition, we also discuss chaotic dynamics involved in the classical limit of nonadiabatic dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A complete low-energy polarized proton-deuteron breakup experiment

An experiment at low energy will be done at the COSY cooler synchrotron and storage ring in order to study the three-particle final states of proton-deuteron scattering reactions measuring a complete set of single and double spin observables over large areas of phase space. The physics objective is to test the predictive power of chiral effective field theory at an energy where convergence is guaranteed and few previous measurements exist. A direct comparison between theoretical predictions and experimental data will be enabled by the use of the so called sampling method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Stellar and gas dynamical model for tidal disruption events in a quiescent galaxy

A detailed model of the tidal disruption events (TDEs) has been constructed using stellar dynamical and gas dynamical inputs that include black hole (BH) mass M{sub •}, specific orbital energy E and angular momentum J, star mass M{sub ⋆} and radius R{sub ⋆,} and the pericenter of the star orbit r{sub p}(E, J, M{sub •}). We solved the steady state Fokker–Planck equation using the standard loss cone theory for the galactic density profile ρ (r) ∝ r{sup −γ} and stellar mass function ξ(m) where m = M{sub ⋆}/M{sub ⊙} and obtained the feeding rate of stars to the BH integrated over the phase space as N-dot {sub t}∝M{sub ∙}{sup β}, where β=−0.3±0.01 for M{sub •} > 10{sup 7}M{sub ⊙} and ∼6.8 × 10{sup −5} Yr{sup −1} for γ = 0.7. We use this to model the in-fall rate of the disrupted debris, M-dot (E,J,m,t), and discuss the conditions for the disk formation, finding that the accretion disk is almost always formed for the fiduciary range of the physical parameters. We also find the conditions under which the disk formed from the tidal debris of a given star with a super Eddington accretion phase. We have simulated the light curve profiles in the relevant optical g band and soft X-rays for both super and sub-Eddington accretion disks as a function of M-dot (E,J,t). Using this, standard cosmological parameters, and mission instrument details, we predict the detectable TDE rates for various forthcoming surveys as a function of γ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Theoretical analysis of the background intensity distribution in X-ray Birefringence Imaging using synchrotron bending-magnet radiation

In the recently developed technique of X-ray Birefringence Imaging, molecular orientational order in anisotropic materials is studied by exploiting the birefringence of linearly polarized X-rays with energy close to an absorption edge of an element in the material. In the experimental setup, a vertically deflecting high-resolution double-crystal monochromator is used upstream from the sample to select the appropriate photon energy, and a horizontally deflecting X-ray polarization analyzer, consisting of a perfect single crystal with a Bragg reflection at Bragg angle of approximately 45°, is placed downstream from the sample to measure the resulting rotation of the X-ray polarization. However, if the experiment is performed on a synchrotron bending-magnet beamline, then the elliptical polarization of the X-rays out of the electron orbit plane affects the shape of the output beam. Also, because the monochromator introduces a correlation between vertical position and photon energy to the X-ray beam, the polarization analyzer does not select the entire beam, but instead selects a diagonal stripe, the slope of which depends on the Bragg angles of the monochromator and the polarization analyzer. In the present work, the final background intensity distribution is calculated analytically because the phase space sampling methods normally used in ray traces are too inefficient for this setup. X-ray Birefringence Imaging data measured at the Diamond Light Source beamline B16 agree well with the theory developed here.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: Adiabatic and non-adiabatic electron-nuclear motion: Quantum and classical dynamics

Using a model for coupled electronic-nuclear motion we investigate the range from negligible to strong non-adiabatic coupling. In the adiabatic case, the quantum dynamics proceeds in a single electronic state, whereas for strong coupling a complete transition between two adiabatic electronic states takes place. It is shown that in all coupling regimes the short-time wave-packet dynamics can be described using ensembles of classical trajectories in the phase space spanned by electronic and nuclear degrees of freedom. We thus provide an example which documents that the quantum concept of non-adiabatic transitions is not necessarily needed if electronic and nuclear motion is treated on the same footing.

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

End-to-end simulation of bunch merging for a muon collider

Muon accelerator beams are commonly produced indirectly through pion decay by interaction of a charged particle beam with a target. Efficient muon capture requires the muons to be first phase-rotated by rf cavities into a train of 21 bunches with much reduced energy spread. Since luminosity is proportional to the square of the number of muons per bunch, it is crucial for a Muon Collider to use relatively few bunches with many muons per bunch. In this paper we will describe a bunch merging scheme that should achieve this goal. We present for the first time a complete end-to-end simulation of a 6D bunch merger for a Muon Collider. The 21 bunches arising from the phase-rotator, after some initial cooling, are merged in longitudinal phase space into seven bunches, which then go through seven paths with different lengths and reach the final collecting "funnel" at the same time. The final single bunch has a transverse and a longitudinal emittance that matches well with the subsequent 6D rectilinear cooling scheme.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

Spin dynamics of large-spin fermions in a harmonic trap

Understanding the collective dynamics in a many-body system has been a central task in condensed matter physics. To achieve this task, we develop a Hartree–Fock theory to study the collective oscillations of spinor Fermi system, motivated by recent experiment on spin-9/2 fermions. We observe an oscillation period shoulder for small rotation angles. Different from previous studies, where the shoulder is found connected to the resonance from periodic to running phase, here the system is always in a running phase in the two-body phase space. This shoulder survives even in the many-body oscillations, which could be tested in the experiments. We also show how these collective oscillations evolve from two- to many-body. Our theory provides an alternative way to understand the collective dynamics in large-spin Fermi systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A Diagnostic Technique for Generalized Synchronization in Systems with a Complex Chaotic Attractor Topology

A diagnostic technique is proposed for the mode of generalized synchronization in systems with a complex chaotic attractor topology based on consideration of tubes of trajectories in the phase space of interacting systems. The method functionality is tested by the numerical simulation of two mutually coupled modified Lorentz systems: one in the chaotic state, and another in the hyperchaotic state. The results are compared to the data obtained by calculating the Lyapunov exponent spectrum and show good agreement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗