Search results

A numerical method that solves the Euler equations for compressible flow is used to study vortex stretching. The particular case simulated is supersonic flow M∞=1.2 α=10 degrees around the twisted and cambered cranked‐and‐cropped TKF delta wing of MBB. This geometry induces multiple leading‐edge vortices in a straining velocity field that brings about flow instabilities but the result is a state of statistical equilibrium. The discretization contains over 600,000 cells and offers sufficient degrees of freedom in the solution to resolve the small‐scale unstable modes that lead to disordered vortex flow.

The large‐scale numerical simulation of fluid flow is described as a discipline within the field of software engineering. As an example of such work, a vortex flowfield is analysed for its essential physical flow features, an appropriate mathematical description is presented (the Euler equations with an artificial viscosity model), a numerical algorithm to solve the mathematical equations is described, and the programming methodology which allows us to attain a very high degree of vectorization on the CYBER 205 is discussed. Four simulated flowfields with vorticity shed from wing edges are computed with up to as many as one million grid points and verify the realism of the simulation model. The computed solutions show all the qualitative features that are expected in these flows. The twisted cranked‐and‐cropped delta case is one where the leading‐edge vortex is highly stretched and unstable, displaying ultimately inviscid large‐scale turbulent‐like phenomena.

The centrifugal instability mechanism of boundary layers over concave surfaces is responsible for the development of quasi-periodic, counter-rotating vortices aligned in a streamwise direction known as Görtler vortices. By distorting the boundary layer structure in both the spanwise and the wall-normal directions, Görtler vortices may modify heat transfer rates. The purpose of this study is to conduct spatial numerical simulation experiments based on a vorticity–velocity formulation of the incompressible Navier–Stokes system of equations to quantify the role of the transition in the heat transfer process.

Experiments are conducted using an in-house, parallel, message-passing code. Compact finite difference approximations and a spectral method are used to approximate spatial derivatives. A fourth-order Runge–Kutta method is adopted for time integration. The Poisson equation is solved using a geometric multigrid method.

Results show that the numerical method can capture the physics of transitional flows over concave geometries. They also show that the heat transfer rates in the late stages of the transition may be greater than those for either laminar or turbulent ones.

The numerical method can be considered as a robust alternative to investigate heat transfer properties in transitional boundary layer flows over concave surfaces.

The purpose of this paper is to apply lattice Boltzmann equation method (LBM) with multiple relaxation time (MRT) model, to investigate lid‐driven flow in a three‐dimensional (3D), rectangular cavity, and compare the results with flow in an equivalent two‐dimensional (2D) cavity.

The second‐order MRT model is implemented in a 3D LBM code. The flow structure in cavities of different aspect ratios (0.25‐4) and Reynolds numbers (0.01‐1000) is investigated. The LBM simulation results are compared with those from numerical solution of Navier‐Stokes (NS) equations and with available experimental data.

The 3D simulations demonstrate that 2D models may predict the flow structure reasonably well at low Reynolds numbers, but significant differences with experimental data appear at high Reynolds numbers. Such discrepancy between 2D and 3D results are attributed to the effect of boundary layers near the side‐walls in transverse direction (in 3D), due to which the vorticity in the core‐region is weakened in general. Secondly, owing to the vortex stretching effect present in 3D flow, the vorticity in the transverse plane intensifies whereas that in the lateral plane decays, with increase in Reynolds number. However, on the symmetry‐plane, the flow structure variation with respect to cavity aspect ratio is found to be qualitatively consistent with results of 2D simulations. Secondary flow vortices whose axis is in the direction of the lid‐motion are observed; these are weak at low Reynolds numbers, but become quite strong at high Reynolds numbers.

The findings will be useful in the study of variety of enclosed fluid flows.

Relevant analyses are presented on the base of the compressible vortex method for simulating the development of two or three co-rotating vortices with different characteristic Mach numbers. The paper aims to discuss this issue.

In addition to having vorticity and dilatation properties, the vortex particles also carry density, enthalpy, and entropy. Taking co-rotating vortices in two-dimensional unsteady compressible flow for an example, truncation of unbounded domains via a nonreflecting boundary condition was considered in order to make the method computationally efficient.

For two identical vortices, the effect of the vortex Mach number on merging process is not evident; if two vortices have the same circulation rather than different radiuses, the vorticity and dilatation fields of the vortex under a vortex Mach number will be absorbed by the vortex under a higher vortex Mach number. For three vortices, if the original arrangement of the vortices is changed, the evolvement of the vorticity and dilatation fields is different.

The paper reveals new mechanism of the three co-rotating vortices by a feasible compressible vortex method.

This paper aims to provide discussions of a numerical method for bubbly flows and the interaction between a vortex ring and a bubble plume.

Small bubbles are released into quiescent water from a cylinder tip. They rise under the buoyant force, forming a plume. A vortex ring is launched vertically upward into the bubble plume. The interactions between the vortex ring and the bubble plume are numerically simulated using a semi-Lagrangian–Lagrangian approach composed of a vortex-in-cell method for the fluid phase and a Lagrangian description of the gas phase.

A vortex ring can transport the bubbles surrounding it over a distance significantly depending on the correlative initial position between the bubbles and the core center. The motion of some bubbles is nearly periodic and gradually extinguishes with time. These bubble trajectories are similar to two-dimensional-helix shapes. The vortex is fragmented into multiple regions with high values of Q, the second invariant of velocity gradient tensor, settling at these regional centers. The entrained bubbles excite a growth rate of the vortex ring's azimuthal instability with a formation of the second- and third-harmonic oscillations of modes of 16 and 24, respectively.

