Local Fluid Flow Research Articles

Pyrite trace metal and sulfur isotopic compositions track metalliferous fluid circulation within the Ordovician/Silurian organic-rich black shales in the eastern Sichuan Basin, southwestern China

Depicting metalliferous fluid flows in sedimentary basins has a remarkable implication for understanding the formation and evolution of organic-rich sediments. The Middle-Upper Yangtze region in South China hosts voluminous gaseous hydrocarbons and MVT-type ZnPb deposits. Natural gases are mainly distributed in the Sichuan Basin, yet ZnPb deposits are found in surrounding regions of the basin. Such a unique distribution pattern implies that the interplay between metalliferous fluids and organic-rich sediments may be extensive in the boundary of these two types of deposits. A typical Ordovician/Silurian (Wufeng/Longmaxi formations) organic-rich black shale outcrop occurs in the eastern boundary of the Sichuan Basin. Pyrites are frequently distributed across this section, providing an ideal opportunity to investigate features of metalliferous fluids and their potential impacts on organic-rich sediments. Pyrites associated with high-angle carbonate veinlets are recognized in the studied area, and this group of pyrites (“Group 2”) commonly display planar-laminated morphologies, moderate δ34S values (0.78 ‰–8.86 ‰), and elevated trace metal contents (Ni, Pb, Mn, Mo, Tl, and REE) than those not associated with carbonate veinlets. These features suggest that this group of pyrites may be precipitated via local metalliferous fluid flows. Besides, pyrites with relatively lower trace metal contents can be further divided into two groups, including a group of euhedral/subhedral pyrites with more depleted δ34S values (−18.06 ‰ – -1.15 ‰; “Group 1”) and a group of planar-laminated/cubic pyrites with enriched δ34S values (10.55 ‰–37.62 ‰; “Group 3”). Pyrites of Group 1 and Group 3 may be formed via bacterial sulfate reduction (BSR) and thermochemical sulfate reduction (TSR), respectively. The discovery of fluid-related, trace-metal-enriched pyrites implies that fluid circulation within organic-rich black shales has the potential to remobilize, transport, and re-deposit trace metals. Besides, metalliferous fluid may also promote organic matter maturation within the Sichuan Basin. The outcomes of this study, combined with previous findings of metalliferous fluid flows in the center of the basin and ZnPb mineralization belts surrounding the basin, imply that a widespread Ediacaran-Palaeozoic fluid circulation system may exist in the Middle-Upper Yangtze region.

Dynamic characteristics of the hydrogen injector-ejector unit in a PEM fuel cell system

The ejector, as a fluid-dynamics controlled passive component, is a promising hydrogen recirculation device in proton exchange membrane (PEM) fuel cell systems. Nevertheless, its performance is significantly affected by dynamic operating conditions and the behavior of the injector. This study investigates the dynamic performance of an ejector integrated with a hydrogen injector in a 100 kW PEM fuel cell system. A dynamic computational fluid dynamics (CFD) model of the integrated injector-ejector unit is developed using a dynamic mesh method and validated through theoretical analysis and experimental data. A thorough analysis of the global performance and local fluid flow characteristics of the injector-ejector unit is conducted. The results show that the hydrogen injector, with a throat diameter of 2.00 mm, exhibits an adjustable linear flow range from 0 to 1.62 g/s under an inlet pressure of 2.0 MPa. Additionally, the primary mass flow rate can be linearly regulated by adjusting the valve gap, while the secondary mass flow rate exhibits asynchronous behavior with the primary flow. The velocity and temperature distributions within the injector-ejector unit undergo significant changes under dynamic conditions. These findings contribute to optimizing the design and control strategies of hydrogen supply components in PEM fuel cell systems.

Normal fault architecture, evolution, and deformation mechanisms in basalts, Húsavik, Iceland: Impact on fluid flow in geothermal reservoirs and seismicity

