Theoretical and Computational Fluid Dynamics最新文献

Resolvent analysis is a powerful tool that can reveal the linear amplification mechanisms between the forcing inputs and the response outputs about a base flow. These mechanisms can be revealed in terms of a pair of forcing and response modes and the associated energy gains (amplification magnitude) at a given frequency. The linear relationship that ties the forcing and the response is represented through the resolvent operator (transfer function), which is constructed through spatially discretizing the linearized Navier–Stokes operator. One of the unique strengths of resolvent analysis is its ability to analyze statistically stationary turbulent flows. In light of the increasing interest in using resolvent analysis to study a variety of flows, we offer this guide in hopes of removing the hurdle for students and researchers to initiate the development of a resolvent analysis code and its applications to their problems of interest. To achieve this goal, we discuss various aspects of resolvent analysis and its role in identifying dominant flow structures about the base flow. The discussion in this paper revolves around the compressible Navier–Stokes equations in the most general manner. We cover essential considerations ranging from selecting the base flow and appropriate energy norms to the intricacies of constructing the linear operator and performing eigenvalue and singular value decompositions. Throughout the paper, we offer details and know-how that may not be available to readers in a collective manner elsewhere. Towards the end of this paper, examples are offered to demonstrate the practical applicability of resolvent analysis, aiming to guide readers through its implementation and inspire further extensions. We invite readers to consider resolvent analysis as a companion for their research endeavors.

We present LinStab2D, an easy-to-use linear stability analysis MATLAB tool capable of handling complex domains, performing temporal and spatial linear stability, and resolvent analysis. We present the theoretical foundations of the code, including the linear stability and resolvent analysis frameworks, finite differences discretization schemes, and the Floquet ansatz. These concepts are explored in five different examples, highlighting and illustrating the different code capabilities, including mesh masking, mapping, imposition of boundary constraints, and the analysis of periodic flows using Cartesian or axisymmetric coordinates. These examples were constructed to be a departure point for studying other flows.

Cluster and void formations are key processes in the dynamics of particle-laden turbulence. In this work, we assess the performance of various neural network models for synthesizing preferential concentration fields of particles in turbulence. A database of direct numerical simulations of homogeneous isotropic two-dimensional turbulence with one-way coupled inertial point particles, is used to train the models using vorticity as the input to predict the particle number density fields. We compare encoder–decoder, U-Net, generative adversarial network (GAN), and diffusion model approaches, and assess the statistical properties of the generated particle number density fields. We find that the GANs are superior in predicting clusters and voids, and therefore result in the best performance. Additionally, we explore a concept of “supersampling”, where neural networks can be used to predict full particle data using only the information of few particles, which yields promising perspectives for reducing the computational cost of expensive DNS computations by avoiding the tracking of millions of particles. We also explore the inverse problem of synthesizing the absolute values of the vorticity fields using the particle number density distribution as the input at different Stokes numbers. Hence, our study also indicates the potential use of neural networks to predict turbulent flow statistics using experimental measurements of inertial particles.

The use of heaving and pitching fins for underwater propulsion of engineering devices poses an attractive outlook given the efficiency and adaptability of natural fish. However, significant knowledge gaps need to be bridged before biologically inspired propulsion is able to operate at competitive performances in a practical setting. One of these relates to the design of structures that can leverage passive deformation and active morphing in order to achieve optimal hydrodynamic performance. To provide insights into the performance improvements associated with passive and active fin deformations, we provide here a systematic numerical investigation in the thrust, power, and efficiency of 2D heaving and pitching fins with imposed curvature variations. The results show that for a given chordline kinematics, the use of curvature can improve thrust by 70% or efficiency by 35% over a rigid fin. Maximum thrust is achieved when the camber variations are synchronized with the maximum heave velocity, increasing the overall magnitude of the force vector while increasing efficiency as well. Maximum efficiency is achieved when camber is applied during the first half of the stroke, tilting the force vector to create thrust earlier in the cycle than a comparable rigid fin. Overall, our results demonstrate that curving fins are consistently able to significantly outperform rigid fins with the same chord line kinematics on both thrust and hydrodynamic efficiency.

The generation mechanism of wall heat flux is one of the fundamental problems in supersonic/hypersonic turbulent boundary layers. A novel heat decomposition formula under the curvilinear coordinate was proposed in this paper. The new formula has wider application scope and can be applied in the configurations with grid deformed. The new formula analyzes the wall heat flux of an interaction between a shock wave and a turbulent boundary layer over a compression corner. The results indicated good performance of the formula in the complex interaction region. The contributions of different energy transport processes were obtained. While the processes by the mean profiles such as molecular stresses and heat conduction, can be ignored, the contributions by the turbulent fluctuations, such as Reynolds stresses and turbulent transfer of heat flux, were greatly increased. Additionally, the pressure work is another factor that affects the wall heat flux. The pressure work in the wall-normal direction is concentrated close to the reattachment point, while the pressure work in the streamwise direction acts primarily in the shear layer and the reattachment point.

abstract

