• Coupled systems mechanics
• /
• 제5권4호
• /
• pp.341-353
• /
• 2016

Manufactures of multi-crystalline silicon ingots by means of the directional solidification system (DSS) is important to the solar photovoltaic (PV) cell industry. The quality of the ingots, including the grain size and morphology, is highly related to the shape of the crystal-melt interface during the crystal growth process. We performed numerical simulations to analyze the thermo-fluid field and the shape of the crystal-melt interface both for steady conditions and transient processes. The steady simulations are first validated and then applied to improve the hot zone design in the furnace. The numerical results reveal that, an additional guiding plate weakens the strength of vortex and improves the desired profile of the crystal-melt interface. Based on the steady solutions at an early stage, detailed transient processes of crystal growth can be simulated. Accuracy of the results is supported by comparing the evolutions of crystal heights with the experimental measurements. The excellent agreements demonstrate the applicability of the present numerical methods in simulating a practical and complex system of directional solidification system.

• 천문학회보
• /
• 제39권1호
• /
• pp.37.2-37.2
• /
• 2014

Gas in disk galaxies interacts nonlinearly with a underlying stellar spiral potential to form galactic spiral shocks. Numerical simulations typically show that these shocks are unstable to the wiggle instability, forming non-axisymmetric structures with high vorticity. While previous studies suggested that the wiggle instability may arise from the Kelvin-Helmholtz instability or orbit crowding of gas elements near the shock, its physical nature remains uncertain. It was even argued that the wiggle instability is of numerical origin, caused by the inability of a numerical code to resolve a shock that is inclined to numerical grids. In this work, we perform a normal-mode linear stability analysis of galactic spiral shocks as a boundary-value problem. We find that the wiggle instability originates physically from the potential vorticity generation at a distorted shock front. As the gas follows galaxy rotation, it periodically passes through multiple shocks, successively increasing its potential vorticity. This sets up a normal-mode that grows exponentially, with a growth rate comparable to the orbital angular frequency. We show that the results of our linear stability analysis are in good agreement with the those of local hydrodynamic simulations of the wiggle instability.

• 천문학회지
• /
• 제36권3호
• /
• pp.111-121
• /
• 2003

In order to explore the cosmic ray acceleration at the cosmological shocks, we have performed numerical simulations of one-dimensional, plane-parallel, cosmic ray (CR) modified shocks with the newly developed CRASH (Cosmic Ray Amr SHock) numerical code. Based on the hypothesis that strong Alfven waves are self-generated by streaming CRs, the Bohm diffusion model for CRs is adopted. The code includes a plasma-physics-based 'injection' model that transfers a small proportion of the thermal proton flux through the shock into low energy CRs for acceleration there. We found that, for strong accretion shocks with Mach numbers greater than 10, CRs can absorb most of shock kinetic energy and the accretion shock speed is reduced up to $20\%$, compared to pure gas dynamic shocks. Although the amount of kinetic energy passed through accretion shocks is small, since they propagate into the low density intergalactic medium, they might possibly provide acceleration sites for ultra-high energy cosmic rays of $E\ll10^{18}eV$. For internal/merger shocks with Mach numbers less than 3, however, the energy transfer to CRs is only about $10-20\%$ and so nonlinear feedback due to the CR pressure is insignificant. Considering that intracluster medium (ICM) can be shocked repeatedly, however, the CRs generated by these weak shocks could be sufficient to explain the observed non-thermal signatures from clusters of galaxies.

• 대한조선학회논문집
• /
• 제55권3호
• /
• pp.265-273
• /
• 2018

