• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.386-390
• /
• 2016

Water aeration is an effective water treatment process, which involves the injection of air or air-water mixture into water treatment reservoir commonly through pipes. The main purpose of water aeration is to maintain healthy levels of dissolved oxygen (DO), which is the most important water quality factor. The pipes' operating conditions are important for controlling the efficiency and effectiveness of aeration process. Many studies have been conducted on two-phase flows in pipes, however, there are a few studies to deal with small s ale in millimeter. The main objective of this study is to perform 2-dimensional two-phase simulations inside various straight pipes using the computational fluid dynamic (CFD) OpenFOAM (Open source Field Operation And Manipulation) tools to examine the influence of flow patterns on bubble size, which is closely related to DO concentration in a water body. The both flow regimes, laminar and turbulence, have been considered in this study. For turbulence, Reynolds-averaged Navier-Stokes (RANS) has been applied. The coalescence and breakage of bubbles caused by random collisions and turbulent eddies, respectively, are considered in this research. Sauter mean bubble diameter and water velocity are compared against experimental data. The simulation results are in good agreement with the experimental measurements.

• Proceedings of the Korea Association of Crystal Growth Conference
• /
• 1996.06a
• /
• pp.179-200
• /
• 1996

The intrinsic instabilities of fluid flow occurred in the melt of the Czochralski crystal growth system Czochralski method, asymmetric flow patterns and temperature profiles in the melt have been studied by many researchers. The idea that the non-symmetric structure of the growing equipment is responsible for the asymmetric profiles is usually accepted at the first time. However further researches revealed that some intrinsic instabilities not related to the non-symmetric equipment structure in the melt could also appear. Ristorcelli had pointed out that there are many possible causes of instabilities in the melt. The instabilities appears because of the coupling effects of fluid flow and temperature profiles in the melt. Among the instabilities, the B nard type instabilities with no or low crucible rotation rates are analyzed by the visualizing experiments using X-ray radiography and the 3-D numerical simulation in this study. The velocity profiles in the Silicon melt at different crucible rotation rates were measured using X-ray radiography method using tungsten tracers in the melt. The results showed that there exits two types of fluid flow mode. One is axisymmetric flow, the other is asymmetric flow. In the axisymmetric flow, the trajectory of the tracers show torus pattern. However, more exact measurement of the axisymmetrc case shows that this flow field has small non-axisymmetric components of the velocity. When fluid flow is asymmetric, the tracers show random motion from the fixed view point. On the other hand, when the observer rotates to the same velocity of the crucible, the trajectory of the tracer show a rotating motion, the center of the motion is not same the center of the melt. The temperature of a point in the melt were measured using thermocouples with different rotating rates. Measured temperatures oscillated. Such kind of oscillations are also measured by the other researchers. The behavior of temperature oscillations were quite different between at low rotations and at high rotations. Above experimental results means that the fluid flow and temperature profiles in the melt is not symmetric, and then the mode of the asymmetric is changed when rotation rates are changed. To compare with these experimental results, the fluid flow and temperature profiles at no rotation and 8 rpm of crucible rotation rates on the same size of crucible is calculated using a 3-dimensional numerical simulation. A finite different method is adopted for this simulation. 50×30×30 grids are used. The numerical simulation also showed that the velocity and flow profiles are changed when rotation rates change. Futhermore, the flow patterns and temperature profiles of both cases are not axisymmetric even though axisymmetric boundary conditions are used. Several cells appear at no rotation. The cells are formed by the unstable vertical temperature profiles (upper region is colder than lower part) beneath the free surface of the melt. When the temperature profile is combined with density difference (Rayleigh-B nard instability) or surface tension difference (Marangoni-B nard instability) on temperature, cell structures are naturally formed. Both sources of instabilities are coupled to the cell structures in the melt of the Czochralski process. With high rotation rates, the shape of the fluid field is changed to another type of asymmetric profile. Because of the velocity profile, isothermal lines on the plane vertical to the centerline change to elliptic. When the velocity profiles are plotted at the rotating view point, two vortices appear at the both sides of centerline. These vortices seem to be the main reason of the tracer behavior shown in the asymmetric velocity experiment. This profile is quite similar to the profiles created by the baroclinic instability on the rotating annulus. The temperature profiles obtained from the numerical calculations and Fourier transforms of it are quite similar to the results of the experiment. bove esults intend that at least two types of intrinsic instabilities can occur in the melt of Czochralski growing systems. Because the instabilities cause temperature fluctuations in the melt and near the crystal-melt interface, some defects may be generated by them. When the crucible size becomes large, the intensity of the instabilities should increase. Therefore, to produce large single crystals with good quality, the behavior of the intrinsic instabilities in the melt as well as the effects of the instabilities on the defects in the ingot should be studied. As one of the cause of the defects in the large diameter Silicon single crystal grown by the

