Two-dimensional Navier-Stokes Research Articles

Simulation of Interfacial Waves in Core-Annular Pipe Flow with a Turbulent Annulus

Abstract Numerical simulations are conducted for the wave initiation, growth, and saturation at the oil-water interface in core-annular flow. The focus is on conditions with a turbulent water annulus, but the laminar water annulus is also considered. The simulation results are compared with lab measurements. The growth rate for the linear instability of different wave lengths in the case of a turbulent water annulus is obtained from two-dimensional axisymmetric Reynolds-Averaged Navier-Stokes simulations with the Launder-Sharma low-Reynolds number k-e model. The latter simulation results provide the most unstable wave length for the turbulent water annulus. Our study also shows the following. The maximum wave growth rate for a turbulent water annulus is significantly higher than for a laminar water annulus. The most unstable wave length in the simulations is about 25% smaller than in the experiments. The wave amplitude for the different wave lengths in the simulations is typically 17% lower than in the experiments.

Approximation of Two-Dimensional Time-Fractional Navier-Stokes Equations involving Atangana-Baleanu Derivative

This article addresses the two analytical methods, i.e., the new iterative transform method (NITM) and the homotopy perturbation transform method (HPTM), along with an Aboodh transform (AT), to approximate the nonlinear system of two-dimensional (2D) time-fractional Navier-Stokes (TFNS) equations. We take the time-fractional derivative in the form of Atangana-Baleanu (AB). The article's suggested examples examine the accuracy and efficacy of the proposed methods, while the graphs demonstrate their potential and effectiveness. The article also provides demonstrations of uniqueness and convergence. The aforementioned techniques are straightforward and support a high rate of convergence, which helps in understanding the dynamics of fractional nonlinear systems.

Reconstruction of dynamic wind turbine wake flow fields from virtual Lidar measurements via physics-informed neural networks

Accurate characterisation of wind turbine wakes is important for the optimal design and operation of wind farms. However, current techniques for full-scale wind measurements are still limited to point characterisation. To address the research challenge in obtaining field characterisation of real-world wind turbine wakes, this work investigates the reconstruction of the dynamic wake flow fields based on a virtual turbine-mounted Lidar and physics-informed neural networks. Specifically, the wake flow field is reconstructed by fusing the sparse measurements with the two-dimensional Navier-Stokes equations without imposing any models for the unsteady wake. Different from supervised machine learning approaches which need the measured values for the quantities of interest in the first place, the proposed method can achieve the prediction of the wind velocity at new locations where there is no measurement available. The reconstruction performance is evaluated via high-fidelity numerical experiments and it is shown that the dynamic wind turbine wake flow fields are predicted accurately, where the main wake features, including the downwind development and crosswind meandering of the wake, are both captured. This work thus paves the way for investigating full-scale in situ wake flow dynamics in real-world wind energy sites.

A MEEVC discretization for two-dimensional incompressible Navier-Stokes equations with general boundary conditions

In this work, we present a mass, energy, enstrophy and vorticity conserving (MEEVC) mixed finite element discretization for two-dimensional incompressible Navier-Stokes equations as an alternative to the original MEEVC scheme proposed in A. Palha and M. Gerritsma (2017) [5]. The present method can incorporate general boundary conditions. Conservation properties are proven. Supportive numerical experiments with both exact and inexact quadrature are provided.

Design, analysis, and testing of a prototype-scale latent heat thermal energy storage (LTES) system

A parametric analysis was performed to design a prototype-scale latent heat thermal energy storage (LTES) system using commercial grade hexahydrate calcium chloride as phase change material (PCM) in a staggered tube array configuration, placed horizontally. The study involved numerical analysis of melting and solidification processes of the PCM within the tube array, considering transient two-dimensional Navier-Stokes equations and a Realizable k-ɛ turbulence model to predict fluid flow and heat transfer. The enthalpy-porosity technique was used to model PCM melting and solidification. The study included a comparison between numerically predicted and experimentally obtained air temperatures at the tube array outlet, as well as PCM temperatures. This comparison revealed an excellent level of agreement. Additionally, the study compared the pressure drop within the array and the average peripheral heat transfer coefficient calculated from the numerical results to the well-established Zukauskas correlation. The resulting design successfully met the predefined performance criteria: achieving a cooling effect of 4 °C for a four-hour duration while maintaining a pressure drop of <100 Pa. The proposed prototype-scale tube array design can be efficiently cooled down and PCM frozen overnight. The energy storage density of the system falls within the range of 22 to 27 kWh/m3.