Faults within layered basaltic sequences significantly influence hydrothermal fluid flow in shallow geothermal reservoirs and potentially during CO2 sequestration and storage. Nevertheless, their characterization regarding fault zone architecture, fluid flow, deformation mechanisms, and seismic potential remains underdeveloped. This study addresses this gap by integrating structural and microstructural observations with X-ray diffraction analyses of exposed normal-transtensional faults associated with the seismically active Húsavík-Flatey Fault in the Tjörnes Fracture Zone, Northern Iceland. Our findings demonstrate that the evolution of basalt-hosted normal-transtensional faults progresses through distinct stages: (1) low-displacement fault propagation from pre-existing cooling joints; (2) fault linkage via dilational jogs; (3) damage zone/fault core growth through brecciation and cataclastic processes; (4) shear localization along sharp slip surfaces; and (5) smearing of volcaniclastic interbeds along the principal fault plane. Evidence of shear localization, truncated clasts, and hydrothermal breccias/veins suggests repeated seismic slip events facilitated by overpressured fluids. Conversely, the presence of clay-rich foliated cataclasite indicates aseismic slips during interseismic periods. Slip along fault jogs, bends, geometric irregularities, and orientation changes causes the dilatant opening of the fault planes and extensional horsetail fractures at fault tips. These structures create main tabular zones for lateral movement of hydrothermal fluids parallel to the fault strike in shallow geothermal reservoirs situated in active extensional-transtensional tectonic settings. In addition, the dilational jogs and the intersection of horsetail veins with the hosting faults may define linear zones of high structural permeability and intense localized fluid flow parallel to the σ2 paleostress orientation and finally mineral precipitation. The results of this study can be utilized to improve models of geothermal fluid flow for enhanced recovery in basaltic reservoirs and assess seismic risk in basaltic faults.

Theoretical model for the elastic properties of cracked fluid-saturated rocks considering the crack connectivity

SUMMARY Cracks are a common rock microstructure and have a large effect on elastic properties during wave propagation. The fluid flow between a crack and its adjacent pore space can cause wave attenuation and dispersion. In this work, we introduce a crack connectivity parameter which is meant to improve the expression of local flow by weighting the contributions of fully connected and isolated cracks. We then update the analytical expression for frequency-dependent moduli by modifying the boundary conditions of the linearized Navier–Stokes equation and mass conservation equation. The proposed model contains the effect of cracks and stiff pores, in which the attenuation and dispersion are determined by squirt-flow and stiff-pore relaxations. The resulting model shows the squirt-flow relaxation frequency depends on not only the crack aspect ratio but also the crack connectivity. However, their contributions are different. The crack connectivity has little effect on the attenuation amplitude of shear modulus, but affects the attenuation amplitude of bulk modulus when multiple sets of cracks exist in the rock. The attenuation frequency band is also affected by the crack connectivity. As the crack connectivity deteriorates, the attenuation peak moves to low frequencies. In addition, by comparing the crack connectivity with the fluid viscosity coefficient, it is observed that the crack connectivity only affects the attenuation frequency band of cracks, whereas the fluid viscosity coefficient affects the attenuation frequency bands of cracks and stiff pores simultaneously. Thus, the introduction of crack connectivity is a supplement to the theoretical model of cracked fluid-saturated rocks. It helps understand the local fluid flow induced by seismic waves and provides a reasonable variation analysis of moduli and attenuation, especially for tight reservoirs.

Steady state regimes in lid driven cavity flows of magnetic fluids in the presence of an external magnetic field

In this work, we study the flow of magnetic fluids in a square top lid driven cavity in the presence of an applied magnetic field. The magnetic field is generated by a steady electric current flowing through a wire placed underneath the cavity. The magnetization of the magnetic fluid is described by an evolution equation that combines magnetization relaxation, advection and interaction with the flow via the vorticity. The governing equations are solved via a vorticity-streamfunction formulation implemented in a finite difference scheme. The main governing physical parameters of the examined cavity flow are the Reynolds number, ranging from 50 to 500, the magnetic pressure coefficient, which measures the relative importance between magnetic and inertial forces, ranging from 0 to 1000, and Péclet number which is O(1) and measures the ratio between the magnetization relaxation time and the flow time scale. We have examined the different mechanisms which cause changes in the local fluid magnetization, i.e., flow vorticity and advection. In particular, we examine how vorticity of the flow intensifies deviations of the magnetization of the fluid with respect to the local equilibrium of magnetization. In addition, we identify a complex topological structure of the flow inside the cavity characterized by the magnetic effects that dominate the generation of secondary structures of vorticity induced by the effect of field gradient, orientation and a non-uniformity in the fluid magnetization inside the cavity.

