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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

What Causes Radiation Belt Enhancements: A Survey of the Van Allen Probes Era

We survey radiation belt enhancement events during the Van Allen Probes era to determine what mechanism is the dominant cause of enhancements and where it is most effective. Two primary mechanisms have been proposed: (1) betatron/Fermi acceleration due to the Earthward radial transport of electrons, which produces monotonic gradients in phase space density (PSD), and (2) “local acceleration” due to gyro/Landau resonant interaction with electromagnetic waves, which produces radially localized growing peaks in PSD. Here, to differentiate between these processes, we examine radial profiles of PSD in adiabatic coordinates using data from the Van Allen Probes and Time History of Events and Macroscale Interactions during Substorms satellites for 80 outer belt enhancement events from October 2012 to April 2017. Here, this study shows that local acceleration is the dominant acceleration mechanism for MeV electrons in the outer belt, with 87% of the enhancement events exhibiting growing peaks. The strong correlation of the location of these with geomagnetic activity further supports this conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Rayleigh-type parametric chemical oscillation

We consider a nonlinear chemical dynamical system of two phase space variables in a stable steady state. When the system is driven by a time-dependent sinusoidal forcing of a suitable scaling parameter at a frequency twice the output frequency and the strength of perturbation exceeds a threshold, the system undergoes sustained Rayleigh-type periodic oscillation, wellknown for parametric oscillation in pipe organs and distinct from the usual forced quasiperiodic oscillation of a damped nonlinear system where the system is oscillatory even in absence of any external forcing. Our theoretical analysis of the parametric chemical oscillation is corroborated by full numerical simulation of two well known models of chemical dynamics, chlorite-iodine-malonic acid and iodine-clock reactions.

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

Kinetic theory of passing energetic ion transport in presence of the resonant interactions with a rotating magnetic island

The enhanced transport of passing energetic ions (PEIs) in presence of the resonant interactions with a rotating magnetic island is investigated within the drift kinetic framework. When the island rotation plays a role in the resonant interaction, we find that the velocities of PEIs satisfy a constraint relation of resonant flux surface in phase space. The resonant flux surfaces overlap with the magnetic flux surfaces in real space. A new transport channel responsible for the PEIs moving across the magnetic flux surfaces, i.e., continuously overlapping, is found. Two kinds of radial motions can be induced by the surface overlapping: one arises from the coupling between the resonance and the collision with the background plasma and the other from not completely overlapping of the two surfaces. The two radial motions and the symmetry-breaking induced radial motion constitute the total radial motion. When the pitch-angle scattering rate is very weak, the surface-shear induced transport is dominant. Only a small increase in the collision rate can significantly influence the total transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

SU-E-T-795: Validations of Dose Calculation Accuracy of Acuros BV in High-Dose-Rate (HDR) Brachytherapy with a Shielded Cylinder Applicator Using Monte Carlo Simulation

Purpose: Acuros BV has become available to perform accurate dose calculations in high-dose-rate (HDR) brachytherapy with phantom heterogeneity considered by solving the Boltzmann transport equation. In this work, we performed validation studies regarding the dose calculation accuracy of Acuros BV in cases with a shielded cylinder applicator using Monte Carlo (MC) simulations. Methods: Fifteen cases were considered in our studies, covering five different diameters of the applicator and three different shielding degrees. For each case, a digital phantom was created in Varian BrachyVision with the cylinder applicator inserted in the middle of a large water phantom. A treatment plan with eight dwell positions was generated for these fifteen cases. Dose calculations were performed with Acuros BV. We then generated a voxelized phantom of the same geometry, and the materials were modeled according to the vendor’s specifications. MC dose calculations were then performed using our in-house developed fast MC dose engine for HDR brachytherapy (gBMC) on a GPU platform, which is able to simulate both photon transport and electron transport in a voxelized geometry. A phase-space file for the Ir-192 HDR source was used as a source model for MC simulations. Results: Satisfactory agreements between the dose distributions calculated by Acuros BV and those calculated by gBMC were observed in all cases. Quantitatively, we computed point-wise dose difference within the region that receives a dose higher than 10% of the reference dose, defined to be the dose at 5mm outward away from the applicator surface. The mean dose difference was ∼0.45%–0.51% and the 95-percentile maximum difference was ∼1.24%–1.47%. Conclusion: Acuros BV is able to accurately perform dose calculations in HDR brachytherapy with a shielded cylinder applicator.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

