Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structural variation”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

The Teton fault, Wyoming: Topographic signature, neotectonics, and mechanisms of deformation

We integrated geophysical and geological methods to evalute the structural evolution of the active Teton normal fault, Wyoming, and its role in the development of the dramatic topography of Teton Range and Jackson Hole. Comparison of variations in surface offsets with the topographic expression of the Teton range crest and drainage divide, and the overall structure of the range, suggests that the effects ofpostglacial faulting cannot be discriminated from the influence of pre-extensional structures and differential; erosion on the footwall topography. In contrast, the effects of multiple scarp-forming normal faulting earthquakes are expressed by the anomalous drainage pattern and westward tilt of the hanging wall, Jackson Hole, toward the Teton fault. Kinematic boundary element fault models suggest that the westward tilt of the valley floor is the product of 110-125 m of displacement on a 45 deg-75 deg E dipping Teton fault in the past 25,000-75,000 years. Comparisons with historic normal faulting earthquake displacements imply that this range of displacement corresponds to 10-50, M greater than 7 scarp-forming earthquakes. A total throw of 2.5 to 3.5 km across the Teton fault is suggested by inverse ray-tracing and forward gravity models. These models also suggest that Laramide age structures have been offset across the Teton fault and obscure its geophysical signature but also continue to influence the structural and topographic expression of the footwall and hanging wall blocks. Paleomagnetic analyses of the approximately 2.0 Ma Huckelberry Ridge Tuff suggest that the overall westward tilt of the Teton Range is a result ofabout 10 deg of west side down tilt across the Teton fault since tuff emplacement. This suggests that much if not all of the throw across the Teton fault has accumulated in the past 2 m.y. Complex demagnetization and rock magnetic behavior and local emplacement of the Huckleberry Ridge Tuff on preexisting topogrpahy preclude determination of the amount or variations in throw along strike of the Teton fault from the paleomagnetic data.

Byrd, John O. D.↗
NASA NTRS1985

Effects of stratospheric lapse rate on numerically simulated thunderstorm cloud top structure

The study of overshooting midlatitude severe thunderstorm tops has led to the detection of some intriguing patterns. This paper provides a summary of an investigation of the effects of stratospheric stability on thunderstorm cloud top structure, taking into account the utilization of a fully three-dimensional (3D) anelastic model discussed by Schlesinger (1975, 1978, 1984). The layout of model experiments is discussed, giving attention to initial environmental soundings, and initial forcing. The time variation of the vertical velocity extreme is examined, and a description of the mature storm structure is presented. Constant-altitude plane views throughout storm depth are considered along with cloud top height fields. With the stratospheric lapse, the warm region is small and breaks down into separate spots with secondary cold spots downshear of them.

Schlesinger, R. E.↗
NASA NTRS2004

Small-scale Interstellar Structure Toward the Open Cluster Chi Persei; Program ID: BO29 (Cycle 2)

The primary objective of this program was to obtain FUSE observations of the interstellar H2 absorption toward six early-type stars in the core of the open star cluster Chi Persei. High resolution optical observations of the interstellar Na I absorption toward these stars and others in the core of Chi Per have revealed a rich variety of small-scale diffuse ISM structure in both the distant Perseus and more local Orion spiral arms. At the 1 arc minute angular scales typically probed by the separations of these stars, this structure corresponds to respective physical length scales of 0.6 and 0.15 pc in the Perseus and Orion arm gas. The principal scientific goal of the FUSE observations was to determine if this small-scale Na I structure toward the core of Chi Per was reflected in the abundance and/or excitation of interstellar H2. In September 2001, we obtained FUSE observations of the Chi Per stars BD +56 563, BD +56 571, BD +56 573, BD +56 574, BD +56 575, and BD +56 578 for this program. The data reveal that the Perseus arm gas exhibits significant small-scale variations in the strength of its molecular hydrogen absorption while the Orion arm H2 absorption is relatively uniform. We presented these results at the January 2003 American Astronomical Society meeting in Seattle, WA (points, Meyer, Lauroesch, Nguyen, and Friedman 2002, Bulletin of the American Astronomical Society, 34,1229). Since Na I is not a dominant ion in H I clouds, the variations detected in the Perseus arm Na I line profiles could be due to spatial differences in either the physical structure or the physical conditions in the gas. Although many of the H2 absorption profiles are too saturated or too weak to be definitive in this regard, it is clear that small-scale variations are present in H2 gas with excitation as high as the J=4 rotational level. In contrast, lower spatial resolution H I 21 cm studies of the Perseus arm gas across the face of the h and Chi Per double cluster have shown smoother, continuous variations on the larger scale. In order to investigate if the global velocity structure of H2 toward h and Chi Per more closely follows the chaotic structure seen in interstellar Na I absorption or the smoother, continuous H I structure, we have expanded our FUSE study with a Cycle 4 program (D902, PI: Points) to survey the H2 absorption toward a number of other stars across the face of the double cluster. We have obtained FUSE spectra of 14 additional h and Chi Per stars through this program and are now combining all of the data into a comprehensive study that will be submitted for publication next year.