On the solvability of a boundary value problem for a two-dimensional system of Navier-Stokes equations in a cone

Due to the fact that the Navier-Stokes equations are involved in the formulation of a large number of interesting problems that are important from an applied point of view, these equations have been the object of attention of mechanics, mathematicians and other scientists for several decades in a row. But despite this, many problems for the Navier-Stokes equation remain unexplored to this day. In this work, we are exploring the solvability of a boundary value problem for a two-dimensional Navier-Stokes system in a non-cylindrical degenerating domain, namely, in a cone with its vertex at the origin. Previously, we studied cases of the linearized Navier-Stokes system or non-degenerating cylindrical domains, so this work is a logical continuation of our previous research in this direction. To the above-mentioned degenerate domain we associate a family of non-degenerate truncated cones, which, in turn, are formed by a oneto-one transformation into cylindrical domains, where for the problem under consideration we established uniform a priori estimates with respect to changes in the index of the domains. Further, using a priori estimates and the Faedo-Galerkin method, we established the existence, uniqueness of solution in Sobolev classes, and its regularity as the smoothness of the given functions increases.

Energy and exergy analysis of a laboratory-scale latent heat thermal energy storage (LTES) using salt-hydrate in a staggered tube arrangement

The energy and exergy analyses were performed for a laboratory-scale latent heat thermal energy storage (LTES) using hexahydrate calcium chloride (CC6) as phase change material (PCM) in a staggered tube array configuration, placed horizontally. The PCM melting and solidification process within the tube array was investigated by performing a numerical analysis using transient two-dimensional Navier-Stokes equations and Realizable k-ɛ turbulence model to predict flow and heat transfer. The enthalpy-porosity technique was applied to model PCM melting and solidification. The accuracy of the numerical model was validated against experimentally obtained data where the numerically predicted and measured air temperature at the tube array outlet and PCM temperatures were compared. Additionally, the pressure drop in the array and the average peripheral heat transfer coefficient calculated from the numerical results were compared to the well-known Zukauskas correlation (Zu). The results show that the melting and solidification process of PCM in the tubes of the array is asymmetric in nature, with the PCM melting taking slightly longer compared to the solidification process. It was also observed that the energy quantity was higher compared to the exergy quantity, as exergy considers entropy generation and accounts for irreversibility within the system.

A comparative study of latent heat thermal energy storage (LTES) system using cylindrical and elliptical tubes in a staggered tube arrangement

Numerical analysis was performed to compare the thermal and hydraulic performance of the elliptical and circular tube geometries in the prototype-scale latent heat thermal energy storage (LTES) system. A staggered tube array configuration was used for the analysis where the tubes were placed horizontally. Commercial grade hexahydrate calcium chloride (CC6) was selected as a phase change material (PCM) due to its phase change occurring at room temperature. The study included a numerical analysis of the melting and solidification processes of the PCM within the tube array for both tube shapes, employing transient two-dimensional Navier-Stokes equations and a Realizable k-ɛ turbulence model to predict fluid flow and heat transfer. The enthalpy-porosity technique was used to model PCM melting and solidification. Numerical predictions indicate that the thermal performance for both tube shapes is almost indistinguishable during the melting and freezing processes. However, due to the aerodynamic shape of the elliptical tube, the pressure drop across the horizontal elliptical tube array is approximately 80% lower compared to that of the circular geometry. This reduced pressure drop implies that less pumping power will be required to achieve the desired flow rate across the tube bank, leading to economic advantages for the horizontal elliptical tube array. The proposed design can be used as a dry cooling technique employing an Air Cooled Condenser (ACC) in a thermal power plant.

Convective thermal losses of long-term underground hot water storage

