Theoretical and Computational Fluid Dynamics最新文献

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.

The effect of non-Boltzmann energy distributions on the free propagation of shock waves through a monoatomic gas is investigated via theory and simulation. First, the non-Boltzmann heat capacity ratio (gamma ), as a key property for describing shock waves, is derived from first principles via microcanonical integration. Second, atomistic molecular dynamics simulations resembling a shock tube setup are used to test the theory. The presented theory provides heat capacity ratios ranging from the well-known (gamma = 5/3) for Boltzmann energy-distributed gas to (gamma rightarrow 1) for delta energy-distributed gas. The molecular dynamics simulations of Boltzmann and non-Boltzmann driven gases suggest that the shock wave propagates about 9% slower through the non-Boltzmann driven gas, while the contact wave appears to be about 4% faster if it trails non-Boltzmann driven gas. The observed slowdown of the shock wave through applying a non-Boltzmann energy distribution was found to be consistent with the classical shock wave equations when applying the non-Boltzmann heat capacity ratio. These fundamental findings provide insights into the behavior of non-Boltzmann gases and might help to improve the understanding of gas dynamical phenomena.

An adjoint-based method is presented for determining manufacturing tolerances for aerodynamic surfaces with natural laminar flow subjected to wavy excrescences. The growth of convective unstable disturbances is computed by solving Euler, boundary layer, and parabolized stability equations. The gradient of the kinetic energy of disturbances in the boundary layer (E) with respect to surface grid points is calculated by solving adjoints of the governing equations. The accuracy of approximations of (Delta E), using gradients obtained from adjoint, is investigated for several waviness heights. It is also shown how second-order derivatives increase the accuracy of approximations of (Delta E) when surface deformations are large. Then, for specific flight conditions, using the steepest ascent and the sequential least squares programming methodologies, the waviness profile with minimum (L2-)norm that causes a specific increase in the maximum value of N- factor, (Delta N), is found. Finally, numerical tests are performed using the NLF(2)-0415 airfoil to specify tolerance levels for (Delta {N}) up to 2.0 for different flight conditions. Most simulations are carried out for a Mach number and angle of attack equal to 0.5 and (1.25^{circ }), respectively, and with Reynolds numbers between (9times 10^6) and (15times 10^6) and for waviness profiles with different ranges of wavelengths. Finally, some additional studies are presented for different angles of attack and Mach numbers to show their effects on the computed tolerances.

Under the low-noise Mach 3 flight conditions for a supersonic passenger aircraft having unswept wings with a thin parabolic airfoil, laminar-turbulent transition is due to amplification of the first mode. Stability of a local self-similar boundary layer over such a wing is investigated both using the (e^{N}) method in the framework of linear stability theory and direct numerical simulation (DNS). It is found that the instability amplitude should reach a maximum over the entire spectral range above the profiles of 2.5% and thicker. The locus of maximum appears at the trailing edge and moves to the leading edge as the profile becomes thicker, while the maximum amplitude decreases. The theoretical findings are supported by DNS of the linear wave packets propagating in the boundary layer. Significance of these results to the design of laminar supersonic wings is discussed.

A low-order physics-based model to simulate the unsteady flow response to airfoils undergoing large-amplitude variations of the camber is presented in this paper. Potential-flow theory adapted for unsteady airfoils and numerical methods using discrete-vortex elements are combined to obtain rapid predictions of flow behavior and force evolution. To elude the inherent restriction of thin-airfoil theory to small flow disturbances, a time-varying chord line is proposed in this work over which to satisfy the appropriate boundary condition, enabling large deformations of the camber line to be modeled. Computational fluid dynamics simulations are performed to assess the accuracy of the low-order model for a wide range of dynamic trailing-edge flap deflections. By allowing the chord line to rotate with trailing-edge deflections, aerodynamic loads predictions are greatly enhanced as compared to the classical approach where the chord line is fixed. This is especially evident for large-amplitude deformations.

Pressure-driven Newtonian fluid flow between grooved and flat surfaces is analysed with no-slip boundary conditions at walls. The effect of corrugation on the fluid flow is investigated using the mesh-free spectral method. The primary aim of the present work is to develop an asymptotic/semi-analytical theory for confined transverse flows to bridge the gap between the limits of thin and thick channels. The secondary aim is to calculate permeability with reference to the effect of wall corrugation (roughness) without the restriction of pattern amplitude. We performed mathematical modelling and evaluated the analytical solution for hydraulic permeability with respect to the flat channel. The Pad(acute{e}) approximate is employed to improve the solution accuracy of an asymptotic model. The results elucidate that permeability always follows a decreasing trend with increasing pattern amplitude using the spectral approach at the long-wave and short-wave limits. The prediction of the spectral model is more accurate than the asymptotic-based model by Stroock et al. (Anal Chem 74(20):5306, 2002) and Pad(acute{e}) approximate, regardless of the grooved depth and wavelength of the channel. The finite-element-based numerical simulation is also used to understand the usefulness of theoretical models. A very low computational time is required using the mesh-free spectral model as compared to the numerical study. The agreement between the present model and the fully resolved numerical results is gratifying. Regarding numerical values, we calculated the relative error for different theoretical models such as an asymptotic model, Pad(acute{e}) approximate, and a mesh-free spectral model. The spectral model always predicts the maximum relative error as less than (3 %), regardless of the large pattern amplitude and wavelength. In addition, the results of the molecular dynamic (MD) simulations by Guo et al. (Phys Rev Fluids 1(7):074102, 2016) and the theoretical model by Wang (Phys Fluids 15(5):1121, 2003) are found to be quantitatively compatible with the predictions of effective slip length from the spectral model in the thick channel limit.

We present a Graphic Processing Units (GPU) implementation of non-Newtonian Hele-Shaw flow that models the displacement of Herschel-Bulkley fluids along narrow eccentric annuli. This flow is characteristic of many long-thin flows that require extensive calculation due to an inherent nonlinearity in the constitutive law. A common method of dealing with such flows is via an augmented Lagrangian algorithm, which is often painfully slow. Here we show that such algorithms, although involving slow iterations, can often be accelerated via parallel implementation on GPUs. Indeed, such algorithms explicitly solve the nonlinear aspects only locally on each mesh cell (or node), which makes them ideal candidates for GPUs. Combined with other advances, the optimized GPU implementation takes (approx 2.5%) of the time of the original algorithm.

The paper presents a three-dimensional numerical study of the acoustic streaming induced by the dissipation of ultrasounds during their propagation in the air. The waves are generated by a circular acoustic source positioned at the center of the left wall of a parallelepipedic cavity. The simulations are performed with the lattice Boltzmann method associated with the D3Q19 multiple relaxation time model. A validation of this model is first performed by comparing the numerical and analytical acoustic intensities along the central axis of the acoustic source. The main objective of this study is to use two different methods to calculate the acoustic streaming flow. The first method is the direct calculation of the mean velocity fields as the mean values of the instantaneous velocities. The second method is an indirect technique, which first calculates the acoustic streaming force and then injects this force into the numerical code to produce the streaming. A comparison between the results obtained by the two methods was carried out and a good agreement was found between them. These different investigations, rather new in three-dimensional configurations, have allowed us to discuss the advantages and limitations of the lattice Boltzmann approach to simulate real situations of wave propagation and acoustic streaming.