• Nuclear Engineering and Technology
• /
• v.52 no.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
• /
• v.53 no.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.06a
• /
• 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.

• Proceedings of the KIEE Conference
• /
• 1998.07a
• /
• 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.

• Membrane Journal
• /
• v.12 no.4
• /
• pp.207-215
• /
• 2002

The hindered diffusion in confined spaces is an important phenomenon to understand in a micro-scale the filtration mechanism determined by the particle motion in membrane pores. Compared to the case of spherical colloids, both the theoretical investigations and the experiments on the hindered diffusion of polyelectrolytes is actually more difficult, due to lots of relevant parameters resulting from the complicated conformational properties of the polyelectrolyte chain. We have successfully performed the Brownian dynamics simulations upon a single polyeiectrolyte confined in a slit-like pore, where a coarse-grained bead-spring model incorporated with Debye-Huckel interaction is properly adopted. For the given sizes of both the polyelectrolyte and the pore width, the hindered diffusion coefficient decreases as the solution ionic concentration decreases. It is evident that a charge effect of the pore wall enhances the hindered diffusion of polyelectrolyte. Simulation results allow us to make sense of the diffusive transport through the micro-pore, which is restricted by the influences of the steric hindrance of polyelectrolytes as well as the electrostatic repulsion between the polyelectrolytes and pore wall.

• Nuclear Engineering and Technology
• /
• v.41 no.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.

• Membrane Journal
• /
• v.9 no.1
• /
• pp.51-62
• /
• 1999

The recovery of trace volatile organic components from water by pervaporation was investigated. Permeation experiments through homogeneous polydimethylsiloxane(PDMS) membrane was carried out and the effect of feed concentrations and membrane thicknesses on the permeation characteristics were investigated. A solution-diffusion model is used to describe the pervaporation transport mechanism. In homogeneous PDMS membrane it appeared that the selectivities of MEK and toluene are constant, and that organic flux has a linear relationship with feed concentration. These results indicate that the coupling effects between organics were negligible. The selectivity of PDMS membranes is invariant with respect to the membrane thickness. The intrinsic membrane permeability of organic components determined by using a solution-diffusion model. Comparing with various composite type membrane, the membrane using PEG treated nonwoven fabric as sublayer showed the best performance in VOC recovery by pervaporation.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.11
• /
• pp.4664-4670
• /
• 2010

Comparison of the classical mobile-immobile zone (MIM) model to the derived model led to several conclusions. If the MIM model is to be applied, the initial concentration in the immobile zone has to be down-scaled by a correction factor that is a function of pore geometry. The MIM model was valid only after sufficiently long time has passed, i.e., only after the diffusion front reaches the deepest pore wall in the immobile zone. The MIM mass-transfer coefficient $\alpha$, was inversely proportional to the square of the pore depth. Also it did not depend on the mobile-zone flow velocity, contrary to the number of laboratory and field observations. The classical MIM model displayed a rapid exponential decay of immobile-zone concentration. Meanwhile at large times, the newly derived model displayed similar exponential decay. This was contrary to the mounting evidence of power-law BTC tails observed in laboratory and field settings.

• The Bulletin of The Korean Astronomical Society
• /
• v.37 no.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
• /
• v.6 no.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.