RESONANCES, CHAOS, AND SHORT-TERM INTERACTIONS AMONG THE INNER URANIAN SATELLITES

The Portia group of Uranian satellites, representing 9 of the planet’s 13 tiny innermost moons, form a densely packed dynamical system. Hubble Space Telescope observations indicate that their orbits have changed significantly over two decades, and long-term numerical integrations show that their orbits are unstable over millions of years. To investigate the dynamical interactions of the Portia group satellites on the decade timescale over which orbital changes have been observed, we have performed a suite of 100–1000 yr N-body numerical integrations for a range of assumed satellite masses, which are at present not tightly constrained by observations. As first noted by Dawson et al. and recently investigated independently by Quillen and (Robert) French, the moons are configured in chains of interlinked first- and second-order eccentric resonances that contribute to chaotic behavior. We explore in detail several of the strongest of these interlinked resonances. The first such chain is a quintet of orbital resonances: Bianca is near a resonance with Cressida, which is itself near a resonance with Desdemona. Desdemona, in turn, is near a resonance with Portia, which is itself near a resonance with Juliet. The five participating resonances are: Cressida and Bianca (16:15), Desdemona and Cressida (47:46), Portia and Desdemona (13:12), and Portia and Juliet (51:49). A second such chain is a set of two interlinked resonances: Cupid and Belinda (58:57) and Belinda and Perdita (44:43). We also report the new identification of a companion set of second-order inclination-type resonances (Cressida and Bianca (32:30), Desdemona and Cressida (94:92), Portia and Desdemona (26:24), Portia and Juliet (51:49), Cupid and Belinda (116:114), and Belinda and Perdita (88:86)), some of which result in quite strongly coupled variations in the inclinations of the interacting satellites. Using a robust formulation of orbital elements that accounts for the oblateness of Uranus, we probe the dynamical interactions among the moons in the time and frequency domains, and also in phase space, using numerical integrations of subsets of the inner moons, for a range of assumed masses. Several of the satellites are near mean-motion resonance with more than one neighbor, and undergo orbital variations at two nearly equal resonant frequencies. Such configurations of two interlinked resonances can result in chaotic behavior, associated with the transition of one resonant argument from circulation to libration. We demonstrate that, even on the short timescales investigated here, the dynamical interactions, onset of chaos, and associated Lyapunov times are highly sensitive to the masses of the interacting satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Nanoparticles and amyloid systems: A fatal encounter?

Nanoparticles (NPs) are used in many products of our daily life, however, there has been concern that they may also be harmful to human health. Recently NPs have been found to accelerate the fibrillation kinetics of amyloid systems. In the past this has been preliminarily attributed to a nucleation effect. Nanoparticle surfaces and interfaces appear to limit the degrees of freedom of amyloid systems (i.e., peptides and proteins) due to a phase space constraint such that rapid cross-beta structures are formed faster than without interface interactions and in turn fibril formation is enhanced significantly. Here we explore if lipid bilayers in the form of liposomes (140nm) also accelerate fibril formation for amyloid systems. We have investigated a fragment NNFGAIL of the Human islet amyloid polypeptide (hIAPP) in contact with 1,2-diphytanoyl-sn-glycero-3-phosphocholine (DPhPC) liposomes in aqueous solution. We found that the lipid bilayer vesicles do accelerate fibril formation in time-resolved off-line detected atomic force microscopy experiments. Characteristic Thioflavine-T fluorescence on the same structures verify that the structures consist of aggregated peptides in a typical cross-β-structure arrangement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Boltzmann equation solver adapted to emergent chemical non-equilibrium

