Search results

This paper, which is concerned with fluid‐structure interaction analysis, is a sequel to our earlier paper which gave an introduction to the numerical treatment of such systems. The paper is divided into five main sections. In the first two, a state‐of‐the‐art review on near‐field and far‐field fluid structure interaction is presented. In attempting to highlight where current research should be directed, only the most widely used computer codes are reviewed in the third section. Conclusions are presented in the fourth section.

The arbitrary Lagrangian—Eulerian (ALE) formulation, which is already well established in the hydrodynamics and fluid‐structure interaction fields, is extended to materials with memory, namely, non‐ linear path‐dependent materials. Previous attempts to treat non‐ linear solid mechanics with the ALE description have, in common, the implicit interpolation technique employed. Obviously, this implies a numerical burden which may be uneconomical and may induce to give up this formulation, particularly in fast‐transient dynamics where explicit algorithms are usually employed. Here, several applications are presented to show that if adequate stress updating techniques are implemented, the ALE formulation could be much more competitive than classical Lagrangian computations when large deformations are present. Moreover, if the ALE technique is interpreted as a simple interpolation enrichment, adequate—in opposition to distorted or locally coarse—meshes are employed. Notice also that impossible computations (or at least very involved numerically) with a Lagrangian code are easily implementable in an ALE analysis. Finally, it is important to observe that the numerical examples shown range from a purely academic test to real engineering simulations. They show the effective applicability of this formulation to non‐linear solid mechanics and, in particular, to impact, coining or forming analysis.

Presents a finite element/volume method for non‐linear aeroelasticity analyses of turbomachinery blades. The method uses an Arbitrary Lagrangian‐Eulerian (ALE) kinematical description of the fluid domain, in which the grid points can be displaced independently of the fluid motion. In addition, it employs an iterative implicit formulation similar to that of the Implicit‐continuous Eulerian (ICE) technique, making it applicable to flows at all speeds. A deforming mesh capability that can move the grid to conform continuously to the instantaneous shape of an aeroelastically deforming body without excessive distortion is also included in the algorithm. The unsteady aerodynamic loads are obtained using inviscid Euler equations. The model for the solid is general and can accommodate any spatial or modal representation of the structure. Determines the flutter stability of the system by studying the aeroelastic time response histories which are obtained by integration of the coupled equations of motion for both the fluid and the structure. Develops and demonstrates in 2D the formulation, which includes several corrections for better numerical stability. The cases studied include NACA64A006 and NACA0012 aerofoils and the EPFL Configuration 4 cascade. Finds the results from the numerical indicate good overall agreement with other published work and hence demonstrates the suitability of an ICED‐ALE formulation for turbomachinery applications.

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.

The simulation of the fluid–solid interaction (FSI) problem is important for both academic studies and engineering applications. However, the numerical approach for simulating the FSI problems is a great challenge owing to the large discrepancy of material properties and inconsistent description of grid motion between the fluid and solid domains. The difficulties will be further increased if there are multiple materials in the fluid region. In these complicated applications, interface reconstruction, multi-material advection and FSI must be all taken into account. This paper aims to present an effective integrated work of multi-material arbitrary Lagrangian Eulerian (MMALE) method, finite element (FE) method and the continuum analogy method to simulate the complex FSI problems involving multi-material flow. The coupled method is used to simulate the three-dimensional CONT test and the blast-plate interaction. The numerical results show good agreement with the benchmark and the experiment data, which indicates that the presented method is effective for solving the complicated FSI problems.

MMALE and FE methods are used to simulate fluid and solid regions, respectively. The interfacial nodes of fluid and solid are required to be coincident in the whole simulation so the interacted force can be easily and accurately calculated. To this end, the continuum analogy method is used in the rezoning phase.

The coupled method is used to simulate the three-dimensional CONT test and the blast-plate interaction. The numerical results show good agreement with the benchmark and the experiment data, which indicates that the presented method is effective for solving the complicated FSI problems.

To the best of the authors’ knowledge, this is the first time that the ALE method, moment of fluid interface reconstruction method, continuum analogy method and the FE method are combined to solve complicated practical problems.

In oil and gas industries, the presence of sand particles in produced oil and natural gas represents a major concern because of the associated erosive wear occurring in various flow passages. Erosion in the tube entrance region of a typical shell and tube heat exchanger is numerically predicted.

