2026-05-08 | | Total: 14
Heat dissipation is critical in modern engineering systems. Polymer additives offer a potential route to improve fluid-based cooling. Here, we study elasticity-enhanced heat transfer in two-dimensional, thermally-stratified Poiseuille flow. At Reynolds numbers, $Re$, $\lesssim 1000$, we observe two types of linearly unstable modes: the recently identified elasticity-induced centre mode (Khalid et al., J. Fluid Mech. 915, 2021) and the classical buoyancy-driven convective mode (Kelly, Adv. Appl. Mech. 31, 35-112, 1994). Direct numerical simulations show that the centre mode develops into a nonlinear `arrowhead' state but yields negligible heat transfer enhancement (typically $\approx 0.03\%$ increase compared to the conductive state). By contrast, polymers can enhance the heat flux associated with the convective mode by up to $1100\%$. The nonlinear convective-mode states take the form of either periodic orbits or travelling waves, and are dominated by hook-like polymer-stress structures that can attach to the walls. The unattached hooks act as `speed bumps' that reduced streamwise velocity and promote wall-normal motion, whereas wall-attached hooks form effective `polymer walls', reorganising the flow into strong counter-rotating rolls and triggering the extreme-enhancement regime. The elasto-buoyant nature of these states is confirmed by perturbation kinetic energy budgets, which show that polymer and buoyancy sustain the states synergistically. The wall-attached hooks enable rapid thermal equilibration but impose a large hydraulic penalty, making them suitable for process streams requiring fast temperature adjustment. Unattached hooks provide a more thermally efficient regime for heat-transport applications. These results highlight the potential of elastic fluids for future heat transfer enhancement technologies.
Recent LLM-based agents have closed substantial portions of the scientific discovery loop in software-only machine-learning research, in chemistry, and in biology. Extending the same loop to high-fidelity physical simulators is harder, because solver completion does not imply physical validity and many failure modes appear only in field-level imagery rather than in solver logs. We present AI CFD Scientist, an open-source AI scientist for computational fluid dynamics (CFD) that, to our knowledge, is the first to span literature-grounded ideation, validated execution, vision-based physics verification, source-code modification, and figure-grounded writing within a single inspectable workflow. Three coupled pathways cover parameter sweeps within a fixed solver, case-local C++ library compilation for new physical models, and open-ended hypothesis search against a reference comparator, all running on OpenFOAM through Foam-Agent. At the center of the framework is a vision-language physics-verification gate that inspects rendered flow fields before any result is accepted, rerun, or written into a manuscript. On five tasks under a shared GPT-5.5 backbone, AI CFD Scientist autonomously discovers a Spalart-Allmaras runtime correction that reduces lower-wall Cf RMSE against DNS by 7.89% on the periodic hill at Reh=5600; under matched LLM cost, two strong general AI-scientist baselines (ARIS, DeepScientist) execute partial CFD workflows but lack the domain-specific validity gates needed to convert runs into defensible scientific claims; and a controlled planted-failure ablation shows that the vision-language gate detects 14 of 16 silent failures missed by solver-level checks. Code, prompts, and run artifacts are released at https://github.com/csml-rpi/cfd-scientist.
We propose a non-intrusive reduced-order modeling framework for parametrized visco-plastic free-surface flows governed by a shallow-water formulation of Herschel--Bulkley fluids. These flows exhibit strong nonlinearities, non-smooth rheology, moving fronts, and yield surfaces, making efficient surrogate modeling particularly challenging. To address this challenge, we employ a tensor-based approach in which the solution manifold is approximated using a low-rank representation obtained via higher-order singular value decomposition of snapshot data over a structured parameter space. The resulting tensorial reduced-order model (TROM) enables rapid online evaluation by directly reconstructing solution trajectories from the compressed representation, thereby avoiding the need to perform time integration of a reduced dynamical system. The proposed non-intrusive framework can be interpreted as an encoder--decoder architecture with a compressed latent representation and efficient multilinear decoding. Numerical experiments demonstrate that the proposed approach accurately captures key flow features, including front propagation, plug and shear regions, and near-stopping dynamics, while achieving substantial computational speedups relative to full-order simulations.
Mixing of miscible liquids is an essential process in multiple industrial settings, usually with the intent to homogenize the product. This seemingly simple process is in fact a complex hydrodynamic problem that has a direct impact on the product quality. In this study, numerical simulations of a stirred tank were performed with a 50/50 ratio of liquids and systematically varied the Reynolds and Richardson numbers. A positive correlation between the mixing time and the Richardson number was observed, as reported in the literature. The influence of the Reynolds number was not as pronounced and clear. Based on the Power, Froude and Richardson numbers, we were able to derive an exponential scaling for the dimensionless mixing time that collapsed all our data onto one master curve.
This paper presents a topology optimization method for designing two-fluid heat exchangers under turbulent conditions using a Darcy flow-based low-fidelity (LF) model. The LF model is calibrated against a high-fidelity (HF) model based on the Reynolds-averaged Navier-Stokes (RANS) equations to increase the accuracy of predictions for fluid flow and heat transfer characteristics. Since the discrepancies between the LF and HF models can be significant, particularly for pressure drops, a multifidelity topology optimization framework is adopted to leverage the strengths of both models. Using the calibrated LF model, we perform topology optimization for various inlet velocities in the boundary conditions and trade-off parameters in the objective function to obtain diverse optimized designs. The optimized designs are then evaluated using the HF model to assess their performance with higher accuracy. The results demonstrate that the optimized designs significantly improve overall heat transfer coefficients while maintaining manageable pressure drops, achieving up to a 22% higher performance evaluation criterion (PEC) compared to a reference design enhanced by conventional twisted tape insertion. The improvements are attributed to the optimized configurations that promote enhanced fluid mixing and increased surface area for heat exchange, yet maintain streamlined flow paths to minimize pressure losses. Overall, the proposed topology optimization method using the Darcy flow-based LF model proves effective in designing high-performance double pipe heat exchangers, showcasing the potential of the multifidelity approach in overcoming the challenges of optimizing heat exchangers under turbulent flow conditions.
Thermal protection remains a critical challenge for oblique detonation engines (ODEs) operating under hypersonic conditions due to the extreme heat release and compact combustor geometry associated with oblique detonation waves (ODWs). In the present study, the effectiveness of film cooling for a kerosene-air ODE combustor is numerically investigated under a flight Mach number of 10 and an altitude of 15 km. Three active cooling strategies are considered, including air film cooling, gaseous-kerosene film cooling, and liquid-kerosene mist cooling. The results show that all cooling strategies preserve stable oblique-detonation propagation and maintain the canonical wave-system structure within the investigated operating range. Air cooling produces stronger disturbances near the initiation region and triple point, resulting in enhanced downstream wave interactions and larger propulsion penalties. In contrast, fuel-based cooling induces milder disturbances and better preserves the global detonation structure. All cooling methods substantially reduce the near-wall thermal load, although their cooling characteristics differ significantly. Gaseous-kerosene film cooling exhibits a spatially periodic near-wall thermal response associated with the discrete cooling hole arrangement, while liquid-kerosene mist cooling produces a smoother near-wall temperature distribution due to enhanced two-phase mixing and phase-change heat absorption. Among the investigated strategies, mist cooling provides the best overall balance between thermal protection and propulsion performance at coolant mass ratios of 1%-3%, whereas gaseous-kerosene film cooling becomes advantageous at higher injection levels due to improved wall coverage continuity. The present results demonstrate the feasibility and potential of fuel-based film cooling for thermal management in hypersonic ODE combustors.
The numerical modeling of hydraulic jumps remains challenging due to complex interactions among free-surface deformation, air entrainment and detrainment, and turbulent bubble transport. Whereas accurate prediction of these flows is essential for the design of hydraulic structures, existing high-fidelity tools require prohibitive computational resources for engineering applications. This study implements a three-phase mixture model based on an Unsteady Reynolds-Averaged Navier Stokes (URANS) framework, to numerically simulate flow and air entrainment across twelve hydraulic jumps with Froude numbers ranging from $1.98$ to $8.48$, representing the first systematic analysis for such a comprehensive range of Froude numbers. The model accurately represents time-averaged velocity fields and air concentration profiles, as well as dynamic features including jump toe oscillation and free-surface deformation, showing good agreement with experimental data from seven facilities. Compared to Improved Delayed Detached Eddy Simulations (IDDES), the proposed approach achieves similar accuracy with approximately 400-fold fewer cells and a 300-fold reduction in computational cost. The investigation shows that the selection of turbulence closure affects the accuracy of the prediction of air entrainment. These findings establish the three-phase mixture approach as a practical engineering tool for hydraulic jump simulation, offering an effective balance of accuracy and computational cost.
Rapidly rotating Rayleigh-Bénard convection admits a class of exact steady single-mode solutions describing high-amplitude convection cells. Using a matched asymptotic analysis in the high-Rayleigh-number limit, we obtain a rigorous characterization of their bulk and boundary-layer structure, yielding explicit scaling laws for the Nusselt and Reynolds numbers, including their dependence on the horizontal wavenumber. We show that, for suitable wavenumbers, these solutions attain the diffusivity-free ultimate scalings frequently assumed for geophysical and astrophysical convection, with additional enhancing logarithmic corrections. This reveals a specific mechanism through which rapidly rotating convection can approach ultimate heat transport via coherent columnar structures with well-defined horizontal scales.
The prediction of aircraft icing is conventionally performed using multishot simulation frameworks that fail to predict the progressive roughening of the ice surface. To understand roughness formation, we investigate droplet impingement on clean and laser-scanned rough ice shapes using a high-fidelity computational framework based on wall-modeled large-eddy simulations and Lagrangian particle tracking. This methodology is validated against experimental data for a NACA 23012 airfoil and a NACA 64A008 swept tail, accurately predicting collection efficiency and supercooled large droplet splashing. The framework is subsequently applied to laser-scanned rime ice geometries to quantify the impact of surface roughness on local impingement distributions. The results reveal that physical roughness induces a highly nonuniform collection efficiency, with droplet impingement intensely concentrated on upstream-faces of roughness elements, creating sheltered shadow zones immediately downstream. While the spanwise-averaged collection efficiency remains remarkably similar to that of an equivalent smooth body, idealized smooth surfaces completely suppress these localized impingement peaks. Ice accretion simulations demonstrate that this localized impingement creates a self-reinforcing feedback loop, actively amplifying existing roughness features over time. These findings provide a direct physical explanation for the formation of characteristic rime ice structures and highlight the critical role of local surface topology in the accretion process.
We analyze the dynamics towards partial thermalization and subsequent cooling in the defocusing two-dimensional nonlinear Schrödinger model, using direct simulations and insights from the wave-kinetic equations (WKE) and a fourth-order differential approximation model (DAM). We show that the evolving WKE spectrum exhibits two distinct similarity ranges--the quasi-thermal core and the ultraviolet tail--whereas in the DAM, an additional range of infrared self-similarity appears. By stretching the quasi-thermal region, the self-similar fronts drive an effective dynamical cooling process towards the formal but ill-defined equilibrium state at vanishing temperature--analogous to an ultraviolet catastrophe in a system of classical waves.
Reconstructing high-fidelity flow fields from low-fidelity observations is a central problem in scientific machine learning, yet recent diffusion and flow-matching models typically rely on iterative sampling, making them costly for latency-sensitive workflows such as ensemble forecasting, real-time visualization, and simulation-in-the-loop inference. We study whether a high-fidelity flow-matching generative model can be compressed into a compact one-step model for fast scientific flow reconstruction. Our approach distills an optimal-transport flow-matching teacher into a one-step consistency model. Low-fidelity observations are incorporated at inference by initializing the generative trajectory from a noised observation along the transport path, allowing an unconditional high-fidelity flow model to perform conditional reconstruction without retraining the teacher. We evaluate this distillation strategy on three fluid benchmarks, Smoke Buoyancy, Turbulent Channel Flow, and Kolmogorov Flow, using coarse-to-fine reconstruction as a controlled testbed at field sizes up to $256 \times 256$. Across these settings, the distilled student retains similar performance of the teacher's model on spectrum metrics, while using roughly half as many parameters and achieving a $12\times$ inference speedup over the flow-matching teacher. Under the same training budget, the distilled student also outperforms a one-step consistency model trained directly from scratch by $23.1\%$ in SSIM, showing that teacher distillation improves training efficiency rather than merely accelerating sampling. These results suggest a promising route for turning future high-capacity scientific generative models into compact reconstruction models that are faster to train, cheaper to run, and easier to deploy.
Cephalopod pulsed-jet locomotion is not a single isolated expulsion event, but a coordinated cycle involving jet expulsion, passive gliding, and mantle refilling. Inspired by this cycle-resolved biological strategy, this paper presents a cephalopod-inspired pulsed-jet robot with a rigid-soft hybrid origami mantle that enables large, actively driven, and geometry-guided body deformation. The proposed mantle integrates rigid folding panels with a compliant silicone framework, allowing a 75% effective cavity-volume reduction during expulsion and reducing the projected cross-sectional drag area by approximately 75.7% in the contracted gliding configuration. Using this platform, we formulate a cycle-resolved framework to separately investigate how expelled volume, glide duration, and refill pathway influence whole-cycle locomotion performance. Experiments show that the robot reaches a peak speed of approximately 0.5 m/s (3.8 BL/s) and an average speed exceeding 0.2 m/s (1.5 BL/s) within the first jetting cycle. The results further demonstrate the roles of high expelled-volume-ratio contraction in speed generation, reduced-drag-area gliding under different glide durations, and mantle-aperture-inspired passive inlet valves in assisting refill. This work provides both a robotic implementation of actively deformable cephalopod-like jet propulsion and a unified experimental platform for studying expulsion-gliding-refilling dynamics in pulsed-jet locomotion.
We propose a mixed quantum-classical hydrodynamic framework to model short-time inertial effects in the non-adiabatic evolution of a quantum solute coupled to a classical polar solvent. Drawing upon the work of Burghardt and Bagchi [Chem. Phys. 329 (2006), 343], we employ the Hamiltonian approach to incorporate consistent backreaction and preserve quantum decoherence beyond standard Ehrenfest dynamics. The solvent is treated as an ideal polar fluid and the quantum solute state is correlated to both the position and molecular orientation coordinates of the liquid. This approach retains essential solute-solvent correlations while significantly reducing the computational complexity of previous approaches. We further incorporate dissipative terms to capture both inertial effects and polarization relaxation. After establishing the general setting for non-local dielectric continua, the Marcus local approximation is integrated into the model thereby extending traditional solvation theory to account for collective fluid sloshing on fast timescales.
Fast surrogate modeling for high-dimensional physical dynamics requires more than low short-term error: useful models must roll out efficiently while preserving the statistical structure of long trajectories. Neural operators provide inexpensive autoregressive forecasts but can drift in turbulent regimes, whereas rolling diffusion and latent generative surrogates can represent stochastic transitions at the cost of multi-step denoising, noise-schedule design, or auxiliary compression models. We propose MeanFlow Long-term Invariant Spatiotemporal Consistency Autoregressive Models (MeLISA), a latent-free autoregressive generative surrogate built on pixel-space MeanFlow. MeLISA defines a blockwise stochastic transition kernel that generates each forecast block with a single model evaluation, avoiding latent encoders and iterative diffusion solvers at inference time. To stabilize long-horizon rollouts, MeLISA combines a Window-Consistency MeanFlow objective that learns conditional spatiotemporal generation from partially observed temporal windows with a Time Increment Consistency loss that constrains multi-lag finite increments and targets temporal-correlation structure. We evaluate MeLISA with compact UNet and scalable DiT backbones on two high-resolution benchmarks, extended 2D Kolmogorov flow at $256 \times 256$ and turbulent channel-flow slice at $192 \times 192$. MeLISA outperforms neural-operator baselines on short-term forecasting accuracy and long-horizon statistical metrics, including energy spectra, turbulent kinetic energy, and mixing-rate-related dynamics, while achieving inference speeds comparable to, and in some cases faster than, neural operators. Compact 3.7-5.7M-parameter variants already deliver strong parameter efficiency, and DiT variants provide a scalable path up to 150M parameters. Overall, MeLISA benefits both rollout efficiency and long-horizon statistical accuracy.