Detailed Velocity Research Articles

Flow state transition induced by emergence of orbiting satellite eddies in two-dimensional turbulent Rayleigh–Bénard convection

We report a numerical investigation of a previously noticed but less explored flow state transition in two-dimensional turbulent Rayleigh–Bénard convection. The simulations are performed in a square domain over a Rayleigh number range of $10^7 \leq Ra \leq 2 \times 10^{11}$ and a Prandtl number range of $0.25 \leq Pr \leq 20$ . The transition is characterized by the emergence of multiple satellite eddies with increasing $Ra$ , which orbit around and interact with the main vortex roll in the system. Consequently, the main roll is squeezed to a smaller size compared with the domain and wanders around in the bulk region irregularly and extensively. This is in sharp contrast to the flow state before the transition, which is featured by a domain-sized circulatory roll with its vortex centre ‘condensed’ near the domain's centre. Detailed velocity field analysis reveals that there exists an abrupt increase in the energy fluctuations of the Fourier modes during the transition. Based on this phase-transition-like signal, the critical condition for the transition is found to follow a scaling relation as $Ra_t \sim Pr^{1.41}$ where $Ra_t$ is the critical Rayleigh number for the transition. This scaling relation is quantitatively explained by a phenomenological model grounded on the bistability behaviour (i.e. spontaneous and stochastic switching between the two flow states) observed at the edge of the transition. The model can also account for the effects of aspect ratio on the transition reported in the literature (van der Poel et al., Phys. Fluids, vol. 24, 2012).

Three-dimensional CFD modeling for hollow fiber gas separation membrane modules based on a dual porous media model

Mathematical modeling is crucial to optimizing the design of the hollow fiber membrane (HFM) module. A Computational Fluid Dynamics (CFD) model enables efficient three-dimensional (3D) simulations of HFM modules, which integrates a dual porous media approach, significantly reducing computational costs. Validation against alternative simulations and experimental data confirms its prediction capability (<7% deviation). Comparative analyses demonstrate the superior reliability of the 3D approach over one-dimensional simulations. Detailed velocity, pressure, and concentration profiles for the shell and tube sides reveal optimization opportunities such as dead zones within the module. Aerometric studies show significant pressure drops (>30%) for smaller fiber diameters (<100μm) on the permeate side. The model's adaptability suggests broader applications in membrane gas separation processes with modifications. This research highlights the efficacy of CFD modeling in enhancing the design and optimization of HFM modules, highlighting the cost savings of the dual porous media approximation.

2D Sedimentary structures at the southeast margin of the Tarim Basin, China, constrained by love wave ambient noise tomography

Summary Imaging the detailed structure of sedimentary basins is crucial for natural resource exploration and essential for better analysis and correction of sediment responses when studying deeper interior earth structures using seismic data. In the Tarim Basin, previous studies on the sedimentary structures are mostly obtained by active-source seismic surveys, which can provide high-resolution underground interface information but are highly costly or environmentally unfriendly. In this paper, ambient noise tomography, an efficient and economical method based on background vibration, is employed to construct the sedimentary velocity structure at the southeast margin of the Tarim Basin. Based on ambient noise data collected from a linear dense short-period seismic array, we extract Love wave signals from T-T component cross-correlation functions (CCFs) and measure Love wave dispersion curves at a period band of ∼0.3–11.5 s. Then, we utilize a one-step direct surface wave tomography method to image a fine 2D sedimentary shear wave velocity structure with a depth reaching 10 km. Our results reveal a clear layered sedimentary structure, the paleo Tadong uplift, and the thrust Cherchen fault. Our study provides reference sedimentary velocity models for the southeastern Tarim Basin, focusing on depths shallower than 10 km. This model is intended to serve as crucial input for studies requiring detailed shallow sedimentary velocity data. Moreover, our research demonstrates that the application of the ambient noise tomography method with dense arrays has great potential for enhancing resource exploration efforts in sedimentary basins.

Biomimetic Wings for Micro Air Vehicles.