Wave-induced fluid flow and reflection/transmission of seismic waves at a fluid/double-porosity thermoelastic medium interface

This study investigates the wave-induced fluid flow and reflection/transmission of seismic waves at the interface of a non-viscous fluid and a double-porosity thermoelastic (DPT) solid. The analytical reflection/transmission coefficients are calculated for two boundary conditions: wholly sealed and open pores using the displacement potentials and Gauss-elimination method. The wave-induced local fluid flow (LFF) is computed analytically using the transmission coefficients of transmitted waves in a DPT medium, and it found that four compressional waves contribute to wave-induced LFF. Further, the energy partitioning between reflected and transmitted waves is also computed. An energy matrix describes the amount of energy transmitted to the DPT medium. The matrix has five diagonal elements representing the five waves’ energy proportions with different properties. The total of all the non-diagonal elements in the matrix indicates the energy involved in the interaction between transmitted waves. A numerical example is considered to perform the computational analysis of distinct propagation characteristics. Finally, the impacts of incident direction, wave frequency, inclusion radius, and pores fluid viscosity on the wave-induced LFF and energy partitioning are investigated graphically. Finally, energy conservation for both kinds of surface pores is found.

Experimental investigation on a microfluidic U-turn channel for heat transfer applications

In the present work, an experimental set-up designed to allow investigations of the carrier fluid local behaviours in terms of velocity and temperature fields is presented. The study is carried out on a U-turn microfluidic channel, to bear a complete view of the thermo-fluidic system. The experiments have been performed by Infrared Thermal Imaging, which allows deriving the thermal field on the microchannel external walls, and inferring the internal thermal profiles through specific transfer functions in laminar and transitional regimes. Maps of thermal transients allow deriving cooling performances, which are used to identify the local thermal efficiency of the microchannel and relate it to the global efficiency. Velocimetry measurements have been conducted with a Micro Particle Image Velocimetry (μPIV) setup at different flow rates. This analysis has been coupled with thermal results to obtain a description of the effects of local fluid flow phenomena in transitional turbulent regime on the global heat transfer efficiency of a U shaped microchannel. The interactions among thermal and flow fields are specifically related to secondary recirculating flows, close to the turns, where high cooling rates and high magnitude of local velocity fluctuations are measured.

Interfacial Effects of Nanostructured Doubly Reentrant Surfaces on the Evolution of Local Concentration and Fluid Flow in an Evaporating Droplet.

Surface modification, such as bioinspired nanostructured doubly reentrant surfaces that have presented superhydrophobic wettability even under low-surface-tension liquid, is a very promising technology for controlling droplet dynamics, heat transfer, and evaporation. In this article, we investigate the interfacial effects of nanostructured doubly reentrant surfaces on the flow behaviors and local concentration evolution during the evaporation of an ethanol/water multicomponent droplet. Using particle image velocimetry (PIV) and novel aggregate-induced emission-based (AIE) techniques, the flow patterns and local concentration distributions on both hydrophobic and nanostructured doubly reentrant surfaces were probed and compared. It is found that in addition to the established Marangoni flow-dominated stage, transition stage, and buoyancy-induced flow-dominated stage, a new transition stage and a rolling stage for the nanostructured doubly reentrant surface are detected in the late evaporation period. Differences in the local concentration distribution evolution occur depending on the hydrophobicity of the surface on which the droplet is placed. For the hydrophobic surface, a nonuniform local concentration distribution exists consistently, with a high water fraction in a shell-shaped region near the liquid-air interface and a secondary concentration gradient within this shell-shaped region. The concentration distribution on the nanostructured doubly reentrant surface evolves in a more complex manner, with a strip-shaped region of high water fraction forming in the intermediate stage and then reorganized by rolling flow in the late stage. Finally, theoretical analysis combining PIV and AIE visualization results reveals that the variations in droplet concentration distributions on surfaces with different hydrophobicities exert a significant impact on evaporative behaviors. These behaviors, in turn, affect the evolution of the local concentration distribution.