• Journal of the Korea Society for Simulation
• /
• v.27 no.2
• /
• pp.35-48
• /
• 2018

The study of air defense against North Korean tactical ballistic missiles (TBM) should consider the rapidly changing battlefield environment. The study for target re-designation for intercept missiles enables effective operation of friendly defensive assets as well as responses to dynamic battlefield. The researches that have been conducted so far do not represent real-time dynamic battlefield situation because the hit probability for the TBM, which plays an important role in the decision making process, is fixed. Therefore, this study proposes a target re-designation algorithm that makes decision based on hit probability which considers real-time field environment. The proposed method contains a trajectory prediction model that predicts the expected trajectory of the TBM from the current position and velocity information by using random forest and moving window. The predicted hit probability can be calculated through the trajectory prediction model and the simulator of the intercept missile, and the calculated hit probability becomes the decision criterion of the target re-designation algorithm for the missile. In the experiment, the validity of the methodology used in the TBM trajectory prediction model was verified and the superiority of using the hit probability through the proposed model in the target re-designation decision making process was validated.

• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.211-216
• /
• 2013

In general, smoothing filters regularize functions by reducing differences between adjacent values. The smoothing filters, therefore, can regularize inverse solutions and produce more accurate subsurface structure when we apply it to full waveform inversion. If we apply a smoothing filter with a constant coefficient to subsurface image or velocity model, it will make layer interfaces and fault structures vague because it does not consider any information of geologic structures and variations of velocity. In this study, we develop a selective smoothing regularization technique, which adapts smoothing coefficients according to inversion iteration, to solve the weakness of smoothing regularization with a constant coefficient. First, we determine appropriate frequencies and analyze the corresponding wavenumber coverage. Then, we define effective maximum wavenumber as 99 percentile of wavenumber spectrum in order to choose smoothing coefficients which can effectively limit the wavenumber coverage. By adapting the chosen smoothing coefficients according to the iteration, we can implement multi-scale full waveform inversion while inverting multi-frequency components simultaneously. Through the successful inversion example on a salt model with high-contrast velocity structures, we can note that our method effectively regularizes the inverse solution. We also verify that our scheme is applicable to field data through the numerical example to the synthetic data containing random noise.

• Journal of Korea Water Resources Association
• /
• v.33 no.2
• /
• pp.181-193
• /
• 2000

This paper presents the flow and the suspended sediment movement over ripples for oscillatory flows. A new numerical model system is developed, and applied to a laboratory experimental condition of regular waves and a fictitious condition of irregular waves. The flow field is obtained from a programme proposed by Kim et. al.(1994), which is a modified version of SOLA based on SMAC scheme. The sub-model solves the continuity and Reynolds momentum equations in the x-z plane. The wave orbital velocities, shear stresses, and pressure are all reasonably reproduced by the model. The model results on the vertical velocity component show good agreement with the measurements. The suspended sediment transport sub-model is newly set up to solve the advection-diffusion equation of suspended sediment using a split method, and involving a special shear entrainment from the whole ripple surface. The calculated suspended sediment concentrations for regular waves show reasonable agreement with measurements at Deltaflume. The model results for random waves show that the suspended sediment concentration is higher than those for regular waves and that the sediment diffuses higher than for regular waves with the significant wave height and the peak wave period of the irregular waves.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.12
• /
• pp.1035-1044
• /
• 2013