In this work, micro air vehicles (MAVs) equipped with bio-inspired wings are investigated experimentally in wind tunnel. The starting point is that insects such as dragonflies, butterflies and locusts have wings with rigid tubular elements (corrugation) connected by flexible parts (profiling). So far, it is important to understand the specific aerodynamic effects of corrugation and profiling as applied to conventional wings for the optimization of low-Reynolds-number aerodynamics. The present study, in comparison to previous investigations on the topic, considers whole MAVs rather than isolated wings. A planform with a low aperture-to-chord ratio is employed in order to investigate the interaction between large tip vortices and the flow over the wing surface at large angles of incidence. Comparisons are made by measuring global aerodynamic loads using force balance, specifically drag and lift, and detailed local velocity fields over wing surfaces, by means of particle image velocimetry (PIV). This type of combined global-local investigation allows describing and relating overall MAV performance to detailed high-resolution flow fields. The results indicate that the combination of wing corrugation and profiling gives effective enhancements in performance, around 50%, in comparison to the classical flat-plate configuration. These results are particularly relevant in the framework of low-aspect-ratio MAVs, undergoing beneficial interactions between tip vortices and large-scale separation.

The Gulf Stream Structure and Meandering Based on the CTD and SADCP Measurements in 1989–1990 and 2014–2015

An review of field studies in the Gulf Stream region carried out by the authors in two periods with a break of 25 years is presented to summarize the results. The studies in the early period included hydrographic surveys in the area of a southern meander of the Gulf Stream (1989) and in the area of dividing of a single jet of the current into separate branches: in the Gulf Stream delta (1990). The second, recent, stage includes on-route surveys with SADCP profiler in 2014–2015 while crossing the meandering Gulf Stream at mid-latitudes to study its detailed high-resolution velocity field structure.

Experimental investigation on aerodynamic noise and flow structures of a vibrissa-shaped cylinder

The noise mitigation effect of bio-inspired geometries has attracted growing attention from both research and industry, such as the vibrissa-shaped cylinder derived from the harbor seal. Experiments were conducted to investigate the far-field noise and the near-field wake of the flow past a vibrissa cylinder, a circular cylinder, and an elliptical cylinder at Re=3.6×104, in the subcritical flow regime. The frequency characteristic of the far-field acoustic pressure and the near-field velocities are analyzed. The mean and fluctuating velocities, dominant flow modes from proper orthogonal decomposition in both vertical and horizontal planes as well as the time-frequency behavior of the dominant flow structures from wavelet transform are also presented to better understand the wake dynamics and the direct relation of these flow structures with the far-field noise. The vibrissa cylinder reduces the overall sound pressure level by 13.2 dB and 8.3 dB compared with the circular and the elliptical cylinders, respectively, with a remarkable attenuation of the tonal peak associated with vortex shedding. From the detailed velocity measurements in multiple wake planes, it is clearly observed that vortex shedding of the vibrissa cylinder is weaker in strength and significantly less coherent in the spanwise direction than the other two cylinder cases, accompanied by more transient changes. The results also reveal the distinct flow behaviors behind the nodal and saddle planes of the vibrissa cylinder, further contributing to this three-dimensional vortex shedding. Consequently, the power spectral density of the tonal peaks associated with the vortex shedding in both near-field velocities and far-field acoustic pressure are attenuated, leading to a lower noise level. Understanding the detailed flow dynamics of the vibrissa cylinder will provide useful insights into more efficient bio-inspired cylinder designs in noise mitigation and wake control.

Wind tunnel experiment on cross-ventilation of generic isolated building with various roof shapes: Impact of roof pitch and eaves

This study examined the effects of roof shape, specifically roof pitch and eaves, on cross-ventilation performance through wind tunnel experiments that measure static pressure on walls and the detailed wind velocity distribution for isolated building models with basic roof geometries. A flat roof, a gable roof, and two types of shed roof with different orientations were compared. The buildings with a flat or shallow roof pitch (15°) showed no significant change in the pressure difference between windward and leeward facades. However, a steeper roof pitch (30°) caused a significant increase in pressure difference, especially for the shed roof that ascends in the wind direction (shed-A roof). This was caused by larger projected building areas creating stronger recirculation zones behind the building, leading to lower pressure on the leeward facade. The shed-A roof with a pitch of 30° had the best ventilation performance among the tested cases. Eaves generally reduced pressure differences, except for the shed-A roof. The effect of eaves was weaker than reported in previous studies. The experimental results obtained in this study can serve as a validation database for a more detailed study of this problem using computational fluid dynamics.

