Search NASA⌕ Search

SEARCH · Search NASA

Results for “migration velocity analysis”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Joint PP and PS plane-wave wave-equation migration-velocity analysis

Conventional joint PP and PS velocity analysis is based on ray tomography. We develop a joint PP and PS wave-equation migration-velocity-analysis method using plane-wave common-image gathers (CIGs) to produce accurate P- and S-wave velocity models. The objective function of our new method consists of three terms: The first and second terms penalize the moveout residuals computed from PP and PS plane-wave CIGs, respectively, and the third term constrains the nonzero relative depth shifts between the PP and PS migration images. The moveout of plane-wave CIGs is automatically picked using a semblance analysis method, and the relative depth shifts between the PP and PS images are automatically computed using dynamic warping or manually picking the depths of certain primary reflectors. The moveout residuals and the relative depth shifts are transformed into weighted image perturbations, and they are then projected into the velocity models to update the P- and S-wave velocity models using the scalar-wave equations and their linearized forms. Here, numerical tests with synthetic and multicomponent field data demonstrate that our method can simultaneously invert for accurate P- and S-wave velocity models for elastic migration.

58 GEOSCIENCES↗
DOE OSTI2020

Efficient 3-D velocity model building using joint inline and crossline plane-wave wave-equation migration velocity analyses

SUMMARY Wave-equation migration velocity analysis (WEMVA) is an image-domain inversion method for velocity model building. Automatic plane-wave WEMVA (PWEMVA) calculates the moveouts of plane-wave common-image gathers (CIGs) by searching a best-fitting parabola with semblance analysis and backprojects residual CIG moveouts into wavefield wave paths with a reflection tomographic kernel. However, 3-D PWEMVA is very computationally expensive because 3-D reflection tomographic inversion requires at least five 3-D reverse-time migrations per iteration and stores two types of source wavefields at model boundaries. We develop a joint inline and crossline PWEMVA method for efficient 3-D velocity model building. We alternatively implement the inline and crossline PWEMVAs with a constraint for each other, in which we iteratively construct the 3-D velocity model update through 1-D spline interpolation of 2-D gradients. The inline and crossline joint inversion is practical since PWEMVA only inverts for low-wavenumber velocity perturbations along wave paths, and the method can take less than 1 per cent of the computational cost of full 3-D PWEMVA. To construct unaliased plane waves for our joint inline and crossline PWEMVA, we develop a 3-D data interpolation method in the frequency–wavenumber (FK) domain to recover regularly and randomly missing traces. The method minimizes the misfit on sufficiently localized data subsets with iterative optimal step lengths and a gradient preconditioner that iteratively selects dominant dips along different azimuths. In numerical experiments, we use a 3-D synthetic seismic data set and a land 3-D field seismic data set acquired at the Farnsworth CO2-EOR (enhanced oil recovery) field to demonstrate the efficacy of our velocity model building and data interpolation methods.

Liu, Xuejian↗
DOE OSTI2021

Migration of inclusions in a matrix due to a spatially varying interface energy

An interfacial energy can be a function of a bulk field such as temperature or electric field. Here we find that in a system with a gradient in temperature or electric potential, the resulting variation in interfacial energy can induce a particle to migrate by either surface or bulk diffusion. For a circular particle under a constant unidirectional gradient in the bulk field, the field dependence of the interfacial energy induces a circular particle to move as a circle at a constant velocity along the direction from higher to lower interfacial energy. A linear stability analysis of this steady state migration suggests that perturbations will damp out as time evolves, and thus under these conditions a migrating circular particle is morphologically stable. Other spatial distributions of interface energy can lead to the distortion of an initially circular shaped particle during migration. A phase field model is developed that captures these distortions and verifies the theoretical results mentioned above.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Structural and Stratigraphic Controls on Methane Hydrate occurrence and distribution: Gulf of Mexico, Walker Ridge 313 and Green Canyon 955: Final Report