Because of the high specific stiffness and strength inherent in the sandwich structure composed of facesheet that resists in-plane loads and a core that resists out-of-plane loads, it is often used for large and light-weighted structures. However, inevitably the increased flexibility allows greater deformation-based disturbances in the structures. Thus, it is necessary to analyze the structural safety. To obtain more accurate analytical results, the input disturbances must more closely simulate real load conditions; to improve accuracy, non-linear elements such as gust effects were considered. In addition, the structural safety was analyzed for the iso-grid sandwich panel structure using fluid-structure interactions. For a more realistic simulation, flow velocity fields, which consider the effects of irregular gust fluctuation, were generated and the coupled field was analyzed by mapping the pressure and displacement.

• The Journal of the Acoustical Society of Korea
• /
• v.34 no.1
• /
• pp.36-45
• /
• 2015

Recently, a lot of studies have been made about the methods used to generate turbulent velocity fields stochastically in order to effectively predict broadband flow noise. Among them, the FRPM (Fast Random Particle Mesh) method which generates turbulence with specific statistical properties using turbulence kinetic energy and dissipation obtained from the steady solution of the RANS (Reynolds Averaged Navier-Stokes) equations has been successfully applied. However, the FRPM method cannot be applied to the flow noise problems involving intrinsic unsteady characteristics such as centrifugal fan. In this paper, to effectively predict the broadband noise generated by centrifugal fan, U-FRPM (unsteady FRPM) method is developed by extending the FRPM method to be combined with the unsteady numerical solutions of the unsteady RANS equations to generate the turbulence considered as broadband noise sources. Firstly, an unsteady flow field is obtained from the unsteady RANS equations through CFD (Computational Fluid Dynamics). Then, noise sources are generated using the U-FRPM method combined with acoustic analogy. Finally, the linear propagation model which is realized through BEM (Boundary Element Method) is combined with the generated sources to predict broadband noise at the listeners' position. The proposed technique is validated to compare its prediction result with the measured data.

• The Journal of Engineering Geology
• /
• v.15 no.3
• /
• pp.245-255
• /
• 2005

Three dimensional Monte-Carlo simulations of non-ergodic transport of a non-reactive solute plume by steady-state groundwater flow under a uniform mean velocity in isotropic heterogeneous aquifers were conducted. The log-normally distributed hydraulic conductivity, K(x), is modeled as a random field. Significant efforts are made to reduce the simulation uncertainties. Ensemble averages of the second spatial moments of the plume, $$lt;S_{ij}'(t',l')$gt;$ and plume centroid variances, $$lt;R_{ij}'(t',l')$gt;$ were simulated with 3200 Monte Carlo runs for three variances of log K, $\omega^2_y1.0,,2.5,$ and 5.0, and three dimensionless lengths of line plume sources ( l=,5 and 10) normal to the mean velocity. The simulated second spatial moment and the plume centroid variance in longitudinal direction fit well to the first order theoretical results while the simulated transverse moments are not fit well with the first order results. The first order theoretical results definitely underestimated the simulated transverse second spatial moments for the aquifers of large u: and small initial plume sources. The ergodic condition for the second spatial moments is far from reaching, and the first order theoretical results of the transverse second spatial moment of the ergodic plume slightly underestimated the simulated moments.

• Journal of Soil and Groundwater Environment
• /
• v.10 no.6
• /
• pp.20-31
• /
• 2005

Three dimensional Monte-Carlo simulations of non-ergodic transport of a lion-reactive solute plume by steady-state groundwater flow under a uniform mean velocity in isotropic heterogeneous aquifers were conducted. The log-normally distributed hydraulic conductivity, K(x), is modeled as a random field. Significant efforts are made to reduce tile simulation uncertainties. Ensemble averages of the second spatial moments of the plume and plume centroid variances were simulated with 1600 Monte Carlo runs for three variances of log K, ${\sigma}_Y^2=0.09,\;0.23$, and 0.46, and three dimensionless lengths of line plume sources normal to the mean velocity. The simulated second spatial moment and the plume centroid variance in longitudinal direction fit well to the first order theoretical results while the simulated transverse moments are generally larger than the first order results. The first order theoretical results significantly underestimated the simulated dimensionless transverse moments for the aquifers of large ${\sigma}_Y^2$ and large dimensionless time. The ergodic condition for the second spatial moments is far from reaching in all cases simulated, and transport In transverse directions may reach ergodic condition much slower than that in longitudinal direction. The evolution of the contaminant transported in a heterogeneous aquifer is not affected by the shape of the initial plume but affected mainly by the degree of the heterogeneity and the size of the initial plume.