Direct numerical and large eddy simulations of transitional flows around studs installed on flat plate and bulbous bow have been performed to investigate an effectiveness of turbulence stimulators on laminar-to-turbulence transition at a very low speed. The flow velocity was determined to be 0.366m/s corresponding to 4 knots of full-scale ship speed when the objective ship was Kriso container ship. The spatial evolution of skin friction coefficient disclosed that a fully development of turbulence was observed behind the second stud installed on flat plate while a rapid transition from laminar to turbulence gave rise to the fully development of turbulence behind the first stud installed on bulbous bow. A comparison of streamwise mean velocity profiles showed that the viscous sublayer and log-layer were in good agreement with previous results although the friction velocity of Smagrosinsky sub-grid scale model was about 10% larger than that of direct numerical simulation. While the turbulence intensities of bulbous bow was similar to those of flat plate in inner region, larger intensities of turbulence were observed in outer region of bulbous bow than those of flat plate.

• 한국안전학회지
• /
• 제25권2호
• /
• pp.12-17
• /
• 2010

Numerical simulations were performed for the prediction of the flame structure during the interaction between hydrogen and hydrocarbon flames. A counterflow flow geometry was introduced to establish the interacting two flames. Methane was used as a representative hydrocarbon fuel in this study. A well-known numerical code for the counterflow flame, OPPDIF, was used for the simulations. The detailed chemistry was adopted to predict the flame structure reasonably. The interaction of two one-dimensional premixed flames established in counterflow burner was investigated with the global strain rate and velocity ratio. It was found that the maximum temperature located near the methane flame surface while the heat release rate of methane was lower than hydrogen flame. The flame thickness become narrow with increasing the velocity ratio while the global strain rate was fixed. The local strain rate and heat release rate at the methane flame surface were correlated with the global strain rate, while those at the hydrogen flame were not correlated with the global strain rate. However, the maximum temperature of the interacting flames was correlated with the global strain rate.

• Computers and Concrete
• /
• 제25권2호
• /
• pp.181-192
• /
• 2020

Nowadays, steel fiber reinforced concretes (SFRCs) are widely used in practical applications. Significant experimental research has thus been carried out to determine the constitutive equations that represent the behavior of SFRCs under multiaxial loadings. However, numerical modelling of SFRCs via FEM has been challenging due to the complexities of the implementation of these constitutive equations. In this study, following the literature, a plasticity model is constructed for the behavior of SFRCs that involves the Willam-Warnke failure surface with the relevant evolution laws and a non-associated flow rule for determining the plastic deformations. For the precise (yet rapid) integration of the constitutive equations, an explicit substepping scheme consisting of yield intersection and drift correction algorithms is employed and thus implemented in ABAQUS via UMAT. The FEM model includes various material parameters that are determined from the experimental data. Three sets of parameters are used in the numerical simulations. While the first set is from the experiments that are conducted in this study on SFRC specimens with various contents of steel fibers, the other two sets are from the experiments reported in the literature. The response of SFRCs under multiaxial compression obtained from various numerical simulations are compared with the experimental data. The good agreement between numerical results and the experimental data indicates that not only the adopted plasticity model represents the behavior of SFRCs very well but also the implemented integration scheme can be employed in practical applications of SFRCs.

• Geomechanics and Engineering
• /
• 제16권5호
• /
• pp.539-545
• /
• 2018

Longitudinal Displacement Profile (LDP) is an appropriate tool for determination of the displacement magnitude of the tunnel walls as a function of the distance to the tunnel face. Some useful formulations for calculation of LDP have been developed based on the monitoring data on site or by 3D numerical simulations. However, the presented equations are only based on the tunnel dimensions and for different quality of rock masses proposed a unique LDP. In the present study, it is tried to present a new formulation, for calculation of LDP, on the basis of Rock mass quality. For this purpose, a comprehensive numerical simulation program was developed to investigate the effect of rock mass quality on the LDP. Results of the numerical modelling were analyzed and the least square technique was used for fitting an appropriate curve on the derived data from the numerical simulations. The proposed formulation in the present study, is a logistic function and the constants of the logistic function were predicted by rock mass quality index (GSI). Results of this study revealed that, the LDP curves of the tunnel surrounded by rock masses with high quality (GSI>60) match together; because the rock mass deformation varies over an elastic range.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1998년도 춘계학술발표회논문집(1)
• /
• pp.149-154
• /
• 1998