The goal of this project was to determine structural and stratigraphic controls on hydrate occurrence and distribution in Green Canyon (GC) 955 and Walker Ridge (WR) 313 blocks using seismic and well data. Gas hydrate was discovered in these blocks in coarse- and fine-grained sediments during the 2009 Joint Industrial project (JIP) Leg 11 drilling expedition. Although the immediate interest of the exploration community is exclusively hydrate which is present in coarse–grained sediments, factors that control hydrate and free gas distribution in the two blocks and whether coarse and fine-grained hydrate-bearing units are related in any manner, formed the core of this research. The project spanned from 10/01/2012 to 07/31/2016. In the project, in both the leased blocks, the interval spanning the gas hydrate stability zone (GHSZ) was characterized using a joint analysis of sparse Ocean Bottom Seismic (OBS) and dense, surface–towed multichannel seismic (MCS) data. The project team had the luxury of calibrating their results with two well logs. Advance processing methods such as depth migration and full-waveform inversion (FWI) were used for seismic data analysis. Hydrate quantification was achieved through interpretation of the FWI velocity field using appropriate rock physics models at both blocks. The seismic modeling/inversion methodology (common to both GC955 and WR313 blocks) was as follows. First, the MCS data were depth migrated using a P-wave velocity (VP) model constructed using inversion of reflection arrival times of a few (four in both cases) key horizons carefully picked in the OBS data to farthest possible offsets. Then, the resolution of the traveltime VP model was improved to wavelength scale by inverting OBS gathers up to the highest frequency possible (21.75 Hz for GC955 and 17.5 for WR313) using FWI. Finally, the hydrate saturation (or the volume fraction) was estimated at the well location assuming one of the other hydrate morphology (filling the primary or the secondary porosity) was extrapolated out from the wells using the FWI VP as a guide. General outcomes were as follows. First and foremost, an imaging methodology using sparse seismic data, which is easily replicable at other sites with similar datasets, has been demonstrated. The end product of this methodology at both the leased blocks is quantitative estimates of hydrate distribution. Second, at both locations there is strong evidence that the base of the GHSZ, which does not appear as a clear Bottom Simulating Reflection (BSR), manifests in the VP perturbations created by FWI, suggesting that FWI is sensitive to subtle compositional changes in shallow sediments and establishes it as a valuable tool for investigations of hydrate-bearing basins. Third, through joint interpretation of the depth migrated image and the FWI VP model, how structure and stratigraphy jointly determine hydrate and free gas distribution in both blocks could be clearly visualized. The joint interpretation also suggests that the coarse and fine grained hydrate-bearing sediments at both leased are connected. Site specific results, in addition to general results, are as follows. At GC955 the overlying fine-grained hydrate-bearing unit could have been sourced from the underlying hydrate coarse-grained channel-levee complex through a chimney feature. The channel-levee system at GC955 is compartmentalized by faults, of which only a few may be impermeable. Although compartmentalized, the channel-levee system in the GC955 as a whole might be in communication except selected zones. At WR313 the overlying fine-grained fracture-filled hydrate unit appears to be sourced from below the GHSZ. The reason that only a particular fine-grained unit has hydrate, despite having lower porosity that the bounding units, could be the presence of secondary porosity (such as those formed from clay dewatering under compaction). In conclusion, the project was a pioneering effort in in joint analysis of OBS and MCS datasets for advancing the knowledge about a hydrate and free–gas system dynamics using advanced processing methods such as FWI and depth migration. Results obtained in this project can greatly advance the tools and techniques used for delineating specific hydrate prospects. Results obtained in this project can also be seamlessly incorporated into other DOE funded project on modeling the potential productivity and commercial viability of hydrate from sand-dominated reservoirs. The OBS and MCS data in this project were acquired in 2012 (after the JIP II drilling) by the USGS and therefore the results are a posteriori. Nonetheless, the seismic inversion workflow established through this project can be used to generate various what-if quantification scenarios even in absence of logs and serve as a valuable tool for guiding drilling operations. Results from this project can augment other DOE sponsored projects on determining the commercial viability of methane production from the Gulf of Mexico.

03 NATURAL GAS↗
DOE OSTI2021

Scale-Dependent Bedform Migration and Deformation in the Physical and Spectral Domains