We present a novel method to solve the spatially homogeneous and isotropic relativistic Boltzmann equation. We employ a basis set of orthogonal polynomials dynamically adapted to allow for emergence of chemical non-equilibrium. Two time dependent parameters characterize the set of orthogonal polynomials, the effective temperature T(t) and phase space occupation factor ϒ(t). In this first paper we address (effectively) massless fermions and derive dynamical equations for T(t) and ϒ(t) such that the zeroth order term of the basis alone captures the particle number density and energy density of each particle distribution. We validate our method and illustrate the reduced computational cost and the ability to easily represent final state chemical non-equilibrium by studying a model problem that is motivated by the physics of the neutrino freeze-out processes in the early Universe, where the essential physical characteristics include reheating from another disappearing particle component (e{sup ±}-annihilation)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On the onset of secondary flow and unsteady solutions through a loosely coiled rectangular duct for large aspect ratio

The study of flows through coiled ducts and channels has attracted considerable attention not only because of their ample applications in Chemical, Mechanical, Civil, Nuclear and Biomechanical engineering but also because of their ample applications in other areas, such as blood flow in the veins and arteries of human and other animals. In this paper, a numerical study is presented for the fully developed two-dimensional flow of viscous incompressible fluid through a loosely coiled rectangular duct of large aspect ratio. Numerical calculations are carried out by using a spectral method, and covering a wide range of the Dean number, Dn, for two types of curvatures of the duct. The main concern of the present study is to find out effects of curvature as well as formation of secondary vortices on unsteady solutions whether the unsteady flow is steady-state, periodic, multi-periodic or chaotic, if Dn is increased. Time evolution calculations as well as their phase spaces are performed with a view to study the non-linear behavior of the unsteady solutions, and it is found that the steady-state flow turns into chaotic flow through various flow instabilities, if Dn is increased no matter what the curvature is. It is found that the unsteady flow is a steady-state solution for small Dn’s and oscillates periodically or non-periodically (chaotic) between two- and twelve-vortex solutions, if Dn is increased. It is also found that the chaotic solution is weak for small Dn’s but strong as Dn becomes large. Axial flow distribution is also investigated and shown in contour plots.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Dynamical modeling of tidal streams

I present a new framework for modeling the dynamics of tidal streams. The framework consists of simple models for the initial action-angle distribution of tidal debris, which can be straightforwardly evolved forward in time. Taking advantage of the essentially one-dimensional nature of tidal streams, the transformation to position-velocity coordinates can be linearized and interpolated near a small number of points along the stream, thus allowing for efficient computations of a stream's properties in observable quantities. I illustrate how to calculate the stream's average location (its 'track') in different coordinate systems, how to quickly estimate the dispersion around its track, and how to draw mock stream data. As a generative model, this framework allows one to compute the full probability distribution function and marginalize over or condition it on certain phase-space dimensions as well as convolve it with observational uncertainties. This will be instrumental in proper data analysis of stream data. In addition to providing a computationally efficient practical tool for modeling the dynamics of tidal streams, the action-angle nature of the framework helps elucidate how the observed width of the stream relates to the velocity dispersion or mass of the progenitor, and how the progenitors of 'orphan' streams could be located. The practical usefulness of the proposed framework crucially depends on the ability to calculate action-angle variables for any orbit in any gravitational potential. A novel method for calculating actions, frequencies, and angles in any static potential using a single orbit integration is described in the Appendix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A Map of the Local Velocity Substructure in the Milky Way Disk

