• 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.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.9 no.2
• /
• pp.109-119
• /
• 2006

This study investigated the effect of artificially enhanced mesozooplankton on the phytoplankton dynamics during fall blooming period using a mesocosm in Jangmok bay located in the Southern Sea of Korea in 2001. The four bags with 2,500 liter seawater containment were directly filled with the ambient water. And then, abundances of mesozooplankton in two experimental bags were treated 6 times higher than those in control bags by towing with net($300{\mu}m$) through the ambient water. Phytoplankton community between control and experimental bags were not significantly different in terms of chlorophyll-a(chl-a) concentration and standing crop (one-way ANOVA, p>0.05) during the study period. Initial high standing crop and chl-a concentration of phytoplankton drastically decreased and remained low until the end of the experiment in all bags. Diatoms, accounting for most of the phytoplankton community, consisted of Skeletonema costatum, Pseudo-nitzschia seriata, Chaetoceros curvisetus, Ch. debilis, Cerataulina pelagica, Thalassiosira pacifica, Cylindrotheca closterium, and Leptocylindrus danicus. Noctiluca scintillans dominated the temporal variation of mesozooplankton abundances, which peaked on Day 10 in the control and experimental bags, while the next dominant copepods showed their peak on Day 7. Shortly after mesozooplankton addition, copepod abundance in the experimental bags was obviously higher than that in the control bags on Day 1, however, it became similar to that in the control bags during the remnant period. It was supported by the higher abundance and length of both ctenophores and hydromedusae in experimental bags relative to the control bags. However, the cascading trophic effect, commonly leading to re-increase of phytoplankton abundance, was not found in the experimental bags, indicating that copepods were not able to control the phytoplankton in the bags based on the low grazing rate of Acartia erythraea. Besides that, rapidly sunken diatoms in the absence of natural turbulence as well as N-limited condition likely contributed the no occurrence of re-increased phytoplankton in the experimental bags.

• Korean Journal of Agricultural and Forest Meteorology
• /
• v.10 no.3
• /
• pp.82-93
• /
• 2008

In order to ensure a standardized data analysis of the eddy covariance measurements, Hong and Kim's quality control program has been updated and used to process eddy covariance data measured at two levels on the main flux tower at Gwangneung site from January to May in 2005. The updated program was allowed to remove outliers automatically for $CO_2$ and latent heat fluxes. The flag system consists of four quality groups(G, D, B and M). During the study period, the missing data were about 25% of the total records. About 60% of the good quality data were obtained after the quality control. The number of record in G group was larger at 40m than at 20m. It is due that the level of 20m was within the roughness sublayer where the presence of the canopy influences directly on the character of the turbulence. About 60% of the bad data were due to low wind speed. Energy balance closure at this site was about 40% during the study period. Large imbalance is attributed partly to the combined effects of the neglected heat storage terms, inaccuracy of ground heat flux and advection due to local wind system near the surface. The analysis of wind direction indicates that the frequent occurrence of positive momentum flux was closely associated with mountain valley wind system at this site. The negative $CO_2$ flux at night was examined in terms of averaging time. The results show that when averaging time is larger than 10min, the magnitude of calculated $CO_2$ fluxes increases rapidly, suggesting that the 30min $CO_2$ flux is influenced severely by the mesoscale motion or nonstationarity. A proper choice of averaging time needs to be considered to get accurate turbulent fluxes during nighttime.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.6
• /
• pp.374-385
• /
• 2013

Sediment particle motions have been numerically simulated over a sinusoidal ripple. Turbulent boundary layer flows are generated by Large Eddy Simulation, and the sediment particle motions are simulated using Lagrangian particle tracking method. Two unsteady flow conditions are used in the experiment by employing two different wave amplitudes while keeping other conditions such as wave period same. As expected, the amount of suspended sediment particles is clearly dependent on the wave amplitude as it is increasing with increasing flow intensity. However, it is also observed that the pattern of suspension may be different as well due to the only different condition caused by wave amplitude. Specially, the time of maximum sediment suspension within the wave period is not coincident between the two cases because sediment suspension is strongly affected by the existence of turbulent eddies that are formed at different times over the ripple between the two cases as well. The role of these turbulent eddies on sediment suspension is important as it is also confirmed in previous researches. However, it is also found the time of these eddies' formation may also dependent on the wave amplitude over rippled beds. Therefore, it has been proved that various flow as well as geometric conditions under waves has to be considered in order to have better understanding on the sediment suspension process over ripples. In addition, it is found that high turbulent energy and strong upward flow velocities occur during the time of eddy formation, which also supports high suspension rate at these time steps. The results indicate that the relationship between the structure of flows and bedforms has to be carefully examined in studying sediment suspension at coastal regions.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.