Corner melting in low Pr metals: A study using lattice Boltzmann method

Evolution of melting boundary in the solid–liquid phase-change problem involves nonplanar geometry due to the presence of Rayleigh–Benard convective cells in the melted region. The transition of flow behavior from steady to oscillatory in the melt zone is particularly an interesting phenomenon to study. In this work, numerical investigation of melting of low Prandtl number materials in a square cavity with two adjacent heated walls has been carried out using total enthalpy-based lattice Boltzmann method. The local thermal and fluid flow behavior within a moving melting front has been investigated. The influence of natural convection in the melt zone has been observed for two different cases: (1) heating from a corner formed by the left and the bottom walls and (2) heating from a corner formed by the top and the right walls. The effect of Rayleigh number in the range of Ra = 102–5 × 107 on the convective flow field is evaluated for a typical parametric value of Stefan number of 0.01 and Prandtl number of 0.025. Results show distinct convection rolls at the melt zone for the cases under investigations. Evolution of flow fields in the melt zone has been described by a set of isotherms and streamlines for both the cases. The interface fronts for both the cases are initially of convex shape, and around Fourier number of 0.5 turn into concave shape, as viewed from the direction of the movement of melting front. At higher Ra, the rates of melting for top-right side and left-bottom side heated cavities are nearly the same above Fourier number of about 1.2. The average heat flux from the walls scales with Ran , where n is a constant which varies between 0.98 and 1.11. In the melt region, the instability pattern of flow and thermal field changes with the change of effective melt area. Transition to periodic flow is observed at Ra of 7.5 × 106 for case 1 and 2.2 × 107 for case 2, respectively.

Exploring therapy transport from implantable medical devices using experimentally informed computational methods.

Implantable medical devices that can facilitate therapy transport to localized sites are being developed for a number of diverse applications, including the treatment of diseases such as diabetes and cancer, and tissue regeneration after myocardial infraction. These implants can take the form of an encapsulation device which encases therapy in the form of drugs, proteins, cells, and bioactive agents, in semi-permeable membranes. Such implants have shown some success but the nature of these devices pose a barrier to the diffusion of vital factors, which is further exacerbated upon implantation due to the foreign body response (FBR). The FBR results in the formation of a dense hypo-permeable fibrous capsule around devices and is a leading cause of failure in many implantable technologies. One potential method for overcoming this diffusion barrier and enhancing therapy transport from the device is to incorporate local fluid flow. In this work, we used experimentally informed inputs to characterize the change in the fibrous capsule over time and quantified how this impacts therapy release from a device using computational methods. Insulin was used as a representative therapy as encapsulation devices for Type 1 diabetes are among the most-well characterised. We then explored how local fluid flow may be used to counteract these diffusion barriers, as well as how a more practical pulsatile flow regimen could be implemented to achieve similar results to continuous fluid flow. The generated model is a versatile tool toward informing future device design through its ability to capture the expected decrease in insulin release over time resulting from the FBR and investigate potential methods to overcome these effects.

Direct assessment of the hydraulic structure of the plate boundary at the toe of the Nankai accretionary prism