Multi-scale bathymetries observed in laboratory channels and natural rivers have posed several challenges to the description of bedform geometry and kinematics. To better quantify scale-dependent migration velocity, high resolution spatio-temporal bed evolution data are analyzed using three independent methods: (i) a bedform tracking method that identifies individual bedforms in longitudinal bed elevation profiles, (ii) a cross-correlation based, image matching technique in consecutive bathymetry scans similar to Particle Image Velocimetry, and (iii) two-dimensional frequency-wavenumber spectra marking a dispersion relation between bedform period and length, allowing estimation of Fourier-based, scale-dependent migration velocities. The comparative studies show that the spectral and bedform tracking methods present small differences in the migration velocity of medium-large bedforms. However, significant deviations occur in the small bedforms, mostly due to bedform sheltering or amplified exposure to the flow, which also contribute to their enhanced scale-dependent deformation. The analysis also shows that some form of scale separation is needed to extract the large and slowly migrating bedforms as they are masked by secondary features when several bedform orders coexist. Bathymetric Image Velocimetry reasonably estimates the averaged migration velocity and exhibits potential to study relations between spanwise bed surface velocity and local bedform kinematics.

58 GEOSCIENCES↗
DOE OSTI2016

Kinematics and dynamics of the East Pacific Rise linked to a stable, deep-mantle upwelling [Kinematics and dynamics of the East Pacific Rise linked to whole mantel convective motions]

Earth’s tectonic plates are generally considered to be driven largely by negative buoyancy associated with subduction of oceanic lithosphere. In this context, mid-ocean ridges (MORs) are passive plate boundaries whose divergence accommodates flow driven by subduction of oceanic slabs at trenches. We show that over the past 80 million years (My), the East Pacific Rise (EPR), Earth’s dominant MOR, has been characterized by limited ridge-perpendicular migration and persistent, asymmetric ridge accretion that are anomalous relative to other MORs. We reconstruct the subduction-related buoyancy fluxes of plates on either side of the EPR. The general expectation is that greater slab pull should correlate with faster plate motion and faster spreading at the EPR. Moreover, asymmetry in slab pull on either side of the EPR should correlate with either ridge migration or enhanced plate velocity in the direction of greater slab pull. Based on our analysis, none of the expected correlations are evident. This implies that other forces significantly contribute to EPR behavior. We explain these observations using mantle flow calculations based on globally integrated buoyancy distributions that require core-mantle boundary heat flux of up to 20 TW. The time-dependent mantle flow predictions yield a long-lived deep-seated upwelling that has its highest radial velocity under the EPR and is inferred to control its observed kinematics. Lastly, the mantle-wide upwelling beneath the EPR drives horizontal components of asthenospheric flows beneath the plates that are similarly asymmetric but faster than the overlying surface plates, thereby contributing to plate motions through viscous tractions in the Pacific region.

58 GEOSCIENCES↗
DOE OSTI2022

Using cloud radar to investigate the effect of rainfall on migratory insect flight

The fate of migrating insects that encounter rainfall in flight is a critical consideration when modelling insect movement, but few field observations of this common phenomenon have ever been collected due to the logistical challenges of witnessing these encounters. Operational cloud radars have been deployed around the world by meteorological agencies to study precipitation physics, and as a byproduct, provide a rich database of insect observations that is freely available to researchers. Although considered unwanted ‘clutter’ by the meteorologists who collect the data, the analysis method presented here enables ecologists to delineate co-occurring signals from insects and raindrops. We present a method that uses image processing techniques on cloud radar velocity spectra to examine the fate of migrating insects when they encounter precipitation. By analysing velocity spectra, we can distinguish flying insects from falling rain and compare the relative density of insects in flight before, during and after the rainfall. We demonstrate the method on a case of insect migration in Oklahoma, USA. Using this method, we show the first reconstructed images of migrating insect layers in flight during rainfall. Our analysis shows that mild to moderate rainfall diminishes the number of insects aloft but does not cause full termination of migratory flight, as has previously been suggested. We hope this technique will spur further investigations of how changing weather conditions impact insect migration, and enable some of the first of such studies in regions of the world that are underrepresented in the literature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The effect of stress on the migration of He gas bubbles under a thermal gradient in Fe by phase-field modeling