A case of underground long-term hot water storage is investigated numerically. The study is based on the unsteady two-dimensional Navier-Stokes equations in Boussinesq approximation applied to a closed cavern with time-dependent temperature boundary conditions on the walls. The problem formulated in a vorticity-stream function statement is solved by finite difference method (FDM) for high values of the Rayleigh number and for the Prandtl number of water. Streamlines, velocity and temperature fields are presented graphically for given moments of time. The evolution of the thermocline thickness in the mid-section of the cavern is slow and illustrates that the hot water zone occupies more than the half of the cavern even after 6 months period. The Nusselt number on the walls shows that the convective thermal losses are small and after certain period of time tend to decrease due to the diminished temperature difference at the walls. The influence of the fluid convection on the thermal losses is evaluated quantitatively, showing high seasonal thermal efficiency of the insulated hot water storage.

A viscous investigation on the hydrodynamic coefficients and wave loads under the interaction of side-by-side cylinders in regular waves

The paper introduces a two-dimensional Reynolds-Averaged Navier-Stokes Computational Fluid Dynamics (RANS CFD) model to investigate the effects of adjacently floating bodies on the hydrodynamic coefficients and wave loads acting on these bodies. The analysis considers two square cylinders with a narrow gap in way of the free water surface. The physical significance of subsections of the hydrodynamic coefficient matrix becomes evident when comparing coefficients obtained from oscillating-oscillating versus oscillating-fixed models. This understanding is significant, particularly in applications where only specific portions of the hydrodynamic coefficient matrix are utilized, such as in ship-ice collision scenarios (Jiang et al., 2023a). The correlation between wave loads and wave elevations reveals that the sway force and roll moment are associated with the disparity in wave elevation at both sides of each cylinder, while heave forces correlate with the averaged wave elevation on both sides of each cylinder. It is concluded that hydrodynamic interactions are sensitive to the displacement of a floater and accordingly can dominate added mass and damping effects.

Stationary and Non-Stationary Periodic Flows Mathematical Modelling using Various Vortex Viscosity Models

Introduction. Mathematical modelling of currents is an urgent research topic in the field of hydrodynamics and oceanography. Despite ongoing research in the field of developing accurate and efficient numerical methods for solving Navier-Stokes equations that take into account vortex viscosity, the problems of accurate prediction and control of turbulence remain unresolved. The influence of nonlinear effects in vortex viscosity models on the accuracy of forecasts and their applicability to various flow conditions also remains relevant. The aim of the study is to study the influence of linearized and quadratic bottom friction and two turbulence models on the numerical solution of stationary and non-stationary periodic flows. Special emphasis is placed on comparing numerical results with analytical solutions within the framework of using various models of bottom friction.Materials and Methods. The computational models used in this study are based on a simplified two-dimensional wave model and full three-dimensional Navier-Stokes equations. The classical model of shallow water motion and the 2D model without taking into account dynamic changes in the geometry of the reservoir surface are derived from a system of equations for a spatially inhomogeneous three-dimensional mathematical model of wave hydrodynamics of a shallow reservoir. Analytical solutions were found by linearization of the equations, which obviously has its limitations. A distinction is made between two types of nonlinear effects – nonlinearities caused by higher-order terms in the equations of motion, i. e. terms of advective acceleration and friction, and nonlinear effects caused by geometric nonlinearities, this is due, for example, to different water depths and reservoir widths, which will be important when modelling a real sea.Results. The results of modeling stationary and non-stationary periodic flows in a schematized rectangular basin using linearized bottom friction are presented. The influence of linearization on the numerical solution is investigated in comparison with analytical profiles using models calculating bottom friction in a quadratic formulation. In combination with quadratic bottom friction, two turbulence models are studied: the constant vortex viscosity and the Prandtl mixing length model. The results obtained as a result of three-dimensional modelling are compared with the results of two-dimensional modeling and analytical solutions averaged in depth.Discussion and Conclusion. New approaches to modelling and studying flows with variable vortex viscosity are proposed, including analysis of the influence of linearization and the use of various turbulence models. For the linearized and quadratic formulations of bottom friction, it is proved that the numerical results for the case of stationary flow show great similarity with analytical solutions, since the surface height is much less than the water depth and advection can be neglected. The numerical results for the unsteady flow also show a good agreement with the theory. Unlike analytical solutions, numerical modelling has minor deviations in the long run. The study of flows, within the framework of using various turbulence models, will make it possible to take into account the influence of nonlinear effects in vortex viscosity models on the accuracy of forecasts and their applicability to various flow conditions. The results obtained make it possible to better understand and describe the physical processes occurring in shallow waters. This opens up new possibilities for applying mathematical modelling to predict and analyze the impact of human activities on the marine environment and to solve other problems in the field of oceanology and geophysics.

