Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PHASE VELOCITY”

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,747 records · Page 97

Self-propagating reactions in Al/Zr multilayers: Anomalous dependence of reaction velocity on bilayer thickness

High temperature, self-propagating reactions are observed in vapor-deposited Al/Zr multilayered foils of overall atomic ratios 3 Al:1 Zr and 2 Al:1 Zr and nanoscale layer thicknesses; however, the reaction velocities do not exhibit the inverse dependence on bilayer thickness that is expected based on changes in the average diffusion distance. Instead, for bilayer thicknesses of 20-30 nm, the velocity is essentially constant at ∼7.7 m/s. We explore several possible explanations for this anomalous behavior, including microstructural factors, changes in the phase evolution, and phase transformations in the reactant layers, but find no conclusive explanations. We determine that the phase evolution during self-propagating reactions in foils with a 3 Al:1 Zr stoichiometry is a rapid transformation from Al/Zr multilayers to the equilibrium intermetallic Al{sub 3}Zr compound with no intermediate crystalline phases. This phase evolution is the same for foils of 90 nm bilayer thicknesses and foils of bilayer thicknesses in the range of 27 nm to 35 nm. Further, for foils with a bilayer thickness of 90 nm and a 3 Al:1 Zr overall chemistry, the propagation front is planar and steady, in contrast to unsteady reaction fronts in foils with 1 Al:1 Zr overall chemistry and similar bilayer thicknesses.

36 MATERIALS SCIENCE↗

Effects of fog droplets on wake vortex decay rate

A simple model for the motion of particles in a laminar line vortex is discussed. The energy required to accelerate a set of these particles was determined and shown to be only a small fraction of the energy content of the vortex flow. It is shown that this energy transfer is unlikely to be sufficient to significantly modify the vortex decay rate. It is further argued that the effect of the particle on the viscous properties of the resulting two phase fluid leads to a slower decay rate than in single phase air flow. However, this conclusion may not necessarily follow for turbulence flows. Results show that the migration of particles to the outer flow results in a redistribution of the velocity profile in the vortex and in a non-uniform two phase viscosity across the core. It is suggested that these effects may accelerate vortex bursting.

Moulden, T. H.↗

Uplink Array Concept Demonstration with the EPOXI Spacecraft

Uplink array technology is currently being developed for NASA's Deep Space Network (DSN) to provide greater range and data throughput for future NASA missions, including manned missions to Mars and exploratory missions to the outer planets, the Kuiper belt, and beyond. Here we describe a novel technique for generating the frequency predicts that are used to compensate for relative Doppler, derived from interpolated earth position and spacecraft ephemerides. The method described here guarantees velocity and range estimates that are consistent with each other, hence one can always be recovered from the other. Experimental results have recently proven that these frequency predicts are accurate enough to maintain the phase of a three element array at the EPOXI spacecraft for three hours. Previous methods derive frequency predicts directly from interpolated relative velocities. However, these velocities were found to be inconsistent with the corresponding spacecraft range, meaning that range could not always be recovered accurately from the velocity predicts, and vice versa. Nevertheless, velocity-based predicts are also capable of maintaining uplink array phase calibration for extended periods, as demonstrated with the EPOXI spacecraft, however with these predicts important range and phase information may be lost. A comparison of the steering-vector method with velocity-based techniques for generating precise frequency predicts specifically for uplink array applications is provided in the following sections.

spacecraft uplink commands↗

Pointing-Vector and Velocity Based Frequency Predicts for Deep-Space Uplink Array Applications

Uplink array technology is currently being developed for NASA's Deep Space Network (DSN) to provide greater range and data throughput for future NASA missions, including manned missions to Mars and exploratory missions to the outer planets, the Kuiper belt, and beyond. Here we describe a novel technique for generating the frequency predicts that are used to compensate for relative Doppler, derived from interpolated earth position and spacecraft ephemerides. The method described here guarantees velocity and range estimates that are consistent with each other, hence one can always be recovered from the other. Experimental results have recently proven that these frequency predicts are accurate enough to maintain the phase of a three element array at the EPOXI spacecraft for three hours. Previous methods derive frequency predicts directly from interpolated relative velocities. However, these velocities were found to be inconsistent with the corresponding spacecraft range, meaning that range could not always be recovered accurately from the velocity predicts, and vice versa. Nevertheless, velocity-based predicts are also capable of maintaining uplink array phase calibration for extended periods, as demonstrated with the EPOXI spacecraft, however with these predicts important range and phase information may be lost. A comparison of the steering-vector method with velocity-based techniques for generating precise frequency predicts specifically for uplink array applications is provided in the following sections.

