• 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

• Advances in aircraft and spacecraft science
• /
• v.5 no.3
• /
• pp.363-383
• /
• 2018

One-dimensional (1D) models of incompressible flows, can be of interest for many applications in which fast resolution times are demanded, such as fluid-structure interaction of flows in compliant pipes and hemodynamics. This work proposes a higher-order 1D theory for the flow-field analysis of incompressible, laminar, and viscous fluids in rigid pipes. This methodology is developed in the domain of the Carrera Unified Formulation (CUF), which was first employed in structural mechanics. In the framework of 1D modelling, CUF allows to express the primary variables (i.e., velocity and pressure fields in the case of incompressible flows) as arbitrary expansions of the generalized unknowns, which are functions of the 1D computational domain coordinate. As a consequence, the governing equations can be expressed in terms of fundamental nuclei, which are invariant of the theory approximation order. Several numerical examples are considered for validating this novel methodology, including simple Poiseuille flows in circular pipes and more complex velocity/pressure profiles of Stokes fluids into non-conventional computational domains. The attention is mainly focused on the use of hierarchical McLaurin polynomials as well as piece-wise nonlocal Lagrange expansions of the generalized unknowns across the pipe section. The preliminary results show the great advantages in terms of computational costs of the proposed method. Furthermore, they provide enough confidence for future extensions to more complex fluid-dynamics problems and fluid-structure interaction analysis.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.51 no.6
• /
• pp.53-62
• /
• 2009

The EFDC (Environmental Fluid Dynamics Code), a numerical model for simulating three-dimensional (3D) flow, transport, and biogeochemical processes in surface water systems including rivers, reservoirs, and estuaries, was applied to assess the effect of sea water and fresh water exchange rates ($Q_e$) on the mixing characteristics of a conservative pollutant (tracer) induced from upstreams and salinity in Saemangeum Lake, Korea. The lake has been closed by a 33 km estuary embankment since last April of 2006, and now seawater enters the lake partially through two sluice gates (Sinsi and Garyuk), which is driving the changes of hydrodynamic and water quality properties of the lake. The EFDC was constructed and calibrated with surveyed bathymetry data and field data including water level, temperature, and salinity in 2008. The model showed good agreement with the field data and adequately replicated the spatial and temporal variations of the variables. The validated model was applied to simulated the tracer and salinity with two different gate operation scenarios: RUN-1 and RUN-2. RUN-1 is the case of real operation condition ($Q_e=25,000,000\;m^3$) of 2008, while RUN-2 assumed full open of Sinsi gate to increase $Q_e$ by $120,000,000\;m^3$. Statistical analysis of the simulation results indicate that mixing characteristics of pollutants from upstream can be significantly affected by the amount of $Q_e$.

• Economic and Environmental Geology
• /
• v.29 no.3
• /
• pp.335-343
• /
• 1996

The dipole-dipol electrical resistivity survey was conducted to investigate the probable contamination of the Han river by leachate from the near-by Nanjido Landfill. The survey line of 3 km was set along the unpaved road toward the Han river. For the convenience of the field work, the survey line was divided into four segments. The complete two-dimensional resistivity section was constructed by connecting the inversion result of each segment. Gravity survey was also carried out along the profile parallel to the resistivity line. Near surface resistivity generally appeared to be of very low value in most part of the survey area and the boundary between the alluvium layer and underlying basement rocks is well discriminated on the resistivity section. These results agree well with those of the preceding Schlumberger depth sounding made at adjacent area by Lee and fun (1995). The variation of thickness of the alluvium layer delineated by gravity anomaly profile also correlates well with the result of the resistivity survey on the qualitative basis. The problem of contamination by leachate from the Nanjido Landfill, where various waste materials have been dumped without any proper treatment facilities, has been remains unsolved yet. Therefore, we present the most probable passages of leachate flow based on the survey results and have briefly discussed about measure for contamination control. Considering the thickness of alluvium and the possible existence of fractured zone, the middle point between 1st and 2nd landfill and the midst of 1st landfill are the most hazardous regions to make leachates flow into the Han river. Since large amounts of leachates are observed from the test wells located on the lines extending from the border between the 1st and 2nd landfill and the middle of the lst landfill, contamination protection barriers are strongly recommended near these regions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.97-106
• /
• 1989

This paper is on the modeling of two-dimensional groundwater flow, which is the first step of the development of Dynamic System Model for groundwater flow and pollutant transport in subsurface porous media. The particular features of the model are its versatility and flexibility to deal with as many real-world problems as possible. Points as well as distributed sources/sinks are included to represent recharges/pumping and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Sources/sinks strength over each element and node, hydraulic head at each Dirichlet boundary node and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution methed for the matrix equation approximating the partial differential equation of groundwater flow. The model also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. The groundwater flow model shall be combined with the model of pollutant transport in subsurface porous media. Then the combined model, with the applications of the Eigenvalue technique and the Dynamic system theory, shall be improved to the Dynamic System Model which can simulate the real groundwater flow and the pollutant transport accurately and effectively for the analyses and predictions.

• Clean Technology
• /
• v.27 no.1
• /
• pp.61-68
• /
• 2021

