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.
The mechanisms, dynamics, and implications of self-sealing and CO 2 resistance in wellbore cements
Here, this study analyzes the dynamics and mechanisms of the interactions of carbonated brine with hydrated-Portland-cement; in particular, the study focuses on self-sealing, a process whereby hydrated-Portland cement reacts with carbonated brine to for silica and calcium carbonate in sufficient quantities to seal the flow pathway. The analysis is based on a comprehensive set of reactive-transport simulations that explore the complex coupled dynamics between the fluid flow and mineral reactions that underlie self-sealing, and it relies heavily on the synthesis of the extensive body of work on wellbore integrity that has been conducted over the past decade. The analysis explores a large chemical and mineralogical diversity and a wide range in physical conditions and flow regimes, attempting to assess the robustness of the analysis. Self-sealing conditions arise over a wide range in cement properties and reservoir conditions. Although some properties and conditions promote a stronger self-sealing response, self-sealing occurs for a wide range of Ca:Si ratios in cement and for various reservoir fluid compositions. Self-sealing conditions move along a wellbore proportional to the flux of the leaking carbonated brine, and the reaction zone spreads out proportional to the fluid velocity, where volumetric flux and velocity are related by porosity (flux = velocity * porosity). However, self-sealing conditions can be maintained in a specific section of a wellbore by controlling the pressure drive and/or effective wellbore permeability, which in turn can limit the flux and velocity of any leaking fluid. Finally, the phases produced by hydrating Portland cement represent a carbonic cement that will react with a carbonated brine to produce end products (calcium carbonate and silica) that can maintain integrity in the presence of carbonic acid. Lastly, the attributes that make hydrated Portland cement phases a carbonic cement are required for self-sealing.
Synchrotron X-ray based particle image velocimetry to measure multiphase streamflow and densitometry
Synchrotron-based radiation techniques are increasingly used in the field of flow condition measurements. As a complicated multiphase flow, cavitation has been investigated in a wide variety of industrial fields, to prevent damages on spillways at large dams or intensify the bacteria eradication rate during water treatments. To control the detrimental effects of cavitating flows, the dynamics of and mechanisms affecting cavitation development shall be visualized and identified in real-time. Due to the large reflection and scattering effects of multiphase flow, standard visible-light imaging technologies cannot penetrate vapor bubbles and provide extremely limited information. In this study, the synchrotron source emitted high-flux X-ray pulses are used to capture flow motion and visualize internal structures with sufficiently high spatial and temporal resolutions. The image processing procedures on raw x-ray images are developed using 2D Fourier transform (FT) and wavelet transform (WT) to implement the contrast enhancement and de-noising. Through the edge detection algorithm, tracking either seeded particles or phase interfaces inside the opaque multiphase flow can be achieved to perform the particle image velocimetry. Further, with decomposing particle-only and particle-eliminated subplots from the X-ray phase-contrast image, the time-resolved velocity and void fraction fields are obtained simultaneously, paving a way for further flow condition analysis such as densitometry.
Frequency–Magnitude Statistics of Laboratory Foreshocks Vary With Shear Velocity, Fault Slip Rate, and Shear Stress
Understanding the nucleation phase of earthquakes is key for advancing earthquake hazard assessment and improving earthquake early warning systems. However, little progress has been made in this area due to a poor understanding of nucleation processes and incomplete seismic and fault zone measurements. The ability to integrate measured fault zone properties with seismic data could significantly improve our understanding of how earthquakes begin and whether there are systematic variations in seismic properties preceding failure. In this work, we use high-resolution laboratory measurements of fault zone properties along with acoustic emission data to document temporal variations of foreshock properties. Our data show that foreshock size increases with shear stress, loading rate, and fault slip rate. We propose that the preseismic fault slip rate and fault zone thickness (i.e., porosity) work in concert to modulate foreshock properties.
Fluid-like shock compression of dilute polymer nanocomposites
The shock Hugoniot of heterogenous mixtures of discrete particles has been experimentally investigated for porous agglomerates and fluids, while the study of full-density solids has been primarily limited to compressed powders and high-dimensional composites. By dispersing ceria nanoparticles in a polyethylene matrix, we are able to examine the hydrodynamic behavior of a nonporous, heterogenous solid in thermal equilibrium during weak shock compression. Phase-driven discontinuities in the Hugoniot particle velocity–shock velocity (u−D) relationship of pure polyethylene are replicated in the nanocomposites but are shifted to lower velocity and to higher pressure with higher particle concentration. The results are explained using an isothermal, two-velocity fluid model under the hydrodynamic approximation. The model, which assumes a theoretical equation-of-state for ceria and either a low-order or high-order fit to the measured polyethylene Hugoniot, reasonably predicts the Hugoniot for two different polyethylene/ceria nanocomposites. Using the model, the mixture Hugoniot is shown to be insensitive to the Hugoniot of the stiffer constituent when the moduli are sufficiently disparate, while dependence on particle density and volume fraction is preserved through fluid-like motion.
Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation
Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.
NUMERICAL SIMULATIONS OF CONVERSION TO ALFVEN WAVES IN SUNSPOTS
We study the conversion of fast magnetoacoustic waves to Alfven waves by means of 2.5D numerical simulations in a sunspot-like magnetic configuration. A fast, essentially acoustic, wave of a given frequency and wave number is generated below the surface and propagates upward through the Alfven/acoustic equipartition layer where it splits into upgoing slow (acoustic) and fast (magnetic) waves. The fast wave quickly reflects off the steep Alfven speed gradient, but around and above this reflection height it partially converts to Alfven waves, depending on the local relative inclinations of the background magnetic field and the wavevector. To measure the efficiency of this conversion to Alfven waves we calculate acoustic and magnetic energy fluxes. The particular amplitude and phase relations between the magnetic field and velocity oscillations help us to demonstrate that the waves produced are indeed Alfven waves. We find that the conversion to Alfven waves is particularly important for strongly inclined fields like those existing in sunspot penumbrae. Equally important is the magnetic field orientation with respect to the vertical plane of wave propagation, which we refer to as 'field azimuth'. For a field azimuth less than 90 Degree-Sign the generated Alfven waves continue upward, but above 90 Degree-Sign downgoing Alfven waves are preferentially produced. This yields negative Alfven energy flux for azimuths between 90 Degree-Sign and 180 Degree-Sign . Alfven energy fluxes may be comparable to or exceed acoustic fluxes, depending upon geometry, though computational exigencies limit their magnitude in our simulations.
The HETDEX pilot survey. V. The physical origin of Lyα emitters probed by near-infrared spectroscopy
We present the results from a Very Large Telescope/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyα emitters (LAEs) at z = 2.1-2.5 in the COSMOS and GOODS-N fields discovered from the Hobby Eberly Telescope Dark Energy Experiment Pilot Survey. We detect rest-frame optical nebular lines (Hα and/or [O III] λ5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas mass fraction, and Lyα velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Lyα relative to the systemic redshift ranging between +85 and +296 km s{sup –1} with a mean of +180 km s{sup –1}. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Lyα velocity offsets show a moderate correlation with the measured SFR (2.5σ), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass (≲1.5σ). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Lyα photons, we find no signature of selective quenching of resonantly scattered Lyα photons. However, we also find no evidence that a clumpy ISM is enhancing the Lyα equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density and less dust, easing the escape of Lyα photons over that in continuum-selected star-forming galaxies.
THE TURBULENCE VELOCITY POWER SPECTRUM OF NEUTRAL HYDROGEN IN THE SMALL MAGELLANIC CLOUD
We present the results of the Velocity Coordinate Spectrum (VCS) technique to calculate the velocity power spectrum of turbulence in the Small Magellanic Cloud (SMC) in 21 cm emission. We present an updated version of the VCS technique that takes into account regular motions, which is an important factor in our SMC VCS analysis. We have obtained a velocity spectral index of −3.85, a cold phase sonic Mach number of 5.6, and an injection scale of 2.3 kpc. The spectral index is steeper than the Kolmogorov index, which is expected for shock-dominated turbulence. The injection scale of 2.3 kpc suggests that HI supershells or tidal interactions with the Large Magellanic Cloud are the dominant drivers of turbulence in this dwarf galaxy. This implies that turbulence may be driven by multiple mechanisms in galaxies and that galaxy–galaxy interactions may play an important role in addition to supernova feedback.
Space-Time Structured Plasma Waves
Electrostatic waves play a critical role in nearly every branch of plasma physics from fusion to advanced accelerators, to astro, solar, and ionospheric physics. The properties of planar electrostatic waves are fully determined by the plasma conditions, such as density, temperature, ionization state, or details of the distribution functions. Furthermore we demonstrate that electrostatic wave packets structured with space-time correlations can have properties that are independent of the plasma conditions. For instance, an appropriately structured electrostatic wave packet can travel at any group velocity, even backward with respect to its phase fronts, while maintaining a localized energy density. These linear, propagation-invariant wave packets can be constructed with or without orbital angular momentum by superposing natural modes of the plasma and can be ponderomotively excited by space-time structured laser pulses like the flying focus.
Earthquake Locations
The location of an earthquake is a classic inverse problem. Arrival times of various seismic phases are combined with the location of the instrument to form the data vector. The inverse operator is applied in an iterative manor until the solution converges and the resultant solution is the hypocentral location of the earthquake in space and time. In this presentation, I will present how seismologists obtain local earthquake hypocenters. I will present how the generalized inverse operator is built using the simple distance = rate x time equation and how the inverse operator is used in an iterative manor to converge on the hypocenter of the earthquake. Once an earthquake location is determined, seismologists can then use the location to infer other information about the medium through which the waves traveled. One possible application of the earthquake locations is a branch of seismology known as tomography. In tomography, source and station locations are used to compute the travel times of seismic phases and these phases are then inverted to obtain a velocity model of the medium. Seismic tomography can be applied in both 2-D and 3-D settings and in this presentation, I will present examples from both. The first example comes from a technique known as ambient noise tomography. Rayleigh waves are obtained by cross correlating the background noise between two stations. Because we know the interstation distance between the two stations and the arrival times of the Rayleigh waves, we can use tomography to build a group velocity map of the study region. I will present results from a study of the Eastern Tennessee Seismic Zone using the EarthScope’s Transportable Array together with different seismic networks in the region, figure 1. The second example uses P and S-waves from 412 earthquakes which occurred during the 1994 Draney Peak swarm to calculate a 3-D velocity model of the area. Preliminary results show a strong velocity contrast at a depth of about 4 km where the majority of the earthquakes occur, figure 2.
Apparatus and method for noninvasive particle detection using doppler spectroscopy
An apparatus and method for noninvasively detecting the presence of solid particulate matter suspended in a fluid flowing through a pipe or an oil and gas wellbore are described. Fluid flowing through a conduit containing the particulate solids is exposed to a fixed frequency (>1 MHz) of ultrasonic vibrations from a transducer attached to the outside of the pipe. The returning Doppler frequency shifted signal derived from the scattering of sound from the moving solid particles is detected by an adjacent transducer. The transmitted signal and the Doppler signal are combined to provide sensitive particulate detection. The magnitude of the signal and the Doppler frequency shift are used to determine the particle size distribution and the velocity of the particles. Measurement of the phase shift between the applied frequency and the detected Doppler shifted may be used to determine the direction of motion of the particles.
Development of integrated core disruptive accident analysis code for FBR - ASTERIA-FBR
The evaluation of consequence at the severe accident is the most important as a safety licensing issue for the reactor core of liquid metal cooled fast breeder reactor (LMFBR), since the LMFBR core is not in an optimum condition from the viewpoint of reactivity. This characteristics might induce a super-prompt criticality due to the core geometry change during the core disruptive accident (CDA). The previous CDA analysis codes have been modeled in plural phases dependent on the mechanism driving a super-prompt criticality. Then, the following event is calculated by connecting different codes. This scheme, however, should introduce uncertainty and/or arbitrary to calculation results. To resolve the issues and obtain the consistent calculation results without arbitrary, JNES is developing the ASTERIA-FBR code for the purpose of providing the cross-check analysis code, which is another required scheme to confirm the validity of the evaluation results prepared by applicants, in the safety licensing procedure of the planned high performance core of Monju. ASTERIA-FBR consists of the three major calculation modules, CONCORD, dynamic-GMVP, and FEMAXI-FBR. CONCORD is a three-dimensional thermal-hydraulics calculation module with multi-phase, multi-component, and multi-velocity field model. Dynamic-GMVP is a space-time neutronics calculation module. FEMAXI-FBR calculates the fuel pellet deformation behavior and fuel pin failure behavior. This paper describes the needs of ASTERIA-FBR development, major module outlines, and the model validation status. (authors)
Stochastic inflation in phase space: is slow roll a stochastic attractor?
An appealing feature of inflationary cosmology is the presence of a phase-space attractor, ''slow roll'', which washes out the dependence on initial field velocities. We investigate the robustness of this property under backreaction from quantum fluctuations using the stochastic inflation formalism in the phase-space approach. A Hamiltonian formulation of stochastic inflation is presented, where it is shown that the coarse-graining procedure—where wavelengths smaller than the Hubble radius are integrated out—preserves the canonical structure of free fields. This means that different sets of canonical variables give rise to the same probability distribution which clarifies the literature with respect to this issue. The role played by the quantum-to-classical transition is also analysed and is shown to constrain the coarse-graining scale. In the case of free fields, we find that quantum diffusion is aligned in phase space with the slow-roll direction. This implies that the classical slow-roll attractor is immune to stochastic effects and thus generalises to a stochastic attractor regardless of initial conditions, with a relaxation time at least as short as in the classical system. For non-test fields or for test fields with non-linear self interactions however, quantum diffusion and the classical slow-roll flow are misaligned. We derive a condition on the coarse-graining scale so that observational corrections from this misalignment are negligible at leading order in slow roll.
A Case Study of Nonresonant Mode 3-s ULF Waves Observed by MMS
The nature of the 3-s ultralow frequency (ULF) wave in the Earth's foreshock region and the associated wave-particle interaction are not yet well understood. In this work, we investigate the 3-s ULF waves using Magnetospheric Multiscale (MMS) observations. By combining the plasma rest frame wave properties obtained from multiple methods with the instability analysis based on the velocity distribution in the linear wave stage, the ULF wave is determined to be due to the ion/ion nonresonant mode instability. The interaction between the wave and ions is analyzed using the phase relationship between the transverse wave fields and ion velocities and using the longitudinal momentum equation. During the stage when ULF waves have sinusoidal waveforms up to |dB|/|B 0 | ~ 3, where dB is the wave magnetic field and B0 is the background magnetic field, the wave electric fields perpendicular to B 0 do negative work to solar wind ions; along B 0 , a longitudinal electric field develops, but the V × B force is stronger and leads to solar wind ion deceleration. During the same wave stage, the backstreaming beam ions gain energy from the transverse wave fields and get deceleration along B 0 by the longitudinal electric field. The ULF wave leads to electron heating, preferentially in the direction perpendicular to the local magnetic field. Secondary waves are generated within the ULF waveforms, including whistler waves near half of the electron cyclotron frequency, high-frequency electrostatic waves, and magnetosonic whistler waves. The work improves the understanding of the nature of 3-s ULF waves and the associated wave-particle interaction.
Measuring the streaming motion in the Milky Way disc with Gaia EDR3+
We map the 3D kinematics of the Galactic disc out to 3.5 kpc from the Sun, and within 0.75 kpc from the mid-plane of the Galaxy. To this end, we combine high-quality astrometry from Gaia EDR3, with heliocentric line-of-sight velocities from Gaia DR2, and spectroscopic surveys including APOGEE, GALAH, and LAMOST. We construct an axisymmetric model for the mean velocity field, and subtract this on a star-by-star basis to obtain the residual velocity field in the Galactocentric components (V Φ , V R , V z ), and V los . The velocity residuals are quantified using the power spectrum, and we find that the peak power (A/[km s –1 ]) in the mid-plane (|z| < 0.25 kpc) is (A Φ , A R , A Z , A los ) = (4.2,8.5,2.6,4.6), at 0.25 < |z|/[kpc] < 0.5, is (A Φ , A R , A Z , A los ) = (4.0,7.9,3.6,5.3), and at 0.5 < |z|/[kpc] < 0.75, is (A Φ , A R , A Z , A los ) = (1.9,6.9,5.2,6.4). Our results provide a sophisticated measurement of the streaming motion in the disc and in the individual components. We find that streaming is most significant in VR, and at all heights (|Z|) probed, but is also non-negligible in other components. Additionally, we find that patterns in velocity field overlap spatially with models for spiral arms in the Galaxy. Our simulations show that phase-mixing of disrupting spiral arms can generate such residuals in the velocity field, where the radial component is dominant, just as in real data. Furthermore, we also find that with time evolution, both the amplitude and physical scale of the residual motion decrease.
The local dark matter phase-space density and impact on WIMP direct detection
We present a new determination of the local dark matter phase-space density. This result is obtained implementing, in the limit of isotropic velocity distribution and spherical symmetry, Eddington's inversion formula, which links univocally the dark matter distribution function to the density profile, and applying, within a Bayesian framework, a Markov Chain Monte Carlo algorithm to sample mass models for the Milky Way against a broad and variegated sample of dynamical constraints. We consider three possible choices for the dark matter density profile, namely the Einasto, NFW and Burkert profiles, finding that the velocity dispersion, which characterizes the width in the distribution, tends to be larger for the Burkert case, while the escape velocity depends very weakly on the profile, with the mean value we obtain being in very good agreement with estimates from stellar kinematics. The derived dark matter phase-space densities differ significantly — most dramatically in the high velocity tails — from the model usually taken as a reference in dark matter detection studies, a Maxwell-Boltzmann distribution with velocity dispersion fixed in terms of the local circular velocity and with a sharp truncation at a given value of the escape velocity. We discuss the impact of astrophysical uncertainties on dark matter scattering rates and direct detection exclusion limits, considering a few sample cases and showing that the most sensitive ones are those for light dark matter particles and for particles scattering inelastically. As a general trend, regardless of the assumed profile, when adopting a self-consistent phase-space density, we find that rates are larger, and hence exclusion limits stronger, than with the standard Maxwell-Boltzmann approximation. Tools for applying our result on the local dark matter phase-space density to other dark matter candidates or experimental setups are provided.
Laminar burning velocity predictions for C 1 and C 2 hydrofluorocarbon refrigerants with air
Due to their high global warming potentials, many existing working fluids are being phased out. Their replacements will often be slightly flammable, and the burning velocity of refrigerant-air mixtures is a metric for ranking their flammability. To allow estimates of the flammability of new blends of agents, predictive tools for the burning velocity are being developed and require a kinetic mechanism. A hydrofluorocarbon (HFC) mechanism was developed 20 years ago to describe hydrocarbon-air flames with added trace amounts of hydrofluorocarbon fire retardants. In the present work, the mechanism has been updated slightly to include new HFC compounds, new rate data. The modified mechanism is used to predict steady, planar, 1D, unstretched burning velocities for mixtures of air with one- and two-carbon saturated HFC compounds R41 (CH 3 F), R32 (CH 2 F 2 ), R161 (C 2 F 5 H), R152 (CH 2 F-CH 2 F), R152a (CH 3 -CHF 2 ), R143 (CH 2 F-CHF 2 ), R143a (CH 3 -CF 3 ), R134 (CHF 2 -CHF 2 ), and R134a (CH 2 F-CF 3 ), for which existing experimental data were available.