2026-11-10 | | Total: 24
The collapse of a drop-impact crater that bypasses the pinch-off singularity produces exceptionally thin, high-velocity jets, often referred to as singular-jets. We perform 2D axisymmetric DNS using Basilisk, employing adaptive mesh refinement to capture the intricate flow features. The numerical results for a drop diameter of $3.64$ mm at ambient pressure match the experiments well under different impact conditions. During the crater collapse, qualitatively different dimple shapes form, including an entrapped bubble beneath the dimple. Quantitatively, an entrapped bubble with a diameter of ~210 μm and a jet velocity of 61.7 m/s for the singular-jet case during cylindrical dimple collapse is obtained from simulations. Subtle morphological changes in crater evolution distinguish the jetting regime, with a critical angle of 109.2 between the crater side wall and its base in the singular-jet case. The simulation provides insight into the internal flow fields, enabling us to delineate the observed phenomenon. The pressure and velocity fields reveal that the high jetting velocity is associated with the formation and rapid collapse of a nearly cylindrical dimple. The resulting pressure distribution drives a strongly upward flow near the dimple base, while the nearly parallel dimple side walls redirect radial liquid motion towards the axis. We extend the simulations into the post-entrainment region and identify another singular-jet near the upper boundary of the entrainment region. The jet reaches a velocity of 63.93 m/s and exhibits morphological characteristics similar to the first singular-jet case. A regime map is constructed from simulations spanning Weber numbers ranging from 100 to 200 and Froude numbers from 45 to 92; the updated lower and upper boundaries for the entrainment region, based on singular-jet formation, are We = 45.5\, Fr^{0.237} and We = 53.5\, Fr^{0.251}, respectively.
Numerical simulation of hypersonic vehicles is an attractive approach to investigate the complex phenomena involved in such technologies. For this, second-order finite volume methods (FVM) are the established method in computational fluid dynamics, with the cost of high dissipation in affordable meshes, thus needing refined ones or adaptive mesh refinements techniques to obtain accurate results. An alternative is to make use of high-order discontinuous Galerkin methods, in particular discontinuous Galerkin spectral element method (DGSEM), once it achieves higher formal order of accuracy on element-local polynomial bases. Therefore the goal of this work is to present the numerical simulations of a two-dimensional scramjet configuration to test and verify the applicability of the DGSEM approach for a complex hypersonic flow where shock-shock interactions occur. Results revealed to be in agreement with reference FVM solvers, thus verifying the code implementation.
Modons, or dipolar vortices, are a form of propagating vortex solution consisting of two counterrotating regions of fluid and moving together under self-advection. They are common features of the upper ocean that can remain stable over long times, thereby transporting fluid over large distances. Here, we present a semi-analytical method for finding modon solutions in a vertically-sheared background flow using the multi-layer quasi-geostrophic framework. This method works for arbitrarily many fluid layers and builds on our previous studies of multi-layer modon solutions. We present three example of modons in a background flow, highlighting the effect of varying relative potential velocity between layers and the impact of depth-dependent potential vorticity gradients arising from surface deformation and velocity shear.
The squeeze film between an electrostatically actuated microbeam can both induce damping and add stiffness to the structure. In addition, the dynamics of the film can interact with the nonlinear response of the microbeam. This paper extends a reduced-order model built on the complex modes of the coupled fluid-structure problem to the computation of periodic responses and their stability for clamped-clamped microbeams. Periodic orbits are continued in the frequency domain, and their stability is determined from their Floquet multipliers. The robustness of the response is assessed through basins of attraction, integrity measures, and dynamic pull-in maps in the plane of the excitation parameters. The model is validated against the measured and computed frequency response of a microswitch. Closed-form expressions are derived for the fluid relaxation modes, together with a correction for strongly coupled films, and define a parameter that characterizes the regime of the film. The approach is applied to different microbeams including cases where fluid modes relax on the timescale of the oscillation. For all simulated devices, a single fluid mode is necessary and sufficient, whereas adding coupled modes cannot replace it, whether the fluid modes respond quasi-statically or carry their own dynamics. Removing the fluid mode displaces both bifurcations that bound the pull-in band, changes the size of the basins of attraction and the position of the attractor within them, misses an abrupt erosion of the safe basin under increasing excitation, and shifts the whole region of dynamic pull-in in the excitation plane. The reduced model, with three real unknowns, makes these global analyses affordable for a fully coupled fluid-structure system.
Variational principles for dissipative processes come in two families: minimum principles at nonequilibrium stationary states (Helmholtz, Kirchhoff, Prigogine) and maximum principles for relaxation in continuum mechanics (Ziegler, Hackl, Fischer). We contribute to the understanding of the thermodynamic extremal principle (TEP) by investigating its extremal properties, and uncover the connection between the two families, whose distinction was recently made precise for diffusion: the extremal character is decided by the regime, which is either relaxation at frozen state, as in a closed system (maximum dissipation), or a driven stationary state in an open system (minimum dissipation), while the underlying Rayleighian functional is convex and minimized in both cases. We further examine the Reynolds transport theorem, a useful tool for fixing the energy rates of continua, and the role of inertia, recently argued to lie beyond the reach of Rayleighian principles. With inertia the driving force must be held fixed during the variation, which is less elegant than in the inertialess case, where it is a state function. Writing the energy rate as rates times driving forces still yields the correct equations, and an extremal statement survives whenever the frozen force does not depend on the varied rate. This holds directly for unsteady Stokes flow, and for Navier-Stokes flow after a change to material coordinates, in which the particle acceleration is the second time derivative of the state and the convective term, nonvariational in the Eulerian description, does not appear. The framework is tested on the Navier-Stokes family (Stokes, unsteady Stokes, and Navier-Stokes flow, in order of increasing complexity) and on dissipative Lagrangian mechanics.
We present SPOD-GPR, a non-intrusive, parametric, spatio-temporal reduced-order model (ROM) that combines Proper Orthogonal Decomposition (POD), Spectral POD (SPOD), and Gaussian Process Regression (GPR). A parametric POD first projects the high-dimensional state onto a low-dimensional subspace. The reduced trajectories are then divided into windowed sections, transformed into the frequency domain, and truncated to the most energetic frequencies. SPOD is then performed independently at each retained frequency, yielding temporally coherent modes and coefficients. The parametric dependence of the POD and SPOD bases is modelled by GPR on a local chart of the Grassmann manifold after Procrustes alignment. The SPOD modal coefficients are learned both as explicit functions of the parameters and through a parameter-conditioned autoregressive model that allows for temporal extrapolation beyond the training window. Gaussian processes are sampled jointly and pathwise via Matheron conditioning, ensuring that each predictive realization carries uncertainty quantification consistent with the posterior distribution. The method is tested on four parametric examples where the parameter values and the prediction horizon lie outside the training data: the two-dimensional anisotropic Kuramoto--Sivashinsky equation in a time-periodic and a travelling-wave regime, each parametrised by one of its two coefficients; flow past a cylinder parametrised by the Reynolds number; and flow past an elliptic obstacle parametrised by its semi-axes. In all cases, SPOD-GPR reduces the normalised mean absolute error compared to a POD-GPR baseline, and improves the calibration of the predictive uncertainty, especially in the fluid flow cases, where POD-GPR is severely overconfident. This opens the way to spectral, uncertainty-aware reduced-order models for many-query and real-time applications, such as digital twins.
Direct numerical simulations examine the phase-dependent turbulence dynamics of pipe flow driven by cardiac-cycle-inspired waveforms. Four cases with and without a low-flow plateau span two Reynolds-number ranges. The response follows four stages of turbulent--turbulent reorganisation. Strong acceleration suppresses transverse fluctuations and vortical activity. During acceleration and early deceleration, streamwise-normal-stress production recovers ahead of pressure--strain redistribution and radial momentum transport, allowing energy to accumulate in elongated streaks. In early deceleration, near-wall high-speed streaks develop predominantly varicose disturbances associated with Q4 sweeps, whereas outer-buffer-layer low-speed streaks develop predominantly sinuous disturbances associated with Q2 ejections. The transition exhibits distinct radial sequences: Q2-based local onset occurs progressively later towards the wall, whereas Q4-based onset is nearly simultaneous across the viscous sublayer and buffer layer, with peaks occurring successively outwards. The two responses overlap between the drag minimum and recovery maximum as redistribution intensifies markedly and multiscale momentum transport recovers rapidly. Turbulence then decays and its scales readjust. Where present, the plateau prolongs this decay and weakens the turbulence inherited by the next acceleration stage. It further delays local instability and lengthens the onset-to-peak intervals of the shear contributions. An unsteady friction decomposition links these phase-dependent dynamics to the drag response, highlighting the role of delayed turbulent shear recovery in sustaining the low-drag interval.
Most turbulent flows are sustained by inertia, which transfers energy across scales and produces characteristic power-law spectra. Yet dilute polymer solutions also exhibit turbulence when inertia is negligible: in this state, known as elastic turbulence, the velocity field remains spatially smooth while both the velocity and the polymer conformation exhibit robust power-law spectra. How a smooth flow, lacking a hierarchy of eddies, generates such multiscale structure has remained unexplained. Here we show that the missing mechanism is elastic frontogenesis, whereby the flow spontaneously organizes polymer stress into thin sheets that are stretched within their plane and compressed across their thickness, causing their amplitude to grow as they become thinner. Because the fronts are sheet-like structures, this amplification-thinning law uniquely determines the polymer conformation spectrum, while the Stokes response transfers the same geometric scaling to the velocity field, reconciling broadband spectra with spatial smoothness. Being geometric rather than model-specific, elastic frontogenesis provides an organizing principle for viscoelastic flows and the viscoelastic analogue of vortex stretching in inertial turbulence.
The velocity gradient tensor is a local quantity, but methods that obtain it by spatial differencing of the velocities of spherical tracers amplify measurement noise. We propose a method that reconstructs the velocity gradient tensor from the orientations of non-spherical tracers such as flakes and rods. Based on the Jeffery-type equation, we show that, for incompressible flow, all components of the velocity gradient tensor are uniquely determined from the orientations of four tracers and their time derivatives in the neighbourhood of a reference point (orientation-based method). The method uses neither tracer positions nor spatial differencing. It therefore requires orientation measurements only near the reference point, without requiring measurements throughout a three-dimensional region, and measurement noise is amplified only through temporal differencing. Consequently, the minimum achievable error scales as $σ^{1/2}$ with the noise level $σ$, compared with $σ^{1/3}$ for the positional PTV-plus-differencing approach (position-based method). Numerical tests further show that the vorticity is reconstructed more accurately than by the position-based method, whereas the strain rate is reconstructed with comparable accuracy. Finally, applying the method to spherical Couette flow, we find that the reconstruction fails where the neighbouring tracers align, which renders the inversion ill-conditioned, and that this degeneracy appears as the characteristic light-and-dark pattern familiar from flow visualization.
When heat and solute both destabilise a rapidly rotating fluid, do their buoyancies add? Codensity models assume so, giving a straight onset diagonal; spherical shells show curved and kinked thresholds. We study the Busse--Simitev rotating annulus, whose normal modes obey an exact complex cubic dispersion relation. Combining the Bilharz criterion with the neutral-curve envelope yields, in the rapid-rotation limit, an explicit swallowtail wave front, in Arnold's sense, that bends the codensity diagonal into smooth and then kinked thresholds as the diffusivity ratio decreases. Its cusp edges lie inside the unstable domain and its exposed kink ends at the swallowtail singularity; these are singularities of wave-number selection, without coalescence of temporal eigenvalues, and the swallowtail diffusivity ratio has analytical small- and large-Prandtl-number limits. At a kink, broad, fast thermal and narrow, slow compositional Rossby waves become neutral together; hidden arcs mark the birth and merger of unstable wave-number bands. The same geometry enters nonlinear selection: mean-scalar feedback factors through neutral-arc normals and scalar responses, so the mean-profile coexistence margin vanishes with the kink angle as the waves merge. Complete cubic coefficients predict robust coexistence at a reference corner and a narrow bistable interval near merger. Exact finite-rotation corners and periodic-cell simulations confirm the coexisting waves and their invasion thresholds.
We present \textit{Hydrograd} v2.0, a significantly improved differentiable 2D shallow-water model written in Python with JAX. The model solves the shallow-water equations on the unstructured meshes of production models such as SRH-2D and HEC-RAS 2D with a well-balanced, positivity-preserving finite-volume scheme that handles wetting and drying. Gradients are tailored to the type of problem. Steady states are found by Newton's method with an exact sparse Jacobian obtained by graph colouring, and their gradients follow from the implicit function theorem at the cost of one linear solve. Unsteady flows are differentiated in reverse mode with checkpointing, or with an implicit time integrator whose gradient costs about two forward runs. All gradients agree with finite differences, to about $10^{-8}$ for steady flows. On the Savannah River case, the steady water surface differs from that of SRH-2D by 4.7\,mm (root mean square). Steady closure learning was 96 times faster than differentiating through time marching and reached a lower loss. An inflow hydrograph was inverted from stage records with a peak error of 2\,\%. In an unsteady roughness calibration, a singular value decomposition analysis of the sensitivity matrix showed that the stage records at the gauges are 70 to 160 times more sensitive to the total resistance of the reach than to its distribution among the roughness zones, which explains why only a gauge layout with gauges in every zone recovered all five roughness values. On one GPU, a mesh of 512,000 cells takes 7.6\,s for a forward run and 51\,s for a gradient.
The Flash evaporation Riemann problem (FeRP) formulated the wave structures of flash evaporation in the vapor-liquid homogeneous equilibrium limit. However, physical flash evaporation involves finite-rate interphase heat and mass transfer, inducing the wave structure to deviate from the equilibrium state. To solve such a non-equilibrium Riemann problem, this paper first derives the governing equation of mechanical equilibrium evolution in a relaxation system, based on which an oscillation-free Homogeneous Relaxation Real Fluid (HRRF) finite volume method is developed. In addition to the Euler equations with mass relaxation, the HRRF method evolves advection equations for the internal energy--pressure affine parameters and the vapor volume fraction, in order to eliminate numerical pressure oscillations inherent in conservative methods while preserving phase-transition-induced pressure changes, and consistently degenerate into a completely pressure-equilibrium method in supercritical flows. The HRRF method is then employed to comprehensively evaluate how thermodynamic relaxation affects the FeRP solution. Results indicate that the non-classical expansion shock of the FeRP loses its traveling-wave characteristics with increasing relaxation time, continuously degenerating into a self-similar rarefaction wave. Compared with the equilibrium limit, an increased mass transfer relaxation time smooths the expansion branch of the RP solution, lowers the intermediate pressure and velocity, and induces a pseudo-splitting of the rarefaction wave as well as a pseudo-advance of the phase transition onset. Meanwhile, the heat transfer relaxation time primarily dictates the temperature difference in the vapor-liquid mixture, thereby modulating the liquid superheat and in turn influencing the mass transfer rate.
Helically stranded cables in aerial towing systems experience oblique inflow, but are often modeled as smooth cylinders, neglecting three-dimensional cross-stream coupling and surface texture. We perform large-eddy simulations and spectral proper orthogonal decomposition for a smooth cylinder and four stranded cables at four inclinations from $22.5^\circ$ to $90^\circ$, all stationary, with $Re \approx 1.0 \times 10^4$. Aerodynamic differences are governed mainly by cross-sectional symmetry and strand-scale groove depth. Force-source decomposition shows that surface pressure contributes 94.9%--98.3% of the mean normal force and dominates the mean lateral force. High-symmetry, quasi-circular envelopes maintain balanced circumferential pressure, whereas deep-grooved, low-symmetry tri-lobed geometries produce asymmetric separation and a persistent lateral pressure imbalance. In the wake, high-symmetry cables sustain spanwise-synchronized Karman shedding at moderate-to-high inclinations, concentrating fluctuation energy in a dominant low-rank mode. Deep grooves disrupt shedding coherence, deflect the recirculation zone laterally, and spread energy across broadband scales. Mean lateral force and lateral-force fluctuations are therefore decoupled: the quasi-circular $(1 \times 19)$ cable has minimal mean lateral force but strong periodic fluctuations, while the deep-grooved $(3 \times 7)$ cable has weak fluctuations but appreciable mean lateral force. These results support construction-dependent cable selection and aerodynamic load modeling.
Blood flow and elastic vessel-wall motion interact with each other: fluid stress drives wall deformation, while the moving wall changes the domain in which blood flows. Given noisy flow images, we study the inverse problem of reconstructing the coupled, unsteady three-dimensional motions. We solve the reconstruction through Gauss--Seidel-type alternating fluid and geometry updates. For the fluid, we use physics-informed networks to recover velocity and pressure by fitting the observations while enforcing the incompressible Navier--Stokes equations on the moving domain, which is described by an arbitrary Lagrangian--Eulerian (ALE) map. For the geometry, a reference-domain correction and sequential registration recover the vessel shape and wall displacement. Interface traction continuity and an elastic momentum equation regularize the registration so that the recovered displacement is also physics-informed. The networks evaluate fields and physical residuals at sampled space--time points without remeshing the moving domain. Experiments on a converging channel and a patient-like aorta demonstrate reconstruction of flow and wall motion from noisy data. Ablation and parameter-sensitivity studies examine the role of structure regularization and imperfect physical coefficients. A further experiment with a three-element Windkessel outlet condition demonstrates how the framework accommodates a more physiologically realistic boundary setting.
Charged amphiphilic surfactants are versatile polymorphic molecules which can form supramolecular structures, such as wormlike micelles. We study the rheology and non-Newtonian fluid dynamics of the cationic surfactant trimethyltetradecylammonium chloride (TTAC) combined with the aromatic counter-ion sodium salicylate (NaSal). In the non-dilute regime, this solution exhibits a linear viscoelastic response consistent with a single-mode Jeffreys model, and in capillary pinch-off experiments, fluid filaments form extended cylindrical threads, indicating strong strain-stiffening. In the dilute regime the linear viscoelasticity and extensional rheology become immeasurably small, but the TTAC/NaSal mixture develops flow-induced structures at high shear rates, as the result of an end-to-end concatenation of micelles produced by shear-induced flow alignment. These flow-induced structures augment both the measured viscosity and the viscous drag in laminar flows at high shear rates (by as much as 50 %). In turbulent Taylor-Couette flow, the flow-induced structures produce a sudden and dramatic destabilizing influence on the stable toroidal Taylor vortices, resulting in upwards of a 200 % enhancement in turbulent drag. However, with increasing Reynolds number, the turbulent flow transitions smoothly to a state of drag reduction, generating up to a 30 % reduction in fluid frictional drag but without the shear-induced degradation that plagues drag reduction strategies based on solutions of high molecular weight polymer, such as polyacrylamide. These results show that surfactants are a unique and effective modifier of fluid frictional drag in both laminar and turbulent flow regimes, capable of producing both significant levels of drag enhancement and drag reduction by tuning the composition and operating conditions.
Targeting a prescribed flow configuration rather than an aerodynamic coefficient is challenging because the desired state may lie outside the natural attractor. We demonstrate cluster-based control in a learned vorticity representation for two-dimensional flow over a NACA0015 airfoil at $Re=500$ and $α=15^\circ$, using a localized streamwise momentum source near the leading edge. An actuation-enriched ensemble is compressed with a variational autoencoder and partitioned into macrostates. We optimize a continuous feedback law interpolating macrostate-dependent amplitudes without an explicit latent dynamical model. For an attached target, the optimized law approaches constant forcing, arrests vortex shedding and maintains the attached state. For a strongly separated target, the law instead establishes attachment, withdraws actuation so that the natural separation transient carries the flow through the target, and restores forcing to recapture attachment. Over the evaluation horizon, the closed-loop trajectory revisits the target approximately every three shedding periods. Companion open-loop calculations isolate the arrest, release and recapture processes.
We study the three-dimensional incompressible Navier--Stokes equation on the cubic torus with a constant odd-viscosity operator given by the vertical derivative of the curl. Arbitrary real, mean-zero, divergence-free initial data and an independent steady force, both in the second-order Sobolev space, admit a unique global strong solution when the magnitude of one fixed odd-viscosity coefficient is at least a norm-uniform sufficient threshold. In the unforced case the solution is smooth at positive times. The qualitative argument adapts the resonant-averaging framework for rapidly rotating fluids to the different odd-viscous dispersion. Weighted lattice counts and signed corrected-energy estimates yield a sufficient bound for the odd-to-ordinary viscosity ratio that is polynomial in the normalized initial-data radius at each fixed normalized force radius, together with a separate bound in a specified force-dominated regime. No optimality is claimed for the polynomial exponents or prefactors; determining the sharp viscosity threshold for the original equation remains open. Direct kinetic and weighted enstrophy barriers also give regularity for every coefficient on explicit bounded data and force regions, including the closed unforced ball of homogeneous second-order Fourier Sobolev radius ninety-seven fortieths times the viscosity and seven closed forced rectangles. We specify how finite certificates enter the analytic argument and distinguish minimization of a sufficient formula from sharpness for the original equation. All three main theorems are formally verified in Lean 4 with Mathlib.
Hierarchical Bayesian models pool measurements, but with the forward model in the likelihood every posterior evaluation costs one forward solve per cell. To reduce this cost, we introduce a simulation-based inference framework in which a triangular transport-map likelihood, trained once per design, serves every cell. We apply it to infer wetting properties from photographs of a stable hole in a thin liquid film. The hole's radius depends, through the nonlinear Young-Laplace equation, on the contact angle and capillary length. One photograph gives rank-one information about them, so we pool many in a hierarchical model whose cells, each a substrate at one volume setting, share the apparatus geometry and the noise scale. A design fixes the number of substrates, the volume settings and the photographs per cell. The map's posteriors agree with direct sampling of the forward model in the wetting properties at a synthetic validation design, and less closely at a larger synthetic scaling design. At the validation design, simulation-based calibration detects no departure for any parameter. There, the map needs several hundred times fewer forward-model evaluations than direct sampling with a gradient-free ensemble sampler for a single posterior, and about five orders of magnitude fewer over a calibration study of 200 posteriors. On real photographs its posterior agrees with direct sampling and is consistent with the measured radii. Trained offline, the amortized map turns new photographs into a posterior with no forward solve. It is thus a major step toward measuring wetting properties in situ without disturbing the liquid.
We derive and numerically validate an odd Darcy law for chiral active fluids in porous media, in which an antisymmetric permeability tensor induces a flow transverse to the applied pressure gradient. We identify two distinct mechanisms leading to this odd response: odd viscosity in the fluid or the nonlinear coupling of edge currents to inertia. We then show how the transverse transport depends on porosity, boundary friction, and chirality.
Unsteady compressible simulations with high-order methods on unstructured meshes require boundaries that let acoustic, entropy and vorticity waves leave the domain without spurious reflection. For the spectral difference (SD) method, the Navier-Stokes characteristic boundary conditions (NSCBC) have only been formulated on quadrilaterals and hexahedra. Their strong-form formalism relies on a line of flux points normal to the boundary face, which does not exist on triangles and tetrahedra. This paper extends the NSCBC to the SD method on simplex elements, including multi-species reacting gases. The characteristic analysis is performed in the physical frame of each flux point of the boundary face, and the normal derivative of the normal flux is evaluated using the Raviart-Thomas divergence operator. The boundary fluxes are recovered from a linear system local to each boundary element, whose matrix depends on the geometry only. The treatment is compared with the tensor-product formulation at a comparable number of degrees of freedom on a series of test cases covering acoustic, vortical, viscous and reacting flows on triangles and tetrahedra. In all cases, the results are in close agreement with those obtained on quadrilaterals and hexahedra, so that the NSCBC become available for the SD method on simplex meshes, which are well suited to complex geometries.
We revisit the protocol known as Exponential Shear (ES), first proposed by Doshi and Dealy in 1986. This is a transient strong flow that can be readily imposed in the laboratory using commercial strain-controlled (separate motor and transducer) rheometers which accurately measure both the (exponentially) time-varying torque and normal force exerted by the sample. Numerous studies over the years have proposed competing definitions of an appropriate material function for reporting data obtained from the exponentially time-varying imposed deformation. We revisit this unresolved challenge and carefully reconsider the exponential material kinematics involved and the proper definition of the corresponding material functions. We show that by understanding the time-varying eigenvalues and eigenvectors of both the shearing kinematics and the evolving state of stress in the sample it is, in fact, possible to use an exponential shear test to quantitatively map out \textit{both} the full flow curve (of shear viscosity vs. shear rate) and the steady planar extensional viscosity curve as a function of deformation rate. These findings are supported by numerical simulations of a number of nonlinear viscoelastic constitutive models, asymptotic analysis, as well as experiments with several polymeric liquids.
Embedded microfluidic cooling is a promising technology for dissipating the intense heat generated by high-performance electronic chips. However, the design and optimization of such coolers impose prohibitive computational costs when using conventional numerical solvers. Here, we present a parameterized physics-informed neural network (PINN) as a deep learning framework for rapidly and accurately predicting the thermohydraulic performance of embedded manifold microchannels. By incorporating four geometric dimensions along with three spatial coordinates as the network inputs, the proposed model enables multiphysics prediction across a continuous design space. Remarkably, our predicted velocity and temperature fields differ from finite-element simulations by relative errors of only 8.1% and 5.6%, respectively, while reducing the computational time and cost by over four orders of magnitude. Key to this framework is a composite strategy combining geometric ansatz embedding, the SOAP optimizer, and a thermal-conductivity-based curriculum learning scheme, which helps overcome the convergence challenges inherent in training PINNs for complex three-dimensional flows. Notably, by leveraging a highly optimized JAX backend, the network achieves convergence on a single consumer-grade GPU in less than four hours. Our approach offers a physics-consistent, real-time surrogate model that can substantially accelerate innovations in advanced microfluidic cooling of future electronics.
The scaling of electrical conductance with the number of nanopores, $N$, in an array is of fundamental importance for applications in membrane separation and ion transport across biological membranes, and has been the subject of longstanding disagreement in the field. Combining finite element calculations with theory, we resolve this disagreement for pores with negligible surface charge by showing that two previously conflicting predictions are each correct in their own limiting reservoir geometry: cuboidal reservoirs commensurate with the interpore spacing, for which the conductance scales linearly with $N$, and infinite reservoirs, for which the scaling is sublinear in $N$ for closely spaced pores. We show the transition between these limits with increasing reservoir size and, for infinite reservoirs, we show that a simple analytical model from electrical contact theory is mathematically equivalent to an existing point-charge capacitance model for nanopore arrays and validate its applicability against finite element and method-of-moments calculations. Using this model, we derive a characteristic coupling length, $\ell_\mathrm{c,95}$, the interpore spacing at which the array conductance reaches $95\%$ of that of the same number of independent pores. $\ell_\mathrm{c,95}$ scales linearly with pore diameter $d$, and grows with $N$ as $\sim\log N$ for linear arrays and $\sim\sqrt{N}$ for square arrays, so that pore coupling persists to surprisingly large separations in large arrays. Our work provides a general framework for predicting nanopore array conductance as a function of pore arrangement and reservoir geometry, of importance to nanofluidics research and the targeted design of functional devices.
Colloid transport and retention in environmental granular porous media underlie groundwater resource protection from, and remediation using, biological and non-biological nano- and micro-particles (colloids). During transport, colloids periodically intercept grain surfaces and translate along them at separation distances across several tens of nm corresponding to the secondary minimum between the opposing forces of longer-range of van der Waals (VDW) attraction and shorter-range electric double layer (EDL) repulsion. Relatively weak diffusive translation of colloids supports their continued translation in secondary minimum association until they are expelled to bulk fluid at rear flow stagnation zones on grain leeward sides. We herein show that this dynamic focuses colloids into a network of fluid streamlines (incomplete mixing lamellae) that enhances their subsequent likelihood of interception with downgradient grains, i.e., increases the so-called collector efficiency (typically denoted by Greek letter eta). Interception focusing of colloids in granular porous media explains well-observed outcomes of colloid transport under environmental conditions, including nonmonotonic nonexponential profiles of retained colloids and the non-Markovian nature of extended tailing of low colloid concentrations following initial elution. Interception focusing is a fundamental process that distinguishes colloids from solutes and renders inapplicable to colloids the assumption that is often invoked for solutes of complete pore scale mixing in granular porous media.