We investigate the controlled K-type breakdown of a flat-plate boundary-layer with highly non-ideal supercritical fluid. Direct numerical simulations are performed at a Mach number of M_∞=0.2 for one subcritical (liquid-like regime) temperature profile and one strongly-stratified transcritical (pseudo-boiling) temperature profile with slightly heated wall. In the subcritical case, the formation of aligned Λ-vortices is delayed compared to the reference ideal-gas case of Sayadi et al. (J. Fluid Mech., vol. 724, 2013, pp. 480–509), with steady longitudinal modes dominating the late-transitional stage. When the wall temperature exceeds the pseudo-critical temperature, streak secondary instabilities lead to the simultaneous development of additional hairpin vortices and near-wall streaky structures near the legs of the primary aligned Λ-vortices. Nonetheless, transition to turbulence is not violent and is significantly delayed compared to the subcritical regime.

This study assesses the technical feasibility of nuclear propulsion for naval vessels by investigating the dynamics of a Very High Temperature Reactor combined with a supercritical carbon dioxide recompression cycle. By applying a dynamic model that includes the reactor, heat exchangers, and turbomachinery, the power dynamics of a nuclear energy conversion system are compared with those of common prime movers on a naval vessel. Results show that the turbine bypass valve, in combination with the dump cooler, enables a power ramp of 90%/min. During this power ramp, the reactor temperature stays within safety limits and experiences temperature variations of less than 25 degrees, the shaft speed remains stable with deviations of less than 0.25% RPM, and turbomachinery performs within design limits regarding temperature, pressure ratio and mass flow rate. However, the current turbine bypass valve design, which maintains a stable reactor output, results in a low overall cycle efficiency at part load. Furthermore, temperature- and pressure gradients of up to ±1.48 °C/s and ±0.38 bar/s occur within the heat exchangers during the power transient, which could affect the integrity of the materials. Further research could focus on a design that limits the thermal integrity concerns within the heat exchanger, and could implement energy storage capabilities to optimize the waste heat of the cycle during part-load.

Supercritical natural circulation loops (NCLs) promise passive cooling for critical systems like nuclear reactors and solar collectors, eliminating the need for mechanical pumps. However, instabilities similar to those seen in two-phase systems can emerge in supercritical NCLs, leading to undesirable oscillatory behaviour, marked by system-wide fluctuations in density, temperature, pressure, and flow rate. This study investigates the stability of NCLs at supercritical pressures (73.7≤p≤110.0bar) using CO₂ in an experimental setup with vertical cooling and vertically adjustable heaters to control convective flow rates and to oppose flow reversal. Oscillations were found to originate in the heater of the NCL, and demonstrated a high sensitivity to the thermodynamic state and proximity to the pseudo-critical line of the system. Increased mass flow rates and added resistance upstream of the heater suppressed the oscillations, while increased pressures and reduced heating rates dampened them. A static model which takes into account the non-ideality of the heat exchangers is introduced to assess the presence of multiple steady states. The system is concluded to be statically stable, and the oscillations are considered to be dynamically induced. In particular, the modulation of the NCL velocity by the traversal of the current oscillations in density is assumed to periodically re-incite non-ideality in the heater. These findings intend to refine our understanding of the stability boundaries in NCLs, to ensure a safer operation of prospective passive cooling and circulation systems employing fluids at supercritical pressure.

FluidsDraskic, M.Westerweel, J.Pecnik, R. display sharp, non-linear variations of thermodynamic properties when they are heated at a supercritical pressure. As such, near-pseudo-critical heat transfer is often characterized by large variations in density, leading to sharp near-wall accelerations or strong stratifications when buoyancy becomes dominant. We study the modulation of heat transfer and turbulence by non-negligible buoyancy in such property-variant flows, for the development of near-pseudo-critical heat exchangers for supercritical energy conversion systems. In particular, a liquid-like, horizontal base flow of carbon dioxide at 88.5 bar and 32.6 ^∘C is considered, which is subjected to a vertical heat flux of up to 12.0 kW/m² at Reynolds numbers of up to Re_Dh≤10.000. Here, optical- and surface temperature measurements are used concurrently to evaluate the flow. Integratced visualizations of the flow field show the onset of strong stratifications with limited heating rates in the near-pseudo-critical region. During unstable stratification, the channel flow is dominated by the upward motion of thermal plumes. When the stratification is stable, any vertical motion and turbulence present in an equivalent neutrally buoyant flow is suppressed. As a result, wall heat is removed more effectively in the unstably stratified configuration than in a forced convective flow, whereas the opposite is true for a stably stratified flow. The difference in the perceived heat transfer between the considered configurations increases as buoyancy becomes more dominant.

