• Ecology and Resilient Infrastructure
• /
• v.8 no.4
• /
• pp.165-178
• /
• 2021

We studied driftwood behaviors including generation and deposition in a tsunami using a numerical simulation. We used an integrated two-dimensional numerical model, which included a driftwood dynamics model. The study area was Sendai, Japan. Observation data collected by Inagaki et al. (2012) were used to verify the simulation results by comparing them with driftwood deposition patterns. A simplified model was developed to consider the threshold of driftwood generation by the drag force of water flows. To consider the volume of driftwood generated, we estimated the total wood number in the study area using Google Earth. Therefore, we simulated more than 13,000 pieces of driftwood that were generated and transported inland from approximately 300,000 trees that were growing in the forest. The final distribution of the driftwood was similar to the observation data. The reproducibility of the generation and deposition patterns of driftwood showed good agreement in terms of longitudinal deposition pattern. In the future, a sensitivity analysis on driftwood parameters, such as the size of the wood, boundary conditions, and grid size, will be implemented to predict the travel patterns of driftwood. Such modeling will be a useful methodology for disaster prediction based on water flow and driftwood.

• Explosives and Blasting
• /
• v.37 no.3
• /
• pp.13-24
• /
• 2019

This paper was prepared to investigate the behavior of fragments in underwater torpedo explosion beneath a frigate or surface ship by using an explicit finite element analysis. In this study, a fluid-structure interaction (FSI) methodology, called the multi-material arbitrary Lagrangian-Eulerian (MM-ALE) approach in LS-DYNA, was employed to obtain the responses of the torpedo fragments and frigate hull to the explosion. The Euler models for the analysis were comprised of air, water, and explosive, while the Lagrange models consisted of the fragment and the hull. The focus of this modeling was to examine whether a worst-case fragment could penetrate the frigate hull located close (4.5 m) to the exploding torpedo. The simulation was performed in two separate steps. At first, with the assumption that the expanding skin of the torpedo had been torn apart by consuming 30% of the explosive energy, the initial velocity of the worst-case fragment was sought based on a well-known experimental result concerning the fragment velocity in underwater bomb explosion. Then, the terminal velocity of the worst-case fragment that is expected to occur before the fragment hit the frigate hull was sought in the second step. Under the given conditions, the possible initial velocities of the worst-case fragment were found to be very fast (400 and 1000 m/s). But, the velocity difference between the fragment and the hull was merely 4 m/s at the instant of collision. This result was likely to be due to both the tremendous drag force exerted by the water and the non-failure condition given to the frigate hull. Anyway, at least under the given conditions, it is thought that the worst-case fragment seldom penetrate the frigate hull because there is no significant velocity difference between them.

• Composites Research
• /
• v.32 no.5
• /
• pp.284-289
• /
• 2019

In the present paper, the dynamic force-moment equilibrium equations, driving power and energy equations are analyzed to formulate the equation for fuel economy(km/liter) equivalent to the driving distance (km) divided by the fuel volume (liter) of the vehicle, a selected model of gasoline powered KIA K3 (1.6v). In addition, the effects of the dynamic parameters such as speed of vehicle (V), vehicle total weight(M), rolling resistance ($C_r$) between tires and road surface, inclined angle of road (${\theta}$), as well as the aerodynamic parameters such as drag coefficient ($C_d$) of vehicle, air density(${\rho}$), cross-sectional area (A) of vehicle, wind speed ($V_w$) have been analyzed. And the possibility of alternative materials such as lightweight metal alloys, fiber reinforced plastic composite materials to replace the conventional steel and casting iron materials and to reduce the weight of the vehicle has been investigated by Ashby's material index method. Through studies, the following results were obtained. The most influencing parameters on the fuel economy at high speed zone (100 km/h) were V, the aerodynamic parameters such as $C_d$, A, ${\rho}$, and $C_r$ and M. While at low speed zone (60 km/h), they are, in magnitude order, dynamic parameters such as V, M, $C_r$ and aerodynamic ones such as $C_d$, A, and ${\rho}$, respectively.

• Journal of the Korean Society of Marine Environment & Safety
• /
• v.28 no.4
• /
• pp.648-656
• /
• 2022