Here a phase-field model is parameterized to study the effect of elastic stresses on the migration of He gas bubbles in Fe under a temperature gradient. Stresses caused by the gas bubble pressure and residual stress in the Fe matrix are considered. The dependence of He bubble migration velocity on the magnitude of the residual stress, average temperature, temperature gradient, and bubble size is measured. In agreement with a theoretical model based on surface diffusion, simulation results demonstrate that He bubbles move towards the high temperature region with velocities in Fe that are orders of magnitude faster than previously reported in UO 2 . It is found that local stresses in the matrix caused by the He bubble have negligible effect on the bubble migration process; however, residual stresses in the Fe matrix, potentially caused by processing or irradiation, can modestly modify bubble kinetics through pressure dependence of the He diffusion coefficients. Compressive residual stress decreases diffusion coefficients for bulk and surface diffusion mechanisms, thus reducing the migration velocity of the gas bubble. In contrast, tensile residual stress increases the diffusion coefficients, resulting in an increase in the gas bubble migration velocity. This pressure dependence is also consistent with a theoretical model. This phase field model lays the foundation for analysis of bubble coalescence-induced fracture in He bubble-containing steels.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Diffraction imaging for seal evaluation using ultra high resolution 3D seismic data

Cap rock is any impermeable or low permeability formation that may trap oil, gas or water, preventing it from migrating to the surface. Evaluation of sealing properties of the cap rock is a critical task for seismic exploration. It is important for oil and gas prospecting as well as for locating reservoirs appropriate for carbon dioxide storage. The latter has been proposed as one solution to global climate change caused by heat-trapping anthropogenic gasses in the atmosphere. A fluid escaped from the reservoir fills voids in the overlying strata thereby causing local changes of acoustic properties. We aim to explore such fluid-saturated areas and thus evaluate sealing properties of underlying cap rocks. We assume that local acoustic impedance changes caused by fluid migration generate seismic diffractions. We isolate diffracted waves from P-Cable seismic data acquired over the area of interest and perform diffraction analysis of the near-surface interval (above 200 m). In the diffraction image obtained, we interpret an extensive system of broad sinuous channels, narrow linear faults, and diffuse fracture networks. In addition, we detect high diffractivity anomalies, which extend vertically and laterally. We note no evidence of these anomalies in a conventional seismic image and associated seismic attribute (similarity and fracture density) volumes. We associate this high diffractivity to a variable pore fluid composition and saturation, which we infer to be due to hydrocarbon migration from the underlying strata. Consequently, we conclude poor seal quality for the underlying cap rocks. During data processing, we test diffraction focusing analysis for velocity model building. In conclusion, this option appears quite valuable for P-Cable data processing because of short offset and limitations of the standard velocity analysis.

58 GEOSCIENCES↗
DOE OSTI2016

Seismic Fracture Characterization Methodologies for Enhanced Geothermal Systems

Executive Summary The overall objective of this work was the development of surface and borehole seismic methodologies using both compressional and shear waves for characterizing faults and fractures in Enhanced Geothermal Systems. We used both surface seismic and vertical seismic profile (VSP) methods. We adapted these methods to the unique conditions encountered in Enhanced Geothermal Systems (EGS) creation. These conditions include geological environments with volcanic cover, highly altered rocks, severe structure, extreme near surface velocity contrasts and lack of distinct velocity contrasts at depth. One of the objectives was the development of methods for identifying more appropriate seismic acquisition parameters for overcoming problems associated with these geological factors. Because temperatures up to 300º C are often encountered in these systems, another objective was the testing of VSP borehole tools capable of operating at depths in excess of 1,000 m and at temperatures in excess of 200º C. A final objective was the development of new processing and interpretation techniques based on scattering and time-frequency analysis, as well as the application of modern seismic migration imaging algorithms to seismic data acquired over geothermal areas. The use of surface seismic reflection data at Brady's Hot Springs was found useful in building a geological model, but only when combined with other extensive geological and geophysical data. The use of fine source and geophone spacing was critical in producing useful images. The surface seismic reflection data gave no information about the internal structure (extent, thickness and filling) of faults and fractures, and modeling suggests that they are unlikely to do so. Time-frequency analysis was applied to these data, but was not found to be significantly useful in their interpretation. Modeling does indicate that VSP and other seismic methods with sensors located at depth in wells will be the most effective seismic tools for getting information on the internal structure of faults and fractures in support of fluid flow pathway management and EGS treatment. Scattered events similar to those expected from faults and fractures are seen in the VSP reported here. Unfortunately, the source offset and well depth coverage do not allow for detailed analysis of these events. This limited coverage also precluded the use of advanced migration and imaging algorithms. More extensive acquisition is needed to support fault and fracture characterization in the geothermal reservoir at Brady's Hot Springs. The VSP was effective in generating interval velocity estimates over the depths covered by the array. Upgoing reflection events are also visible in the VSP results at locations corresponding to reflection events in the surface seismic. Overall, the high temperature rated fiber optic sensors used in the VSP produced useful results. Modeling has been found useful in the interpretation of both surface reflection seismic and VSP data. It has helped identify possible near surface scattering in the surface seismic data. It has highlighted potential scattering events from deeper faults in the VSP data. Inclusion of more detailed fault and fracture specific stiffness parameters are needed to fully interpret fault and fracture scattered events for flow properties (Pyrak-Nolte and Morris, 2000, Zhu and Snieder, 2002). Shear wave methods were applied in both the surface seismic reflection and VSP work. They were not found to be effective in the Brady's Hot Springs area. This was due to the extreme attenuation of shear waves in the near surface at Brady's. This does not imply that they will be ineffective in general. In geothermal areas where good shear waves can be recorded, modeling suggests they should be very useful for characterizing faults and fractures.

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