Full waveform inversion of ocean bottom seismometer data from the oceanic Pacific Plate of the Japan trench

SUMMARY The geometrical and seismic structure of the Pacific oceanic Plate of the Japan Trench is essential for the understanding of earthquake activity in the area. Ocean Bottom Seismometer (OBS) data can be used, via ray-based tomography, to obtain estimates of properties such as crust thickness and structure, hydration and depth to the Moho boundary. The spatial resolution of these properties can be substantially improved by using the full waveform inversion (FWI) method. Most OBS data in this area are acquired with a sparse receiver spacing of 5–6 km, whereas FWI is assumed to work best with denser (1–2 km) receiver spacing. We show that FWI can be adapted to sparsely sampled data with better resolution than traveltime tomography. Using a 500 km long OBS longitudinal profile from the Japan Trench we obtain a detailed velocity structure of the crust, a better definition of the Moho boundary, a well-defined low-velocity layer in the lower crust and a clear spatial definition of areas with velocity inversions.

Enhancing Accuracy In 3D Defocusing Particle Tracking Using 3D-Printed Ramp

Three-dimensional defocusing particle tracking velocimetry (3D DPTV) is an emerging tool for investigating fluid motion within microfluidic devices, where observation is typically limited to a single direction. This technique enables the reconstruction of 3D flow fields by analyzing the defocused particle images. The accuracy of 3D DPTV significantly relies on the collection of a proper calibration image stack containing particle images with precisely known defocus distances. One major challenge is to insert a microscale calibration target inside the microfluidic devices due to the small length scales of microfluidic devices. While the alternative approaches such as using tilted glass slides or conducting z-scans of particles fixed at the bottom wall, face several limitations. Here, we present the application of a two-photon polymerization (2PP) based method to manufacture a microscale reference ramp as the calibration target. We demonstrate the fabrication of a precise reference ramp structure that spans the desired calibration range to enhance accuracy in 3D defocusing particle tracking. We experimentally and computationally compare the result with the commonly employed methods to obtain the calibration image stack. The performance of the 2PP-based calibration ramps is compared with the z-scanning method using the General Particle tracking (GPTV). The results demonstrate that the z-scanning method underestimates the defocused distances compared to the ground truth provided by the ramp approach, with the error increasing as the defocused distance grows. The integration of microscale calibration targets using 2PP fabrication method offers a significant advancement in the reliability and accuracy of 3D DPTV measurements in microfluidic systems. By enabling precise, device-specific calibration, this approach can enhance the fidelity of flow field reconstructions, benefiting a wide range of microfluidic applications that rely on detailed velocity field information.

Subsurface velocity structures at the Egyptian seismological network stations retrieved by diffuse field assumption for Earthquakes

Retrieving detailed subsurface velocity structures down to the seismic bedrock at any given site is a crucial step to delineate the site amplification factors accurately and precisely. The present research work contributes first new estimations for detailed velocity structures down to the seismic bedrock beneath the Egyptian National Seismological Network (ENSN) stations, which are distributed in Egypt nationwide. We used the diffuse field assumption for earthquakes to reproduce the horizontal to vertical spectral ratios (EHVSR) at these stations. We accepted waveform database of 424 earthquakes recorded at 75 ENSN stations. After achieving the inverted subsurface velocity structures, we establish site-specific frequency-depth regression and map the VS30 and seismic bedrock depth beneath Egypt. Because of comparability regarding the seismic site class of B and C, the regression coefficients of the newly established frequency-depth regression exhibit similarity with those achieved from previous regressions in Japan. Furthermore, we observe modest consistency between the seismic bedrock depths and the various geologic features, particularly agreement between basin-shape seismic bedrock depths and the existence of Cretaceous and Jurassic extensional basins. Our findings suggest that the Precambrian basement rocks can be interpreted as the seismic bedrock in Egypt. One of the most significant obstacles in the present work is the low-dense distribution of ENSN stations nationwide. However, the achieved results raise new questions and challenges regarding precise and accurate future estimations for site amplification calculations at different localities in Egypt, particularly for urban planning by engineers.