Selective catalytic reduction (SCR) is widely used as a method of removing nitrogen oxide in large-capacity thermal power generation systems. Uniform mixing of the injected ammonia and the inlet flue gas is very important to the performance of the denitrification reduction process in the catalyst bed. In the present study, a computational analysis technique was applied to the ammonia injection system design process of a denitrification facility. The applied model is the denitrification facility of an 800 MW class coal-fired power plant currently in operation. The flow field to be solved ranges from the inlet of the ammonia injection system to the end of the catalyst bed. The flow was analyzed in the two-dimensional domain assuming incompressible. The steady-state turbulent flow was solved with the commercial software named ANSYS-Fluent. The nozzle arrangement gap and injection flow rate in the ammonia injection system were chosen as the design parameters. A total of four (4) cases were simulated and compared. The root mean square of the NH3/NO molar ratio at the inlet of the catalyst layer was chosen as the optimization parameter and the design of the experiment was used as the base of the optimization algorithm. The case where the nozzle pitch and flow rate were adjusted at the same time was the best in terms of flow uniformity.

• Applied Biological Chemistry
• /
• v.43 no.2
• /
• pp.116-123
• /
• 2000

A numerical model of three-dimensional soil water distribution for drip irrigation management under cropped conditions was developed using Richards equation in Cartesian coordinates. The model accounts for both seasonal and diurnal changes in evaporation and transpiration, and the growth of plant root and the shape of root zone. Solutions were numerically approximated using the Crank-Nicolson implicit finite difference technique on the block-centered grid system and the Gauss-Seidel elimination in tandem. The model was tested under several conditions to allow the flow rates and configurations of drip emitters vary. In general, simulation results agreed well with experimental results and were as follows. The velocity of soil-water flow decreased drastically with distance from the drip source, and the rate of expansion of the wetted zone decreased rapidly during irrigation. The wetting front of wetted zone from a surface drip emitter traveled farther in vertical direction than in horizontal direction. Under this experimental weather condition, water use efficiency of a drip-irrigated apple field was greatest for 4-drip-emitter system buried at 25 cm, resulting from 10% increase in transpiration but 20% reduction in soil evaporation compared to those for surface 1-drip emitter system. Soil moisture retention curve obtained using disk tension infiltrometer showed significant difference from the curve obtained with pressure plate extractor.

• International Journal of Aeronautical and Space Sciences
• /
• v.15 no.4
• /
• pp.366-373
• /
• 2014

A conceptual design of a three-surface strike fighter was studied and stealth performance was taken into account to enhance survivability and battle effectiveness. CATIA was used to design the aircraft's three-dimensional prototype model and the weapon carriage arrangement was also studied. The aircraft's RCS characteristics and distributions under X, S, C, and L bands were simulated using the RCSPlus software, which is based on the PO method. Pressure and velocity distributions of the flow field were also simulated using CFD. A turbulence model was based on standard $k-{\varepsilon}$ function and N-S functions were used during the CFD computation. Lift coefficients, drag coefficients, and lift-to-drag ratio were obtained by aerodynamic simulation. The results showed that: (1) the average value of head-on RCS between ${\pm}30^{\circ}$ is below -3.197 dBsm, and (2) the lift coefficient is 0.34674, the drag coefficient is 0.04275, and the lift-to-drag ratio is 8.11087 when the attack angle is $2.5^{\circ}$.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.56 no.1
• /
• pp.81-88
• /
• 2014

Accurate wind data is essential for predicting airborne spread of virus. OpenFOAM was used for computational fluid dynamics (CFD) simulation procedure which is under GNU GPL (General Public License). Using complex terrain, DEM (Digital Elevation Map) that was prepared from GIS information covering a research site is converted to a three dimensional surface mesh that is composed by quad and full hexahedral space meshes. Around this surface mesh, an extended computational domain volume was designed. Atmospheric flow boundary conditions were used at inlet and roughness height and was considered at terrain by using rough wall function. Two different wind conditions that was relatively stable during certain periods were compared in 3 different locations for validating the accuracy of the CFD computed solution. The result shows about 10 % of difference between the calculated result and measured data. This procedure can simulate a prediction of time-series data for airborne virus spread that can be used to make a web-based forecasting system of airborne virus spread.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.7 no.3 s.18
• /
• pp.21-29
• /
• 2004

Numerical analysis of aerodynamic heating for KPSAM is performed using aerodynamic heating model suitable to KPSAM, which has complex flow field resulting from the spike attached to the dome, such as large separation area and the strong shock/boundary layer interaction region around reattachment point on the dome. The aerodynamic heating model is validated and modified through the comparison between the flight test measurement and the thermal analysis results. TFD temperature sensors are installed on the dome to measure surface temperature during the flight. Computation results, obtained from the heat transfer analysis on the sensors, agree well with flight test data. The aerodynamic heating model provides heat transfer rate into surface as a boundary condition of unsteady 1D/axisymmetric thermal analysis on the missile structure. The axisymmetric thermal analysis using FLUENT is more versatile than the 1D analysis and can be applied to the heating problem related with complex structures and multi-dimensional heat transfer problems such as prediction of temperature rise at contact surface of different materials.