uplink arraying↗

Investigation of the feasibility of optical diagnostic measurements at the exit of the SSME

Under Contract NAS8-36861 sponsored by NASA Marshall Space Flight Center, the United Technologies Research Center is conducting an investigation of the feasibility of remote optical diagnostics to measure temperature, species concentration and velocity at the exit of the Space Shuttle Main Engine (SSME). This is a two phase study consisting of a conceptual design phase followed by a laboratory experimental investigation. The first task of the conceptual design studies is to screen and evaluate the techniques which can be used for the measurements. The second task is to select the most promising technique or techniques, if as expected, more than one type of measurement must be used to measure all the flow variables of interest. The third task is to examine in detail analytically the capabilities and limitations of the selected technique(s). The results of this study are described in the section of this report entitled Conceptual Design Investigations. The conceptual design studies identified spontaneous Raman scattering and photodissociative flow-tagging for measurements respectively of gas temperature and major species concentration and for velocity. These techniques and others that were considered are described in the section describing the conceptual design. The objective of the second phase of investigations was to investigate experimentally the techniques identified in the first phase. The first task of the experimental feasibility study is to design and assemble laboratory scale experimental apparatus to evaluate the best approaches for SSME exit optical diagnostics for temperature, species concentrations and velocity, as selected in the Phase I conceptual design study. The second task is to evaluate performance, investigate limitations, and establish actual diagnostic capabilities, accuracies and precision for the selected optical systems. The third task is to evaluate design requirements and system trade-offs of conceptual instruments. Spontaneous Raman scattering excited by a KrF excimer laser pulse was investigated for SSME exit plane temperature and major species concentration measurements. The relative concentrations of molecular hydrogen and water vapor would be determined by measuring the integrated Q-branch scattering signals through narrow bandpass filters in front of photomultipliers. The temperature would be determined by comparing the signal from a single hydrogen rotational Raman line to the total hydrogen Q-branch signal. The rotational Raman line would be isolated by a monochromator and detected with a PMT.

Shirley, John A.↗

An Advanced N -body Model for Interacting Multiple Stellar Systems

We construct an advanced model for interacting multiple stellar systems in which we compute all trajectories with a numerical N -body integrator, namely the Bulirsch–Stoer from the SWIFT package. We can then derive various observables: astrometric positions, radial velocities, minima timings (TTVs), eclipse durations, interferometric visibilities, closure phases, synthetic spectra, spectral energy distribution, and even complete light curves. We use a modified version of the Wilson–Devinney code for the latter, in which the instantaneous true phase and inclination of the eclipsing binary are governed by the N -body integration. If all of these types of observations are at one’s disposal, a joint χ {sup 2} metric and an optimization algorithm (a simplex or simulated annealing) allow one to search for a global minimum and construct very robust models of stellar systems. At the same time, our N -body model is free from artifacts that may arise if mutual gravitational interactions among all components are not self-consistently accounted for. Finally, we present a number of examples showing dynamical effects that can be studied with our code and we discuss how systematic errors may affect the results (and how to prevent this from happening).

79 ASTRONOMY AND ASTROPHYSICS↗

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis of toroidal phasing of resonant magnetic perturbation effects on edge transport in the DIII-D tokamak