Stress drops of intermediate-depth intraslab earthquakes beneath Tohoku, northern Japan

We calculated stress drops for 2875 small intraslab earthquakes at intermediate depths beneath Tohoku, Japan. We applied an S-coda-wave spectral ratio method to almost 900,000 event pairs. Detailed velocity values for the oceanic crust (OC) were adopted from previous observational studies. The median stress drops in the OC are about half those in the oceanic mantle (OM). The median stress drop for earthquakes in the OC decreases from depths of 70 to 120 km and increases from 120 to 170 km. Our preferred interpretation is that the rigidity in the OC decreases and then increases with depth due to combined effects of the dehydration associated with the eclogite formation and the increasing temperature with depth. These depth variations are consistent with results of a similar study beneath Hokkaido. The median stress drops in the oceanic plate beneath Tohoku are generally smaller than those beneath Hokkaido. Previous studies imaging the seismic structure at shallow depths and b-value analyses of intraslab earthquakes indicate that the near-trench region of the oceanic plate off Tohoku is more hydrated than that off Hokkaido. Taken together, these results suggest that differences in the degree of hydration of the oceanic plate in the near-trench regions could produce the different behaviors of stress drops of intermediate-depth earthquakes observed in Tohoku and Hokkaido.Graphical

An experimental investigation of constrained melting of a phase change material (PCM) in circular geometries, part I: Velocity characterization

This study is focused on the detailed characterization of the transient velocity fields of a phase change material (PCM) enclosed in a geometry with a circular cross-section subjected to uniform external heating. Particle image velocimetry (PIV) was used to study the velocity behavior related to the buoyancy-driven convection. The PCM was initially subcooled, and the walls suddenly heated to exceed the melting temperature. Experiments were conducted using Rubitherm RT26 in a quasi-2D circular enclosure at three applied heating conditions. Detailed velocity fields show strong symmetric convective vortices at the sides of the PCM chamber, and a three-dimensional recirculation was detected near the bottom. These structures strongly influence the melting behavior. The results suggest self-similarity in flow behavior between heating conditions. Detailed transient temperature fields and an extended discussion of the associated heat transfer behavior are presented in the companion paper (Part II).

Investigation of fluid flow during flow boiling inside a horizontal rectangular channel with single-sided heating using particle image velocimetry

Subcooled flow boiling is a highly efficient cooling systems for thermal management systems. This study explores the intricate dynamics of subcooled flow boiling within a horizontal channel, investigating the impact of vapor generation on liquid-phase velocity using Particle Image Velocimetry (PIV) and advanced image processing techniques. Four mass flow rates ranging from 5–20 g/s with subcooled inlet conditions are investigated in a rectangular channel with single-sided heating. Three regions of interest along the heated channel are investigated for instantaneous PIV analysis. The PIV system captures detailed velocity profiles, illustrating the impact of varying mass flow rates and heat flux levels on flow behavior. Vapor masking techniques are introduced to enhance the precision of PIV data by mitigating interference from the vapor phase. Results demonstrate the influence of vapor bubbles on flow resistance, revealing non-uniform velocity distributions and turbulence near the liquid–vapor interface. The study emphasizes the critical role of inertia and buoyancy forces in shaping the velocity profiles. Moreover, the investigation sheds light on the effects of flow rates on the interfacial behaviors, hinting at a transition point between 10 and 15 g/s. In summary, this research contributes valuable insights into the nuanced dynamics of flow boiling, laying the foundation for future studies on turbulence, heat transfer, and phase-change phenomena in two-phase thermal management systems.

Uppermost Mantle Pn Velocity and Anisotropy Structures Beneath the Sakhalin–Kuril–Kamchatka Region