Meyer, David M.↗
DOE OSTI2022

Quantum algorithm for electronic band structures with local tight-binding orbitals

While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum “toolbox”. With this in mind, we consider the quantum calculation of a periodic system’s single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Milky Way tomography with APOGEE: intrinsic density distribution and structure of mono-abundance populations

ABSTRACT The spatial distribution of mono-abundance populations (MAPs, selected in [Fe/H] and [Mg/Fe]) reflect the chemical and structural evolution in a galaxy and impose strong constraints on galaxy formation models. In this paper, we use APOGEE data to derive the intrinsic density distribution of MAPs in the Milky Way, after carefully considering the survey selection function. We find that a single exponential profile is not a sufficient description of the Milky Way’s disc. Both the individual MAPs and the integrated disc exhibit a broken radial density distribution; densities are relatively constant with radius in the inner Galaxy and rapidly decrease beyond the break radius. We fit the intrinsic density distribution as a function of radius and vertical height with a 2D density model that considers both a broken radial profile and radial variation of scale height (i.e. flaring). There is a large variety of structural parameters between different MAPs, indicative of strong structure evolution of the Milky Way. One surprising result is that high-α MAPs show the strongest flaring. The young, solar-abundance MAPs present the shortest scale height and least flaring, suggesting recent and ongoing star formation confined to the disc plane. Finally we derive the intrinsic density distribution and corresponding structural parameters of the chemically defined thin and thick discs. The chemical thick and thin discs have local surface mass densities of 5.62 ± 0.08 and 15.69 ± 0.32 M⊙pc−2, respectively, suggesting a massive thick disc with a local surface mass density ratio between thick to thin disc of 36 per cent.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2022

Validation and Interpretation of a Three-dimensional Configuration of a Magnetic Cloud Flux Rope

One strong magnetic cloud (MC)with a magnetic field magnitude reaching∼40 nT at 1 au during 2012 June16–17 is examined in association with a preexisting magnetic flux rope (MFR)identified on the Sun. The MC is characterized by a quasi-three-dimensional(3D) flux rope model based on in situ measurements from the Wind spacecraft. The contents of the magnetic flux and other parameters are quantified. In addition, a correlative study with the corresponding measurements of the same structure crossed by the Venus Express (VEX)spacecraft at a heliocentric distance of 0.7 au and with an angular separation of∼6°in longitude is performed to validate the MC modeling results. The spatial variation between the Wind and VEX magnetic field measurements is attributed to the3D configuration of the structure appearing as a knotted bundle of flux. A comparison of the magnetic flux contents between the MC and the preexisting MFR on the Sun indicates that the 3D reconnection process accompanying anM1.9flare may correspond to the magnetic reconnection between the field lines of the preexisting MFR rooted in the opposite polarity foot points. Such a process reduces the amount of the axial magnetic flux in the erupted flux rope, by approximately 50%, in this case.

Qiang Hu↗
DOE OSTI2025

Correlating energy density induced residual stress, porosity, and mechanical property variations in directed energy deposition using neutron diffraction and imaging techniques

In a directed energy deposition (DED) process, the energy density deposited into the materials drives the structure and properties. The energy density is influenced by several printing parameters, such as laser power, hatch spacing, scanning speed, etc. In this study, 316L stainless steels samples were fabricated with varying hatch spacing aiming to induce variations in sample microstructure and properties. A combination of techniques including engineering neutron diffraction, neutron computed tomography, and tensile testing, were employed to correlate the structure and macroscopic properties. Porosity was predominantly observed at the bottom of the deposited materials and was effectively reduced by decreased hatch spacing. Correlatively, compressive residual stress was observed at the bottom of the specimens, while internal stress is largely determined by the hatch spacing. The residual stress in general decreases as the hatch spacing decreases. The maximum ultimate tensile strength (UTS) was found to increase from 517.4 to 528.6 MPa with a strain increased from 0.847 to 0.897 as the hatch spacing decreased from 0.45 to 0.42 mm. However, further increase in energy density by reducing the hatch spacing resulted in a significant decrease in UTS (as low as 452 MPa). This work provides new insights from bulk non-destructive techniques into correlating energy density induced residual stress, porosity, and mechanical property variations in directed energy deposition.