We confirm, quantify, and provide a table of the coherent velocity substructure of the Milky Way disk within 2 kpc of the Sun toward the Galactic anticenter, with a 0.2 kpc resolution. We use the radial velocities of ∼340,000 F-type stars obtained with the Guoshoujing Telescope (also known as the Large Sky Area Multi-Object Fiber Spectroscopic Telescope, LAMOST), and proper motions derived from the PPMXL catalog. The PPMXL proper motions have been corrected to remove systematic errors by subtracting the average proper motions of galaxies and QSOs that have been confirmed in the LAMOST spectroscopic survey, and that are within 2.°5 of the star’s position. We provide the resulting table of systematic offsets derived from the PPMXL proper motion measurements of extragalactic objects identified in the LAMOST spectroscopic survey. Using the corrected phase-space stellar sample, we find statistically significant deviations in the bulk disk velocity of 20 km s{sup −1} or more in the three-dimensional velocities of Galactic disk stars. The bulk velocity varies significantly over length scales of half a kiloparsec or less. The rotation velocity of the disk increases by 20 km s{sup −1} from the Sun’s position to 1.5 kpc outside the solar circle. Disk stars in the second quadrant, within 1 kpc of the Sun, are moving radially toward the Galactic center and vertically toward a point a few tenths of a kiloparsec above the Galactic plane; looking down on the disk, the stars appear to move in a circular streaming motion with a radius of the order of 1 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Parallel Fokker–Planck-DSMC algorithm for rarefied gas flow simulation in complex domains at all Knudsen numbers

A major challenge for the conventional Direct Simulation Monte Carlo (DSMC) technique lies in the fact that its computational cost becomes prohibitive in the near continuum regime, where the Knudsen number (Kn)—characterizing the degree of rarefaction—becomes small. In contrast, the Fokker–Planck (FP) based particle Monte Carlo scheme allows for computationally efficient simulations of rarefied gas flows in the low and intermediate Kn regime. The Fokker–Planck collision operator—instead of performing binary collisions employed by the DSMC method—integrates continuous stochastic processes for the phase space evolution in time. This allows for time step and grid cell sizes larger than the respective collisional scales required by DSMC. Dynamically switching between the FP and the DSMC collision operators in each computational cell is the basis of the combined FP-DSMC method, which has been proven successful in simulating flows covering the whole Kn range. Until recently, this algorithm had only been applied to two-dimensional test cases. In this contribution, we present the first general purpose implementation of the combined FP-DSMC method. Utilizing both shared- and distributed-memory parallelization, this implementation provides the capability for simulations involving many particles and complex geometries by exploiting state of the art computer cluster technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Exploring Microphase Separation in Semi-Fluorinated Diblock Copolymers: A Combined Experimental and Modeling Investigation

We report the combined experimental and theoretical study of the bulk self-assembly behavior of polystyrene-blockpoly( 2,3,4,5,6-pentafluorostyrene) diblock copolymers. These block copolymers were designed to create highly antagonistic blocks (with a high Flory−Huggins interaction parameter, χ) with minimum disruption to the molecular construct (i.e., only replacing five hydrogen atoms with five fluorine atoms). A large library of diblock copolymers (41 samples) was synthesized by reversible addition− fragmentation chain transfer (RAFT) polymerization to map out a major portion of the phase space. All block copolymers exhibited narrow molecular weight distributions with dispersity (D) values between 1.07 and 1.32, and subsequent thermal annealing revealed phase separation into well-defined nanoscale morphologies depending on their molecular composition, as determined from small-angle X-ray scattering and transmission electron microscopy analyses, with an experimental phase diagram being constructed. The χ value at 25 °C for this block copolymer was estimated to be 0.2 using strong segregation theory, based on trends in phase-separated domain spacing and interfacial width. When applying theoretical approaches, the majority of the domain spacing data trends were captured by a coil−coil diblock copolymer model; however, a better fit to the data for samples with shorter fluorinated blocks was obtained with a rod−coil model, indicating that the chains in these fluorinated blocks likely have a higher inherent stiffness and were thus rod-like. This observation demonstrates that, due to the very high value of χ, a transition from coil−coil to rod−coil behavior can be obtained purely by reducing the length of the stiffer of the two blocks and without varying temperature or the chemical composition of the polymers. Here, this work showcases the presence of strong microphase separation within AB diblock copolymers despite the relatively similar chemical composition of the constituent “A” and “B” units, with a clear transition from rod−coil to coil−coil segregation behavior.