SUMMARY The Nankai Trough is a locus of slow slip, low-frequency earthquakes, and large magnitude (Mw > 8) earthquakes. It is usually assumed that high pore pressure contributes substantially to earthquake dynamics. Hence, a thorough understanding of the hydraulic regime of the Nankai accretionary prism is needed to understand this diversity of behaviours. We focus on the toe of the accretionary prism by studying data from Hole C0024A, part of the Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) project, that intersected the décollement at 813 m below seafloor (mbsf) about 3 km away from the trench. We contribute to this understanding by innovatively integrating drilling and logging data to derive high-resolution hydraulic profiles along the borehole. A quantitative re-analysis of the variation in the downhole annular pressure monitored during drilling show localized fluid flow from the formation to the borehole in excess of 0.05 m3 s−1, especially in the damage zones at the footwall of the décollement. To validate the fluid flow profile, pore pressure was estimated independently from empirical relationships between pore pressure, porosity and P-wave velocity, obtained from consolidation experiments and Eaton-type methods based on drilling or sonic velocity data. The formation fluids are becoming significantly overpressurized with depth in the few hundred metres above décollement. The hydraulic profile suggests that the core of the décollement acts as a barrier inhibiting upward fluid convection, whereas the damage zone acts as an efficient longitudinal channel able to diffuse high pressure from the deeper part of the subduction.

Acoustically Actuated Flow in Microrobots Powered by Axisymmetric Resonant Bubbles

Bubble‐driven microsystems inspired by the therapeutic use of microbubbles under clinical ultrasound actuation offer innovative remote manipulation in biological settings. Ultrasound‐powered microrobots, benefiting from a distribution of vibrating microbubbles (1–100 μm in diameters), exhibit localized fluid flow for self‐propulsion. Microbubbles also contribute to the design of acoustic metamaterials, where distributed actuation of bubbles enables exotic material properties (e.g., negative refractive index). Herein, a metamaterial‐inspired microrobot design, which exploits the streaming induced by the collective oscillation of microbubble arrays (> 10), is reported. Such a large distribution of bubbles offers twofold advantages: a mesoscale microrobot allowing ease of handling and motion at a considerably low acoustic driving pressure of 50 kPa. In this work, the oscillatory amplitude of a single bubble that acts as a unit cell of the metamaterial‐based microrobot is first characterized. Thereon, this amplitude is related to the flow generated by the collective vibration of all the bubbles close to the theoretically predicted resonant frequency of the bubbles. Finally, the hovering motion of the microrobot induced by the streaming and flow‐assisted debris clearance is demonstrated. Such auxiliary functionalities of the microrobot can be useful for applications like contactless sample extraction in inaccessible biological environments.

Rupture and Afterslip Controlled by Spontaneous Local Fluid Flow in Crustal Rock

AbstractShear rupture and fault slip in crystalline rocks like granite produce large dilation, impacting the spatiotemporal evolution of fluid pressure in the crust during the seismic cycle. To explore how fluid pressure variations are coupled to rock deformation and fault slip, we conducted laboratory experiments under upper crustal conditions while monitoring acoustic emissions and in situ fluid pressure. Our results show two separate faulting stages: initial rupture propagation, associated with large dilatancy and stabilized by local fluid pressure drops, followed by sliding on the newly formed fault, promoted by local fluid pressure recharge from the fault walls. This latter stage had not been previously recognised and can be understood as fluid‐induced afterslip, co‐located with the main rupture patch. Upscaling our laboratory results to the natural scale, we expect that spontaneous fault zone recharge could be responsible for early afterslip in locally dilating regions of major crustal faults, independently from large‐scale fluid flow patterns.

Effect of Multiscale Cracks on Seismic Wave Propagation in Tight Sandstones

AbstractSeismic wave propagation is affected by wave‐induced local fluid flow between stiff pores and multiscale fractures. To investigate this phenomenon, forced vibration (1–100 Hz) and ultrasonic (106 Hz) measurements are performed on two tight sandstones with complex pore geometry in dry and water‐saturated scenarios. Porosity, permeability, and ultrasonic velocities were also measured at different differential pressures. The results indicate that the nonlinear behavior of these properties is strongly influenced by the presence of cracks, and the correlations between the permeability/ultrasonic velocities and porosity are different. A wave propagation model is then developed in which penny‐shaped cracks are inclusions introduced stepwise into a porous medium to describe the wave anelasticity in a wide frequency band at different pressures. As a result, the model provides good agreement with the measured P‐wave velocity dispersion, and the pore aspect ratio spectrum and crack radii are determined. We then compare the estimated crack radii and pore size distributions from nuclear magnetic resonance spectroscopy. Published data of a tight sandstone and a low porosity sandstone in the frequency range (2–200, 106) and (1–3,000) Hz at different differential pressures are also analyzed to validate the model. The aspect ratios, volume fractions, and radii of pores/cracks are used to describe the measured permeability. The present work can provide new insights into the geophysical properties of reservoir rocks with complex pore geometry.