316L stainless steel↗
NASA NTRS1991

Mantle viscosity - A comparison of models from postglacial rebound and from the geoid, plate driving forces, and advected heat flux

Models of the radial variation of effective viscosity inferred from the earth's response to surface loads associated with Pleistocene deglaciation are compared to structures inferred from models of geodynamic phenomena associated with convection: the geoid, plate-driving forces, and advected heat flux. While observations of the earth's response to surface loads do not have sufficient resolution to justify more than two viscous layers, adequately matching the observed long-wavelength geoid anomalies associated with density contrasts in the lower mantle (inferred from seismic tomography) and in the upper mantle (inferred from a model of subducted slabs) requires more structure. It is possible to explain the geoid, observed plate velocities, the advected heat flux in the lower mantle, and relative sea-level variations in oceanic regions, all with a mantle with a high-viscosity/elastic lid, an asthenospheric channel of 2 x 10 exp 19 Pa s from 100 to 400-km depth, a 6 x 10 exp 20 Pa s transition zone, and a lower mantle of 6 x 10 exp 21 Pa s. The uplift history of Australia, Fennoscandia, and Laurentia can be explained with an asthenospheric viscosity less than a factor of 10 higher. Lateral variations in lower mantle viscosity are not required. Transient creep appears to be unimportant for the recent response-to-surface loads from Pleistocene deglaciation.

Hager, Bradford H.↗
DOE OSTI2022

Tuning Metal–Organic Framework (MOF) Topology by Regulating Ligand and Secondary Building Unit (SBU) Geometry: Structures Built on 8-Connected M 6 (M = Zr, Y) Clusters and a Flexible Tetracarboxylate for Propane-Selective Propane/Propylene Separation

Topology evolution originating from variations of linker and SBU (Secondary Building Unit) geometries could largely enrich the chemistry of metal–organic frameworks (MOFs). Here we report the synthesis and characterization of three MOF structures built on the same organic linker, N,N,N',N'-Tetrakis(4-carboxyphenyl)-1,4-phenylenediamine (tcppda) and similar 8-connected M 6 (M = Zr or Y) clusters. The three compounds, HIAM-402, HIAM-403, and HIAM-311, feature 4,8-connected sqc, scu, and flu topology, respectively. Detailed structural analysis revealed that different geometries of the inorganic M 6 SBUs as well as the organic linker have led to the formation of distinct MOF nets. In particular, HIAM-402 features exceptional framework stability and high porosity and acts as a propane-selective adsorbent for the discrimination of propane and propylene. Its balanced adsorption selectivity (S propane/propylene = 1.43) and capacity (Q propane = 133.3 cm 3 /g, 298 K and 1 bar) endow it with the capability of separating propane and propylene mixtures and one-step production of highly pure propylene (purity >99.9%), as validated by column breakthrough measurements, with the presence of moisture or propyne. As a result, ab initio calculations further confirm that the propane-selective behavior of HIAM-402 is a result of its higher binding energy toward propane compared to that of propylene.

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

SCOMAP-XD : atomistic deuterium contrast matching for small-angle neutron scattering in biology

The contrast-variation method in small-angle neutron scattering (SANS) is a uniquely powerful technique for determining the structure of individual components in biomolecular systems containing regions of different neutron scattering length density ρ . By altering the ρ of the target solute and the solvent through judicious incorporation of deuterium, the scattering of desired solute features can be highlighted. Most contrast-variation methods focus on highlighting specific bulk solute elements, but not on how the scattering at specific scattering vectors q , which are associated with specific structural distances, changes with contrast. Indeed, many systems exhibit q -dependent contrast effects. Here, a method is presented for calculating both bulk contrast-match points and q -dependent contrast using 3D models with explicit solute and solvent atoms and SASSENA , an explicit-atom SANS calculator. The method calculates the bulk contrast-match points within 2.4% solvent D 2 O accuracy for test protein–nucleic acid and lipid nanodisc systems. The method incorporates a general model for the incorporation of deuterium at non-exchangeable sites that was derived by performing mass spectrometry on green fluorescent protein. The method also decomposes the scattering profile into its component parts and identifies structural features that change with contrast. The method is readily applicable to a variety of systems, will expand the understanding of q -dependent contrast matching and will aid in the optimization of next-generation neutron scattering experiments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Evolutionary lineage explains trait variation among 75 coexisting grass species