AbstractIn this study, we used the Pn tomography method to obtain detailed velocity and anisotropy structures of the uppermost mantle beneath Sakhalin–Kuril–Kamchatka region for improving the understanding of plate subduction, arc–arc collision, and volcanism. We found low Pn velocities beneath volcanoes and areas characterized by pronounced tectonic activity and high Pn velocities with strong anisotropy in the subducting plate. Low Pn velocity anomalies beneath southern Sakhalin connected the low‐velocity anomalies in the mantle and crust, indicating the ascent of fluid or melt, and may provide a magmatic source for the rear‐arc Rishiri volcano. In the absence of plate subduction, low‐velocity anomalies and north–south Pn anisotropy manifested beneath northern Kamchatka, revealing the lateral propagation of mantle flow beneath this northern region. We suggest that the eastern boundary of the slab window at the Kamchatka–Aleutian junction is likely to be located near the Komandorsky Islands.

Faster determination of an initial velocity model for full-waveform inversion based on simulated annealing

Full-waveform inversion (FWI) is a prominent method commonly used to create detailed velocity models of the subsurface. However, because it relies on gradient methods, it suffers from the limitation of getting trapped in local minima. To avoid this problem, FWI needs to start from an initial velocity model that lies in the same region of convexity as the global minimum. Global optimization methods, such as simulated annealing, can be used to find such an initial velocity model. However, the iterative process of simulated annealing entails high computational cost, first because of the large number of iterations needed to explore the variables space, and second because of the large number of simulations needed to adjust the simulation parameters. Although the first issue is a common concern of previous work, the second issue is usually neglected, relying on trial and error. We introduce an approach to tackle the issues and then apply it within a simulated annealing framework we create. Our experiments using the Marmousi data yielded promising results when compared with previous work. We find that our approach almost eliminated the computational effort to fine-tune several simulation parameters. In addition, the number of iterations needed to explore the variables space is reduced by two orders of magnitude. FWI is able to find a detailed velocity model with high quality when using the initial velocity model generated by the method.

The Magmatic Patterns Formed by the Interaction of the Hainan Mantle Plume and Lei–Qiong Crust Revealed through Seismic Ambient Noise Imaging

Magmatism on continental lithospheres induced by mantle plumes is more complex compared to oceanic intraplate volcanism owing to the heterogeneous nature of continental crustal and lithospheric structures. Substantial evidence points to the deep-oriented Hainan mantle plume beneath the Lei–Qiong region, the southernmost of the South China block. In this study, we present a detailed shear wave velocity model of the crust and uppermost mantle in the Lei–Qiong volcanic region, derived from 3-year seismic data (2016–2018) from 34 stations and the use of the ambient noise tomography method. An evident columnar low-velocity anomaly was imaged in the crust and uppermost mantle beneath the Wushi Sag (WSS), Beibu Gulf, potentially suggesting that the center of either one branch or the entirety of the Hainan mantle plume impacts the crust here. This low-velocity anomaly is overlaid by a local Moho deepening, indicative of underplating beneath the existing crust. The Maanling–Leihuling Volcanic Field (MLVF) in northern Hainan Island, previously considered the center of the hotspot, does not exhibit such distinct velocity anomalies. Instead, subtle lower crustal anomalies beneath the MLVF are linked with the upper mantle low-velocity zone beneath the WSS. Additionally, the high-conductivity bodies beneath the MLVF indicate lateral magma transport. Earthquake swarms and deep-seated seismic events beneath the WSS further support the presence of magmatic processes. This study indicates that in the Lei–Qiong region, the interaction of the continental crust with the mantle plume centered in the WSS results in magma exhibiting both vertical ascent and lateral migration, leading to a dual low-velocity shear wave pattern in the upper crust, which significantly influences the surface volcanic activity.

A lower crust shear zone facilitates delamination and continental subduction under the Apennines