As environmental regulations such as the International Maritime Organization (IMO)'s strategy to reduce greenhouse gases(GHG) are strengthened, technology development such as eco-friendly ships and alternative fuels is expanding. As part of this, ship propulsion technology using energy reduction and wind propulsion technology is emerging, especially in shipping companies and shipbuilders. By securing wind propulsion technology and introducing empirical research into shipbuilding and shipping, a high value-added market using eco-friendly technology can be created. Moreover, by reducing the fuel consumption rate of operating ships, GHG can be reduced by 6-8%. Rotor Sail (RS) technology is to generate a hydrodynamic lift in the vertical direction of the cylinder when the circular cylinder rotates at a constant speed and passes through the fluid. This is called the Magnus effect, and this study attempted to propose a plan to increase propulsion efficiency through a numerical analysis study on turbulence flow characteristics around RS, a wind power assistance propulsion system installed on a ship. Therefore, C_L and C_D values according to SR and AR changes were derived as parameters that affect the aerodynamic force of the RS, and the flow characteristics around the rotor sail were compared according to EP application.

• Wind and Structures
• /
• v.35 no.1
• /
• pp.1-20
• /
• 2022

This research work presents an experimental study's outcomes to reveal the impact of an O-ring on the flow control over a sphere placed in a turbulent boundary layer. The investigation is performed quantitatively and qualitatively using particle image velocimetry (PIV) and dye visualization. The sphere model having a diamater of 42.5 mm is located in a turbulent boundary layer flow over a smooth plate for gap ratios of 0≤G/D≤1.5 at Reynolds number of 5 × 10³. Flow characteristics, including patterns of instantaneous vorticity, streaklines, time-averaged streamlines, velocity vectors, velocity fluctuations, Reynolds stress correlations, and turbulence kinetic energy (), are compared and discussed for a naked sphere and spheres having O-rings. The boundary layer velocity gradient and proximity of the sphere to the flat plate profoundly influence the flow dynamics. At proximity ratios of G/D=0.1 and 0.25, a wall jet is formed between lower side of the sphere and flat plate, and velocity fluctuations increase in regions close to the wall. At G/D=0.25, the jet flow also induces local flow separations on the flat plate. At higher proximity ratios, the velocity gradient of the boundary layer causes asymmetries in the mean flow characteristics and turbulence values in the wake region. It is observed that the O-ring with various placement angles (𝜃) on the sphere has a considerable alteration in the flow structure and turbulence statistics on the wake. At lower placement angles, where the O-ring is closer to the forward stagnation point of the sphere, the flow control performance of the O-ring is limited; however, its impact on the flow separation becomes pronounced as it is moved away from the forward stagnation point. At G/D=1.50 for O-ring diameters of 4.7 (2 mm) and 7 (3 mm) percent of the sphere diameter, the -ring exhibits remarkable flow control at 𝜃=50° and 𝜃=55° before laminar flow separation occurrence on the sphere surface, respectively. This conclusion is yielded from narrowed wakes and reductions in turbulence statistics compared to the naked sphere model. The O-ring with a diameter of 3 mm and placement angle of 50° exhibits the most effective flow control. It decreases, in sequence, streamwise velocity fluctuations and length of wake recovery region by 45% and 40%, respectively, which can be evaluated as source of decrement in drag force.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.36 no.2
• /
• pp.93-103
• /
• 2023

Warships widely spread numerous chaffs using a blast, which form chaff clouds that create false radar cross-sections to deceive enemy radars. In this study, we established a numerical framework based on a one-way coupling of computational fluid dynamics and discrete element method to simulate the spatiotemporal distribution of chaff clouds for warships in the air. Using the framework, we investigated the effects of wind, initial chaff cartridge angle, and blast pressure on the distribution of chaff clouds. We observed three phases for the chaff cloud diffusion: radial diffusion by the explosion, omnidirectional diffusion by turbulence and collision, and gravity-induced diffusion by the difference in the fall speed. The wind moved the average position of the chaff clouds, and the diffusion due to drag force did not occur. The direction of radial diffusion by the explosion depended on the initial angle of the cartridge, and a more vertical angle led to a wider distribution of the chaffs. As the blast pressure increased, the chaff clouds spread out more widely, but the distribution difference in the direction of gravity was not significant.

• Advances in nano research
• /
• v.16 no.4
• /
• pp.325-340
• /
• 2024