Summary Evolutionary history plays a key role driving patterns of trait variation across plant species. For scaling and modeling purposes, grass species are typically organized into C 3 vs C 4 plant functional types (PFTs). Plant functional type groupings may obscure important functional differences among species. Rather, grouping grasses by evolutionary lineage may better represent grass functional diversity. We measured 11 structural and physiological traits in situ from 75 grass species within the North American tallgrass prairie. We tested whether traits differed significantly among photosynthetic pathways or lineages (tribe) in annual and perennial grass species. Critically, we found evidence that grass traits varied among lineages, including independent origins of C 4 photosynthesis. Using a rigorous model selection approach, tribe was included in the top models for five of nine traits for perennial species. Tribes were separable in a multivariate and phylogenetically controlled analysis of traits, owing to coordination of important structural and ecophysiological characteristics. Our findings suggest grouping grass species by photosynthetic pathway overlooks variation in several functional traits, particularly for C 4 species. These results indicate that further assessment of lineage‐based differences at other sites and across other grass species distributions may improve representation of C 4 species in trait comparison analyses and modeling investigations.

Donnelly, Ryan C.↗
NASA NTRS2004

Density, Velocity and Ionization Structure in Accretion-Disc Winds

This was a project to exploit the unique capabilities of FUSE to monitor variations in the wind- formed spectral lines of the luminous, low-inclination, cataclysmic variables(CV) -- RW Sex. (The original proposal contained two additional objects but these were not approved.) These observations were intended to allow us to determine the relative roles of density and ionization state changes in the outflow and to search for spectroscopic signatures of stochastic small-scale structure and shocked gas. By monitoring the temporal behavior of blue-ward extended absorption lines with a wide range of ionization potentials and excitation energies, we proposed to track the changing physical conditions in the outflow. We planned to use a new Monte Carlo code to calculate the ionization structure of and radiative transfer through the CV wind. The analysis therefore was intended to establish the wind geometry, kinematics and ionization state, both in a time-averaged sense and as a function of time.

Sonneborn, George↗
NASA NTRS2009

Structured DC Electric Fields With and Without Associated Plasma Density Gradients Observed with the C/NOFS Satellite