Physical properties and structure of the lithosphere are the first step to constrain the evolution of mountain belts. Here we show detailed shear wave velocity profiles of the lithosphere in the Apennines that clarify a controversial aspect of continental subduction: the intricate mechanism of crust delamination from the downgoing plate. From the analysis of complete and dense teleseismic Receiver Function data set, we find that the delamination of the continental lithosphere is favored by the development of a low seismic shear wave velocity zone in the middle-lower crust. We observe a double Moho below the external portions of the present mountain range, suggesting the progressive formation of the shallow interface. The delamination edge is located in the forearc, far eastward than expected, implying that the re-equilibration of the thermal unbalance, generated by the mantle substitution, may last 10-7 Myr.

In-situ material damping measurements using the crosshole seismic method

Crosshole seismic testing is commonly used to determine detailed and reliable wave velocity profiles. Material damping measurements can also be performed in crosshole seismic testing, but these measurements require a window function to maintain similar amplitude shapes over the frequency range of interest for the calculation. In-situ seismic testing, including downhole and seismic cone penetration testing (SCPT), also requires a window function for the material damping calculations. However, the length of the window function has the largest influence on the material damping calculation. Therefore, the half-power bandwidth method is introduced to determine the material damping ratio without the window function in crosshole testing by replicating the unconfined, free-free, resonant column (Fr-Fr) test. To demonstrate the influence in the length of the window function, the first cycle of the signal is used in the material damping calculation. The half- power bandwidth method is verified using synthetic signals and field data. Two sets of crosshole and downhole tests were performed, one at a backfill test pad and the second at a Hornsby Bend research site operated by the University of Texas at Austin. The in-situ material damping ratio calculated from these two sets of crosshole tests using the half-power bandwidth method are compare with the Spectral Ratio Slope (SRS) method applied to the downhole testing at these sites. The material damping calculated from the half-power bandwidth method using the full signal results in a reliable and precise damping value compared to the SRS method applied to the downhole data.

Validating AU Microscopii d with Transit Timing Variations

AU Mic is a young (22 Myr), nearby exoplanetary system that exhibits excess transit timing variations (TTVs) that cannot be accounted for by the two known transiting planets nor stellar activity. We present the statistical “validation” of the tentative planet AU Mic d (even though there are examples of “confirmed” planets with ambiguous orbital periods). We add 18 new transits and nine midpoint times in an updated TTV analysis to prior work. We perform the joint modeling of transit light curves using EXOFASTv2 and extract the transit midpoint times. Next, we construct an O−C diagram and use Exo-Striker to model the TTVs. We generate TTV log-likelihood periodograms to explore possible solutions for d’s period, then follow those up with detailed TTV and radial velocity Markov Chain Monte Carlo modeling and stability tests. We find several candidate periods for AU Mic d, all of which are near resonances with AU Mic b and c of varying order. Based on our model comparisons, the most-favored orbital period of AU Mic d is 12.73596 ± 0.00793 days (T C,d = 2458340.55781 ± 0.11641 BJD), which puts the three planets near 4:6:9 mean-motion resonance. The mass for d is 1.053 ± 0.511 M ⊕, making this planet Earth-like in mass. If confirmed, AU Mic d would be the first known Earth-mass planet orbiting a young star and would provide a valuable opportunity in probing a young terrestrial planet’s atmosphere. Additional TTV observations of the AU Mic system are needed to further constrain the planetary masses, search for possible transits of AU Mic d, and detect possible additional planets beyond AU Mic c.

Effect of leading-edge geometry on flow beneath a simulated ice block

ABSTRACTThe present article conducted detailed velocity measurements beneath simulated ice blocks with different leading-edge geometries (round, rectangular, upward and downward triangular). The results examined flow separation at the leading-edge, vortex generation, and subsequent vortex propagation. The instantaneous velocity field depicts an unsteady flow dominated by large-scale vortices, with the Kelvin–Helmholtz type instability dominating the shear layer interface. The mode of vortex generation and propagation was influenced by the geometry of the leading edge. These vortices were dominant for the rectangular and upward triangular configurations. Propagation of these vortices creates low-pressure zones beneath the simulated ice block, which can affect the ice block stability. For all ice blocks, the mean flow accelerated, due to flow separation, and this can result in fluctuations in the dynamic pressure field. These events can lead to greater under-turning moments, as well as the interfacial melting of the ice.