Influence of Liquid Splitting Behavior at Intersections on Infiltration Dynamics in an Unsaturated Fracture Network

AbstractUnderstanding the macro‐scale flow characteristics in the fractured vadose zone is of great importance for subsurface hydrological and environmental applications. Here we develop an idealized fracture network model composed of a series of linked intersections, aiming to reveal the roles of local fluid flow, storage and splitting behaviors at intersections in controlling macroscopic unsaturated flow in a fracture network. By setting different local flow rules, unsaturated flow dynamics and structure in networks are systematically investigated. Numerical simulations suggest that avalanche infiltration mode, that is, sudden release of a large amount of water from the network, emerges spontaneously, and is modulated by the local splitting behavior. In terms of water infiltration structure, we find that uneven liquid splitting at intersections inevitably leads to preferential flow in networks. The preferential paths strongly rely on the competition between gravitational and capillary forces, relating to network structure and water saturation. When flow is dominated by gravity, preferential paths extend along large‐aperture routes; conversely, when flow is dominated by capillarity, the small‐aperture paths are preferred. In terms of infiltration dynamics, we show that the power spectral density of the water saturation time series in the network follows a power law with an exponent of −2 for all simulations with different structural parameters and local flow rules, suggesting a universal self‐organized criticality behavior for unsaturated flow in fractured rocks. The improved understanding from this study may shed new light on the complex flow dynamics for unsaturated flow in fractured media.

Natural fractures in ultra-deep reservoirs of China: A review

In recent years, China has made significant progress in exploring oil and gas in sandstone and carbonate rocks with burial depths of more than 6000 m in the Tarim, Sichuan, and Junggar Basins. These ultra-deep reservoirs have widely developed natural fractures that are important storage spaces and primary flow pathways influencing hydrocarbon enrichment and production. For ultra-deep tight sandstone, fractures caused by tectonic, diagenetic, tectonic-diagenetic, and overpressure processes are widespread. In the foreland thrust belt, which is divided into five fracture domains, fault-related folds affect fracture distribution in ultra-deep sandstones. For ultra-deep carbonate rocks, strike-slip faults control scale, structural style, and linkage evolution of fractures. Fractures formed over multiple periods. Many fractures in Cretaceous sandstones formed recently whereas some of those in carbonate rocks formed in the Paleozoic when strata were buried less than 5000 m. With the increase in burial depth, fracture porosity (and, probably) connectivity is reduced by mineral deposits. Fracture-forming processes include overpressure, present in situ burial and tectonic loading. Other variables modifying fractures include when fractures formed (timing), regional and local structural position, and scale and fluid flow leading to heterogeneity in whether fractures are mineral filled, open, or enhanced by dissolution. Our review provides insights into the origins, characteristics, evolution, and effectiveness of fractures in ultra-deep sandstone and carbonate rocks.

Semi-analytical simulation of seismic waves generated by a point source in 1-D layered half-space of double-porosity media

