2026-08-20 | | Total: 8
The near-field shock structures of underexpanded sonic jets exiting into a subsonic coflow are investigated over a range of nozzle pressure ratio (NPR) and coflow-to-nozzle-exit velocity ratio ($U_c$), representative of a propulsive nozzle in subsonic flight. Time-averaged statistics from fully-resolved axisymmetric simulations and inviscid method-of-characteristics (MOC) analysis are used to understand how coflow alters the shock-cell structures, in particular the Mach-disk formation. It is well established that increasing NPR transitions the centerline reflection from regular (characterized by oblique shocks) to Mach reflection (characterized by a near-normal Mach disk). We find that coflow has the opposite influence: a strong coflow shrinks the Mach disk until it vanishes, reverting Mach reflection to regular reflection, so the NPR for this transition increases with $U_c$. This effect has previously been attributed to a reduction in the jet-boundary inclination at the nozzle lip, which confines the lip Prandtl-Meyer fan to a smaller angle and weakens the embedded shock. We show instead that this inclination is determined by the non-uniform pressure the coflow imposes along the jet boundary, which is the primary driver of the shock-structure transitions in coflowing jets. The non-uniform pressure weakens the boundary-reflected compression waves and orients them at shallower angles, so the embedded shock reflects regularly or fails to form. A simulation-informed MOC analysis with this non-uniform pressure boundary condition reproduces the transition behavior with increasing coflow. Coflow also lengthens the first shock cell linearly, which is accurately estimated by a simple correction to Prandtl classical shock-cell length scaling.
In this work, particle-resolved direct numerical simulations are performed to investigate flow past the five Platonic solids, which represent a progression in particle sphericity with an increasing number of faces. The simulations cover particle Reynolds numbers in the range 0.1 <= Re_p <= 300 and multiple particle orientations relative to the incoming flow. Based on the numerical data, new correlations are developed for the drag, lift, and torque coefficients. The proposed drag correlation explicitly accounts for both Reynolds number and particle orientation, whereas the lift and torque coefficients are represented by orientation-dependent trigonometric and exponential basis functions whose coefficients vary with Reynolds number. The simulations are conducted using the immersed boundary method, and the resulting drag correlation accurately reproduces the numerical data. The lift and torque correlations capture the principal trends observed in the numerical simulations, including the strong dependence on particle orientation. The proposed correlations provide a computationally efficient framework for incorporating orientation-dependent hydrodynamic forces and torques into Euler--Lagrange and point-particle simulations, enabling a more realistic representation and predictions of non-spherical particle transport in multiphase flows.
We experimentally and theoretically study the hydrodynamic analog of the classical brachistochrone problem: the {\it time-} and {\it energy-minimizing} paths for a spherical particle rolling down an incline within a viscous fluid. We show that in the presence of viscous dissipation, the paths of minima diverge from the classical cycloid, into curves of opposing curvature for time and energy, and are characterized by an effective dimensionless parameter, $St_p$, representing the ratio of the particle's viscous response time scale to its gravitational time scale. Using a generalized variational framework, we show that the fastest path reduces to {nearly straight ramps}, however, beginning and terminating in localized cycloids of curvature, $κ_c \sim St_p^{-2}$. Remarkably, the path of fastest descent on a given energy budget requires navigating a non-monotonic path ({\it``S-shaped''}) with an interior point of inflection. Our findings reveal a unification of temporal and energetic optimality for transport through dissipative media, and expand the celebrated brachistochrone solutions to the hydrodynamic regime.
Wind turbines operating downstream of others in a farm are routinely exposed to waked inflow, with reduced mean velocity and elevated turbulence driving power deficits and additional structural fatigue. The direct effect of wakes on blade-level structural loading remains under-explored experimentally, owing partly to the sparse spatial coverage of conventional point-based strain sensors. Here, we present a wind-tunnel study of wake-induced blade dynamics using two $1\,\mathrm{m}$-diameter turbine models, in which one blade of a downstream turbine ($WT_2$) is instrumented with distributed Rayleigh-backscattering fibre-optic strain sensors, providing spatially continuous strain measurement across the blade span. By changing the relative position of the upstream turbine ($WT_1$) to the downstream, waked turbine $WT_2$ across the streamwise and spanwise extent, we map power output, spanwise strain, and accumulated representative fatigue relevant loading across the wake profile. Full wake impingement suppresses blade loading through the associated velocity deficit, while partial wake overlap generates the strongest load intermittency and highest relative fatigue relevant loading, despite an intermediate power recovery. A combined performance-to-loading metric shows this partial-wake regime offers the least favourable trade-off between energy yield and structural loading. These results show that minimising partial-wake exposure, not only mean velocity deficits, should be a design consideration for wind-farm layout and turbine spacing.
Granular flows are ubiquitous in natural and industrial systems, yet their complex dynamics remain difficult to characterize. For inverse problems involving unknown inlet, outlet, and wall boundary conditions, where CFD simulations are challenging, reconstructing complete flow fields from sparse observations constitutes a challenging inverse problem. In this study, a physics-informed neural network framework driven by both physical mechanisms and measurement data is developed to reconstruct the steady-state full-field distribution of granular flows in a pipe. The proposed approach integrates sparse measurement data with governing equations and constitutive relations and is trained using high-fidelity datasets generated by CFD solutions of a continuum model. The framework incorporates a dimensionless loss formulation, physics-informed initialization, dynamic global weighting, and a locally weighted granular temperature data-loss strategy. These treatments enable accurate reconstruction of the complete flow-field evolution. This work establishes a robust methodological framework for flow-field reconstruction in complex granular flow systems.
Predicting flow-induced sound from moving and deforming bodies is computationally demanding because the near-field hydrodynamics and the far-field acoustics require substantially different spatial resolutions and domain extents. A fully parallel hybrid framework is developed to address this disparity by coupling an incompressible Navier-Stokes solver to an acoustic perturbation equation (APE) solver on independently generated, non-conforming Cartesian grids. A sharp-interface ghost-cell immersed boundary method, with radial-basis-function reconstruction, imposes the boundary conditions for complex moving geometries on both grids. The converged flow field supplies the acoustic source through a one-way, precomputed parallel interpolation operator. This arrangement confines the flow grid to the body and wake while allowing the acoustic grid to extend independently into the far field. The framework is validated for Gaussian-pulse propagation, pulse scattering by a rigid cylinder, tonal sound from flow past a cylinder, and radiation from a traveling wavy foil. The predicted waveforms, wavelengths, pressure amplitudes, and radiation patterns agree closely with analytical solutions and published reference data. Applications to eel and Jack fish locomotion, a four-eel school, a manta ray, and a harbor seal further demonstrate the treatment of realistic three-dimensional morphologies, large boundary deformation, and multiple interacting swimmers. The results resolve morphology-dependent acoustic signatures and interference-driven changes in far-field directivity without requiring the flow grid to span the acoustic far field.
Liquid bridges suspended between two parallel horizontal cylinders are studied using experiments, reduced-order mathematical modelling and numerical simulations. Both non-electrified and electrified configurations are considered, with the cylinders acting as electrodes between which a potential difference can be applied. The initial focus is on equilibrium bridge shapes, while the dynamics is examined only in the non-electrified case. In the experiments, transformer-oil bridges are investigated and modelled as perfect dielectrics. The electric field counteracts the gravitational effects and causes the bridges to rise and become flatter, with shapes that are well described by Young--Laplace-type equations containing non-local electric-field contributions evaluated using a boundary-element method. Bifurcation diagrams of these solutions are constructed by pseudo-arclength continuation to characterise the dependence of bridge shapes on liquid volume and electric-field strength. In the absence of an electric field, the transient relaxation to equilibrium is analysed using a reduced-order model developed using Onsager's variational principle, with comparison to direct numerical simulations. The steady states predicted by the reduced-order model agree closely with those of the full formulation over a wide parameter range, and the dynamics is well captured in the overdamped regime.
This study compares three reduced-order modeling (ROM) approaches for flow simulations of generalized Newtonian fluids described by the Carreau rheological model. All three methods rely on offline snapshot generation in the rheological parameter space using the full-order model (FOM), followed by a proper orthogonal decomposition (POD) of the snapshot matrix to obtain a reduced basis, but they differ in how they reconstruct the solution for new parameter values in the online phase. The three ROM approaches examined are: (i) intrusive Galerkin projection onto the reduced basis with full operator reassembly (ROM-FULL), (ii) intrusive hyper-reduced Galerkin projection using the discrete empirical interpolation method with GappyPOD for the nonlinear term (ROM-DEIM), and (iii) a non-intrusive interpolation approach using radial basis function interpolation (ROM-RBF). We demonstrate these three ROM approaches on two benchmark flows: a lid-driven cavity and a sphere settling in a closed container, spanning boundary-driven and force-driven flows. ROM-FULL achieves the highest accuracy but requires reassembling the full-order nonlinear operator during the online phase, whereas ROM-RBF is fully non-intrusive, and its accuracy is closely tied to data availability and deteriorates outside the training data range. ROM-DEIM offers a balance between efficiency and accuracy, even when data are sparse. The results provide guidelines for selecting an appropriate ROM strategy based on solver accessibility, computational efficiency, and desired accuracy.