RAFT polymerization↗
DOE OSTI2019

Two body relativistic wave equations

The relativistic quantum mechanics of two interacting particles is considered. We first present a covariant formulation of kinematics and of reduced phase space, giving a short outline of the classical results. We then quantize the systems for the scalar–scalar, fermion–scalar and fermion–fermion cases. We study the spectrum and the spherical waves solutions of the free case. The interaction with central scalar and vector potentials is introduced and the explicit equations are deduced. The one particle and the non relativistic limits are recovered and the general lines for the solution of the boundary value problems are given. We make a numerical analysis of the first two cases with Coulomb interaction. For the two fermions we largely revisit the model we had previously derived in order to uniformize the description for all the three cases. In order to give a complete review we report in the Appendix some of the most interesting results obtained for atomic and mesonic systems with Coulomb and Cornell potential interactions respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

The SLUGGS Survey: A Catalog of Over 4000 Globular Cluster Radial Velocities in 27 Nearby Early-type Galaxies

Here, we present positions and radial velocities for over 4000 globular clusters (GCs) in 27 nearby early-type galaxies from the SLUGGS survey. The SLUGGS survey is designed to be representative of elliptical and lenticular galaxies in the stellar mass range 10 < log M{sub ∗}/M {sub ⊙} < 11.7. The data have been obtained over many years, mostly using the very stable multi-object spectrograph DEIMOS on the Keck II 10 m telescope. Radial velocities are measured using the calcium triplet lines, with a velocity accuracy of ±10–15 km s{sup −1}. We use phase space diagrams (i.e., velocity–position diagrams) to identify contaminants such as foreground stars and background galaxies, and to show that the contribution of GCs from neighboring galaxies is generally insignificant. Likely ultra-compact dwarfs are tabulated separately. We find that the mean velocity of the GC system is close to that of the host galaxy systemic velocity, indicating that the GC system is in overall dynamical equilibrium within the galaxy potential. We also find that the GC system velocity dispersion scales with host galaxy stellar mass, in a similar manner to the Faber–Jackson relation for the stellar velocity dispersion. Publication of these GC radial velocity catalogs should enable further studies in many areas, such as GC system substructure, kinematics, and host galaxy mass measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Supermassive black holes in galactic nuclei with tidal disruption of stars. II. Axisymmetric nuclei

The tidal disruption (TD) of stars by supermassive central black holes from dense rotating star clusters is modeled by high-accuracy direct N-body simulations. As in a previous paper on spherical star clusters, we study the time evolution of the stellar tidal disruption rate and the origin of tidally disrupted stars, which are now accorded to several classes of orbits that only occur in axisymmetric systems (short-axis tube and saucer orbits). Compared with that in spherical systems, we found a higher TD rate in axisymmetric systems. The enhancement can be explained by an enlarged loss cone in phase space that stems from the fact that the total angular momentum J is not conserved. As in the case of spherical systems, the distribution of the last apocenter distance of tidally accreted stars peaks at the classical critical radius. However, the angular distribution of the origin of the accreted stars reveals interesting features. Inside the influence radius of the supermassive black hole the angular distribution of disrupted stars has a conspicuous bimodal structure with a local minimum near the equatorial plane. Outside of the influence radius this dependence is weak. We show that the bimodal structure of orbital parameters can be explained by the presence of two families of regular orbits, namely short-axis tube and saucer orbits. Also, we present the consequences of our results for the loss cone in axisymmetric galactic nuclei.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Long-term Dynamical Stability in the Outer Solar System. I. The Regular and Chaotic Evolution of the 34 Largest Trans-Neptunian Objects