Inertial migration of deformable droplets in a microchannel

The microfluidic inertial effect is an effective way of focusing and sorting droplets suspended in a carrier fluid in microchannels. To understand the flow dynamics of microscale droplet migration, we conduct numerical simulations on the droplet motion and deformation in a straight microchannel. The results are compared with preliminary experiments and theoretical analysis. In contrast to most existing literature, the present simulations are three-dimensional and full length in the streamwise direction and consider the confinement effects for a rectangular cross section. To thoroughly examine the effect of the velocity distribution, the release positions of single droplets are varied in a quarter of the channel cross section based on the geometrical symmetries. The migration dynamics and equilibrium positions of the droplets are obtained for different fluid velocities and droplet sizes. Droplets with diameters larger than half of the channel height migrate to the centerline in the height direction and two equilibrium positions are observed between the centerline and the wall in the width direction. In addition to the well-known Segré-Silberberg equilibrium positions, new equilibrium positions closer to the centerline are observed. This finding is validated by preliminary experiments that are designed to introduce droplets at different initial lateral positions. Small droplets also migrate to two equilibrium positions in the quarter of the channel cross section, but the coordinates in the width direction are between the centerline and the wall. The equilibrium positions move toward the centerlines with increasing Reynolds number due to increasing deformations of the droplets. The distributions of the lift forces, angular velocities, and the deformation parameters of droplets along the two confinement direction are investigated in detail. Comparisons are made with theoretical predictions to determine the fundamentals of droplet migration in microchannels. In addition, existence of the inner equilibrium position is linked to the quartic velocity distribution in the width direction through a simple model for the slip angular velocities of droplets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Hydrogeologic Modeling at the Sylvania Corning FUSRAP Site - 13419

