• Ocean Systems Engineering
• /
• v.12 no.1
• /
• pp.63-86
• /
• 2022

A refined projection-based purely Lagrangian meshfree method is presented towards reliable numerical analysis of fluid flow interactions with saturated/unsaturated porous media of uniform/spatially-varying porosities. The governing equations are reformulated on the basis of two-phase mixture theory with incorporation of volume fraction. These principal equations of mixture are discretized in the context of Incompressible SPH (Smoothed Particle Hydrodynamics) method. Associated with the consideration of governing equations of mixture, a new term arises in the source term of PPE (Poisson Pressure Equation), resulting in modified source term. The linear and nonlinear force terms are included in momentum equation to represent the resistance from porous media. Volume increase of fluid particles are taken into consideration on account of the presence of porous media, and hence multi-resolution ISPH framework is also incorporated. The stability and accuracy of the proposed method are thoroughly examined by reproducing several numerical examples including the interactions between fluid flow and saturated/unsaturated porous media of uniform/spatially-varying porosities. The method shows continuous pressure field, smooth variations of particle volumes and regular distributions of particles at the interface between fluid and porous media.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.4
• /
• pp.398-407
• /
• 2007

A smoothed particle hydrodynamics (SPH) method is applied for simulating two-dimensional free-surface problems. The SPH method based on the Lagrangian formulation provides realistic flow motions with violent surface deformation, fragmentation and reunification. In this study, the effect of computational parameters in SPH simulation is explored through two-dimensional dam-breaking and sloshing problem. The parameters to be considered are the speed of sound, the frequency of density re-initialization, the number of particle and smoothing length. Through a series of numerical test. detailed information was obtained about how SPH solution can be more stabilized and improved by adjusting computational parameters. Finally, some numerical simulations for various fluid flow problem were carried out based on the parameters chosen through the sensitivity study.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.1
• /
• pp.26-33
• /
• 2017

In this study, the surface marker method, one of the particle tracking methods, used to track the free surface is extended to cover the more general cases easily including the collision and separation of the free surface. In surface particle method to redistribute particles effectively using the grid system, the free surface is composed of the sum of quadrilaterals having four curves where fixed markers are placed at ends of each curve. Fixed markers are used to know how curves are connected to each other. The position of fixed markers can move as the free surface deforms but all fixed markers cannot be deleted during all time of simulation to keep informations of curve connection. In the case of the collision or separtion of the free surface where several curves can be intersected disorderly, severe difficulties can occur to define newly states of curve connection. In this study, the adaptable surface parTicle method without fixed markers is introduced. Intersection markers instead of the fixed markers are used to define quadrilaterals. The position of the intersection markers is defined to be the intersection point between the free surface and the edge of the grid and it can be added or deleted during the time of simulation to allow more flexibilities. To verify numerical schemes, two flow cases are simulated and the numerical results are compared with other's one and shown to be valid.

• Journal of Drive and Control
• /
• v.21 no.1
• /
• pp.8-15
• /
• 2024

This thesis investigates the jetting characteristics of an aerosol jet printing (AJP) process as a function of design and operating conditions. The governing equations of the AJP system are derived for experimentation and analysis. To understand the characteristics of the AJP system, it analyzes the jetting characteristics as a function of the flow rate of the carrier gas and the sheath gas, and the variation of the linewidth with the nozzle exit size based on particle tracking. The thesis focuses on computational fluid dynamics (CFD), which is a computer simulation. The particle tracking results obtained by CFD were analyzed using MATLAB. CFD analytical models can be analyzed in environments with different conditions and consider more specific situations than mathematical computational models. The validity of the CFD analysis is shown by comparing the experimental results with the CFD analysis.

• Journal of Drive and Control
• /
• v.21 no.2
• /
• pp.8-14
• /
• 2024

This paper investigates the jetting characteristics of an aerosol jet printing (AJP) process as a function of design and operating conditions. The governing equations of the AJP system are derived for experimentation and analysis. To understand the characteristics of the AJP system, this thesis analyzes the jetting characteristics as a function of the flow rate of the carrier gas and the sheath gas, and the variation of the linewidth with the nozzle exit size based on particle tracking. This thesis focuses on computational fluid dynamics (CFD), which is a computer simulation. The particle tracking results obtained by CFD were analyzed using MATLAB. CFD analytical models can be analyzed in environments with different conditions and consider more specific situations than mathematical computational models. The validity of the CFD analysis is shown by comparing the experimental results with the CFD analysis.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• Wind and Structures
• /
• v.13 no.4
• /
• pp.347-362
• /
• 2010

Researches on the Reynolds number effect on bridge decks have made slow progress due to the complicated nature of the subject. Heretofore, few studies on this topic have been made. In this paper, aerostatic coefficients, Strouhal number ($S_t$), pressure distribution and Reynolds number ($R_e$) of Great Belt East Bridge and Sutong Bridge were investigated based on deterministic vortex method (DVM). In this method, Particle Strength Exchange (PSE) was chosen to implement the simulation of the flow around bluff body and to analyze the micro-mechanism of the aerostatic loading and Reynolds number effect. Compared with the results obtained from wind tunnel tests, reliability of numerical simulation can be proved. Numerical results also showed that the Reynolds number effect on aerostatic coefficients and Strouhal number of the two bridges can not be neglected. In the range of the Reynolds number from $10^5$ to $10^6$, it has great effect on the Strouhal number of Sutong Bridge, while the St is difficult to obtain from wind tunnel tests in this range.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• Journal of the Korean Society of Visualization
• /
• v.9 no.4
• /
• pp.47-53
• /
• 2011

A water driven ejector loop was designed and constructed for air absorption. The used ejector was horizontally installed in the loop and annular water jet at the throat entrained air through the circular pipe placed at the center of the ejector. Wide range of water flow rate was provided using two kinds of pumps in the loop. The tested range of water flow rate was 100${\ell}$ /min to 1,000 ${\ell}$/min. Two-phase flow inside the ejector loop was simulated by CFD analysis. Homogeneous particle model was used for void fraction prediction. Water and air flow rates and pressure drop through the ejector were automatically recorded by using the LabView based data acquisition system. Flow characteristics and air bubble velocity field downstream of the ejector were investigated by two-phase flow visualization and PIV measurement based on bubble shadow images. Overall performance of the two-phase ejector predicted by the CFD simulation agrees well with that of the experiment.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.7 no.6
• /
• pp.1034-1043
• /
• 2015

This is the continuation of a numerical study on vortex shedding from a blunt trailing-edge of a hydrofoil. In our previous work (Lee et al., 2015), numerical schemes for efficient computations were successfully implemented; i.e. multiple domains, the approximation of domain boundary conditions using cubic spline functions, and particle-based domain decomposition for better load balancing. In this study, numerical results through a hybrid particle-mesh method which adopts the Vortex-In-Cell (VIC) method and the Brinkman penalization model are further rigorously validated through comparison to experimental data at the Reynolds number of $2{\times}10^6$. The effects of changes in numerical parameters are also explored herein. We find that the present numerical method enables us to reasonably simulate vortex shedding phenomenon, as well as turbulent wakes of a hydrofoil.