DC electric field observations and associated plasma drifts gathered with the Vector Electric Field Investigation on the Air Force Communication/Navigation Outage Forecasting System (C/NOFS) satellite typically reveal considerable variation at large scales (approximately 100's of km), in both daytime and nighttime cases, with enhanced structures usually confined to the nightside. Although such electric field structures are typically associated with plasma density depletions and structures, as observed by the Planar Langmuir Probe on C/NOFS, what is surprising is the number of cases in which large amplitude, structured DC electric fields are observed without a significant plasma density counterpart structure, including their appearance at times when the ambient plasma density appears relatively quiescent. We investigate the relationship of such structured DC electric fields and the ambient plasma density in the C/NOFS satellite measurements observed thus far, taking into account both plasma density depletions and enhancements. We investigate the mapping of the electric fields along magnetic field lines from distant altitudes and latitudes to locations where the density structures, which presumably formed the original seat of the electric fields, are no longer discernible in the observations. In some cases, the electric field structures and spectral characteristics appear to mimic those associated with equatorial spread-F processes, providing important clues to their origins. We examine altitude, seasonal, and longitudinal effects in an effort to establish the origin of such structured DC electric fields observed both with, and without, associated plasma density gradients

Pfaff, R.↗
NASA NTRS1983

Saturn's E, G, and F rings - Modulated by the plasma sheet?

Saturn's broad E ring, the narrow G ring, and the structured and apparently time-variable F ring(s) contain many micron and submicron-sized particles, which make up the 'visible' component. These rings (or ring systems) are in direct contact with magnetospheric plasma. Fluctuations in the plasma density and/or mean energy, due to magnetospheric and solar wind processes, may induce stochastic charge variations on the dust particles, which in turn lead to an orbit perturbation and spatial diffusion. In addition, Coulomb drag forces may be important, in particular for the E ring. The possibility that electromagnetic effects may play a role in determining the F ring structure and its possible time variations is critically examined. Sputtering of micron-sized dust particles in the E ring by magnetospheric ions yields lifetimes of 100 to 10,000 years. This effect as well as the plasma induced transport processes require an active source for the E ring, probably Enceladus.

Morfill, G. E.↗
NASA NTRS2008

Characterizing ISS Charging Environments with On-Board Ionospheric Plasma Measurements

Charging of the International Space Station (ISS) is dominated by interactions of the biased United States (US) 160 volt solar arrays with the relatively high density, low temperature plasma environment in low Earth orbit. Conducting surfaces on the vehicle structure charge negative relative to the ambient plasma environment because ISS structure is grounded to the negative end of the US solar arrays. Transient charging peaks reaching potentials of some tens of volts negative controlled by photovoltaic array current collection typically occur at orbital sunrise and sunset as well as near orbital noon. In addition, surface potentials across the vehicle structure vary due to an induced v x B (dot) L voltage generated by the high speed motion of the conducting structure across the Earth's magnetic field. Induced voltages in low Earth orbit are typically only approx.0.4 volts/meter but the approx.100 meter scale dimensions of the ISS yield maximum induced potential variations ofapprox.40 volts across the vehicle. Induced voltages are variable due to the orientation of the vehicle structure and orbital velocity vector with respect to the orientation of the Earth's magnetic field along the ISS orbit. In order to address the need to better understand the ISS spacecraft potential and plasma environments, NASA funded development and construction of the Floating Potential Measurement Unit (FPMU) which was deployed on an ISS starboard truss arm in August 2006. The suite of FPMU instruments includes two Langmuir probes, a plasma impedance probe, and a potential probe for use in in-situ monitoring of electron temperatures and densities and the vehicle potential relative to the plasma environment. This presentation will describe the use of the FPMU to better characterize interactions of the ISS with the space environment, changes in ISS charging as the vehicle configuration is modified during ISS construction, and contributions of FPMU vehicle potential and plasma environment measurements to investigations of on-orbit anomalies in ISS systems.

Minow, Jospeh I.↗
DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗
NASA NTRS1987

Overview of NASA PTA propfan flight test program

During the last several years high-speed propellers have made the transition from a wind tunnel curiosity to a very likely near-term, fuel-efficient propulsion system that could revolutionize the subsonic commercial air transport industry. A key ingredient in this remarkable progress is the advanced turboprop industry. Working together, NASA and industry have developed and flight tested two propeller propulsion systems to provide answers to key technical questions and concerns. An industry team is currently developing a third propeller propulsion system for flight testing late this year. The progress of one of the NASA-industry flight test programs, called the Propfan Test Assessment (PTA) Program is reported. In PTA, a 9 foot diameter propfan was installed on the left wing of a Gulfstream GII executive jet and is undergoing extensive flight testing at Dobbins Air Force Base to evaluate propfan structural integrity, near and far field noise, and cabin interior noise characteristics. This research testing includes variations in propeller tip speed and power loading, nacelle tilt angle, and aircraft Mach number and altitude. As a result, extensive parametric data will be obtained to verify and improve computer codes for predicting propeller structural aeroelastic, aerodynamic, and acoustic characteristics. Over 600 measurements are being recorded for each of approximately 600 flight test conditions.

Graber, Edwin J.↗
DOE OSTI2023

ESCARGOT: Mapping Vertical Phase Spiral Characteristics Throughout the Real and Simulated Milky Way

Abstract The recent discovery of a spiral pattern in the vertical kinematic structure in the solar neighborhood provides a prime opportunity to study nonequilibrium dynamics in the Milky Way from local stellar kinematics. Furthermore, results from simulations indicate that even in a limited volume, differences in stellar orbital histories allow us to trace variations in the initial perturbation across large regions of the disk. We present ESCARGOT , a novel algorithm for studying these variations in both simulated and observed data sets. ESCARGOT automatically extracts key quantities from the structure of a given phase spiral, including the time since perturbation and the perturbation mode. We test ESCARGOT on simulated data and show that it is capable of accurately recovering information about the time since the perturbation occurred as well as subtle differences in phase spiral morphology due to stellar locations in the disk at the time of perturbation. We apply ESCARGOT to kinematic data from data release 3 of the Gaia mission in bins of guiding radius. We show that similar structural differences in morphology occur in the Gaia phase spirals as a function of stellar orbital history. These results indicate that the phase spirals are the product of a complex dynamical response in the disk with large-scale coupling between different regions of phase space.

Darragh-Ford, Elise (ORCID:0000000288005652)↗