A semi-Lagrangian–Lagrangian approach is applied to simulate the interactions between a vortex ring and a bubble plume. The simulations provide the detail features of the interactions.

This study aims to provide discussions of the numerical method and the bubbly flow characteristics of an annular bubble plume.

The bubbles, released from the annulus located at the bottom of the domain, rise owing to buoyant force. These released bubbles have diameters of 0.15–0.25 mm and satisfy the bubble flow rate of 4.1 mm³/s. The evolution of the three-dimensional annular bubble plume is numerically simulated using the semi-Lagrangian–Lagrangian (semi-L–L) approach. The approach is composed of a vortex-in-cell method for the liquid phase and a Lagrangian description of the gas phase.

First, a new phenomenon of fluid dynamics was discovered. The bubbly flow enters a transition state with the meandering motion of the bubble plume after the early stable stage. A vortex structure in the form of vortex rings is formed because of the inhomogeneous bubble distribution and the fluid-surface effects. The vortex structure of the flow deforms as three-dimensionality appears in the flow before the flow fully develops. Second, the superior abilities of the semi-L–L approach to analyze the vortex structure of the flow and supply physical details of bubble dynamics were demonstrated in this investigation.

The semi-L–L approach is applied to the simulation of the gas–liquid two-phase flows.

This paper aims to study the flow and heat transfer inside a wavy enclosure filled with Cu-water nanofluid under magnetic field effect by parallel implemented meshfree approach.

The simulation has been carried out for a two-dimensional model with steady, laminar and incompressible flow of the nanofluid filled inside wavy enclosure in which one of the walls is sinusoidal such that the amplitude (A = 0.15) and number of undulations (n = 2) are fixed. A uniform magnetic field B₀ has been applied at an inclination angle γ. The governing equations for the transport phenomena have been solved numerically by implementing element-free Galerkin method (EFGM) with the sequential as well as parallel approach. The effect of various parameters, namely, nanoparticle volume fraction (φ), Rayleigh number (Ra), Hartmann number (Ha) and magnetic field inclination angle (γ) has been studied on the natural convection flow of nanofluid.

The results are obtained in terms of average Nusselt number calculated at the cold wavy wall, streamlines and isotherms. It has been observed that the increasing value of Rayleigh number results in increased heat transfer rate while the Hartmann number retards the fluid motion. On the other hand, the magnetic field inclination angle gives rise to the heat transfer rate up to its critical value. Above this value, the heat transfer rate starts to decrease.

The implementation of the magnetic field and its inclination has provided very interesting results on heat and fluid flow which can be used in the drug delivery where nanofluids are used in many physiological problems. Another important novelty of the paper is that meshfree method (EFGM) has been used here because the domain is irregular. The results have been found to be very satisfactory. In addition, parallelization of the scheme (which has not been implemented earlier in such problems) improves the computational efficiency.

The purpose of this study is to design numerical simulations of bubbly flow around a cylinder to better understand the characteristics of flow around a rigid obstacle.

The bubbly flow around a circular cylinder was numerically simulated using a semi-Lagrangian–Lagrangian method composed of a vortex-in-cell method for the liquid phase and a Lagrangian description of the gas phase. Additionally, a penalization method was applied to account for the cylinder inside the flow. The slip condition of the bubbles on the cylinder’s surface was enforced, and the outflow conditions were applied to the liquid flow at the far field.

The simulation clarified the characteristics of a bubbly flow around a circular cylinder. The bubbles were shown to move around and separate from both sides of the cylinder, because of entrainment by the liquid shear layers. Once the bubbly flow fully developed, the bubbles distributed into groups and were dispersed downstream of the cylinder. A three-dimensional vortex structure of various scales was also shown to form downstream, whereas a quasi-stable two-dimensional vortex structure was observed upstream. Overall, the proposed method captured the characteristics of a bubbly flow around a cylinder well.

A semi-Lagrangian–Lagrangian approach was applied to simulate a bubbly flow around a circular cylinder. The simulations provided the detail features of these flow phenomena.

This paper aims to provide improvements to the newest version of the k‐ ω turbulence model of Wilcox for convective heat transfer prediction in turbulent axisymmetric jets impinging onto a flat plate.

Improvements to the heat transfer prediction in the impingement zone are obtained using the stagnation flow parameter of Goldberg and the vortex stretching parameter of Wilcox. The third invariant of the strain rate tensor in the form of Shih et al. and the blending function of Menter are applied in order make negligible the influence of the impingement modifications in the benchmark flows for turbulence models. Further, it is demonstrated that for two‐dimensional jets impinging onto a flat plate the stagnation region Nusselt number predicted by the original k‐ ω model is in good agreement with direct numerical simulation (DNS) and experimental data. Also for two‐dimensional jets, the proposed modification is deactivated.

The proposed modification has been applied to improve the convective heat transfer predictions in the stagnation flow regions of axisymmetric jets impinging onto a flat plate with nozzle‐plate distances H/D = 2, 6, 10 and Reynolds numbers Re = 23,000 and 70,000. Comparison of the predicted and experimental mean and fluctuating velocity profiles is performed. The heat transfer rates along a flat plate are compared to experimental data. Significant improvements are obtained with respect to the original k‐ ω model.

The proposed modification is simple and can be added to the k‐ ω model without causing stability problems in the computations.