A comparative groundwater hydrogeologic modeling analysis is presented herein to simulate potential contaminant migration pathways in a sole source aquifer in Nassau County, Long Island, New York. The source of contamination is related to historical operations at the Sylvania Corning Plant ('Site'), a 9.49- acre facility located at 70, 100 and 140 Cantiague Rock Road, Town of Oyster Bay in the westernmost portion of Hicksville, Long Island. The Site had historically been utilized as a nuclear materials manufacturing facility (e.g., cores, slug, and fuel elements) for reactors used in both research and electric power generation in early 1950's until late 1960's. The Site is contaminated with various volatile organic and inorganic compounds, as well as radionuclides. The major contaminants of concern at the Site are tetrachloroethene (PCE), trichloroethene (TCE), nickel, uranium, and thorium. These compounds are present in soil and groundwater underlying the Site and have migrated off-site. The Site is currently being investigated as part of the Formerly Utilized Sites Remedial Action Program (FUSRAP). The main objective of the current study is to simulate the complex hydrogeologic features in the region, such as numerous current and historic production well fields; large, localized recharge basins; and, multiple aquifers, and to assess potential contaminant migration pathways originating from the Site. For this purpose, the focus of attention was given to the underlying Magothy formation, which has been impacted by the contaminants of concern. This aquifer provides more than 90% of potable water supply in the region. Nassau and Suffolk Counties jointly developed a three-dimensional regional groundwater flow model to help understand the factors affecting groundwater flow regime in the region, to determine adequate water supply for public consumption, to investigate salt water intrusion in localized areas, to evaluate the impacts of regional pumping activity, and to better understand the contaminant transport and fate mechanisms through the underlying aquifers. This regional model, developed for the N.Y. State Department of Environmental Conservation (NYSDEC) by Camp Dresser and McKee (CDM), uses the finite element model DYNFLOW developed by CDM, Cambridge, Massachusetts. The coarseness of the regional model, however, could not adequately capture the hydrogeologic heterogeneity of the aquifer. Specifically, the regional model did not adequately capture the interbedded nature of the Magothy aquifer and, as such, simulated particles tended to track down-gradient from the Site in relatively straight lines while the movement of groundwater in such a heterogeneous aquifer is expected to proceed along a more tortuous path. This paper presents a qualitative comparison of site-specific groundwater flow modeling results with results obtained from the regional model. In order to assess the potential contaminant migration pathways, a particle tracking method was employed. Available site-specific and regional hydraulic conductivity data measured in-situ with respect to depth and location were incorporated into the T-PROG module in GMS model to define statistical variation to better represent the actual stratigraphy and layer heterogeneity. The groundwater flow characteristics in the Magothy aquifer were simulated with the stochastic hydraulic conductivity variation as opposed to constant values as employed in the regional model. Contaminant sources and their exact locations have been fully delineated at the Site during the Remedial Investigation (RI) phase of the project. Contaminant migration pathways originating from these source locations at the Site are qualitatively traced within the sole source aquifer utilizing particles introduced at source locations. Contaminant transport mechanism modeled in the current study is based on pure advection (i.e., plug flow) and mechanical dispersion while molecular diffusion effects are neglected due to relatively high groundwater velocities encountered in the aquifer. In addition, fate of contaminants is ignored hereby to simulate the worst-case scenario, which considers the contaminants of concern as tracer-like compounds for modeling purposes. The results of the modeling analysis are qualitatively compared with the County's regional model, and patterns of contaminant migration in the region are presented. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Cell migration and organization in three-dimensional in vitro culture driven by stiffness gradient: Cell Migration Driven by Stiffness Gradient

Durotaxis, a phenomenon that cells move according to changes in stiffness of the extra cellular matrix, has emerged as a crucial parameter controlling cell migration behavior. The current study provides a simple method to generate three-dimensional continuous stiffness variations without changing other physical characteristics of the extra cellular environment. Using Finite Element simulations, the stiffness and the stiffness gradient variations are evaluated quantitatively, leading to an analysis of the dependence of cell migration behavior on the substrate stiffness parameters. We tested various cell lines on several 3-D environments. In this work, the durotaxis results show that the cell migration velocity does not have any consistency with the stiffness of the substrate, rather it is more related to the stiffness gradient of the substrate. This finding suggests a new mechanism underlying the durotaxis phenomenon, highlighting the importance of the substrate stiffness gradient, rather than the stiffness itself.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Computational study of inertial migration of prolate particles in a straight rectangular channel

Inertial migration of spherical particles has been investigated extensively using experiments, theory, and computational modeling. Yet, a systematic investigation of the effect of particle shape on inertial migration is still lacking. Herein, we numerically mapped the migration dynamics of a prolate particle in a straight rectangular microchannel using smoothed particle hydrodynamics at moderate Reynolds number flows. After validation, we applied our model to 2:1 and 3:1 shape aspect ratio particles at multiple confinement ratios. Their effects on the final focusing position, rotational behavior, and transitional dynamics were studied. In addition to the commonly reported tumbling motion, for the first time, we identified a new logrolling behavior of a prolate ellipsoidal particle in the confined channel. This new behavior occurs when the confinement ratio is above an approximate threshold value of K = 0.72. Our microfluidic experiments using cell aggregates with similar shape aspect ratio and confinement ratio confirmed this new predicted logrolling motion. In this study, we also found that the same particle can undergo different rotational modes, including kayaking behavior, depending on its initial cross-sectional position and orientation. Furthermore, we examined the migration speed, angular velocity, and rotation period as well as their dependence on both particle shape aspect ratio and confinement ratio. Our findings are especially relevant to the applications where particle shape and alignment are used for sorting and analysis, such as the use of barcoded particles for biochemical assays through optical reading, or the shape-based enrichment of microalgae, bacteria, and chromosomes.

