Theoretical and Computational Fluid Dynamics最新文献

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.

In open flow simulations, the dispersion characteristics of disturbances near synthetic boundaries can lead to unphysical boundary scattering interactions that contaminate the resolved flow upstream by propagating numerical artifacts back into the domain interior. This issue is exacerbated in flows influenced by real or apparent body forces, which can significantly disrupt the normal stress balance along outflow boundaries and generate spurious pressure disturbances. To address this problem, this paper develops a zero-parameter, physics-based outflow boundary condition (BC) designed to minimize pressure scattering from body forces and pseudo-forces and enhance transparency of the artificial boundary. This “balanced outflow BC” is then compared against other common BCs from the literature using example axisymmetric and three-dimensional open swirling flow computations. Due to centrifugal and Coriolis forces, swirling flows are known to be particularly challenging to simulate in open geometries, as these apparent forces induce non-trivial hydrostatic stress distributions along artificial boundaries that cause scattering issues. In this context, the balanced outflow BC is shown to correspond to a geostrophic hydrostatic stress correction that balances the induced pressure gradients. Unlike the alternatives, the balanced outflow BC yields accurate results in truncated domains for both linear and nonlinear computations without requiring assumptions about wave characteristics along the boundary.