2024-10-30 | | Total: 5
A posteriori analysis based upon a recently proposed non-dissipative large-eddy simulation framework for transcritical wall-bounded turbulence has been carried out. Due to the complexities arisen in such flows, the discretization requires kinetic-energy- and pressure-equilibrium-preservation schemes to yield stable and non-dissipative scale-resolving simulations. On the basis of this framework, the objectives are to (i) compute wall-resolved and wall-modeled large-eddy simulations of a high-pressure transcritical turbulent channel flow, and (ii) assess the thermofluid performance with respect to a direct numerical simulation. In this regard, three different subgrid-scale stress tensor models have been considered, together with models for the unresolved scales of the filtered pressure transport equation and equation of state. In terms of wall-modeling, models based on the "standard law of the wall" and velocity-temperature coupled approaches have been assessed. The results show that for the wall-resolved approaches, the subgrid-scale stress tensors examined slightly deviate from the time-averaged velocity and temperature reference profiles. In terms of bulk performance, it has been found that the Nusselt number and skin-friction coefficient are relatively well captured at the cold and hot walls, respectively. While heat transfer phenomena are fairly well reproduced, particularly the Prandtl and heat flux trends along wall-normal direction. Additionally, the wall-modeled strategies improve the recovery of first-order statistics, although they do not attain the profile in the log-law region. However, they enhance the wall metrics and the heat flux prediction. It is, thus, concluded that dedicated efforts by the research community are needed to improve the prediction accuracy of existing subgrid-scale and wall models for wall-bounded transcritical turbulence.
The paper proposes a physically consistent numerical discretization approach for simulating viscous compressible multicomponent flows. It has two main contributions. First, a contact discontinuity (and material interface) detector is developed. In those regions of contact discontinuities, the THINC (Tangent of Hyperbola for INterface Capturing) approach is used for reconstructing appropriate variables (phasic densities). For other flow regions, the variables are reconstructed using the Monotonicity-preserving (MP) scheme (or Weighted essentially non-oscillatory scheme (WENO)). For reconstruction in the characteristic space, the THINC approach is used only for the contact (or entropy) wave and volume fractions and for the reconstruction of primitive variables, the THINC approach is used for phasic densities and volume fractions only, offering an effective solution for reducing dissipation errors near contact discontinuities. The second contribution is the development of an algorithm that uses a central reconstruction scheme for the tangential velocities, as they are continuous across material interfaces in viscous flows. In this regard, the Ducros sensor (a shock detector that cannot detect material interfaces) is employed to compute the tangential velocities using a central scheme across material interfaces. Using the central scheme does not produce any oscillations at the material interface. The proposed approach is thoroughly validated with several benchmark test cases for compressible multicomponent flows, highlighting its advantages. The numerical results of the benchmark tests show that the proposed method captured the material interface sharply compared to existing methods.
The presence of salt in seawater significantly affects the melt rate and morphological evolution of ice. This study investigates the melting process of a vertical cylinder in saline water using a combination of laboratory experiments and direct numerical simulations. The two-dimensional (2D) direct numerical simulations and 3D experiments achieve thermal Rayleigh numbers up to $\text{Ra}_{T}= \mathcal{O}\left(10^{9}\right)$ and saline Rayleigh numbers up to $\text{Ra}_{S}=\mathcal{O}\left(10^{12}\right)$. Some 3D simulations of the vertical ice cylinder are conducted at $\text{Ra}_{T}= \mathcal{O}\left(10^{5}\right)$ to confirm that the results in 2D simulations are qualitatively similar to those in 3D simulations. The mean melt rate exhibits a non-monotonic relationship with ambient salinity, initially decreasing towards the point that thermal and saline effects compensate each other, and then increasing again. Based on the ambient salinity, the flow can be categorized into three regimes: temperature-driven flow, salinity-driven flow, and thermal-saline competing flow. In the temperature-driven and competing flow regimes, we find that the mean melt rate follows a $\text{Ra}_{T}^{1/4}$ scaling. In contrast, in the salinity-driven flow regime, we see a transition from a $\text{Ra}_{T}^{1/4}$ to a $\text{Ra}_{T}^{1/3}$ scaling. Additionally, the mean melt rate follows a $\text{Ra}_{S}^{1/3}$ scaling in this regime. The ice cylinder develops distinct morphologies in different flow regimes. In the thermal-saline competing flow regime, distinctive scallop (dimpled) patterns emerge along the ice cylinder due to the competition between thermal buoyancy and saline buoyancy. We observe these scallop patterns to migrate downwards over time, due to local differences in the melt rate, for which we provide a qualitative explanation.
Tropical cyclones (TCs) are powerful, natural phenomena that can severely impact populations and infrastructure. Enhancing our understanding of the mechanisms driving their intensification is crucial for mitigating these impacts. To this end, researchers are pushing the boundaries of TC simulation resolution down to scales of just a few meters. However, higher resolution simulations come with significant computational challenges, increasing both time and energy costs. Adaptive mesh refinement (AMR) is a technique widely used in computational fluid dynamics but has seen limited application in atmospheric simulations. This study explores the use of h-adaptive grids using the spectral element discretization technique to accelerate TC simulations while allowing very high resolutions in certain parts of the domain. By applying AMR to a rapidly intensifying TC test case, we demonstrate that AMR can replicate the results of uniform grid simulations in terms of mean and local wind speed maxima while dramatically reducing computational costs. We show that AMR can speed up dry TC simulations by a factor of 2 - 13 for the set of tested refinement criteria. Additionally, we show that TC intensity changes as resolution is increased and that AMR can deliver high-resolution simulations at the cost of coarser static simulations. Our findings indicate that AMR and spectral element methods are promising tools for enhancing TC simulations.
Computational physics simulation can be a powerful tool to accelerate industry deployment of new scientific technologies. However, it must address the challenge of computationally tractable, moderately accurate prediction at large industry scales, and training a model without data at such large scales. A recently proposed component reduced order modeling (CROM) tackles this challenge by combining reduced order modeling (ROM) with discontinuous Galerkin domain decomposition (DG-DD). While it can build a component ROM at small scales that can be assembled into a large scale system, its application is limited to linear physics equations. In this work, we extend CROM to nonlinear steady Navier-Stokes flow equation. Nonlinear advection term is evaluated via tensorial approach or empirical quadrature procedure. Application to flow past an array of objects at moderate Reynolds number demonstrates $\sim23.7$ times faster solutions with a relative error of $\sim 2.3\%$, even at scales $256$ times larger than the original problem.