Direct numerical simulations of long-range infrasound propagation: Implications for source spectra estimation.

The evolution of observed dominant frequencies from a high-intensity infrasonic pulse with receiver range and stratospheric temperature is investigated using direct numerical simulations of the two-dimensional unsteady compressible Navier-Stokes equations. There is a high level of uncertainty in estimating source dominant frequencies based on received signals at sparse points on the ground. Nonlinear propagation effects in the ground-level thermospheric arrivals are found to significantly alter dominant frequency measurements compared to stratospheric arrivals with smaller amplitude sources. With a larger amplitude source, variations in observations are minimized as a result of nonlinear effects being ubiquitous across all atmospheric components of received signals but have a greater offset to the source dominant frequency. An approach to determine the source dominant frequency and minimize atmospheric variability is presented by calculating a source-to-receiver spectral transfer function averaged across the atmospheric states. This method reduces atmospheric variability in source frequency estimates within the pseudo-linear propagation regime and the average error to the known source frequency with a large amplitude source. The reduction of errors in source frequency estimates demonstrates the feasibility of using remote infrasound measurements as an indicator of source frequency and, in turn, the explosive yield of clandestine nuclear weapon test explosions.

Effect of confinement on the propagation patterns of lean hydrogen–air flames

The effect of confinement on the propagation of lean hydrogen–air flames in a Hele-Shaw cell, formed by two parallel plates separated by a small distance, is studied numerically. The influence of momentum loss on the flame structure and propagation rates is analyzed by examining two limits: (i) the limit of very small distance between the plates, for which a quasi-three-dimensional (quasi-3D) Darcy’s law-based approximation is used, and (ii) the limit of unconfined geometry, for which the full set of two-dimensional (2D) Navier-Stokes equations is considered. The influence of heat losses is also included. The chemistry is described by a one-step reduced kinetics with a simple model for transport properties. The results demonstrate that momentum loss primarily increases the flame surface area by elongating the cellular structures in a finger-like form, with a relatively small enhancement of the reactivity at the flame front, consistent with a mechanism of hydrodynamic nature. This elongation leads to approximately a 50% increase in flame speed compared to the absence of confinement. Conversely, heat losses are observed to flatten the large-scale continuous flame front and, if significant enough, can lead to the formation of isolated propagating flame cells of the two-headed or, ultimately, one-headed type. The present results are in good agreement with recent experiments.

Numerical comparison of two-dimensional Navier-Stokes flows on the whole plane and the periodic domain

Numerical comparison of two-dimensional Navier-Stokes flows on the whole plane and the periodic domain

Effect of discrete heating on the key parameters and entropy generation in a corrugated enclosure filled with hybrid nanofluid

This research work aims to study the effect of discrete heating on natural convective heat transfer, fluid flow, and entropy generation within a corrugated enclosure filled with water-based molybdenum disulfide (MoS2)–silicon dioxide (SiO2) hybrid nanofluid. In the study, various heights of the discrete heaters in terms of aspect ratios ( AR = 0.2 , 0.4 , 0.6 , 0.8 ) are considered for various Rayleigh numbers ( Ra = 10 3 − 10 6 ) with varying nanoparticle volume fractions ( ϕ = 0 % − 4 % ). Stream function-vorticity formulations of two-dimensional (2D) Navier-Stokes equations and energy equation are used as the governing equations, and a higher-order compact (HOC) finite difference scheme is employed to discretize those equations. The results obtained from the current developed code are first validated with the existing numerical and experimental data. Finally, the present computed results are studied, both qualitatively and quantitatively, in terms of streamlines, isotherms, entropy generations, average Nusselt numbers, and Bejan numbers. At Ra = 10 6 and the aspect ratio of heater length AR = 0.2 , a different and interesting flow and heat transfer phenomenon is observed. The results show that an increase in the aspect ratio of heater length from 0.2 to 0.8 decreases the average Nusselt number and total entropy generation by 57.49% and 13.76%, respectively, for Ra = 10 6 and ϕ = 4 % . In addition, the effect of nanoparticle concentration on total entropy generation increases up to 73.40% when the aspect ratio of heater length increases to 0.8.

