• The Journal of the Acoustical Society of Korea
• /
• v.39 no.2
• /
• pp.77-86
• /
• 2020

In this study, the Eulerian/Lagrangian one-way coupling method is proposed to predict flow noise due to Blade-Tip Vortex Cavitation (BTVC). The proposed method consists of four sequential steps: flow field simulation using Computational Fluid Dynamics (CFD) techniques, reconstruction of wing-tip vortex using vortex model, generation of BTVC using bubble dynamics model and acoustic wave prediction using the acoustic analogy. Because the CFD prediction of tip vortex structure generally suffers from severe under-prediction of its strength along the steamwise direction due to the intrinsic numerical damping of CFD schemes and excessive turbulence intensity, the wing-tip vortex along the freestream direction is regenerated by using the vortex modeling. Then, the bubble dynamics model based on the Rayleigh-Plesset equation was employed to simulate the generation and variation of BTVC. Finally, the flow noise due to BTVC is predicted by modeling each of spherical bubbles as a monople source whose strength is proportional to the rate of time-variation of bubble volume. The validity of the proposed numerical methods is confirmed by comparing the predicted results with the measured data.

• Advances in Computational Design
• /
• v.5 no.4
• /
• pp.417-443
• /
• 2020

Tunnels have been an integral part of human civilization. Due to complexity in its design and structure, the stability of underground structures under extreme loading conditions has utmost importance. Increased terrorism and geo-political conflicts have forced the engineers and researchers to study the response of underground structures, especially tunnels under blast loading. The present study has been carried out to seek the response of tunnel structures under blast load using the finite element technique. The tunnel has been considered in quartzite rock of northern India. The Mohr-Coulomb constitutive model has been adopted for the elastoplastic behaviour of rock. The rock model surrounding the tunnel has dimensions of 30 m x 30 m x 35 m. Both unlined and lined (concrete) tunnel has been studied. Concrete Damage Plasticity model has been considered for the concrete lining. Four different parameters (i.e., tunnel diameter, liners thickness, overburden depth and mass of explosive) have been varied to observe the behaviour under different condition. To carry out blast analysis, Coupled-Eulerian-Lagrangian (CEL) modelling has been adopted for modelling of TNT (Trinitrotoluene) and enclosed air. JWL (Jones-Wilkins-Lee) model has been considered for TNT explosive modelling. The paper concludes that deformations in lined tunnels follow a logarithmic pattern while in unlined tunnels an exponential pattern has been observed. The stability of the tunnel has increased with an increase in overburden depth in both lined and unlined tunnels. Furthermore, the tunnel lining thickness also has a significant effect on the stability of the tunnel, but in smaller diameter tunnel, the increase in tunnel lining thickness has not much significance. The deformations in the rock tunnel have been decreased with an increase in the diameter of the tunnel.

• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.80-88
• /
• 1997

A new wall impaction model for diesel spray is described in this paper. The gas phase is modelled in terms of the Eulerian continuum conservation equations of mass, momentum, energy and fuel vapour fraction. The liquid phase is modelled following the discrete droplet model approach. The droplet parcel contains many thousands of drops assumed to have the same size, temperature and velocity components. The droplet parcel equations of trajectory, momentum, mass and energy are written in Lagrangian form. The new drop-wall interaction model is proposed, which is based on experimental investigations on individual drops, and it is applied for the general non-orthogonal grid. The model is then assessed through comparison with experiments over a wide range of test conditions of sprays. The results are in good agreement with experimental data.

• Journal of Mechanical Science and Technology
• /
• v.14 no.7
• /
• pp.796-806
• /
• 2000

To understand the transient behavior of droplets after impingement in a diesel engine, a numerical model for diesel sprays impinging on a flat wall is newly developed by the proposition of several mathematical formulae to determine the post-impingement characteristics of droplets. The new model consists of three representative regimes such as rebound, deposition and splash. The gas phase is modeled in terms of the Eulerian conservation equations, and the dispersed phase is calculated using a discrete droplet model. To validate the new model, the calculated results are compared with several experimental data. The results show that the new model is generally in good agreement with the experimental data. Therefore, it is thought that the new model is acceptable for the prediction of transient behavior of wall sprays.

• Journal of Advanced Marine Engineering and Technology
• /
• v.26 no.2
• /
• pp.181-191
• /
• 2002

Spray calculation has been studied to understand the behavior of the spray in a combustion chamber But the spray dispersion has not been predicted properly in a high velocity injection spray or a wall impaction spray. In this study the extended grazing collision model is applied to improve the problem. The gas phase is modelled by the Eulerian continuum conservation equations of mass, momentum, energy and fuel vapour fraction. The liquid phase is modelled following the discrete droplet model approach in Lagrangian form. The droplet distributions, penetration, width and gas flows are compared for the cases with or without extended model. The extended collision model makes the results better.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.23-36
• /
• 2015

