Fluid Dynamics
See recent articles
Showing new listings for Wednesday, 25 March 2026
- [1] arXiv:2603.22556 [pdf, html, other]
-
Title: Vibrissa inspired geometries enhance sensitivity of wake-induced vibrationsSubjects: Fluid Dynamics (physics.flu-dyn)
We report on experiments designed to characterize the vortex-induced vibration (VIV) and wake-induced vibration (WIV) experienced by bluff bodies immersed in both steady and unsteady flows. Using a real-time Cyber-Physical System (CPS) we systematically prescribe the virtual mass, spring constant, and damping of elastically mounted models. This allows us to characterize the forces and displacements of the free vibration of a circular cylinder, elliptical cylinder, and a seal whisker inspired vibrissa model with undulating elliptical geometry. In a free flow, the circular cylinder exhibits high VIV, while the reduced aspect ratio objects have minimal vibration across all structural frequencies. When a flow disturbance of a pitching and heaving hydrofoil is introduced, the reduced aspect ratio objects are excited by WIV with highest amplitude oscillations occurring when structural frequency of the test object matches wake frequency of the upstream foil. To further understand the benefits of an undulated geometry over a classic elliptical cylinder, we assess the nonlinear fluid damping experienced by each test object by comparing experimental data to quadratic drag and Van der Pol damping models. Our results show that the amplitude dependent Van der Pol damping model better describes the physical system for both test objects by capturing the suppression of large amplitude WIV, but recovering small amplitude VIV. However, the strength of the fitted Van der Pol damping coefficient is greater for the elliptical cylinder than the vibrissa. We find the vibrissa experiences lower damping than the elliptical cylinder across all tested structural frequencies, indicating how the vibrissa geometry may serve as a higher sensitivity sensor.
- [2] arXiv:2603.22737 [pdf, html, other]
-
Title: Modification of the k-omega0 model for roughnessSubjects: Fluid Dynamics (physics.flu-dyn)
Surface roughness plays a substantial role in many flows for which Reynolds averaged prediction is needed. The transformation used in the k-omega0 model is extended to rough surfaces by adding an effective origin. The log-layer offset is computed as a function of this effective origin, thereby creating a correspondence between effective origin and equivalent sandgrain roughness. A formula is derived for the virtual origin of the fully rough log law. It is shown how the present model is consistent with the fully rough limit.
- [3] arXiv:2603.22902 [pdf, html, other]
-
Title: A diffuse-interface model for N-phase flows with liquid-solid phase changeComments: 19 pages, 13 figuresSubjects: Fluid Dynamics (physics.flu-dyn)
In this work, we first propose a diffuse interface model for simulating N phase flows with solid liquid phase change. In this model, a phase field approach is adopted to capture multiphase fluid interfaces, and an enthalpy based formulation is used to describe the phase change. The volume changes resulting from density differences during phase change are incorporated by introducing a source term into the continuity equation. The method also satisfies the reduction consistent property, allowing it to rigorously degenerate to both the conservative phase field method for N phase flows and the classical enthalpy method for solid liquid phase change. Then a coupled lattice Boltzmann (LB) method is developed to solve this diffuse interface model. Some numerical tests, including film freezing, single droplet freezing, and compound droplet freezing are performed, and the results are in good agreement with the analytical solutions and data reported in the previous works. Furthermore, the proposed method is applied to study freezing dynamics of complex systems with insoluble impurities, capturing the interaction between the advancing freezing front and embedded impurities. It is found that the proposed diffuse interface method is accurate and efficient for studying N phase systems with phase change.
- [4] arXiv:2603.22927 [pdf, html, other]
-
Title: Direct Numerical Simulation of MILD Combustion: Mixing and Autoignition from Non-Premixed StreamsComments: 8 pages for the main paper and 6 pages for the supplementary material, 7 figures for the main paper and 6 figures for the supplementary materialSubjects: Fluid Dynamics (physics.flu-dyn)
Moderate or intense low-oxygen dilution (MILD) combustion is achieved by strongly diluting and preheating the reactants through mixing with hot combustion products before ignition. To better understand how fuel/air/product mixing and interaction govern MILD combustion dynamics, a novel direct numerical simulation (DNS) dataset of a temporally evolving three-stream mixing layer consisting of fuel, air, and hot combustion products has been performed. In this configuration, both fuel-air and air-hot products mixing processes are considered with varying time scales, through four carefully designed DNS cases, to assess how their combined interaction controls ignition under MILD conditions. It is observed that the cases with higher dilution levels fall within the MILD combustion regime, whereas those with lower dilution correspond to non-MILD conditions. The results show that, as long as MILD conditions are observed, ignition is mainly driven by mixing with hot products. Flame index (FI) combined with chemical explosive mode analysis (CEMA) further identifies the local combustion mode: in MILD cases, ignition occurs predominantly through a premixed-autoignition mode, while in non-MILD scenarios, the premixed-deflagrative contribution to the heat release rate is more substantial. Conditional analysis of scalar dissipation rates shows that the combustion modes in MILD conditions are sensitive to mixing by both the fuel and hot products, whereas the combustion modes in non-MILD conditions are mainly influenced by the mixing of the fuel with the surrounding gases.
- [5] arXiv:2603.23100 [pdf, html, other]
-
Title: A numerical study on the coefficient of restitution of wet collisionsComments: Submitted to Powder TechnologySubjects: Fluid Dynamics (physics.flu-dyn); Computational Physics (physics.comp-ph)
Using smoothed particle hydrodynamics (SPH) simulations, we investigate the coefficient of restitution (COR) in wet collisions and identify a scaling law governing its behavior. The simulations employ an updated-Lagrangian, mesh-free framework that is validated against experimental measurements. We neglect surface tension effects since the impact conditions correspond to a moderate-to-high Weber number regime. The COR is found to depend on the Stokes number and a dimensionless film thickness defined as the ratio of the liquid film thickness to the diameter of the impacting solid bead. Two distinct regimes are observed, each characterized by different power-law exponents.
- [6] arXiv:2603.23181 [pdf, html, other]
-
Title: Steady cone-jet mode of electrospray for single-cell depositionSubjects: Fluid Dynamics (physics.flu-dyn)
We propose using the electrospray cone-jet mode operated near its minimum-flow-rate stability limit for single-cell deposition. Because the jet is much thinner than the cells themselves, individual cells can be clearly visualized or detected during deposition. At such low flow rates, individual cells can be placed at distinct, user-defined locations, even at relatively high cell concentrations. In this sense, our approach provides a spatial resolution at the scale of a single cell. We demonstrate the method's capabilities by depositing cells onto a millimeter-scale droplet of a standard cell-culture medium. Cell viability assays indicate that many cells maintain membrane integrity after exposure to the electrosprayed liquid, suggesting that most damage is reversible.
New submissions (showing 6 of 6 entries)
- [7] arXiv:2603.22827 (cross-list from physics.comp-ph) [pdf, html, other]
-
Title: Wafer-to-Wafer Bonding: Part: I -- The Coupled Physics Problem and the 2D Finite Element ImplementationSubjects: Computational Physics (physics.comp-ph); Applied Physics (physics.app-ph); Fluid Dynamics (physics.flu-dyn)
Wafer-to-wafer (WxW) bonding is a key enabler for three-dimensional integration, including hybrid bonding for fine-pitch Cu-Cu interconnects. During bonding, wafer deformation and the air entrapped between the wafers interact through a strongly coupled, time-dependent fluid-structure interaction (FSI) that can produce non-intuitive bonding dynamics and process sensitivities. This paper develops a mathematically consistent reduced-order model for WxW bonding by deriving a Kirchhoff-Love plate equation for wafer bending from three-dimensional linear elasticity and coupling it to a Reynolds lubrication equation for the inter-wafer air film. The resulting nonlinear plate-Reynolds system is discretized and solved monolithically in the high-performance FEniCSx framework using a $C^0$ interior-penalty formulation for the fourth-order plate operator, standard continuous Galerkin discretization for the pressure field, implicit time integration, and a Newton solver with automatic differentiation. Simulations reproduce experimentally reported probe-displacement histories for multiple initial gaps and verify force equilibrium at the bond front, where the Reynolds pressure acts as an effective contact reaction. Parametric studies reveal nonlinear, and in some cases non-monotonic, sensitivities of bonding-front kinetics to the initial gap, air viscosity, and interfacial energy, providing actionable trends for process optimization.
- [8] arXiv:2603.23371 (cross-list from math-ph) [pdf, other]
-
Title: Interfacial instability as a trigger for dryout inception in two-phase CO2 flowSubjects: Mathematical Physics (math-ph); Fluid Dynamics (physics.flu-dyn)
Progress in particle physic leads to increasing in detector luminosity and a consequent increasing overheating induced by Joule effect. An effective cooling strategy is the exploitation of CO\textsubscript{2} heat latency in phase-change. An additional challenge, relevant to detectors for High Energy Particles, is the consequent geometrical constrain due to the limited space avialable for the cooling system within the detector arrangement, leading to the implementation of cooling system by means of millichannels. In this context, at relative high vapour quality the liquid phase exhibits annular flow, anticipating the dryout. Dryout is a critical condition where the heat transfer coefficient dramatically drops and dangerous temperature levels can be reached, potentially leading to catastrophic consequences. Experimental evidences reveal that its behavior in two-phase annular flows differs from conventional refrigerants and the fundamental inception-mechanism is not yet understood. This study aims at investigating the key new idea whereby dryout inception is triggered by instability of the liquid-vapour interface. A mathematical model for two-phase annular flow is presented and the stability of the interface between the two fluids is studied through the linear theory. The stability analysis reduces to solving a coupled forth-order differential eigenvalue problem that is treated numerically with an in-house code based on the Chebyshev-$\tau$ method. Numerical investigations identify a critical value for the vapour quality, named $x_{dry}$, that leads to interface instability. The resulting predictions on $x_{dry}$ are confirmed by experimental data collected from two independent experimental campaigns, validating the hypothesis that dryout inception is governed by interfacial instabilities.
Cross submissions (showing 2 of 2 entries)
- [9] arXiv:2503.13659 (replaced) [pdf, other]
-
Title: Scaling Laws Governing Droplet Spreading and Merging Dynamics on Solid Surfaces: A Molecular Simulation StudyComments: 15 FiguresSubjects: Fluid Dynamics (physics.flu-dyn)
This study employs molecular dynamics simulations to investigate droplet dynamics when a stationary droplet on a solid surface is struck by another droplet of similar size from above. The focus is on the jumping behavior of the merged droplet and the associated energy conversion. The process is primarily governed by the amount of energy converted into kinetic energy after dissipation. At high impact velocities, the energy conversion efficiency becomes constant, with only about 1% lost due to surface adhesion, an effect that diminishes with increasing velocity. Factors such as impact velocity, droplet size, surface texture, and wettability significantly influence the jumping velocity. Scaling laws are developed for the maximum spreading time, spreading factor, and restitution coefficient based on the Weber (We) and Reynolds (Re) numbers, which differ from those for single droplet impacts. On superhydrophobic surfaces, the spreading time is approximated as three times droplet radius to impact velocity, and its dimensionless form scales linearly with We 0.31. The general scaling laws for the spreading factor for velocity dependent and velocity independent spreading regimes are developed.
- [10] arXiv:2509.09315 (replaced) [pdf, html, other]
-
Title: On the role of water activity on the formation of a protein-rich coffee ring in an evaporating multicomponent dropSubjects: Fluid Dynamics (physics.flu-dyn); Soft Condensed Matter (cond-mat.soft)
The coffee-ring effect is a universal feature of evaporating sessile droplets with pinned contact line, wherein solutes or particles are advected to the droplet's edge due to evaporation-driven flows. While existing models have successfully described this phenomenon in particle-laden droplets, they often assume that the evaporative flux, and thus hydrodynamics, are decoupled from solute transport. This assumption breaks down in complex fluids, such as protein or polymeric solutions, where the solute can influence evaporation through changes in water activity. Here, we investigate model respiratory droplets primarily composed of water, salt, and a type of the glycoprotein mucin. Using fluorescence microscopy, we observe the formation of a well-defined protein ring at the droplet edge as water evaporates. The growth and morphology of this ring exhibit a strong dependence on ambient relative humidity ($H_r$), revealing dynamics that existing models cannot capture. Specifically, we find that protein accumulation at the edge is governed by the feedback between local solute concentration and evaporation rate. To account for this, we develop a minimal theoretical model based on the lubrication approximation, incorporating the coupling between hydrodynamics and solute transport through the evaporation rate. Our framework reproduces key features of the experimental observations and suggests a physical basis for the $H_r$-dependent stability and infectivity of respiratory droplets containing viruses.
- [11] arXiv:2509.25983 (replaced) [pdf, html, other]
-
Title: Rare-event detection in a backward-facing-step flow using live optical-flow velocimetry: observation of an upstream jet burstSubjects: Fluid Dynamics (physics.flu-dyn)
Rare and extreme events in turbulent flows play a critical role in transport, mixing and transition, yet are notoriously difficult to capture experimentally. Here we report, to our knowledge, the first direct experimental detection of an upstream-directed jet burst in a backward-facing step (BFS) flow at $Re_h=2100$, using long-duration Live Optical Flow Velocimetry (L-OFV). Continuous monitoring over 1.5 h enabled a data-driven definition of extremes as rare velocity probes excursions deep into the observed distribution's tails; in practice, large negative events ($u: Z < -6$, $v: Z < -5$ at $(x,y) = (2h,h / 2)$, where $|Z| > > 0$ stands for large deviations from the mean value) triggered the live capture of surrounding velocity fields. The recording is triggered when the probes surpass the defined threshold, using live analysis of the velocity fields. The detected event features a jet-like intrusion into the recirculation region initiated by the collapse of a merged Kelvin-Helmholtz vortex and sustained by counter-rotating vortices, and is accompanied with heavy-tailed probe statistics and simultaneous amplification of fluctuating kinetic energy and enstrophy. While a single event was recorded, underscoring its rarity, the results establish L-OFV as a viable platform for rare-event detection in separated shear layers and document a previously unreported mechanism of upstream jet bursting in BFS flow.
- [12] arXiv:2512.09140 (replaced) [pdf, html, other]
-
Title: Simulating surfactant effects in phase-transforming fluidsSubjects: Fluid Dynamics (physics.flu-dyn); Computational Physics (physics.comp-ph)
Surfactants are critical in natural processes and engineering, but measuring their concentrations in non-equilibrium conditions and in the presence of flow is difficult. Therefore, computational methods are a key tool for improving our understanding. Predicting the effect of surfactants on liquid-vapor transformations is particularly challenging due to (1) simultaneous mass transfer, non-equilibrium thermodynamics and Marangoni stresses, and (2) the phenomenological assumptions underlying many liquid-vapor phase-change models. Starting from the Navier-Stokes-Korteweg equations, a first-principles approach to liquid-vapor phase transformations, we developed a model of liquid-vapor flows with surfactants. We performed simulations of bubbles under equilibrium and liquid-vapor interface oscillations to demonstrate that the model successfully reproduces surfactant-mediated reductions in surface tension. We also investigated the mechanisms whereby surfactant affects bubble coalescence and condensation. Overall, this work provides a new framework for studying the effect of surfactants on liquid-vapor transformations and suggests multiple areas for future research, including the impact of complex surface chemistries on flow around bubbles and the acoustic response of bubbles with surfactants.
- [13] arXiv:2512.12628 (replaced) [pdf, html, other]
-
Title: Integrating Fourier Neural Operator with Diffusion Model for Autoregressive Predictions of Three-dimensional TurbulenceSubjects: Fluid Dynamics (physics.flu-dyn)
Accurately autoregressive prediction of three-dimensional (3D) turbulence has been one of the most challenging problems for machine learning approaches. Diffusion models have demonstrated high accuracy in predicting two-dimensional (2D) turbulence, but their applications in 3D turbulence are relatively limited. To achieve reliable autoregressive predictions of 3D turbulence, we propose the DiAFNO model which integrates the implicit adaptive Fourier neural operator (IAFNO) with diffusion model. IAFNO can effectively capture the global frequency and structural features, which is crucial for global consistent reconstructions of the denoising process in diffusion models. Furthermore, based on conditional generation from diffusion models, we design an autoregressive framework in DiAFNO to achieve long-term stable predictions of 3D turbulence. The proposed DiAFNO model is systematically tested with fixed hyperparameters in several types of 3D turbulence, including forced homogeneous isotropic turbulence (HIT) at Taylor Reynolds number $Re_{\lambda}\approx100$, decaying HIT at initial Taylor Reynolds number at $Re_{\lambda}\approx100$ and turbulent channel flow at friction Reynolds numbers $Re_{\tau}\approx395$ and $Re_{\tau}\approx590$. The results in the a posteriori tests demonstrate that DiAFNO exhibits a significantly higher accuracy in terms of the velocity spectra, the root-mean-square (RMS) values of both velocity and vorticity, and Reynolds stresses, as compared to the elucidated diffusion model (EDM) and the traditional large-eddy simulation (LES) using dynamic Smagorinsky model (DSM). Meanwhile, the well-trained DiAFNO is faster than LES with the DSM.
- [14] arXiv:2601.19351 (replaced) [pdf, html, other]
-
Title: Physics-Informed Transformer operator for the prediction of three-dimensional turbulenceSubjects: Fluid Dynamics (physics.flu-dyn)
Data-driven turbulence prediction methods often face challenges related to data dependency and lack of physical interpretability. In this paper, we propose a physics-informed Transformer operator (PITO) and its implicit variant (PIITO) for predicting three-dimensional (3D) turbulence, which are developed based on the vision Transformer (ViT) architecture with an appropriate patch size. Given the current flow field, the Transformer operator computes its prediction for the next time step. By embedding the large-eddy simulation (LES) equations into the loss function, PITO and PIITO can learn solution operators without using labeled data. Furthermore, PITO can automatically learn the subgrid scale (SGS) coefficient using a single set of flow data during training. Both PITO and PIITO exhibit excellent stability and accuracy on the predictions of various statistical properties and flow structures for the situation of long-term extrapolation exceeding 25 times the training horizon in decaying homogeneous isotropic turbulence (HIT), and outperform the physics-informed Fourier neural operator (PIFNO). Furthermore, PITO exhibits a remarkable accuracy on the predictions of forced HIT where PIFNO fails. Notably, PITO and PIITO reduce GPU memory consumption by 79.5\% and 91.3\% while requiring only 31.5\% and 3.1\% of the parameters, respectively, compared to PIFNO. Moreover, both PITO and PIITO models are much faster compared to traditional LES method.
- [15] arXiv:2603.07151 (replaced) [pdf, html, other]
-
Title: Optimize discrete loss with finite-difference physics constraint and time-stepping for PDE solvingSubjects: Fluid Dynamics (physics.flu-dyn)
Computational Fluid Dynamics (CFD) is an important approach for analyzing flow phenomena and predicting engineering-relevant quantities. The governing physics is formulated as partial differential equations(PDEs) and solved numerically on computational grids. Physics-informed neural networks(PINNs) have emerged as a popular optimization-based approach for solving PDEs, but they often suffer from ill-conditioned objectives and the high cost of automatic differentiation. Optimization-based discretizations such as ODIL mitigate several PINN drawbacks by optimizing discrete variables directly, yet accuracy and efficiency remain limited on body-fitted geometries and for time-dependent problems. This paper proposes FDTO, a finite-difference time-stepping loss-optimization solver that defines physics losses from discrete residuals. FDTO couples curvilinear coordinate transforms with body-fitted structured grids and decomposes long-horizon evolution into sequential, well-conditioned subproblems consistent with time marching. The method is primarily evaluated on incompressible Navier-Stokes flows, including lid-driven cavity benchmarks, external airfoil aerodynamics (lift/drag consistency), and a cylinder case on a multi-block structured mesh with cross-block coherent solutions. Additional validations on diffusion and flow-mixing problems further demonstrate generality. Compared with representative PINN-based solvers, FDTO reduces GPU memory by about 82.6% on the lid-driven cavity case and achieves 3-5 times lower relative error on the flow-mixing problem. These results indicate that FDTO enables accurate, stable, and memory-efficient discrete-loss optimization for incompressible-flow solutions, while remaining applicable to other PDE models.
- [16] arXiv:2603.18913 (replaced) [pdf, html, other]
-
Title: Geometric Dynamics of Turbulence: A Minimal Oscillator Structure from Non-local ClosureComments: Second revised version. Includes Supplementary MaterialSubjects: Fluid Dynamics (physics.flu-dyn); Chaotic Dynamics (nlin.CD); Classical Physics (physics.class-ph)
Turbulence remains one of the central open problems in classical physics, largely due to the absence of a closed dynamical description of the Reynolds stress. Existing approaches typically rely either on local constitutive assumptions or on high-dimensional statistical representations, without identifying a minimal set of dynamical variables governing the cascade response. Here we show that the non-local stress response implied by the Navier-Stokes equations admits a systematic reduction onto a low-dimensional anisotropic sector of the turbulent cascade. This reduction leads to a minimal dynamical system with the structure of a damped oscillator, arising from the coupling between the leading angular mode and its nonlinear transfer to higher-order sectors. Within this framework, classical turbulent behaviors -- including inertial-range scaling, shear-driven transport, and wall-bounded logarithmic profiles -- emerge as different realizations of the same underlying dynamical structure. Universal quantities such as the Kolmogorov constant and the von Kármán constant appear as leading-order consequences of internal consistency conditions applied across homogeneous and shear-driven regimes. These results suggest that turbulence admits a minimal dynamical backbone governed by non-local cascade response, providing a unified perspective that connects spectral transfer, anisotropy, and mean-flow interaction within a single reduced framework.
- [17] arXiv:2507.08338 (replaced) [pdf, html, other]
-
Title: Discontinuity-aware KAN-based physics-informed neural networksSubjects: Computational Physics (physics.comp-ph); Fluid Dynamics (physics.flu-dyn)
Physics-informed neural networks (PINNs) have proven to be a promising method for the rapid solving of partial differential equations (PDEs) in both forward and inverse problems. However, due to the smoothness assumption of functions approximated by general neural networks, PINNs are prone to spectral bias and numerical instability and suffer from reduced accuracy when solving PDEs with sharp spatial transitions or fast temporal evolution. To address this limitation, a discontinuity-aware physics-informed neural network (DPINN) method is proposed. It incorporates an adaptive Fourier-feature embedding layer to mitigate spectral bias and capture steep gradients, a discontinuity-aware network that generalizes the Kolmogorov representation theorem to the discontinuous regime for the modeling of shock-wave properties, mesh transformation to accelerate convergence across complex geometries, and learnable local artificial viscosity to stabilize the algorithm near discontinuities. In numerical experiments regarding the inviscid Burgers' equation, Riemann problems, and transonic and supersonic airfoil flows, DPINN demonstrated superior accuracy in capturing discontinuities compared to existing methods.
- [18] arXiv:2509.24180 (replaced) [pdf, html, other]
-
Title: Design of model Boger fluids with systematically controlled viscoelastic propertiesComments: 16 pages and 12 figures (including the supplementary material), in press at Journal of RheologySubjects: Soft Condensed Matter (cond-mat.soft); Fluid Dynamics (physics.flu-dyn)
The subject of viscoelastic flow phenomena is crucial to many areas of engineering and the physical sciences. Although much of our understanding of viscoelastic flow features stems from carefully designed experiments, preparation of model viscoelastic fluids remains a challenge; for example, fabricating a series of fluids with different fluid shear moduli $G_0$, but with an identical relaxation time $\tau$, is nontrivial. In this work, we harness the non-ideality of nearly constant-viscosity elastic fluids, commonly known as `Boger fluids', made with polyisobutylene, to develop an experimental methodology that produces a set of fluids with desired viscoelastic properties, specifically, $G_0$, $\tau$, and the first normal stress difference coefficient $\psi_1$. Through a linear algebraic relation between the rheological properties of interest ($G_0$, $\tau$, $\psi_1$) and the fluid compositions in terms of polymer concentration $c$, molecular weight $M_w$, and solvent viscosity $\eta_s$, we developed a `design equation' that takes $G_0$, $\tau$, $\psi_1$ as inputs and calculates values for $c$, $M_w$, $\eta_s$ as outputs. Using this method, fabrication of dilute viscoelastic fluids whose rheological properties are \textit{a priori} known can be achieved.
- [19] arXiv:2603.02306 (replaced) [pdf, html, other]
-
Title: Bridging the Prandtl number gap: 3D simulations of thermohaline convection in astrophysical regimesComments: 11 pages (7 ignoring appendices/references), 3 figures; accepted for publication in ApJL; minor text revisions compared to v1; code and data at DOI: https://doi.org/10.5281/zenodo.19187749Subjects: Solar and Stellar Astrophysics (astro-ph.SR); Fluid Dynamics (physics.flu-dyn)
Thermohaline convection (also known as fingering convection or thermohaline mixing) occurs in stellar radiation zones where a sufficient inversion of the mean molecular weight is present. This process mixes chemicals radially and occurs in a variety of stars, including near the luminosity bump on the red giant branch and potentially in polluted white dwarfs. Previous efforts to characterize this process using 3D simulations have been restricted to regimes far from actual stars: The Prandtl number $\Pr$--the ratio of the kinematic viscosity to thermal diffusivity--assumes values as low as $10^{-6}$ in stars, but 3D simulations have been restricted to $\Pr \gtrsim 10^{-2}$. For this reason, disagreements between observations and simulations are routinely dismissed as stemming from this $\Pr$ gap. This letter bridges this gap and demonstrates that 3D simulations of thermohaline convection can be performed in stellar parameter regimes. Using a suite of simulations spanning previously studied regimes with $\Pr \gtrsim 10^{-2}$ down to $\Pr = 10^{-6}$, we demonstrate that the chemical mixing model of Brown, Garaud, & Stellmach (2013) remains consistent with 3D simulations across both regimes. Therefore, tensions between this model and observations cannot be dismissed as resulting from a $\Pr$ gap, and must be resolved by considering additional physics.