2026-03-23 | | Total: 5
Localized cavitation in liquids and soft tissues, typically initiated by the rarefaction phase of high-amplitude ultrasound waves, is leveraged in several biomedical applications such as ablation techniques and drug delivery with vaporizing agents. However, safety considerations aimed at avoiding unwanted bubble activity outside the targeted region pose a limit to the maximum allowed peak rarefaction pressure, which on the other hand can hinder the therapeutic efficacy of these techniques. This study shows that a purely compressive shock wave can generate localized, negative pressure and initiate cavitation inside a sub-millimetric perfluorohexane droplet, without requiring any externally applied rarefaction wave. The Gouy phase shift is identified as the physical mechanism responsible for the conversion of positive pressure into tension during shock focusing, and its occurrence is demonstrated through numerical simulations and direct experimental measurements. Comparison of the regions affected by cavitation, visualized \emph{in-situ} by means of high-speed x-ray phase-contrast imaging, with prediction from Classical Nucleation Theory suggests homogeneous nucleation as the underlying mechanism behind bubble formation. The presented findings offer valuable insights into the physics of shock wave propagation which can inspire the development of novel acoustic driving strategies for cavitation generation, facilitating the reduction of negative pressures outside the target region and improving the safety and precision of biomedical treatments.
Large-eddy simulations (LES) require closures for filtered production rates because the resolved fields do not contain all correlations that govern chemical source terms. We develop a graph neural network (GNN) that predicts filtered species production rates on non-uniform meshes from inputs of filtered mass fractions and temperature. Direct numerical simulations of turbulent premixed hydrogen-methane jet flames with hydrogen fractions of 10%, 50%, and 80% provide the dataset. All fields are Favre filtered with the filter width matched to the operating mesh, and learning is performed on subdomain graphs constructed from mesh-point connectivity. A compact set of reactants, intermediates, and products is used, and their filtered production rates form the targets. The model is trained on 10% and 80% blends and evaluated on the unseen 50% blend to test cross-composition generalization. The GNN is compared against an unclosed reference that evaluates rates at the filtered state, and a convolutional neural network baseline that requires remeshing. Across in-distribution and out-of-distribution cases, the GNN yields lower errors and closer statistical agreement with the reference data. Furthermore, the model demonstrates robust generalization across varying filter widths without retraining, maintaining bounded errors at coarser spatial resolutions. A backward facing step configuration further confirms prediction efficacy on a practically relevant geometry. These results highlight the capability of GNNs as robust data-driven closure models for LES on complex meshes.
Turbulent heat and freshwater transport at ice-ocean interfaces controls glacier and iceberg melt rates, yet the underlying physics remains poorly constrained. Parameterizations that assume shear boundary layer scaling are commonly used, which neglects meltwater buoyancy-driven convective processes. Using Direct Numerical Simulations with realistic salt diffusivity, which is critical for representing the thin solutal boundary layer (deltaS ~ 0.4 mm) and resulting convective instabilities, we investigate ice-ocean boundary layer turbulence across varying temperature, salinity, stratification, external velocity, and interfacial slope angles. Our simulations agree with laboratory measurements of melt rate and interfacial temperature. In the absence of external flows, we find no transition from buoyancy-controlled to shear-controlled regimes and convection is important even at near-horizontal slopes. External shear becomes significant only when it is strong enough to thin the thermal and solutal boundary layers, which starts influence melting substantially above background flow speeds of 5 cm/s. Understanding how shear and convection compete to determine the ice-ocean diffusive boundary layer enables accurate melt rate predictions across the parameter space relevant to ice shelves and marine-terminating glaciers.
Turbulent separated flows in axisymmetric expansions can sustain fundamentally different transport organization despite nearly identical mean topology. Using stereo-PIV and time-resolved planar PIV, we compare abrupt $90^\circ$ (step) and gradual $45^\circ$ (wedge) axisymmetric expansions at step height Reynolds numbers of 25000 and 35000. Despite similar reattachment lengths, near-separation turbulence differs, with the wedge exhibiting higher turbulent kinetic energy over a broader shear layer, while the step confines production to a thinner region near the corner, where a secondary vortex weakens momentum and fluctuations. The spatial spectra reveal a pronounced spectral hump in the out-of-plane velocity fluctuations near separation. This feature is consistently observed across all cases and reflects the expansion effects on the redistribution of fluctuation energy associated with the interaction between the separating shear layer and the recirculating flow. Temporal spectra show no geometry-specific dominant frequencies, and space-time correlations indicate similar normalized convection velocities across all cases. The primary effect of geometry, therefore, does not lie in the characteristic scales or transport speeds, but in the spatial organization and persistence of coherence. The step cases exhibit stronger spectral concentration, longer local integral time scales, and a broader distribution of space-time correlations in convection velocities associated with momentum-depleted return flow. Finite-time Lyapunov exponent (FTLE) fields confirm that these differences extend to material transport, as the wedge produces larger and less fragmented deformation regions, while the step yields a more segmented pattern that persists downstream of reattachment.
Porous materials -- natural or engineered -- often exhibit dual pore-network structures that govern processes such as mineral exploration and hydrocarbon recovery from tight shales. Double porosity/permeability (DPP) mathematical models describe incompressible fluid flow through two interacting pore networks with inter-network mass exchange. Despite significant advances in numerical methods, there remains a need for computational frameworks that enable rapid forecasting, data assimilation, and reliable inverse analysis. To address this, we present a physics-informed neural network (PINN) framework for forward and inverse modeling of DPP systems. The proposed approach encodes the governing equations in mixed form, along with boundary conditions, directly into the loss function, with adaptive weighting strategies to balance their contributions. Key features of the framework include adaptive weight tuning, dynamic collocation point selection, and the use of shared trunk neural architectures to efficiently capture the coupled behavior of the dual pore networks. It is inherently mesh-free, making it well-suited for complex geometries typical of porous media. It accurately captures discontinuities in solution fields across layered domains without introducing spurious oscillations commonly observed in classical finite element formulations. Importantly, the framework is well-suited for inverse analysis, enabling robust parameter identification in scenarios where key physical quantities -- such as the mass transfer coefficient in DPP models -- are difficult to measure directly. In addition, a systematic convergence analysis is provided to rigorously assess the stability, accuracy, and reliability of the method. The effectiveness and computational advantages of the approach are demonstrated through a series of representative numerical experiments.