• Nuclear Engineering and Technology
• /
• 제52권4호
• /
• pp.819-826
• /
• 2020

The method of the Laplace transform has been used to obtain an analytical solution of the three-dimensional steady state advection diffusion equation for the airborne radionuclide release from any nuclear installation such as the power reactor in an area source. The present treatment takes into account the removal of the pollutants through the nuclear reaction. We assume that the pollutants are emitted as a constant rate from the area source. This physical consideration is achieved by assuming that the vertical eddy diffusivity coefficient should be a constant. The prevailing wind speed is a constant in 𝑥- direction and a linear function of the vertical height z. The present model calculations are compared with the other models and the available data of the atmospheric dispersion experiments that were carried out in the nuclear power plant of Angra dos Reis (Brazil). The results show that the present treatment performs well as the analytical dispersion model and there is a good agreement between the values computed by our model and the observed data.

• Nuclear Engineering and Technology
• /
• 제53권3호
• /
• pp.833-841
• /
• 2021

During the postulated severe accident of nuclear reactor, eutectic reaction leads to low-temperature melting of fuel cladding and early failure of core structure. In order to model eutectic melting with the moving particle semi-implicit (MPS) method, the eutectic reaction model is developed to simulate the eutectic reaction phenomenon. The coupling of mass diffusion and phase diagram is applied to calculate the eutectic reaction with the uniform temperature. A heat transfer formula is proposed based on the phase diagram to handle the heat release or absorption during the process of eutectic reaction, and it can combine with mass diffusion and phase diagram to describe the eutectic reaction with temperature variation. The heat transfer formula is verified by the one-dimensional melting simulations and the predicted interface position agrees well with the theoretical solution. In order to verify the eutectic reaction models, the eutectic reaction of uranium and iron in two semi-infinite domains is simulated, and the profile of solid thickness decrease over time follows the parabolic law. The modified MPS method is applied to calculate Transient Reactor Test Facility (TREAT) experiment, the penetration rate in the simulations are agreeable with the experiment results. In addition, a hypothetical case based on the TREAT experiment is also conducted to validate the eutectic reaction with temperature variation, the results present continuity with the simulations of TREAT experiment. Thus the improved method is proved to be capable of simulating the eutectic reaction in the severe accident.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 1997년도 제15회 KOSCO SYMPOSIUM 논문집
• /
• pp.113-122
• /
• 1997

The conditional moment closure formulation considering the molecular and turbulent diffusion is derived. A simplified solution procedure is proposed to reduce the computational burden due to the increased dimensionality of the conditionally averaged variables. A conditionally averaged variable is expressed as a linear weighted average of the two extremes, 'no reaction' and 'equilibrium' states. The modified elliptic-type conditional moment closure formulation is implemented to simulate a two dimensional nonpremixed mixing layer reacting flow. Results show good agreement for the conditional averages of the species concentration in Bilger et al.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1998년도 하계학술대회 논문집 A
• /
• pp.113-115
• /
• 1998

This paper presents the indirect boundary integral equation method(BIEM) to analyze 3D moving conductor problem. Instead of an artificial upwind algothm, the proposed method uses a fundamental Green's function which is a particular solution of diffusion equation. Therefore, this method yields a stable and accurate solution regardless of the Peclet number. The indirect BIEM is compared with 3D upwind FEM for a numerical model which has analytic solutions.

• 멤브레인
• /
• 제12권4호
• /
• pp.207-215
• /
• 2002

한정된 미세공간에서의 제한확산(hindered diffusion)은 멤브레인 기공(pore)에서 입자들의 운동에 의해 결정되는 여과 메카니즘을 매우 미세한 수준에서 이해하는데 중요한 현상이다. 구형(spherical) 콜로이드 입자에 비해 보다 복잡한 형태(conformation)인 고분자사슬 구조를 갖는 다가전해질(polyelectrolyte)의 제한확산 거동에는 다양한 인자들이 관련되어 있기 때문에, 이론 접근은 물론 실험적 접근도 한층 어려운 것이 사실이다. 본 연구에서는, 슬릿형 미세기공에 한정되어 있는 단일한 다가전해질(single polyelectrolyte)에 coarse-grained bead spring model과 먼거리(long-range) 정전상호작용(electrostatic interaction)인 Debye-Huckel potential을 적용하여 분자시뮬레이션 기법인 브라운 동력학 모사를 수행하였다. 기공과 다가전해질 사슬(Polyelectrolyte chain)의 주어진 크기에서, 용액의 전해질 이온농도가 감소함에 따른 사슬의 신장(extension)효과는 제한확산계수를 감소시켰고, 기공 벽면의 하전성은 제한확산계수를 더욱 감소시켰다. 이는, 다가전해질 사슬(polyelectrolyte chain)의 입체적 장애(steric hindrance)와 함께 정전반발력이 미세기공에서의 확산이동을 억제함을 의미한다.

• Nuclear Engineering and Technology
• /
• 제41권7호
• /
• pp.885-892
• /
• 2009