There are several novel uses for dispersing many nanoparticles into a conventional fluid, including dynamic sealing, damping, heat dissipation, microfluidics, and more. Therefore, melting heat and mass transfer characteristics of a 3-D MHD Hybrid Nanofluid flow over a rotating disc with presenting dufour and soret effects are assessed numerically in this study. In this instance, we investigated both ferric sulfate and molybdenum disulfide as nanoparticles suspended within base fluid water. The governing partial differential equations are transformed into linked higher-order non-linear ordinary differential equations by the local similarity transformation. The collection of these deduced equations is then resolved using a Chebyshev spectral collocation-based algorithm built into the Mathematica software. To demonstrate how different instances of hybrid/ nanofluid are impacted by changes in temperature, velocity, and the distribution of nanoparticle concentration, examples of graphical and numerical data are given. For many values of the material parameters, the computational findings are shown. Simulations conducted for different physical parameters in the model show that adding hybrid nanoparticle to the fluid mixture increases heat transfer in comparison to simple nanofluids. It has been identified that hybrid nanoparticles, as opposed to single-type nanoparticles, need to be taken into consideration to create an effective thermal system. Furthermore, porosity lowers the velocities of simple and hybrid nanofluids in both cases. Additionally, results show that the drag force from skin friction causes the nanoparticle fluid to travel more slowly than the hybrid nanoparticle fluid. The findings also demonstrate that suction factors like magnetic and porosity parameters, as well as nanoparticles, raise the skin friction coefficient. Furthermore, It indicates that the outcomes from different flow scenarios correlate and are in strong agreement with the findings from the published literature. Bar chart depictions are altered by changes in flow rates. Moreover, the results confirm doctors' views to prescribe hybrid nanoparticle and particle nanoparticle contents for achalasia patients and also those who suffer from esophageal stricture and tumors. The results of this study can also be applied to the energy generated by the melting disc surface, which has a variety of industrial uses. These include, but are not limited to, the preparation of semiconductor materials, the solidification of magma, the melting of permafrost, and the refreezing of frozen land.

• Journal of Astronomy and Space Sciences
• /
• v.25 no.4
• /
• pp.405-414
• /
• 2008

To better understand the physical processes that control the high-latitude lower thermospheric dynamics, we analyze the divergence and vorticity of the high-latitude neutral wind field in the lower thermosphere during the southern summertime for different IMF conditions. For this study the National Center for Atmospheric Research Thermosphere-Ionosphere Electrodynamics General Circulation Model (NCAR-TIEG CM) is used. The analysis of the large-scale vorticity and divergence provides basic understanding flow configurations to help elucidate the momentum sources that ulti-mately determine the total wind field in the lower polar thermosphere and provides insight into the relative strengths of the different sources of momentum responsible for driving winds. The mean neutral wind pattern in the high-latitude lower thermosphere is dominated by rotational flow, imparted primarily through the ion drag force, rather than by divergent flow, imparted primarily through Joule and solar heating. The difference vorticity, obtained by subtracting values with zero IMF from those with non-zero IMF, in the high-latitude lower thermosphere is much larger than the difference divergence for all IMF conditions, indicating that a larger response of the thermospheric wind system to enhancement in the momentum input generating the rotational motion with elevated IMF than the corresponding energy input generating the divergent motion. the difference vorticity in the high-latitude lower thermosphere depends on the direction of the IMF. The difference vorticity for negative and positive $B_y$ shows positive and negative, respectively, at higher magnetic latitudes than $-70^{\circ}$. For negative $B_z$, the difference vorticities have positive in the dusk sector and negative in the dawn sector. The difference vorticities for positive $B_z$ have opposite sign. Negative IMF $B_z$ has a stronger effect on the vorticity than does positive $B_z$.

• Journal of Korean Society of Environmental Engineers
• /
• v.32 no.2
• /
• pp.165-174
• /
• 2010

The numerical modeling of a coal gasification reaction occurring in an entrained flow coal gasifier is presented in this study. The purposes of this study are to develop a reliable evaluation method of coal gasifier not only for the basic design but also further system operation optimization using a CFD(Computational Fluid Dynamics) method. The coal gasification reaction consists of a series of reaction processes such as water evaporation, coal devolatilization, heterogeneous char reactions, and coal-off gaseous reaction in two-phase, turbulent and radiation participating media. Both numerical and experimental studies are made for the 1.0 ton/day entrained flow coal gasifier installed in the Korea Institute of Energy Research (KIER). The comprehensive computer program in this study is made basically using commercial CFD program by implementing several subroutines necessary for gasification process, which include Eddy-Breakup model together with the harmonic mean approach for turbulent reaction. Further Lagrangian approach in particle trajectory is adopted with the consideration of turbulent effect caused by the non-linearity of drag force, etc. The program developed is successfully evaluated against experimental data such as profiles of temperature and gaseous species concentration together with the cold gas efficiency. Further intensive investigation has been made in terms of the size distribution of pulverized coal particle, the slurry concentration, and the design parameters of gasifier. These parameters considered in this study are compared and evaluated each other through the calculated syngas production rate and cold gas efficiency, appearing to directly affect gasification performance. Considering the complexity of entrained coal gasification, even if the results of this study looks physically reasonable and consistent in parametric study, more efforts of elaborating modeling together with the systematic evaluation against experimental data are necessary for the development of an reliable design tool using CFD method.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.40 no.3
• /
• pp.196-205
• /
• 2004