The Monin–Obukhov similarity theory (MOST) is a cornerstone of atmospheric science for describing turbulence in stable boundary layers. Extending MOST to stably stratified turbulent channel flows, however, is non-trivial due to confinement by solid walls. In this study, we investigate the applicability of MOST in closed channels and identify where and to what extent the theory remains valid. A key finding is that the ratio of the half-channel height to the Obukhov length serves as a governing parameter for identifying distinct flow regions and determining their corresponding mean velocity scaling. Hence, we propose a relation to estimate this ratio directly from the governing input parameters: the friction Reynolds and friction Richardson numbers (Re_τ and Ri_τ). The framework is tested against a series of direct numerical simulations across a range of Re_τ and Ri_τ. The reconstructed velocity profiles enable accurate prediction of the skin-friction coefficient crucial for quantifying pressure losses in stratified flows in engineering applications.

The objective of this work is to investigate the unexplored laminar-to-turbulent transition of a heated flat-plate boundary layer with a fluid at supercritical pressure. Two temperature ranges are considered: a subcritical case, where the fluid remains entirely in the liquid-like regime, and a transcritical case, where the pseudo-critical (Widom) line is crossed and pseudo-boiling occurs. Fully compressible direct numerical simulations are used to study (i) the linear and nonlinear instabilities, (ii) the breakdown to turbulence, and (iii) the fully developed turbulent boundary layer. In the transcritical regime, two-dimensional forcing generates not only a train of billow-like structures around the Widom line, resembling Kelvin–Helmholtz instability, but also near-wall travelling regions of flow reversal. These spanwise-oriented billows dominate the early nonlinear stage. When high-amplitude subharmonic three-dimensional forcing is applied, staggered Λ-vortices emerge more abruptly than in the subcritical case. However, unlike the classic H-type breakdown under zero pressure gradient observed in ideal-gas and subcritical regimes, the H-type breakdown is triggered by strong shear layers caused by flow reversals – similar to that observed in adverse pressure gradient boundary layers. Without oblique wave forcing, transition is only slightly delayed and follows a naturally selected fundamental breakdown (K-type) scenario. Hence in the transcritical regime, it is possible to trigger nonlinearities and achieve transition to turbulence relatively early using only a single two-dimensional wave that strongly amplifies background noise. In the fully turbulent region, we demonstrate that variable-property scaling accurately predicts turbulent skin-friction and heat-transfer coefficients.

This chapter focuses on modeling compressible wall-bounded turbulent flows. We first make the distinction between variable-property effects and intrinsic compressibility effects. Variable-property effects are associated with changes in density, viscosity, etc., primarily due to heat transfer, while intrinsic compressibility effects are associated with volume changes induced by pressure fluctuations. Both phenomena significantly impact turbulence and thus complicate turbulence modeling. The focus is on adapting Reynolds-averaged Navier-Stokes (RANS) models to account for these effects, which can also be applied to wall-modeled large eddy simulations (WMLES). Canonical flows are examined to illustrate the necessary modifications needed to incorporate compressibility effects into general RANS models. These improvements are applicable to many turbulence models and for complex flow geometries.

In this work we propose a novel method to ensure important entropy inequalities are satisfied semi-discretely when constructing reduced order models (ROMs) on nonlinear reduced manifolds. We are in particular interested in ROMs of systems of nonlinear hyperbolic conservation laws. The so-called entropy stability property endows the semi-discrete ROMs with physically admissible behaviour. The method generalizes earlier results on entropy-stable ROMs constructed on linear spaces. The ROM works by evaluating the projected system on a well-chosen approximation of the state that ensures entropy stability. To ensure accuracy of the ROM after this approximation we locally enrich the tangent space of the reduced manifold with important quantities. Using numerical experiments on some well-known equations (the inviscid Burgers equation, shallow water equations and compressible Euler equations) we show the improved structure-preserving properties of our ROM compared to standard approaches and that our approximations have minimal impact on the accuracy of the ROM. We additionally generalize the recently proposed polynomial reduced manifolds to rational polynomial manifolds and show that this leads to an increase in accuracy for our experiments.