The erosion rates are obtained for different flow rates and particle sizes assuming low particle concentration. The erosion prediction is based on using a mathematical model for simulating the fluid velocity field and another model for simulating the motion of solid particles. The fluid velocity (continuous phase) model is based on the solution of the time‐averaged governing equations of 3D turbulent flow while the particle‐tracking model is based on the solution of the governing equation of each particle motion taking into consideration the viscous and gravity forces as well as the effect of particle rebound behavior.

The results show that the location and number of eroded tubes depend mainly on the particle size and velocity magnitude at the header inlet. The rate of erosion depends exponentially on the velocity. The particle size shows negligible effect on the erosion rate at high velocity values and the large‐size particles show less erosion rates compared to the small‐size particles at low values of inlet flow velocities.

In oil and gas industries, the presence of sand particles in produced oil and natural gas represents a major concern because of the associated erosive wear occurring in various flow passages. The results indicate that erosion in shell and tube heat exchanger can be minimized through the control of velocity inlet to the header.

This study is concerned with the direct numerical simulation (DNS) of a turbulent channel flow by an improved vortex in cell (VIC) method. The paper aims to discuss these issues.

First, two improvements for VIC method are proposed to heighten the numerical accuracy and efficiency. A discretization method employing a staggered grid is presented to ensure the consistency among the discretized equations as well as to prevent the numerical oscillation of the solution. A correction method for vorticity is also proposed to compute the vorticity field satisfying the solenoidal condition. Second, the DNS for a turbulent channel flow is conducted by the improved VIC method. The Reynolds number based on the friction velocity and the channel half width is 180.

It is highlighted that the simulated turbulence statistics, such as the mean velocity, the Reynolds shear stress and the budget of the mean enstrophy, agree well with the existing DNS results. It is also shown that the organized flow structures in the near-wall region, such as the streaks and the streamwise vortices, are favourably captured. These demonstrate the high applicability of the improved VIC method to the DNS for wall turbulent flows.

This study enables the VIC method to perform the DNS for wall turbulent flows.

This paper aims to validate an axial turbine modified nozzle design, looking for a reduction of the nozzle erosion process during operation in power plants.

The approach taken is numerical simulation using the computational fluid dynamics (CFD) tool, comparing original and proposed/modified nozzle designs.

The paper provides information about how to achieve a solution of the turbine operational problem (abrasive wear) by an analysis of flow patterns under a variety of conditions.

It does not give a detailed interpretation of flow behaviour due to the lack of validation data.

A very useful flow simulation methodology that can be used in industry is provided.

The proposed design modification of an axial turbine nozzle with the aid of CFD simulation has not been performed yet. This paper investigates the possibility of nozzle erosion reduction by modifying local flow patterns.

The internal combustion engine operated on gaseous fuels shows great potential in terms of integration of the renewable and traditional sources for an effective solution for clean energy production challenge. Different fuel mixtures that can be used to power the engine are characterized by various combustion properties, which can affect its overall efficiency. The purpose of this paper is to provide reasonable answer, how the operation condition can change due to different fuel, without enormous cost of prototyping processes using physical models a digital model can be seen as promising technique.

Presented work discusses the application, and extensive description of two commercial codes Ansys Fluent and Forte for modeling stationary engine fueled by compressed natural gas (CNG) and biogas. To check the model accuracy, all carried out numerical results were compared against experimental data collected at in-house test rig of single cylinder four stroke engine. The impacts of tested gaseous fuel on the engine working conditions and emission levels were investigated.

Carried out simulations showed good agreement with experimental data for investigated cases. Application on numerical models give possibility to visualize flame front propagation and pollutant formation for tested fuels. The biogas fuel has shown the impaired early flame phase, which led to longer combustion, lower efficiency, power output, repeatability and in some cases higher HC and carbon monoxide (CO) emissions as a result of combustion during the exhaust stroke. Looking at the CO formation it was observed that it instantly accrue with flame front propagation as a result of methane oxidation, while for NOx formation revers effect was seen.

The application of new approach for modeling combustion process in stationary engines fueled by CNG and alternative biogas fuel has been discussed. The cons and pros of the Forte code in terms of its application for engine prototaping process has been discussed.