Theoretical and Computational Fluid Dynamics最新文献

Shear-induced droplet formation is important in many industrial applications, primarily focusing on droplet sizes and pinch-off frequency. We propose a one-dimensional mathematical model that describes the effect of shear forces on the droplet interface evolution. The aim of this paper is to simulate paraffin wax droplets in a co-flowing fluid using the proposed model to estimate the droplet volume rate for different flow velocities. Thus, the study focuses only on the dripping regime. This one-dimensional model has a single parameter that arises from the force balance on the interface. This parameter is related to the shear layer thickness and hence influenced by the change in quantities like velocity, viscosity, and surface tension. The correlation describing the dependence of the parameter on these quantities using non-dimensional numbers is presented. The model is then cross-validated with the previous computational and experimental data. We use PETSc, an open-source solver toolkit, to implement our model using a mixed finite element discretization. We present the simulation results for liquid paraffin wax under fast-moving airflow with a range of velocities.

Employing direct numerical simulations, we investigate water and water-glycerol (85 wt%) droplets ((sim )25 µL) moving on smooth surfaces, with contact angles of around 90(^{circ }), at varying inclinations. Our focus is on elucidating the relative contribution of local viscous forces in the wedge and bulk regions in droplets to the total viscous force. We observe that, for fast-moving droplets, both regions contribute comparably, while the contribution of the wedge region dominates in slow-moving cases. Comparisons with existing estimates reveal the inadequacy of previous predictions in capturing the contributions of wedge and bulk viscous forces in fast-moving droplets. Furthermore, we demonstrate that droplets with identical velocities can exhibit disparate viscous forces due to variations in internal fluid dynamics.

We study generalised quasilinear (GQL) approximations applied to turbulent plane Couette flow. The GQL framework is explored in conjunction with a Galerkin reduced-order model (ROM) recently developed by Cavalieri and Nogueira (Phys Rev Fluids 7:102601, 2022), which considers controllability modes of the linearised Navier–Stokes system as basis functions, representing coherent structures in the flow. The velocity field is decomposed into two groups: one composed by high-controllability modes and the other by low-controllability modes. The former group is solved with the full nonlinear equations, whereas the equations for the latter are linearised. We also consider a new GQL framework wherein the linearised equations for the low-controllability modes are driven by nonlinear interactions of modes in the first group, which are characterised by large-scale coherent structures. It is shown that GQL-ROMs successfully recover the statistics of the full model with relatively high controllability thresholds and sparser nonlinear operators. Driven GQL-ROMs were found to converge more rapidly than standard GQL approximations, providing accurate description of the statistics with a larger number of linearised modes. This indicates that the forcing of linearised flow structures by large-scale coherent structures is an important feature of turbulence dynamics that should be considered in GQL models. The results presented here reveal that further model reductions are attainable with GQL-ROMs, which can be valuable to extend these models to larger Reynolds numbers.

This study presents a physics-based, low-order model for the trailing edge (TE) noise generated by an airfoil at low angle of attack. The approach employs incompressible resolvent analysis of the mean flow to extract relevant spanwise-coherent structures in the transitional boundary layer and near wake. These structures are integrated into Curle’s solution to Lighthill’s acoustic analogy to obtain the scattered acoustic field. The model has the advantage of predicting surface pressure fluctuations from first principles, avoiding reliance on empirical models, but with a free amplitude set by simulation data. The model is evaluated for the transitional flow ((text {Re} = 5e4)) around a NACA0012 airfoil at 3 deg angle of attack, which features TE noise with multiple tones. The mean flow is obtained from a compressible large eddy simulation, and spectral proper orthogonal decomposition (SPOD) is employed to extract the main hydrodynamic and acoustic features of the flow. Comparisons between resolvent and SPOD demonstrate that the physics-based model accurately captures the leading coherent structures at the main tones’ frequencies, resulting in a good agreement of the reconstructed acoustic power with that of the SPOD (within 4 dB). Discrepancies are observed at high frequencies, likely linked to nonlinearities that are not considered in the resolvent analysis. The model’s directivity aligns well with the data at low Helmholtz numbers, but it fails at high frequencies where the back-scattered pressure plays a significant role in directivity. This modeling approach opens the way for efficient optimization of airfoil shapes in combination with low-fidelity mean flow solvers to reduce TE noise.

Modal decomposition techniques are important tools for the analysis of unsteady flows and, in order to provide meaningful insights with respect to coherent structures and their characteristic frequencies, the modes must possess a robust spatial support. In this context, although widely used, methods based on singular value decomposition (SVD) may produce modes that are difficult to interpret when applied to problems dominated by intermittent and transient events. Fortunately, specific modal decomposition techniques have been recently developed to analyze such problems, but a proper comparison between them is still lacking from the literature. Therefore, this work compares two recent methods: the fast adaptive multivariate empirical mode decomposition (FA-MVEMD) and the multiresolution dynamic mode decomposition (mrDMD). These techniques are employed here for the study of flow databases involving transient and intermittent dynamics. Specifically, the investigated problems include an SD7003 airfoil subjected to deep dynamic stall conditions, and a steady NACA0012 airfoil operating at a transitional Reynolds number. In the former case, the methods are employed to investigate the onset and evolution of the dynamic stall vortex (DSV), while in the latter case, intermittent vortex pairing is analyzed. We show that the combination of a multidimensional EMD with the Hilbert transform provides modes with superior spatial support when compared to the mrDMD, also allowing the characterization of instantaneous frequencies of coherent structures. Moreover, the EMD also condenses a larger amount of information within a single intrinsic mode function (IMF).