We present a massively parallel GPU-accelerated solver for direct numerical simulations of transitional and turbulent flat-plate boundary layers and channel flows involving fluids in non-ideal thermodynamic states. While several high-fidelity solvers are currently available as open source, all of them are restricted to the ideal-gas region. In contrast, the CUBic Equation of state Navier-Stokes solver (CUBENS) can accurately model and simulate the non-ideal thermodynamics of single-phase compressible fluids in the vicinity of the vapor-liquid saturation line or the thermodynamic critical point. By employing high-order finite-difference schemes and convective terms in split, kinetic-energy-, and entropy-preserving form, the solver is numerically stable, and robust with minimal numerical dissipation, enabling it to capture the steep variations of non-ideal thermodynamic properties. For cost-effective high-fidelity simulations, in addition to MPI parallelization, CUBENS is GPU-accelerated using OpenACC directives for computation offloading, and asynchronous GPU-aware MPI for efficient GPU-GPU communication. Moreover, CUBENS is compatible with both NVIDIA and AMD GPU architectures, achieving significant performance results while ensuring energy-efficient simulations. For instance, using 64 NVIDIA A100 GPUs compared to 8192 CPUs at the same computational cost results in a speedup of approximately 130×. In multi-node and multi-GPU configurations ranging from 2 to 128 compute nodes (8 to 512 GPUs), a strong scaling efficiency of around 52% and a weak scaling efficiency of 0.88 with 1024³ points per GPU, corresponding to approximately 5 billion degrees of freedom, are achieved. The CUBENS solver is validated against selected cases from the literature, covering transitional to turbulent ideal and non-ideal flows up to the transonic regime. In particular, we demonstrate the solver's suitability and applicability for direct numerical simulations of transitional boundary layers with fluids at supercritical pressure and with buoyancy effects. The development of this high-fidelity solver offers the potential for future fundamental research in non-ideal compressible fluid dynamics. Program summary: Program Title: CUBic Equation of state Navier-Stokes (CUBENS) CPC Library link to program files: https://doi.org/10.17632/6jfy758gyv.1 Developer's repository link: https://github.com/pcboldini/CUBENS Licensing provisions: MIT Programming language: Fortran 90, OpenACC, MPI, Python, MATLAB Nature of problem: This code solves the three-dimensional Navier-Stokes equations for non-ideal gas flows in a Cartesian domain, applicable to boundary layers and channels. Solution method: This code uses high-order central finite-differences with split-convective form, preserving kinetic energy and entropy (KEEP) and pressure-equilibrium-preserving (PEP) property, for spatial discretization. The time advancement is performed with a third-order Total Variation Diminishing low-storage Runge-Kutta scheme. Flow non-ideality is accounted for by cubic equations of state and complex transport-properties models. Alongside MPI parallelization, the solver is GPU-accelerated using OpenACC for computation offloading and CPU-GPU data transfer, along with GPU-aware MPI for GPU-GPU communication.

Heating in industrial processes is responsible for approximately 13% of greenhouse gas emissions in Europe. Switching from fossil-fuel based boilers to heat pumps can help mitigate the effect of global warming. The present work proposes novel high-temperature transcritical heat pump cycles targeted at heating air with a mass flow rate of 10 kg/s up to 200 °C for spray drying processes. Four low-GWP refrigerants, R1233zd(E), R1336mzz(Z), n-Butane, and Ammonia are considered as the candidate working fluids. The pressure ratio of the compressor is optimized to achieve a maximum coefficient of performance (COP) for the four working fluids. A shell & tube heat exchanger is considered as the gas cooler. Using a generalized version of the ϵ-NTU method, the gas cooler is sized and a second law analysis is conducted. Striking a balance between the first- and second-law performance and size of the gas cooler, the R1233zd(E) transcritical heat pump cycle with a COP of 3.6 is judged to be the most promising option.