A previously developed numerical-multilayer modeling approach for systems of governing equations is extended so that unwanted terms, resulting from vertical variations in certain background parameters, can be removed from the dispersion-relation polynomial associated with the system. The new approach is applied to linearized anelastic and compressible systems of governing equations for gravity waves including molecular viscosity and thermal diffusion. The ability to remove unwanted terms from the dispersion-relation polynomial is crucial for solving the governing equations when realistic background parameters, such as horizontal velocity and temperature, with strong vertical gradients, are included. With the unwanted terms removed, previously studied dispersion-relation polynomials, for which methods for defining upgoing and downgoing vertical wavenumber roots already exist, are obtained. The new methods are applied to a comprehensive set of medium-scale time-wavepacket examples, with realistic background parameters, lower boundary conditions at 30 km altitude, and modeled wavefields extending up to 500 km altitude. Results from the compressible and anelastic model versions are compared, with compressible governing-equation solutions understood as the more physically accurate of the two. The new methods provide significantly less computationally expensive alternatives to nonlinear time-step methods, which makes them useful for comprehensive studies of the behavior of viscous/diffusive gravity waves and also for large studies of cases based on observational data. Additionally, they generalize previously existing Fourier methods that have been applied to inviscid problems while providing a theoretical framework for the study of viscous/diffusive gravity waves.

Graphic abstract

In this paper we present a numerical scheme based on spectral collocation methods to investigate the flow of a piezo-viscous fluid, i.e., a fluid in which the rheological parameters depend on the pressure. In particular, we consider an incompressible Navier–Stokes fluid with pressure dependent viscosity flowing in: (i) a two-dimensional non-symmetric planar channel; (ii) a three-dimensional axisymmetric non-straight conduit. For both cases we impose the Navier slip boundary conditions that can be reduced to the classical no-slip condition for a proper choice of the slip parameter. We assume that the dependence of the viscosity on the pressure is of exponential type (Barus law), even though the model can be replaced by any other viscosity function. We write the mathematical problem (stress based formulation) and discretize the governing equations through a spectral collocation scheme. The advantage of this numerical procedure, which to the authors’ knowledge has never been used before for this class of fluids, lies in in the ease of implementation and in the accuracy of the solution. To validate our model we compare the numerical solution with the one that can be obtained in the case of small aspect ratio, i.e., the leading order lubrication solution. We perform some numerical simulation to investigate the effects of the pressure-dependent viscosity on the flow. We consider different wall functions to gain insight also on the role played by the channel/duct geometry. In both cases (i), (ii) we find that the increase of the coefficient appearing in the viscosity function results in a global reduction of the flow, as physically expected.

The effects of inlet Mach number on the unsteadiness of shock-boundary layer interactions (SBLIs) over curved surfaces are investigated for a supersonic turbine cascade using wall-resolved large eddy simulations. Three inlet Mach numbers, 1.85, 2.00, and 2.15 are considered at a chord-based Reynolds number 395,000. The curved walls of the airfoils impact the SBLIs due to the state of the incoming boundary layers and local pressure gradients. On the suction side, due to the convex wall, the boundary layer entering the SBLI evolves under a favorable pressure gradient and bulk dilatation. On the other hand, the concave wall on the pressure side imposes an adverse pressure gradient and bulk compression. Variations in the inlet Mach number induce different shock impingement locations, enhancing these effects. A detailed characterization of the suction side boundary layers indicates that a higher Mach number leads to larger shape factors, favoring separation and larger bubbles, while the reverse holds for the pressure side. A time-frequency analysis reveals the presence of intermittent events in the separated flow occurring predominantly at low-frequencies on the suction side and at mid-frequencies on the pressure side. Increasing the inlet Mach number leads to an increase in the time scales of the intermittent events on the suction side, which are associated with instants when high-speed streaks penetrate the bubble, causing local flow reattachment and bubble contractions. Instantaneous flow visualizations show the presence of streamwise vortices developing on the turbulent boundary layers on both airfoil sides and along the bubbles. These vortices influence the formation of the large-scale longitudinal structures in the boundary layers, affecting the mass imbalance inside the separation bubbles.

We conduct a comprehensive analysis of two data assimilation methods: the first utilizes the discrete adjoint approach with a correction applied to the production term of the turbulence transport equation, preserving the Boussinesq approximation. The second is a state observer method that implements a correction in the momentum equations alongside a turbulence model, both applied to fluid dynamics simulations. We investigate the impact of varying computational mesh resolutions and experimental data resolutions on the performance of these methods within the context of a periodic hill test case. Our findings reveal the distinct strengths and limitations of both methods, which successfully assimilate data to improve the accuracy of a RANS simulation. The performance of the variational model correction method is independent of input data and computational mesh resolutions. The state observer method, on the other hand, is sensitive to the resolution of the input data and CFD mesh.

Two-dimensional free-surface flow past a submerged rectangular disturbance in an open channel is considered. The forced Korteweg–de Vries model of Binder et al. (Theor Comput Fluid Dyn 20:125–144, 2006) is modified to examine the effect of varying obstacle length and height on the response of the free-surface. For a given obstacle height and flow rate in the subcritical flow regime an analysis of the steady solutions in the phase plane of the problem determines a countably infinite set of discrete obstacle lengths for which there are no waves downstream of the obstacle. A rich structure of nonlinear behaviour is also found as the height of the obstacle approaches critical values in the steady problem. The stability of the steady solutions is investigated numerically in the time-dependent problem with a pseudospectral method.