2024-12-02 | | Total: 15
The simulation of high Reynolds number (Re) separated turbulent flows faces significant problems for decades: large eddy simulation (LES) is computationally too expensive, and Reynolds-averaged Navier-Stokes (RANS) methods and hybrid RANS-LES methods often provide unreliable results. This has serious consequences, we are currently unable to reliably predict very high Re regimes, which hampers applications and our understanding of turbulence structures. The paper reports the advantages of a strict mathematical approach, continuous eddy simulation (CES), to derive partially resolving turbulence models. In contrast to popular hybrid RANS-LES, this minimal error approach includes a dynamic modification of the turbulence model in response to the actual flow resolution: the model can increase (decrease) its contribution to the simulation in dependence of a low (high) flow resolution. This property is the essential requirement to seamlessly cover RANS and LES regimes. The CES modeling approach offers essential advantages regarding its functionality: basically, it is independent of a variety of simulation settings applied in popular hybrid RANS-LES to improve the model performance. In addition, the CES computational cost can be below the cost of other hybrid RANS-LES and LES by orders of magnitude. Essential simulation performance advantages of CES simulations are described here with respect to three complex flow applications: periodic hill flows at high Reynolds number, the NASA wall-mounted hump flow, and the Bachalo & Johnson axisymmetric transonic bump flow.
The discovery of the law of the wall, the log-law including the von Karman constant, is seen to be one of the biggest accomplishments of fluid mechanics. However, after more than ninety years there is still a controversial debate about the validity and universality of the law of the wall. Clarity about this question matters: in absence of alternatives, a reliable and universal theory involving the law of the wall is needed to provide essential guideline for the validation of theory, computational methods, and experimental studies of very high Reynolds number (Re) flows. The paper presents an analysis of concepts used to derive controversial conclusions. It is shown that nonuniversality is a consequence of simplified modeling concepts, which leads to unrealizable models. On the other hand, realizability implies universality: models in consistency with physical requirements do not need to be adjusted to different flows. There are essential advantages of a universal law of the wall: it enables the design of accurate turbulence models and it provides a bridge between finite Re observations and asymptotic structural theories of turbulence.
This study investigates conjugate heat transfer analysis of a rocket nozzle featuring a charring ablative thermal protection system. The study is carried out by weakly coupling a commercial computational fluid dynamics (CFD) flow solver Fluent with an in-house material thermal response solver. The coupling of the two solvers is carried out by exchanging the boundary conditions at the fluid-solid interface using a non-iterative approach. Validation of the numerical framework is performed. A blowing correlation is used to model the pyrolysis gas behavior at the interface. Results highlight a significant surface heating and ablation particularly at the nozzle throat.
Gravity currents are a ubiquitous density driven flow occurring in both the natural environment and in industry. They include: seafloor turbidity currents, primary vectors of sediment, nutrient and pollutant transport; cold fronts; and hazardous gas spills. However, while the energetics are critical for their evolution and particle suspension, they are included in system scale models only crudely, so we cannot yet predict and explain the dynamics and run-out of such real-world flows. Herein, a novel depth-averaged framework is developed to capture the evolution of volume, concentration, momentum, and turbulent kinetic energy from direct integrals of the full governing equations. For the first time, we show the connection between the vertical profiles, the evolution of the depth-averaged flow, and the energetics. The viscous dissipation of mean-flow energy near the bed makes a leading order contribution, and an energetic approach to entrainment captures detrainment of fluid through particle settling. These observations allow a reconsideration of particle suspension, advancing over 50-years of research. We find that the new formulation can describe the full evolution of a shallow dilute current, with the accuracy depending primarily on closures for the profiles and source terms. Critically, this enables accurate and computationally efficient hazard risk analysis and earth surface modelling.
Analytic bounds on the projection of the global attractor of 2D incompressible turbulence in the palinstrophy--enstrophy plane [Dascaliuc, Foias, and Jolly 2005, 2010] are observed to vastly overestimate the values obtained from numerical simulations. This is due to the lack of a good estimate for the inner product $(\cal{B}(u,u),A^2u)$ of the advection term and the biLaplacian. Sobolev inequalities like Ladyzhenskaya or Agmon's inequalities yield an upper bound that we show is not sharp. In fact, for statistically isotropic turbulence, the expected value of $(\cal{B}(u,u),A^2u)$ is zero. The implications for estimates on the behaviour of the global attractor are discussed.
We present and explore a new shock-capturing particle hydrodynamics approach. We enhance a commonly used SPH discretization with Roe's approximate Riemann solver where we use slope-limited linear reconstruction in the dissipative terms of the solver. All gradients of the method are calculated with linearly reproducing kernels that are constructed to enforce the two lowest order consistency relations so that constant and linear functions can be recovered to machine precision. We scrutinize the method with a set of challenging 3D benchmark problems ranging from shocks over instabilities to Schulz-Rinne-type vorticity creating shocks. In all of them we find excellent agreement with analytic/reference solutions.
Mesh generation is essential for accurate and efficient computational fluid dynamics simulations. To resolve critical features in the flow, adaptive mesh refinement (AMR) is routinely employed in certain regions of the computational domain, where gradients or error estimates of the solution are often considered as the refining criteria. In many scenarios, however, these indicators can lead to unnecessary refinement over a large region, making the process a matter of trial and error and resulting in slow convergence of the computation. To this end, we propose a heuristic strategy that employs the residuals of the governing partial differential equations (PDEs) as a novel criterion to adaptively guide the mesh refining process. In particular, we leverage on the physics-informed neural networks (PINNs) to integrate imprecise data obtained on a coarse mesh and the governing PDEs. Once trained, PINNs are capable of identifying regions of highest residuals of the Navier-Stokes/Euler equations and suggesting new potential vertices for the coarse mesh cells. Moreover, we put forth two schemes to maintain the quality of the refined mesh through the strategic insertion of vertices and the implementation of Delaunay triangulation. By applying the residuals-guided AMR to address a multitude of typical incompressible/compressible flow problems and comparing the outcomes with those of gradient-based methods, we illustrate that the former effectively attains a favorable balance between the computational accuracy and cost.
An automatically differentiable, high-order non-oscillatory finite volume shallow water dynamic core has been constructed on a cubed sphere grid. This dynamic core has four advantageous properties: high order accuracy, essential non-oscillation, mass conservation, and scalability. Besides, the code development is based on PyTorch, enabling the model to run seamlessly on both CPU and GPU, while naturally possessing the capability of automatic differentiation. We named the new dynamic core as High Order Prediction Environment (HOPE). The spatial reconstruction is based on the two-dimensional tensor product polynomial and the genuine two-dimensional Weighted Essentially Non-Oscillatory (WENO) scheme. A novel panel boundary approach ensures that the accuracy can reach arbitrary order. The algorithm has very high degree of compatibility with GPU architecture, allowing the computational overhead to be mitigated through the utilization of GPU. The LMARS (Low Mach number Approximate Riemann Solver) scheme is adopted as Riemann solvers to determine fluxes on the Gaussian points on edges. Flux across the interface between each cell edge is computed using linear Gaussian quadrature, and the tendency of prognostic values is obtained by summing all fluxes around the cell. This shallow water dynamic core demonstrates outstanding performance in ideal shallow water test cases. In the steady-state geostrophic flow, the 11th order scheme reduces errors to nearly round-off error even on coarse grids (2 degrees). Furthermore, the Rossby-Haurwitz wave preserves its shape for over 90 days. To test the non-oscillation property, we designed a cylinder dam break case. The WENO approach effectively suppresses non-physical oscillation, and the genuine two-dimensional reconstruction exhibits significantly better isotropy than the dimension-by-dimension scheme.
Solving the Bhatnagar-Gross-Krook (BGK) equation with a stochastic particle approach enables efficient and flexible simulations of flows in the transition regime, between continuum and free molecular flow. However, the usual first-order operator splitting between particle movement and relaxation imposes restrictions on the time step, causing the computational cost to increase with the gas density. The Crank-Nicolson stochastic particle BGK (CN-SPBGK) method is introduced here as an advanced particle-based kinetic solver designed for multi-scale gas flow simulations. This method integrates the BGK equation with second-order accuracy across all Knudsen number regimes without requiring additional parameters, while asymptotically preserving the Navier-Stokes flux in the continuum regime. Comparisons with pre-existing particle BGK methods are conducted on several test cases, with CN-SPBGK demonstrating more consistent convergence and accuracy.
Specializing in high-speed testing, Holloman High-Speed Test Track (HHSTT) uses water braking to stop vehicles on the test track. This method takes advantage of the higher density of water, compared to air, to increase braking capability through momentum exchange by increasing the water content in that section at the end of the track. By studying water braking using computational fluid dynamics (CFD), the forces acting on tracked vehicles can be approximated and prepared before actual testing through numerical simulations. In this study, emphasis will be placed on the brake component of the tracked sled, which is responsible for interacting with water to brake. By discretizing a volume space around our brake, we accelerate the water and air to simulate the brake coupling relatively. The multiphase flow model uses the governing equations of the gas and liquid phases with the finite volume method to perform 3D simulations. By adjusting the air and water inlet velocity, it is possible to simulate HHSTT sled tests at various operating speeds.
We consider traveling waves on a surface of an ideal fluid of finite depth. The equation describing Stokes waves in conformal variables formulation are referred to as the Babenko equation. We use a Newton-Conjugate-Gradient method to compute Stokes waves for a range of conformal depths from deep to shallow water. In deep water, we compute eigenvalues of the linearized Babenko equation with Fourier-Floquet-Hill method. The secondary bifurcation points that correspond to double period bifurcations of the Stokes waves are identified on the family of waves. In shallow water, we find solutions that have broad troughs and sharp crests, and which resemble cnoidal or soliton-like solution profiles of the Korteweg-de Vries equation. Regardless of depth, we find that these solutions form a $2\pi/3$ angle at the crest in the limit of large steepness.
A hybrid sharp-interface immersed-boundary/front-tracking (IB/FT) method is developed for interface-resolved simulation of evaporating droplets in incompressible multiphase flows. A one-field formulation is used to solve the flow, species mass fraction and energy equations in the entire computational domain with appropriate jump conditions at the interface. An image point and ghost cell methodology is coupled with a front-tracking method to achieve an overall second order spatial accuracy for the mass fraction boundary condition on the droplet surface. The immersed-boundary method is also extended to simulate mass transfer from a solid sphere in a convective environment. The numerical method is first validated for the standard benchmark cases and the results are found to be in good agreement with analytical solutions. The method is shown to be overall second order accurate in space. Employing a moving reference frame methodology, the method is then applied to simulate evaporation of a deformable droplet in a convective environment and the results are compared with the existing evaporation models widely used in spray combustion simulations.
Introduction: the Navier-Stokes equations are essential in fluid dynamics, describing the motion of fluids like liquids and gases. Solving these equations, especially in complex flows and high-Reynolds-number regimes, is a significant challenge. Numerical simulations provide some insights, but often under restrictive assumptions that limit applicability. Recent geometric and algebraic methods have emerged, focusing on the equations' structure, yet questions about the uniqueness and stability of weak solutions persist. Objective: this paper aims to reformulate the Navier-Stokes equations in covariant form and develop new equations that facilitate the search for potential solutions, emphasizing symmetries. Geometric Approach: a covariant formulation of the Navier-Stokes equations is presented, applying a Fourier transform on a bounded manifold and seeking smoothness and viable solutions through convergence of manifold elements. Transformations: the study examines transformations between manifolds, investigating symmetries and interpretations related to homeomorphisms, isometries, and diffeomorphisms, including inertial frames of reference. Discussion and Conclusions: this study introduces a geometric reformulation of the Navier-Stokes equations, proposing new equations to enhance convergence and smoothness of solutions. It presents a novel class of solutions and transformations, with significant interdisciplinary connections. Further simulations, experimental validation, and ongoing development are essential to broaden the applicability of these equations and their solutions.
Recent clinical results indicate that aberrant erythrocyte aggregation in hematological disorders is accompanied by endothelial damage and glycocalyx disruption, but the underlying biophysical mechanisms remain unclear. This study uses direct computational modeling to explore how red blood cell (RBC) aggregation impacts shear stress in small blood vessels, highlighting the increased risk of vascular damage. RBC aggregation creates a heterogeneous distribution, leading to variations in the cell-free layer thickness and fluctuating wall shear stress, especially near vessel walls. This effect aligns with experimental findings on endothelial disruption linked to RBC clustering near the wall, potentially reducing the protective glycocalyx layer. The power spectral density analysis of wall shear stress fluctuations reveals that, with RBC aggregation, there is a distinct peak near frequency f = 0.04, indicating increased fluctuations due to aggregated RBC clusters traveling close to the vessel wall. The presence of aberrant cells in blood disorders, modeled here by sickle cells, further amplifies these effects, as aggregation-enhanced margination drives sickle cells closer to vessel walls, exacerbating shear stress fluctuations and increasing the likelihood of vascular injury and inflammation. Simulations show that curved vascular geometry, with curvature accentuating RBC clustering near vessel walls, intensifies aggregation-induced wall shear stress fluctuations and increases the risk of vascular damage, particularly in sickle cell disease where sickle cells marginate closer to the wall.
Intermittency as it occurs in fast dynamos in the MHD framework is evaluated through the examination of relations between normalized moments at third order (skewness S) and fourth order (kurtosis K) for both the velocity and magnetic field, and for their local dissipations. As investigated by several authors in various physical contexts such as fusion plasmas (Krommes (2008)), climate evolution (Sura & Sardeshmukh (2008)), fluid turbulence or rotating stratified flows (Pouquet et al. (2023)), approximate parabolic K(S) ~ alpha laws emerge whose origin may be related to the applicability of intermittency models to their dynamics. The results analyzed herein are obtained through direct numerical simulations of MHD flows for both Taylor-Green and Beltrami ABC forcing at moderate Reynolds numbers, and for up to 314000 turn-over times. We observe for the dissipation 0.2 < alpha <3.0, an evaluation that varies with the field, the forcing, and when filtering for high-skewness intermittent structures. When using the She & Leveque (1994) intermittency model, one can compute alpha analytically; we then find alpha = 2.5, clearly differing from a (strict) parabolic scaling, a result consistent with the numerical data.