SUMMARYIn this study, we present a semi-analytical method to simulate the propagation of seismic waves in horizontally stratified double-porosity media. We solve the governing equations in the frequency–wavenumber domain and then compute the time–space domain solutions by Hankel and fast Fourier transforms. We conduct numerical simulations to study the properties of the seismic waves in the double-porosity media. The results show the existence of three P waves, one fast P wave and two slow P waves. The two slow P waves are highly attenuated and can be observed by assuming a low fluid viscosity. By comparing the fast P wave in a single-porosity medium and that in a double-porosity medium, we find that the double-porosity model predicts much higher attenuation of the fast P wave due to the local fluid flow between the background medium and the inclusions. We also study the parameters affecting the attenuation, such as the radius of the spherical inclusions, the viscosity of the pore fluid, the permeability of the background medium and the porosity of the inclusions. Finally, we present a cross-well survey model to investigate the seismic responses in the double-porosity media and compare them with the responses in the single-porosity media. We find that both the amplitude and velocity of the fast P wave decrease with the increase of the volume ratio of the inclusions.

Shape Memory Soft Robotics with Yield Stress Fluids

Biological movement is a source of inspiration for designing soft robots that use fluidic actuation for adaptive gripping and locomotion. While many biological systems use networks of non‐Newtonian fluid for movement, to date, most soft robots use Newtonian fluids or pneumatics. Herein, yield stress fluids to manufacture and operate soft devices are exploited, particularly to create soft actuators that exhibit shape memory. Our soft robots are fabricated through embedded 3D printing where the suspension media is a yield stress fluid. Moreover, this complex fluid is encapsulated and used as the hydraulic transmission fluid. Diagnostic designs are developed to characterize the force and shape memory of the yield stress fluid, and the findings are used to create a gripper common in modern soft robotic applications. The diagnostic devices have deformable reservoirs that demonstrate force response, flow behavior, and deformation profiles dependent on the yield stress features of the transmission fluid. The actuation using the yield stress fluid from the retained suspension media creates avenues for partial shape retention and unconventional expansion from localized fluid flow. Looking toward the future of soft robotics, these fabrication and operational approaches using yield stress fluids can provide greater tunability for applications requiring nonlinear actuation and shape memory.

Ground Surface Temperature History Since the Last Glacial Maximum in Northeast Asia: Reconstructions From the Borehole Geotherms of the International Continental Scientific Drilling Program

AbstractPast ground surface temperature (GST), one of the important aspects of paleoclimate reconstructions, can be inverted from borehole temperature measurements. Here, we report continuous 6,100‐m temperature logs in the International Continental Scientific Drilling Program SK‐2e. We inverted the past GST changes from upper borehole temperature logging (<600 m). Below this depth, localized fluid flow masks the paleoclimate record. Inversions yield an approximately 2 K GST rise since 0.1–0.6 Kyr BP and an approximately 10 K rise since 20 Kyr BP. Assuming a ±5 K influence from the deep groundwater flow, the inverted temperature rise has varied between 8 and 12 K since 20 Kyr BP, which is consistent with previous reports since the Last Glacial Maximum. Our results emphasize the potential of borehole heat‐flow profiles as a record of climate changes and the importance of climate correction for heat‐flow determinations.

Surface waves at a fluid/double-porosity medium interface

SUMMARYWe consider surface-wave propagations at an interface separating a fluid layer and a double-porosity medium embedded with cracks. The theory is based on a generalization of the Biot-Rayleigh model from spherical cavities to penny-shaped cracks randomly embedded into a host medium, where mesoscopic local fluid flow (LFF) plays an important role. We derive closed-form dispersion equations of surface waves, based on potentials and suitable boundary conditions (BCs), to obtain the phase velocity and attenuation by using numerical iterations. Two special cases are considered by letting the thickness of the fluid (water) layer to be zero and infinity. We obtain pseudo-Rayleigh and pseudo-Stoneley waves for zero and infinite thickness and high-order surface modes for finite nonzero thickness. Numerical examples confirm that the LFF affects the propagation at low frequencies, causing strong attenuation, whereas the impact of BCs is mainly observed at high frequencies, due to the propagation of slow wave modes. The crack density mainly affects the level of attenuation, whereas the aspect ratio the location of the relaxation peak. The fundamental mode undergoes a significant velocity dispersion, whose location moves to low frequencies as the thickness increases. In all cases, there also exist two slower surface modes that resemble the two slow body waves, only present for sealed BCs.