As far as an opening device of fishing gears is concerned, applications of a kite are under development around the world. The typical examples are found in the opening device of the stow net on anchor and the buoyancy material of the trawl. While the stow net on anchor has proved its capability for the past 20 years, the trawl has not been wildly used since it has been first introduced for the commercial use only without sufficient studies and thus has revealed many drawbacks. Therefore, the fundamental hydrodynamics of the kite itself need to ne studied further. Models of plate and canvas kite were deployed in the circulating water tank for the mechanical test. For this situation lift and drag tests were performed considering a change in the shape of objects, which resulted in a different aspect ratio of rectangle and trapezoid. The results obtained from the above approaches are summarized as follows, where aspect ratio, attack angle, lift coefficient and maximum lift coefficient are denoted as A, B, $C_L$ and $C_{Lmax}$ respectively : 1. Given the rectangular plate, $C_{Lmax}$ was produced as 1.46${\sim}$1.54 with A${\leq}$1 and 40$^{\circ}$${\leq}$B${\leq}$42$^{\circ}$. And when A${\geq}$1.5 and 20$^{\circ}$${\leq}$B${\leq}$22$^{\circ}$, $C_{Lmax}$ was 10.7${\sim}$1.11. Given the rectangular canvas, $C_{Lmax}$ was 1.75${\sim}$1.91 with A${\leq}$1 and 32$^{\circ}$${\leq}$B${\leq}$40$^{\circ}$. And when A${\geq}$1.5 and 18$^{\circ}$${\leq}$B${\leq}$22$^{\circ}$, $C_{Lmax}$ was 1.24${\sim}$1.40. Given the trapezoid kite, $C_{Lmax}$ was produced as 1.65${\sim}$1.89 with A${\leq}$1.5 and 34$^{\circ}$${\leq}$B${\leq}$44$^{\circ}$. And when A=2 and B=14${\sim}$48, $C_L$ was around 1. Given the inverted trapezoid kite, $C_{Lmax}$ was 1.57${\sim}$1.74 with A${\leq}$1.5 and 24$^{\circ}$${\leq}$B${\leq}$36$^{\circ}$. And when A=2, $C_{Lmax}$ was 1.21 with B=18$^{\circ}$. 2. For a model with A=1/2, an increase in B caused an increase in $C_L$ until $C_L$ has reached the maximum. Then there was a tendency of a gradual decrease in the value of $C_L$ and in particular, the rectangular kite showed a more rapid decrease. For a model with A=2/3, the tendency of $C_L$ was similar to the case of a model with A=1/2 but the tendency was a more rapid decrease than those of the previous models. For a model with A=1, and increase in B caused an increase in $C_L$ until $C_L$ has reached the maximum. Soon after the tendency of $C_L$ decreased dramatically. For a model with A=1.5, the tendency of $C_L$ as a function of B was various. For a model with A=2, the tendency of $C_L$ as a function of B was almost the same in the rectangular and trapezoid model. There was no considerable change in the models with 20$^{\circ}$${\leq}$B${\leq}$50$^{\circ}$. 3. The tendency of kite model's $C_L$ in accordance with increase of B was increased rapidly than plate models until $C_L$ has reached the maximum. Then $C_L$ in the kite model was decreased dramatically but in the plate model was decreased gradually. The value of $C_{Lmax}$ in the kite model was higher than that of the plate model, and the kite model's attack angel at $C_{Lmax}$ was smaller than the plate model's. 4. In the relationship between aspect ratio and lift force, the attack angle which had the maximum lift coefficient was large at the small aspect ratio models, At the large aspect ratio models, the attack angle was small. 5. There was camber vertex in the position in which the fluid pressure was generated, and the rectangular & trapezoid canvas had larger value of camber vertex when the aspect ratio was high, while the inverted trapezoid canvas was versa. 6. All canvas kite had larger camber ratio when the aspect ratio was high, and the rectangular & trapezoid canvas had larger one when the attack angle was high.