Resonant Magnetic Perturbation (RMP) fields produced by external control coils are considered a viable option for the suppression of Edge Localized Modes in present and future tokamaks. In DIII-D, the RMPs are generated by six pairs of I-coils, each spanning 60° in toroidal angle, with the currents flowing in opposite directions in adjacent pairs of I-coils. Reversal of the currents in all I-coils, which produces a 60° toroidal shift in the RMP field configuration, generates uniquely different edge pedestal profiles of the density, temperature, and rotation velocities, implying different effects on the related edge transport phenomena caused by the difference in toroidal phase of the I-coil currents. The diffusive and non-diffusive transport effects of this RMP toroidal phase reversal are analyzed by comparing the ion and electron heat diffusivities, angular momentum transport frequencies, ion diffusion coefficients, and the particle pinch velocities interpreted from the measured profiles for the two phases of the I-coil currents.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Substructure in the stellar halo near the Sun: I. Data-driven clustering in integrals-of-motion space

Context. Merger debris is expected to populate the stellar haloes of galaxies. In the case of the Milky Way, this debris should be apparent as clumps in a space defined by the orbital integrals of motion of the stars. Aims. Our aim is to develop a data-driven and statistics-based method for finding these clumps in integrals-of-motion space for nearby halo stars and to evaluate their significance robustly. Methods. We used data from Gaia EDR3, extended with radial velocities from ground-based spectroscopic surveys, to construct a sample of halo stars within 2.5 kpc from the Sun. We applied a hierarchical clustering method that makes exhaustive use of the single linkage algorithm in three-dimensional space defined by the commonly used integrals of motion energy E, together with two components of the angular momentum, L z and L ⊥ . To evaluate the statistical significance of the clusters, we compared the density within an ellipsoidal region centred on the cluster to that of random sets with similar global dynamical properties. By selecting the signal at the location of their maximum statistical significance in the hierarchical tree, we extracted a set of significant unique clusters. By describing these clusters with ellipsoids, we estimated the proximity of a star to the cluster centre using the Mahalanobis distance. Additionally, we applied the HDBSCAN clustering algorithm in velocity space to each cluster to extract subgroups representing debris with different orbital phases. Results. Our procedure identifies 67 highly significant clusters (> 3σ), containing 12% of the sources in our halo set, and 232 subgroups or individual streams in velocity space. In total, 13.8% of the stars in our data set can be confidently associated with a significant cluster based on their Mahalanobis distance. Inspection of the hierarchical tree describing our data set reveals a complex web of relations between the significant clusters, suggesting that they can be tentatively grouped into at least six main large structures, many of which can be associated with previously identified halo substructures, and a number of independent substructures. This preliminary conclusion is further explored in a companion paper, in which we also characterise the substructures in terms of their stellar populations. Conclusions. Our method allows us to systematically detect kinematic substructures in the Galactic stellar halo with a data-driven and interpretable algorithm. The list of the clusters and the associated star catalogue are provided in two tables available at the CDS.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of H2O on Upper Mantle Phase Transitions in MgSiO3: is the Seismic X-discontinuity an Indicator of Mantle Water Content

The mantle X-discontinuity, usually assigned to positive seismic velocity reflectors in the 260-330 km depth range, has proved difficult to explain in terms of a single mineralogical phase transformation in part because of its depth variability. The coesite to stishovite transition of SiO{sub 2} matches deeper X-discontinuity depths but requires 5-10% free silica in the mantle to match observed impedance contrast. The orthoenstatite (OEn) to high-pressure clinoenstatite (HPCen) transformation of MgSiO{sub 3} also broadly coincides with depths of the X but requires chemically depleted and orthoenstatite-rich lithology at 300 km depth in order to match observed seismic impedance contrast. On the basis of high-pressure infrared spectroscopy, X-ray diffraction, and Raman spectroscopy, we show that 1300 ppm variation of H{sub 2}O content in MgSiO{sub 3} can displace the transition of low-pressure clinoenstatite (LPCen) to HPCen by up to 2 GPa, similar to previous quench experiments on the OEn to HPCen phase transition, where about 30-45 km (1.0-1.5 GPa) of deflection could occur per 0.1 wt% H{sub 2}O. If the mantle X-discontinuity results from pyroxene transitions in a depleted harzburgite layer, because of the strong influence of minor amounts of water on the transformation boundary, the depth of the mantle X-discontinuity could serve as a potentially sensitive indicator of water content in the uppermantle.