This manuscript will discuss a numerical method where the six equations of two-phase flow, the solid heat conduction equations, and the two equations that describe neutron diffusion and precursor concentration are solved together in a tightly coupled, nonlinear fashion for a simplified model of a nuclear reactor core. This approach has two important advantages. The first advantage is a higher level of accuracy. Because the equations are solved together in a single nonlinear system, the solution is more accurate than the traditional "operator split" approach where the two-phase flow equations are solved first, the heat conduction is solved second and the neutron diffusion is solved third, limiting the temporal accuracy to $1^{st}$ order because the nonlinear coupling between the physics is handled explicitly. The second advantage of the method described in this manuscript is that the time step control in the fully implicit system can be based on the timescale of the solution rather than a stability-based time step restriction like the material Courant limit required of operator-split methods. In this work, a pilot code was used which employs this tightly coupled, fully implicit method to simulate a reactor core. Results are presented from a simulated control rod movement which show $2^{nd}$ order accuracy in time. Also described in this paper is a simulated rod ejection demonstrating how the fastest timescale of the problem can change between the state variables of neutronics, conduction and two-phase flow during the course of a transient.

• 멤브레인
• /
• 제9권1호
• /
• pp.51-62
• /
• 1999

본 연구는 투과증발법을 이용한 미량의 휘발성 유기용제 회수에 관한 것으로서 polydimethylsiloxane (PDMS) 균질막을 이용하여 혼합물의 농도, 막의 두께에 따른 투과 특성을 살펴보고, 이에 따라 용해-확산 모델을 이용하여 해석하였다. PDMS 균질막을 통한 MEK , 톨루엔 혼합수용액의 투과실험에서 각 물질의 선택도는 공급액의 농도와 다른 물질의 농도에 상관없이 일정한 값을 나타내었으며, 투과유량도 공급액 농도와 선형의 비례관계을 보였다. 따라서 본 연구에서 각 물질간의 간섭효과는 나타나지 않았다. 균질막의 두께를 변화시키며 행한 투과 실험에서 막 두께에 다른 선택도의 변화는 관철되지 않았다. 이에 따라서 용해-확산모델을 이용하여 각 물질의 PDMS막에 대한 투과계수를 구하였다. 막의 기능성과 실용성을 향상시키기 위한 복합막을 제조하여 행한 투과 실험을 한 결과 PEG 처리한 부직포위에 PDMS를 도포한 막이 선택도나 용질투과유량면에서 가장 우수한 성능을 보였다.

• 한국산학기술학회논문지
• /
• 제11권11호
• /
• pp.4664-4670
• /
• 2010

본 연구에서는 기존 MIM Zone model과 새로이 유도한 모델과 비교를 통해 몇가지 결론을 도출하였다. MIM model 이 적용되면 immobile 영역에서의 초기농도는 기공의 형태에 의해 달라지는 보정계수에 의해 실제농도보다 저평가되었으며 이는 오직 확산이 기공의 가장 깊은 부분까지 시행된 이후에 유효함을 확인하였다. 물질이동계수, $\alpha$는 기공의 깊이에 따라 반비례하며, 유동 구역의 유속에는 종속되어지지 않는다. 기존의 MIM model 은 확산주도영역의 농도가 급속하게 감소하는 현상을 보여주는데 새로이 유도된 모델의 경우도 충분한 시간이 경과한후 비슷한 현상을 보였으며 이는 기존의 여러 실험에서 관찰된 power-law BTC 의 상반되는 결과를 보여준다.

• 천문학회보
• /
• 제37권2호
• /
• pp.134.1-134.1
• /
• 2012

The radiation belt structure can be approximately reproduced by a form of diffusion equation, which takes into account the radial diffusion process as well as those in pitch angle and energy. The solution of the equation depends on several factors including initial and boundary conditions, diffusion coefficients, and plasmapause location. In this paper, we have attempted to determine a set of approximate functions for the energetic electron fluxes near the outer edge of the outer belt in terms of solar wind variable. We used the electron flux data from SST onboard the THEMIS spacecraft and determined its correlation with solar wind conditions in a systematic way. The functions were determined separately for different energy channels from ~30 keV up to 719 keV. Our determination of these functions allows us to predict the radial boundary condition for the electron flux, which can be implemented in a forecast model.

• Bulletin of the Korean Chemical Society
• /
• 제6권3호
• /
• pp.128-131
• /
• 1985

A series of copolymer membranes of 2-hydroxyethylmethacrylate (HEMA) with selected hydrophobic monomers were prepared without crosslinking agents. The equilibrium water content, the partition coefficient, and the permeability of the solutes such as urea, methylurea, 1,3-di-methylurea, and acetamide via these membranes were measured. The partition coefficient data show that as the hydrophobicity of solutes increased, the partition of solutes were dictated by hydrophobic interaction between solute and polymer matrix. Diffusion coefficients obtained in these experiments decrease as the water content of polymer membrane decreases. This decrease is blunt as the excess heat capacities, ${\phi}C^0_p$ (excess) in aqueous solution at infinite dilution of solute increases. To investigate the relationship between water content and diffusion coefficient, the results of the diffusion experiments were examined in light of a free-volume model of diffusive transport. The remarkable increase of urea mobility in the polymer network containing relatively larger bulk water can be considered as water structure breaking effect.