The competitiveness of hydrogen as a sustainable energy carrier depends greatly on its transportation and storage costs. Liquefying hydrogen offers advantages such as enhanced purity, versatility, and higher density, yet current industrial liquefaction processes face efficiency and cost challenges. Although various large-scale and efficient liquefaction concepts exist in the literature, they often overlook the economic and technical viability of such plants. Here, we addresses this issue by establishing a framework for modeling a large-scale hydrogen liquefaction concept and conducting both technical and economic assessments, with a specific focus on 125 tonnes per day (TPD) high-pressure hydrogen Claude-cycle concept. The technical analysis involves preliminary designs of key process components, while the economic assessment utilizes Aspen Process Economic Analyzer. Our findings indicate that at an electricity price of €0.1/kWh, the Claude-cycle liquefier concept yields a specific liquefaction cost (SLC) of €1.55/kg_LH2. A sensitivity analysis was performed, which shows that electricity price has a significant influence on the economics. Further investigation on the compressors design shows that incorporating high-speed centrifugal compressors could reduce the SLC by 5.42% and potentially more. Scaling up to 250 and 500 TPD reveals further cost improvements, while cost projections indicate substantial declines as the technology matures. Ultimately, this paper presents novel cost-scaling and experience curves of hydrogen liquefaction technology, demonstrating the compelling economic viability of integrating large-scale hydrogen liquefaction into sustainable energy infrastructure.

Experimental investigations into the characterization of vortices in hyperbolic funnels have shown efficient aeration properties. Certain regimes of vortices have been observed to exhibit high gas dissolution rates. This phenomenon has prompted inquiries into the underlying physical mechanisms at both micro and macroscopic scales. The present study employs computational fluid dynamics to numerically analyze the flow field organization inside these vortices, aiming to elucidate the observed high gas transfer rates. Transient simulations are performed on a three-dimensional radially structured hexahedral mesh, utilizing a multiphase Euler-Euler approach-based volume of fluid method for modeling, along with shear stress transport turbulence modeling based on k − ω equations with curvature correction. The evaluation of the two vortex regimes was conducted in terms of hydraulic retention time, water volume in the reactor, air-water interfacial area, and bulk mixing. Instabilities resembling Taylor vortices observed in Taylor-Couette flow systems emerge in the secondary flow field of these vortical structures, facilitating turbulent mixing. A qualitative analysis of the strength of these instabilities in terms of average vorticity per unit mass of water explains the high gas transfer efficiency. Despite high gas transfer rates, water exiting the funnel remains undersaturated under given operating conditions due to the short hydraulic retention time.

In the region close to the thermodynamic critical point and in the proximity of the pseudoboiling (Widom) line, strong property variations substantially alter the growth of modal instabilities, as revealed in Ren et al. [J. Fluid Mech. 871, 831 (2019)0022-112010.1017/jfm.2019.348]. Here, we study nonmodal disturbances in the spatial framework using an eigenvector decomposition of the linearized Navier-Stokes equations under the assumption of locally parallel flow. To account for nonideality, a new energy norm is derived. Several heat transfer scenarios at supercritical pressure are investigated, which are of practical relevance in technical applications. The boundary layers with the fluid at supercritical pressure are heated or cooled by prescribing the wall and free-stream temperatures so that the temperature profile is either entirely subcritical (liquidlike), supercritical (gaslike), or transcritical (across the Widom line). The free-stream Mach number is set to 10-3. In the nontranscritical regimes, the resulting streamwise-independent streaks originate from the lift-up effect. Wall cooling enhances the energy amplification for both subcritical and supercritical regimes. When the temperature profile is increased beyond the Widom line, a strong suboptimal growth is observed over very short streamwise distances due to the Orr mechanism. Due to the additional presence of transcritical Mode II, the optimal energy growth at large distances is found to arise from an interplay between lift-up and Orr mechanism. As a result, optimal disturbances are streamwise-modulated streaks with strong thermal components and with a propagation angle inversely proportional to the local Reynolds number. The nonmodal growth is put in perspective with modal growth by means of an N-factor comparison. In the nontranscritical regimes, modal stability dominates regardless of a wall-temperature variation. In contrast, in the transcritical regime, nonmodal N factors are found to resemble the imposition of an adverse pressure gradient in the ideal-gas regime. When cooling beyond the Widom line, optimal growth is greatly enhanced, yet strong inviscid instability prevails. When heating beyond the Widom line, optimal growth could be sufficiently large to favor transition, particularly with a high free-stream turbulence level.