58 GEOSCIENCES↗

Cosmological simulations of Milky Way-sized galaxies

We introduce a new set of eight Milky Way-sized cosmological simulations performed using the AMR code Adaptive Refinement Tree + Hydrodynamics in a ΛCDM cosmology. The set of zoom-in simulations covers present-day virial masses that range from 8.3×10{sup 11} M{sub ⊙} to 1.56×10{sup 12} M{sub ⊙} and is carried out with our simple but effective prescriptions for deterministic star formation (SF) and “explosive” stellar feedback. The work focuses on showing the goodness of the simulated set of “field” Milky Way-sized galaxies. To this end, we compare some of the predicted physical quantities with the corresponding observed ones. Our results are as follows. (a) In agreement with some previous works, we found curves of circular velocity that are flat or slightly peaked. (b) All simulated galaxies with a significant disk component are consistent with the observed Tully–Fisher, radius–mass, and cold gas–stellar mass correlations of late-type galaxies. (c) The disk-dominated galaxies have stellar specific angular momenta in agreement with those of late-type galaxies, with values around 10{sup 3} km s{sup −1} kpc{sup −1}. (d) The SF rates at z = 0 of all runs but one are comparable to those estimated for the star-forming galaxies. (e) The two most spheroid-dominated galaxies formed in halos with late active merger histories and late bursts of SF, but the other run that ends also dominated by a spheroid never had major mergers. (f) The simulated galaxies lie in the semi-empirical correlation of stellar to halo mass of local central galaxies, and those that end up as disk-dominated evolve mostly along the low-mass branch of this correlation. Moreover, the growth histories of baryonic and stellar mass of these galaxies are proportional to their growth histories of halo mass over the last 6.5–10 Gyr. (g) Within the virial radii of the simulations, ≈25%–50% of the baryons are missed; the amount of gas in the halo is similar to the amount in stars in the galaxy, and most of this gas is in the warm–hot phase. (h) The profiles of z∼0 vertical gas velocity dispersion, σ{sub z}(r), are nearly flat and can be mostly explained by the kinetic energy injected by stars. The average values of σ{sub z} increase at higher redshifts, following roughly the shape of the SF history.

79 ASTRONOMY AND ASTROPHYSICS↗

THE M31 VELOCITY VECTOR. III. FUTURE MILKY WAY M31-M33 ORBITAL EVOLUTION, MERGING, AND FATE OF THE SUN

We study the future orbital evolution and merging of the Milky Way (MW)-M31-M33 system, using a combination of collisionless N-body simulations and semi-analytic orbit integrations. Monte Carlo simulations are used to explore the consequences of varying all relevant initial phase-space and mass parameters within their observational uncertainties. The observed M31 transverse velocity from Papers I and II implies that the MW and M31 will merge t = 5.86{sup +1.61}{sub -0.72} Gyr from now. The first pericenter occurs at t = 3.87{sup +0.42}{sub -0.32} Gyr, at a pericenter distance of r = 31.0{sup +38.0}{sub -19.8} kpc. In 41% of Monte Carlo orbits, M31 makes a direct hit with the MW, defined here as a first-pericenter distance less than 25 kpc. For the M31-M33 system, the first-pericenter time and distance are t = 0.85{sup +0.18}{sub -0.13} Gyr and r = 80.8{sup +42.2}{sub -31.7} kpc. By the time M31 gets to its first pericenter with the MW, M33 is close to its second pericenter with M31. For the MW-M33 system, the first-pericenter time and distance are t = 3.70{sup +0.74}{sub -0.46} Gyr and r = 176.0{sup +239.0}{sub -136.9} kpc. The most likely outcome is for the MW and M31 to merge first, with M33 settling onto an orbit around them that may decay toward a merger later. However, there is a 9% probability that M33 makes a direct hit with the MW at its first pericenter, before M31 gets to or collides with the MW. Also, there is a 7% probability that M33 gets ejected from the Local Group, temporarily or permanently. The radial mass profile of the MW-M31 merger remnant is significantly more extended than the original profiles of either the MW or M31, and suggests that the merger remnant will resemble an elliptical galaxy. The Sun will most likely ({approx}85% probability) end up at a larger radius from the center of the MW-M31 merger remnant than its current distance from the MW center, possibly further than 50 kpc ({approx}10% probability). There is a {approx}20% probability that the Sun will at some time in the next 10 Gyr find itself moving through M33 (within 10 kpc), but while dynamically still bound to the MW-M31 merger remnant. The arrival and possible collision of M31 (and possibly M33) with the MW is the next major cosmic event affecting the environment of our Sun and solar system that can be predicted with some certainty.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization of a Laminate Flat Plate Diffusion Flame in Microgravity using PIV, Visible and CH Emissions