Stability of planar exterior stationary flows with suction

We consider the two-dimensional Navier-Stokes system in a domain exterior to a disk. The system admits a stationary solution with critical decay O(|x|−1) written as a linear combination of the pure rotating flow and the flux carrier. We prove its nonlinear stability in large time for initial disturbances in L2 under smallness conditions, assuming that there is suction across the boundary, namely that the sign of coefficients of the flux carrier is negative. This result partially solves an open problem in the literature.

Numerical simulations of two-dimensional incompressible Navier-Stokes equations by the backward substitution projection method

The backward substitution method is a newly developed meshless method that has been used for the simulation of many problems in science and engineering with high accuracy and efficiency. In this paper, we explore the feasibility of employing the backward substitution method for simulating two-dimensional incompressible flows. Two non-increment pressure correction projection methods are considered to decompose the original velocity and pressure coupling system into two boundary value problems of intermediate velocity and pressure. Then, two boundary value problems are solved by the backward substitution method in each time iteration step. Five numerical examples are provided to demonstrate the accuracy, computational efficiency and convergence of the method. Comparisons with some existing meshless methods verify the method's advantages and potential applications to engineering.

Non-orthogonal multiple-relaxation-time lattice Boltzmann method for vorticity-streamfunction formulation

In this work, we propose a non-orthogonal multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) for the two-dimensional (2D) incompressible Navier-Stokes equations in vorticity-stream function formulation. In this method, the vorticity and the streamfunction equations are calculated by two separate LB equations, which are constructed by using the non-orthogonal transformation matrix. Most importantly, to eliminate the residual term that may appear in the recovered vorticity equation, the convection term is handled as a source term in the LB equation rather than being directly recovered by the evolution equation. In addition, due to the convection term in our method is computed with a local scheme, it retains the parallel feature of the standard LBM. We then test the accuracy and the stability of the present LBM by considering various benchmark examples. The results show that the present method has a second-order convergence rate in space. Our numerical results indicate that the present LBM is a good candidate for solving streamfunction-vorticity formulation.

Numerical simulation of acoustic streaming in standing waves

This study presents a numerical investigation of acoustic streaming motion (of the Rayleigh type) in a compressible gas inside two-dimensional rectangular enclosures. To numerically study the effects of the sound field intensity on the formation process of streaming structures, we propose to discretize the fully compressible form of the two-dimensional Navier-Stokes equations using a high-order (formally greater or equal to fourth-order) accurate numerical scheme in both space and time. The proposed numerical solver utilizes high-order compact schemes along each spatial dimension combined with a filtering procedure when it is necessary. Acoustic standing waves are excited inside the enclosures and the resulting acoustic streaming patterns are investigated for low and high-intensity waves, in both linear and nonlinear regimes following closely the work of Daru et al. (2013). An extended investigation indicates that without the incorporation of an appropriate filter, the application of the high-order compact schemes in case of fast streaming results in spurious oscillations which inhibit their applicability. Following the recent relevant literature, the numerical simulations performed demonstrate the ability of the proposed numerical approach to reproduce efficiently and robustly the transitions from regular acoustic streaming to irregular streaming and the relevant phenomena confirming results that have been previously presented.

Large deviations for the two-time-scale stochastic convective Brinkman-Forchheimer equations

The convective Brinkman-Forchheimer (CBF) equations are employed to characterize the motion of incompressible fluid in a saturated porous medium. This work investigates the small noise asymptotic of two-time-scale stochastic CBF equations in two and three dimensional bounded domains. More precisely, we establish a Wentzell-Freidlin type large deviation principle for stochastic partial differential equations that have slow and fast time-scales. The slow component is the stochastic CBF equations in two or three dimensions perturbed by a small multiplicative Gaussian noise, while the fast component is a stochastic reaction-diffusion equation with damping. The results are obtained by using a variational method (based on weak convergence approach) developed by Budhiraja and Dupuis, Khasminkii's time discretization approach and stopping time arguments. In particular, the findings from this study are also applicable to two-dimensional stochastic Navier-Stokes equations.