We carried out an extensive analysis of the stability of the outer solar system, making use of the frequency analysis technique over short-term integrations of nearly 100,000 test particles, as well as a statistical analysis of 200 1 Gyr long numerical simulations, which consider the mutual perturbations of the giant planets and the 34 largest trans-Neptunian objects (we have called all 34 objects “dwarf planets,” DPs, even if probably only the largest of them are true DPs). From the frequency analysis, we produced statistical diffusion maps for a wide region of the a–e phase-space plane; we also present the average diffusion time for orbits as a function of perihelion. We later turned our attention to the 34 DPs, making an individualized analysis for each of them and producing a first approximation of their future stability. From the 200 distinct realizations of the orbital evolution of the 34 DPs, we classified the sample into three categories, including 17 stable, 11 unstable, and 6 resonant objects; we also found that, statistically, two objects from the sample will leave the trans-Neptunian region within the next gigayear, most likely being ejected from the solar system, but with a nonnegligible probability of going inside the orbit of Neptune, either to collide with a giant planet or even falling to the inner solar system, where our simulations are no longer able to resolve their continuous evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Dynamical transitions in large systems of mean field-coupled Landau-Stuart oscillators: Extensive chaos and cluster states

In this paper, we study dynamical systems in which a large number N of identical Landau-Stuart oscillators are globally coupled via a mean-field. Previously, it has been observed that this type of system can exhibit a variety of different dynamical behaviors. These behaviors include time periodic cluster states in which each oscillator is in one of a small number of groups for which all oscillators in each group have the same state which is different from group to group, as well as a behavior in which all oscillators have different states and the macroscopic dynamics of the mean field is chaotic. We argue that this second type of behavior is “extensive” in the sense that the chaotic attractor in the full phase space of the system has a fractal dimension that scales linearly with N and that the number of positive Lyapunov exponents of the attractor also scales linearly with N. An important focus of this paper is the transition between cluster states and extensive chaos as the system is subjected to slow adiabatic parameter change. We observe discontinuous transitions between the cluster states (which correspond to low dimensional dynamics) and the extensively chaotic states. Furthermore, examining the cluster state, as the system approaches the discontinuous transition to extensive chaos, we find that the oscillator population distribution between the clusters continually evolves so that the cluster state is always marginally stable. This behavior is used to reveal the mechanism of the discontinuous transition. We also apply the Kaplan-Yorke formula to study the fractal structure of the extensively chaotic attractors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

OmegaWINGS: The First Complete Census of Post-starburst Galaxies in Clusters in the Local Universe

Galaxies that abruptly interrupt their star formation in <1.5 Gyr present recognizable features in their spectra (no emission and Hδ in absorption) and are called post-starburst (PSB) galaxies. By studying their stellar population properties and their location within the clusters, we obtain valuable insights on the physical processes responsible for star formation quenching. We present the first complete characterization of PSB galaxies in clusters at 0.04<z<0.07, based on WINGS and OmegaWINGS data, and contrast their properties to those of passive (PAS) and emission-line (EML) galaxies. For V<20, PSBs represent 7.2 ± 0.2% of cluster galaxies within 1.2 virial radii. Their incidence slightly increases from the outskirts toward the cluster center and from the least toward the most luminous and massive clusters, defined in terms of X-ray luminosity and velocity dispersion. The phase-space analysis and velocity-dispersion profile suggest that PSBs represent a combination of galaxies with different accretion histories. Moreover, PSBs with the strongest Hδ are consistent with being recently accreted. PSBs have stellar masses, magnitudes, colors, and morphologies intermediate between PAS and EML galaxies, typical of a population in transition from being star-forming to passive. Comparing the fraction of PSBs to the fraction of galaxies in transition on longer timescales, we estimate that the short-timescale star formation quenching channel contributes two times more than the long timescale one to the growth of the passive population. Processes like ram-pressure stripping and galaxy–galaxy interactions are more efficient than strangulation in affecting star formation.

79 ASTRONOMY AND ASTROPHYSICS↗