Motivated by fire safety concerns and the advent of long-term micro-gravity facilities, a cooperative program has been developed to study the mechanisms and material properties that control flow assisted (co-current) flame spread. This program has used as a common fire scenario a reacting steady-state boundary layer. Preliminary studies explored the aerodynamics of a reacting boundary layer by simulating a condensed fuel by means of a gas burner. Stability curves for ethane air flames were obtained and different burning regimes were identified. An important feature of this study was the independent identification of the different mechanisms leading to the instability of the flow. It was observed that fuel injection velocity and thermal expansion independently contributed to the separation of the flow at the leading edge of the burner. The occurrence of separation resulted in complex three-dimensional flow patterns that have a dominant effect on critical fire safety parameters such as the stand-off distance and flame length. This work was extended to a solid fuel (PMMA) leading to a Sounding Rocket experiment (Mini-Texus-6). The solid phase showed similar flow patterns, mostly present at low flow velocities (<100 mm/s) but the results clearly demonstrated that the thermal balance at the pyrolyzing fuel surface is the dominant mechanism that controls both stand-off distance and flame length. This thermal balance could be described in a global manner by means of a total mass transfer or "B" number. This "B" number incorporates surface re-radiation, radiative feedback and in-depth heat conduction as first prescribed by Emmons. The mass transfer number becomes the single parameter that determines the evolution of these fire safety variables (flame length, stand-off distance) and therefore can be used as a ranking criterion to assess the flammability of materials. The particular configuration is representative of the NASA upward flame spread test (Test 1) therefore this approach can be used in the interpretation of the results obtained from this test. Nevertheless, complete validation of this approach has not been fully achieved due, mainly because all the measurements necessary to compare with the theoretical predictions have not been obtained. Following these studies two different directions have been taken. The first attempts to elucidate the details of the gas phase combustion reaction and the associated flow field by means of quantitative and qualitative measurements. The second approach, a more practical one, is to apply this methodology to the assessment of material flammability. The former is currently being conducted with a gas burner because it allows for easier control and longer experimentation time. The results obtained so far will be presented in more detail. The latter is a new program therefore only a brief summary of the objectives will be presented.

Joulain, P.↗

Magma intrusion near Volcan Tancitaro: Evidence from seismic analysis

Between May and June 2006, an earthquake swarm occurred near Volcan Tancítaro in Mexico, which was recorded by a temporary seismic deployment known as the MARS network. We located ~1000 events from this seismic swarm. Previous earthquake swarms in the area were reported in the years 1997, 1999 and 2000. We relocate and analyze the evolution and properties of the 2006 earthquake swarm, employing a waveform cross-correlation-based phase repicking technique. Hypocenters from 911 events were located and divided into eighteen families having a correlation coefficient at or above 0.75. 90% of the earthquakes provide at least sixteen phase picks. We used the single-event location code Hypo71 and the P-wave velocity model used by the Jalisco Seismic and Accelerometer Network to improve hypocenters based on the correlation-adjusted phase arrival times. We relocated 121 earthquakes, which show clearly two clusters, between 9–10 km and 3–4 km depth respectively. The average location error estimates are <1 km epicentrally, and <2 km in depth, for the largest event in each cluster. Depths of seismicity migrate upward from 16 to 3.5 km and exhibit a NE-SW trend. The swarm first migrated toward Paricutin Volcano but by mid-June began propagating back toward Volcán Tancítaro. In addition to its persistence, noteworthy aspects of this swarm include a quasi-exponential increase in the rate of activity within the first 15 days; a b-value of 1.47; a jug-shaped hypocenter distribution; a shoaling rate of ~5 km/month within the deeper cluster, and a composite focal mechanism solution indicating largely reverse faulting. As a result, these features of the swarm suggest a magmatic source elevating the crustal strain beneath Volcan Tancítaro.