The $2^{nd}$ generation ice accretion analysis program has been developed and validated for various icing conditions. The essential feature of the $2^{nd}$ generation code lies in its capability of handling general 3-D geometry and improved accuracy. The entire velocity fields are obtained based on Navier-Stokes equations in order to take the massively separated flow field into account. Unlike $1^{st}$ generation code, the droplet trajectories are calculated using Eulerian approach, which is adopted to yield appropriate collection efficiency even in the shadow region. For improved thermodynamic analysis on the surfaces, water film model and modified Messinger model are newly included in the present analysis. The ice shape for a given time step is obtained by considering the exact amount of ice accreted on the surface. Each module of the icing analysis code has been seamlessly integrated on the OpenFOAM platform. The developed code was validated against available experimental data for 2D airfoils and 3D DLR-F4. Due to the lack of experimental data, the computed results of DLR-F4 were compared with those obtained from FENSAP-ICE, which is state-of-the-art 3D icing analysis code. It was clearly shown that the present code produces comparable results to those of FENSAP-ICE, in terms of prediction accuracy and the capability of handling general 3-D geometries.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.39 no.11
• /
• pp.861-869
• /
• 2015

In this study, wall to bed heat transfer and hydrodynamic characteristics in a conical fluidized bed combustor was investigated using computational fluid dynamics method. A two-fluid Eulerian-Eulerian model was used with applying the kinetic theory for granular flow(KTGF). The effects of the two drag models, Gidaspow and the Syamlal-O'Brien model, different inlet velocities($1.4U_{mf}{\sim}4U_{mf}$) and different particle sizes on the hydrodynamics and heat transfer were studied. The results showed that the hydrodynamic characteristics such as bed expansion ratio and pressure drop were not affected significantly by the drag models. But the heat transfer coefficient was different for the two drag models, especially at lower gas inlet velocities and small particle sizes.

• Nuclear Engineering and Technology
• /
• v.48 no.4
• /
• pp.859-869
• /
• 2016

Component-scale modeling of boiling is predominantly based on the Eulerian-Eulerian two-fluid approach. Within this framework, wall boiling is accounted for via the Rensselaer Polytechnic Institute (RPI) model and, within this model, the bubble is characterized using three main parameters: departure diameter (D), nucleation site density (N), and departure frequency (f). Typically, the magnitudes of these three parameters are obtained from empirical correlations. However, in recent years, efforts have been directed toward mechanistic modeling of the boiling process. Of the three parameters mentioned above, the departure diameter (D) is least affected by the intrinsic uncertainties of the nucleate boiling process. This feature, along with its prominence within the RPI boiling model, has made it the primary candidate for mechanistic modeling ventures. Mechanistic modeling of D is mostly carried out through solving of force balance equations on the bubble. Forces incorporated in these equations are formulated as functions of the radius of the bubble and have been developed for, and applied to, low-pressure conditions only. Conversely, for high-pressure conditions, no mechanistic information is available regarding the growth rates of bubbles and the forces acting on them. In this study, we use direct numerical simulation coupled with an interface tracking method to simulate bubble growth under high (up to 45 bar) pressure, to obtain the kind of mechanistic information required for an RPI-type approach. In this study, we compare the resulting bubble growth rate curves with predictions made with existing experimental data.

• Structural Engineering and Mechanics
• /
• v.32 no.2
• /
• pp.283-300
• /
• 2009

The characteristic response of a structure to blast load may be divided into two distinctive phases, namely the direct blast response during which the shock wave effect and localized damage take place, and the post-blast phase whereby progressive collapse may occur. A reliable post-blast analysis depends on a sound understanding of the direct blast effect. Because of the complex loading environment and the stress wave effects, the analysis on the direct effect often necessitates a high fidelity numerical model with coupled fluid (air) and solid subdomains. In such a modelling framework, an appropriate representation of the blast load and the high nonlinearity of the material response is a key to a reliable outcome. This paper presents a series of calibration study on these two important modelling considerations in a coupled Eulerian-Lagrangian framework using a hydrocode. The calibration of the simulated blast load is carried out for both free air and internal explosions. The simulation of the extreme dynamic response of concrete components is achieved using an advanced concrete damage model in conjunction with an element erosion scheme. Validation simulations are conducted for two representative scenarios; one involves a concrete slab under internal blast, and the other with a RC column under air blast, with a particular focus on the simulation sensitivity to the mesh size and the erosion criterion.