A subcritical reactor driven by a linear proton accelerator has been considered as a nuclear waste incinerator at Korea Atomic Energy Research Institute(KAERI). Since the multiplication factor of a subcritical reactor is less than unity, to compensate exponentially decreasing fission neutrons from spallation reactions are essentially required for operating the reactor in its steady state. furthermore, the profile of accelerator beam currents is very important in controlling a subcritical reactor, because the reactor power varies in accordance of the profile of external neutrons. We have developed a code system to find numerical solutions of reactor kinetics equations, which are the simplest dynamic model for controlling reactors. In a due course of our previous numerical study of point kinetics equations for critical reactors, however, we learned that the same code system can be used in studying dynamic behavior of the subcritical reactor. Our major motivation of this paper is to investigate responses of subcritical reactors for small changes in thermal hydraulic parameters. Building a thermal hydraulic model for the subcritical reactor dynamics, we performed numerical simulations for dynamic responses of the reactor based on point kinetics equations with a source term. Linearizing a set of coupled differential equations for reactor responses, we focus our research interest on dynamic responses of the reactor to variations of the thermal hydraulic parameters in transient phases.

• Computers and Concrete
• /
• 제22권4호
• /
• pp.395-405
• /
• 2018

In this paper, the performance of post-installed adhesive bonded anchor embedded in concrete is assessed using numerical simulations. This study aims at studying the influence of parameters on the performance of a chemically bonded anchorage system. Non-linear finite element modelling and simulations are carried out by properly using the material properties and phenomenon. Materials parameters such as characteristic length, fracture energy, damage criteria, tension retention and crack width of concrete and interface characteristics are carefully assigned so as to obtain a most realistic behaviour of the chemical anchor system. The peak strength of two different anchor systems obtained from present numerical studies is validated against experimental results. Furthermore, validated numerical models are used to study the load transferring mechanism and damage progression characteristics of various anchors systems where strength of concrete, strength of epoxy, and geometry and disposition of anchors are the parameters. The process of development of strain in concrete adjacent to the anchor and energy dissipated during the course of damage progression are analysed. Results show that the performance of the considered anchorage system is, though a combined effect of material and geometric parameters, but a clear distinction could be made on the parameters to achieve a desired performance based on strength, slip, strain development or dissipated energy. Inspite the increase in anchor capacity with increase in concrete strength, it brings some undesirable performance as well. Furthermore, the pullout capacity of the chemical anchor system increases with a decrease in disparity among the strength of concrete and epoxy.

• Structural Engineering and Mechanics
• /
• 제68권3호
• /
• pp.335-344
• /
• 2018

The execution of an experiment is a complex affair. It includes the preparation of test specimens, the measurement process itself and also the evaluation of the experiment as such. Financial requirements can differ significantly. In contrast, the cost of numerical simulations can be negligible, but what is the credibility of a simulated experiment? Discussions frequently arise concerning the methodology used in simulations, and particularly over the geometric model used. Simplification, rounding or the complete omission of details are frequent reasons for differences that occur between simulation results and the results of executed experiments. However, the creation of a very complex geometry, perhaps all the way down to the resolution of the very structure of the material, can be complicated. The subject of the article is therefore a means of creating the material structure of concrete contained in a test specimen. Because a complex approach is taken right from the very start of the numerical simulation, maximum agreement with experimental results can be achieved. With regard to the automation of the process described, countless material structures can be generated and randomly produced samples simulated in this way. Subsequently, a certain degree of randomness can be observed in the results obtained, e.g., the shape of the failure - just as is the case with experiments. The first part of the article presents a description of a complex approach to the creation of a geometry representing real concrete test specimens. The second part presents a practical application in which the numerical simulation of the compressive testing of concrete is executed using the generated geometry.