The stability of the flow past a circular cylinder in the presence of a wavy ground is investigated numerically in this paper. The wavy ground consists of two complete waves with a wavelength of 4D and an amplitude of 0.5D, where D is the cylinder diameter. The vertical distance between the cylinder and the ground is varied, and four different cases are considered. The stability analysis shows that the critical Reynolds number increases for cases close to the ground when compared to the flow past a cylinder away from the ground. The maximum critical Reynolds number is obtained when the cylinder is located in front of the waves. The wavy ground adds layers of clockwise (negative) vorticity due to flow separation from the wave peak, to the oscillating Kármán vortex. This negative vorticity from the wave peak also cancels part of the positive (counterclockwise) vorticity shed from the bottom half of the cylinder. In addition, the negative vorticity from the wave peak strengthens the clockwise (negative) vorticity shed from the top half of the cylinder. These interactions combined with the ground effect skewed the flow away from the ground. The base flow is skewed upward for all the near-ground cases. However, this skew is larger when the cylinder is located over the wavy ground. The vortex shedding frequency is also altered due to the presence of the waves. The main eigenmode found for plain flow past a cylinder appears to become suppressed for cases closer to the ground. Limited particle image velocimetry experiments are reported which corroborate the finding from the stability analysis.

An inviscid vortex shedding model is numerically extended to simulate falling flat plates. The body and vortices separated from the edge of the body are described by vortex sheets. The vortex shedding model has computational limitations when the angle of incidence is small and the free vortex sheet approaches the body closely. These problems are overcome by using numerical procedures such as a method for a near-singular integral and the suppression of vortex shedding at the plate edge. The model is applied to a falling plate of flow regimes of various Froude numbers. For (text {Fr}=0.5), the plate develops large-scale side-to-side oscillations. In the case of (text {Fr}=1), the plate motion is a combination of side-to-side oscillations and tumbling and is identified as a chaotic type. For (text {Fr}=1.5), the plate develops to autorotating motion. Comparisons with previous experimental results show good agreement for the falling pattern. The dependence of change in the vortex structure on the Froude number and its relation with the plate motion is also examined.

A theoretical and experimental study was conducted to investigate the effect of injection angle on surface waves. Linear stability theory was utilized to obtain the analytical relation. In the experimental study, high-speed photography and shadowgraph techniques were used. Image processing codes were developed to extract information from photos. The results obtained from the theoretical relation were validated with the experimental results at different injection angles. In addition, at the injection angle of 90({^circ }), the theoretical results were evaluated with the experimental results of other researchers. This evaluation showed that the theory results were in good agreement with the experimental data. The proper orthogonal decomposition (POD) and the power spectra density (PSD) analysis were also used to investigate the effect of the injection angle on the flow structures. The results obtained from the linear stability were used to determine the maximum waves’ growth rate, and a relation was presented for the breakup length of the liquid jet at different injection angles. The breakup length results were compared with theory and published experimental data. The presented relation is more consistent with experimental data than other theories due to considering the nature of waves. The results showed that the instability of the liquid jet is influenced by three forces: inertial, surface tension, and aerodynamic. Therefore, Rayleigh–Taylor, Kelvin–Helmholtz, Rayleigh–Plateau, and azimuthal instabilities occur in the process. Decreasing the injection angle changes the nature of waves and shifts from Rayleigh–Taylor to Kelvin–Helmholtz. That reduces the wavelength and increases the growth rate of the waves. Axial waves have a significant impact on the physics of the waves and influence parameters. If axial waves are not formed, the growth rate of the waves is independent of the injection angle. An increase in the gas Weber number causes a change in the type of dominant waves and a greater instability of the liquid jet. In contrast, an increase in the liquid Weber number causes an enhancement in the resistance of the liquid jet against the transverse flow without changing the type of the dominant waves. Decreasing the density ratio reduces the effect of Rayleigh–Taylor waves and strengthens the Kelvin–Helmholtz waves. It causes two trends to be observed for the growth rate of waves at low spray angles, while one trend occurs at high spray angles.

The focus here is on a thin solid body passing through a channel flow and interacting with the flow. Unsteady two-dimensional interactive properties from modelling, analysis and computation are presented along with comparisons. These include the effects of a finite dilation or constriction, as the body travels through, and the effects of a continuing expansion of the vessel. Finite-time clashing of the body with the channel walls is investigated as well as the means to avoid clashing. Sustained oscillations are found to be possible. Wake properties behind the body are obtained, and broad agreement in trends between full-system and reduced-system responses is found for increased body mass.

Recently, Rim (Ocean Engng 239:711, 2021; J Ocean Engng Mar Energy 9:41-51, 2023 ) suggested an exact DtN artificial boundary condition to study the three-dimensional wave diffraction by stationary bodies. This paper is concerned with three-dimensional linear interaction between a submerged oscillating body with arbitrary shape and the regular water wave with finite depth. An exact Dirichlet-to-Neumann (DtN) boundary condition on a virtual cylindrical surface is derived, where the virtual surface is chosen so as to enclose the body and extract an interior subdomain with finite volume from the horizontally unbounded water domain. The DtN boundary condition is then applied to solve the interaction between the body and the linear wave in the interior subdomain by using boundary integral equation. Based on verification of the present model for a submerged vertical cylinder, the model is extended to the case of a submerged chamfer box with fillet radius in order to study 6-DoF oscillatory motion of the body under the free surface wave.