58 GEOSCIENCES↗

Acoustoelastic characterisation of materials

Recent developments in the use of nonlinear acoustoelastic techniques for the characterization of materials are reviewed, with particular attention given to the quantitative modeling of the acoustoelastic properties of multiphase alloys in terms of material composition. Topics discussed include general theory of acoustoelasticity, models of two-phase systems, measurement techniques, acoustic nonlinearity parameters and material properties, stress-temperature dependence of sound velocity, and low field magnetoacoustics. The discussion is illustrated by results obtained for steels and aluminum alloys with second phase precipitates.

Cantrell, J. H.↗

Testing approximate predictions of displacements of cosmological dark matter halos

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at z =0,0.5,1, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods, and this is true even for mesh sizes as large as the inter-particle distance. This test sets an upper limit on the ability of these methods to reproduce the clustering of halos, for the cases when these objects are reconstructed at the object-by-object level.

79 ASTRONOMY AND ASTROPHYSICS↗

Pore-scale investigation on stress-dependent characteristics of granular packs and the impact of pore deformation on fluid distribution

Understanding the effect of changing stress conditions on multiphase flow in porous media is of fundamental importance for many subsurface activities including enhanced oil recovery, water drawdown from aquifers, soil confinement, and geologic carbon storage. Geomechanical properties of complex porous systems are dynamically linked to flow conditions, but their feedback relationship is often oversimplified due to the difficulty of representing pore-scale stress deformation and multiphase flow characteristics in high fidelity. In this work, we performed pore-scale experiments of single- and multiphase flow through bead packs at different confining pressure conditions to elucidate compaction-dependent characteristics of granular packs and their impact on fluid flow. A series of drainage and imbibition cycles were conducted on a water-wet, soda-lime glass bead pack under varying confining stress conditions. Simultaneously, X-ray micro-CT was used to visualize and quantify the degree of deformation and fluid distribution corresponding with each stress condition and injection cycle. Micro-CT images were segmented using a gradient-based method to identify fluids (e.g., oil and water), and solid phase redistribution throughout the different experimental stages. Changes in porosity, tortuosity, and specific surface area were quantified as a function of applied confining pressure. Results demonstrate varying degrees of sensitivity of these properties to confining pressure, which suggests that caution must be taken when considering scalability of these properties for practical modeling purposes. Changes in capillary number with confining pressure are attributed to the increase in pore velocity as a result of pore contraction. Furthermore, this increase in pore velocity was found to have a marginal impact on average phase trapping at different confining pressures.

58 GEOSCIENCES↗

Using principal component analysis to distinguish different dynamic phases in superconducting vortex matter

Vortices in type-II superconductors driven over random disorder are known to exhibit a remarkable variety of distinct nonequilibrium dynamical phases that arise owing to the competition between vortex-vortex interactions, the quenched disorder, and the drive. These include pinned states, elastic flows, plastic or disordered flows, and dynamically reordered moving crystal or moving smectic states. The plastic flow phases can be particularly difficult to characterize since the flows are strongly disordered. Here, we perform principal component analysis (PCA) on the positions and velocities of vortex matter moving over random disorder for different disorder strengths and drives. We find that PCA can distinguish the known dynamic phases as well as or better than previous measures based on transport signatures or topological defect densities. In addition, PCA recognizes distinct plastic flow regimes, a slowly changing channel flow and a moving amorphous fluid flow, that do not produce distinct signatures in the standard measurements. In conclusion, our results suggest that this position and velocity-based PCA approach could be used to characterize dynamic phases in a broader class of systems that exhibit depinning and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