42 ENGINEERING↗
DOE OSTI2020

A phase field study of the thermal migration of gas bubbles in UO 2 nuclear fuel under temperature gradient

Phase field models are developed to study the gas bubble migration in uranium dioxide nuclear fuel in which a large temperature gradient exists during the operation. In this work, thermal diffusion mechanism for nanosized gas bubbles and vapor transport process for micron-sized gas bubbles are considered, respectively. In both cases, gas bubbles migrate to the high-temperature area. Due to the velocity difference between leading and trailing edges of the gas bubbles, nanosized gas bubbles are elongated along the temperature gradient direction when thermal diffusion is dominated. Micron-sized gas bubbles are either compressed along temperature gradient direction to form lenticular shape bubbles or elongated along temperature gradient direction, depending on the location of the gas bubbles within the fuel pellet. Initial gas bubble radius has no significant effect on the gas bubble migration velocity for both thermal diffusion and vapor transport mechanisms. We notice that the shape change of the gas bubble due to vapor transport mechanism has no significant effect on the migration velocity. Furthermore, the center cavity formation is also captured by our model which is due to the migration and accumulation of lenticular gas bubbles at the center of the fuel pellet. The modeling results compare well with experimental observations and theoretical analysis in the literature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Ion velocity effect governs damage annealing process in defective KTaO 3

Effects of electronic to nuclear energy losses (S e /S n ) ratio on damage evolution in defective KTaO 3 have been investigated by irradiating pre-damaged single crystal KTaO 3 with intermediate energy O ions (6 MeV, 8 MeV and 12 MeV) at 300 K. By exploring these processes in pre-damaged KTaO 3 containing a fractional disorder level of 0.35, the results demonstrate the occurrence of a precursory stage of damage production before the onset of damage annealing process in defective KTaO 3 that decreases with O ion energy. The observed ionization-induced annealing process by ion channeling analysis has been further mirrored by high resolution transmission electron microscopy analysis. In addition, the reduction of disorder level is accompanied by the broadening of the disorder profiles to greater depth with increasing ion fluence, and enhanced migration is observed with decreasing O ion energy. Since S e (~3.0 keV nm –1 ) is nearly constant for all 3 ion energies across the pre-damaged depth, the difference in behavior is due to the so-called 'velocity effect': the lower ion velocity below the Bragg peak yields a confined spread of the electron cascade and hence an increased energy deposition density. Here, the inelastic thermal spike calculation has further confirmed the existence of a velocity effect, not previously reported in KTaO 3 or very scarcely reported in other materials for which the existence of ionization-induced annealing has been reported. In other words, understanding of ionization-induced annealing has been advanced by pointing out that ion velocity effect governs the healing of pre-existing defects, which may have significant implication for the creation of new functionalities in KTaO 3 through atomic-level control of microstructural modifications, but may not be limited to KTaO 3 .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

A Highly Eccentric Warm Jupiter Orbiting TIC 237913194

The orbital parameters of warm Jupiters serve as a record of their formation history, providing constraints on formation scenarios for giant planets on close and intermediate orbits. Here, we report the discovery of TIC 237913194b, detected in full-frame images from Sectors 1 and 2 of the Transiting Exoplanet Survey Satellite (TESS), ground-based photometry (Chilean–Hungarian Automated Telescope, Las Cumbres Observatory Global Telescope), and Fiber-fed Extended Range Optical Spectrograph radial velocity time series. We constrain its mass to M P = 1.942 -0.091 +0.091 M J and its radius to R P = 1.117 -0.047 +0.054 R J , implying a bulk density similar to Neptune’s. It orbits a G-type star (M ∗ = 1.026 -0.055 +0.057 M ⊙ , V = 12.1 mag) with a period of 15.17 days on one of the most eccentric orbits of all known warm giants (e ≈ 0.58). This extreme dynamical state points to a past interaction with an additional, undetected massive companion. A tidal evolution analysis showed a large tidal dissipation timescale, suggesting that the planet is not a progenitor for a hot Jupiter caught during its high-eccentricity migration. TIC 237913194b further represents an attractive opportunity to study the energy deposition and redistribution in the atmosphere of a warm Jupiter with high eccentricity.

79 ASTRONOMY AND ASTROPHYSICS↗