Some kinds of break in the reactor coolant system may cause the coolant to exit rapidly from the failure site,which leads to the loss of coolant accident(LOCA).In this paper,a stress analysis of an AP1000 reactor cont...Some kinds of break in the reactor coolant system may cause the coolant to exit rapidly from the failure site,which leads to the loss of coolant accident(LOCA).In this paper,a stress analysis of an AP1000 reactor containment is performed in an LOCA,with the passive containment cooling system(PCCS) being available and not available for cooling the wall's containment.The variations in the mechanical properties of the wall's containment,including elastic modulus,strength,and stress,are analyzed using the ABAQUS code.A general two-phase model is applied for modeling thermal-hydraulic behavior inside the containment.Obtained pressure and temperature from thermal-hydraulic models are considered as boundary conditions of the ABAQUS code to obtain distributions of temperature and stress across steel shell of the containment in the accident.The results indicate that if the PCCS fails,the peak pressure inside the containment exceeds the design value.However,the stress would still be lower than the yield stress value,and no risk would threaten the integrity of the containment.展开更多
The highest thermal-hydraulic pressure in the containment occurs when reactor coolant in the first loop and steam in the secondary loop discharge simultaneously,and when the maximum amount of energy from reactor unit ...The highest thermal-hydraulic pressure in the containment occurs when reactor coolant in the first loop and steam in the secondary loop discharge simultaneously,and when the maximum amount of energy from reactor unit enters to containment volume. In this paper, we investigate temperature and pressure variations in the VVER 1000 containment compartments owing to concurrent break in the pipelines of the primary and secondary loops. A two-phase, multicellular model is applied in the presence of non-condensable gases. Convection and conduction through the main heat structures inside the containment are also considered. The predicted results agree well with available data. Maximum values of pressure and temperature in the containment are then calculated and compared to the design values. If LOCA and MSLB occur simultaneously, the maximum pressure would exceed the design value and integrity of the containment would be threatened.展开更多
The present paper discusses entropy generation in fully developed turbulent flows through a subchannel,arranged in square and triangle arrays. Entropy generation is due to contribution of both heat transfer and pressu...The present paper discusses entropy generation in fully developed turbulent flows through a subchannel,arranged in square and triangle arrays. Entropy generation is due to contribution of both heat transfer and pressure drop. Our main objective is to study the effect of key parameters such as spacer grid, fuel rod power distribution,Reynolds number Re, dimensionless heat power ω, lengthto-fuel-diameter ratio λ, and pitch-to-diameter ratio ξ on subchannel entropy generation. The analysis explicitly shows the contribution of heat transfer and pressure drop to the total entropy generation. An analytical formulation is introduced to total entropy generation for situations with uniform and sinusoidal rod power distribution. It is concluded that power distribution affects entropy generation.A smoother power profile leads to less entropy generation.The entropy generation of square rod array bundles is more efficient than that of triangular rod arrays, and spacer grids generate more entropy.展开更多
Optical features of a semiconductor–dielectric photonic crystal are studied theoretically. Alternating layers of micrometer sized SiO2/In Sb slabs are considered as building blocks of the proposed ideal crystal. By i...Optical features of a semiconductor–dielectric photonic crystal are studied theoretically. Alternating layers of micrometer sized SiO2/In Sb slabs are considered as building blocks of the proposed ideal crystal. By inserting additional layers and disrupting the regularity, two more defective crystals are also proposed. Photonic band structure of the ideal crystal and its dependence on the structural parameters are explored at the first step. Transmittance of the defective crystals and its changes with the thicknesses of the layers are studied. After extracting the optimum values for the thicknesses of the unit cells of the crystals, the optical response of the proposed structures at different temperatures and incident angles are investigated. Changes of the defect layers’ induced mode(s) are discussed by taking into consideration of the temperature dependence of the In Sb layer permittivity. The results clearly reflect the high potential of the proposed crystals to be used at high temperature terahertz technology as a promising alternative to their electronic counterparts.展开更多
The drift-flux model has a practical importance in two-phase flow analysis.In this study,a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh,leading to the d...The drift-flux model has a practical importance in two-phase flow analysis.In this study,a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh,leading to the development of a fully implicit discretization method.The main advantage of the fully implicit method is its unconditional stability.Newton's scheme is a popular method of choice for the solution of a nonlinear system of equations arising from fully implicit discretization of field equations.However,the lack of convergence robustness and the construction of Jacobian matrix have created several difficulties for the researchers.In this paper,a fully implicit model is developed based on the SIMPLE algorithm for two-phase flow simulations.The drawbacks of Newton's method are avoided in the developed model.Different limiter functions are considered,and the stabilized method is developed under steady and transient conditions.The results obtained by the numerical modeling are in good agreement with the experimental data.As expected,the results prove that the developed model is not restricted by any stability limit.展开更多
Measurement of dose distribution in patients during radiotherapy is impossible. The Monte Carlo simulation is an alternative method for dose calculations. In routine radiotherapy, the source-to-surface distance(SSD)me...Measurement of dose distribution in patients during radiotherapy is impossible. The Monte Carlo simulation is an alternative method for dose calculations. In routine radiotherapy, the source-to-surface distance(SSD)method is not practical for an isocentric unit because it requires numerous values of tissue–air ratios and inverse square law. Therefore, this method is time consuming. In this paper, the curves of relative depth doses were obtained for three different SSDs using the MCNP4C Monte Carlo simulation and approximated with a single curve called calibration curve. This curve was compared to the curve obtained by published data, differing in approximately 5% in the worst case. It was also observed that the obtained results were more accurate for distances between-5 and 10 cm from source-to-axis distance.展开更多
This paper deals with the investigation of the power absorption in helicon plasma excited through a half-helix antenna driven at 13.56 MHz. The simulations were carried out by means of a code,HELIC. They were carried ...This paper deals with the investigation of the power absorption in helicon plasma excited through a half-helix antenna driven at 13.56 MHz. The simulations were carried out by means of a code,HELIC. They were carried out by taking into account different inhomogeneous radial density profiles and for a wide range of plasma densities, from 10^(11) cm^(-3) to 10^(13) cm^(-3). The magnetic field was 200, 400, 600 and 1000 G. A three-parameter function was used for generating various density profiles with different volume gradients, edge gradients and density widths. The density profile had a large effect on the efficient Trivelpiece–Gould(TG) and helicon mode excitation and antenna coupling to the plasma. The fraction of power deposition via the TG mode was extremely dependent on the plasma density near the plasma boundary. Interestingly, the obtained efficient parallel helicon wavelength was close to the anticipated value for Gaussian radial density profile.Power deposition was considerably asymmetric when the n/B_0 ratio was more than a specific value for a determined density width. The longitudinal power absorption was symmetric at approximately n_0 =10^(11) cm^(-3), irrespective of the magnetic field supposed. The asymmetry became more pronounced when the plasma density was 10^(12) cm^(-3). The ratio of density width to the magnetic field was an important parameter in the power coupling. At high magnetic fields, the maximum of the power absorption was reached at higher plasma density widths. There was at least one combination of the plasma density, magnetic field and density width for which the RF power deposition at both side of the tube reached its maximum value.展开更多
Multiphase flow meters as the potential alternatives to separation and metering techniques have been in rapid development since 1980 s.Before its field operation,the instrument should be calibrated in a standard test-...Multiphase flow meters as the potential alternatives to separation and metering techniques have been in rapid development since 1980 s.Before its field operation,the instrument should be calibrated in a standard test-facility.In spite of the known medium and large scale facilities all over the world,we developed a laboratory scale instrument for component fraction measurements.It has a two-phase flow homogenizer loop with the clamp-on potential of the meters to provide a regime independent measurement.It is capable of delivering a complete homogenization by γ-ray densitometer.With an error of±5%in component fraction measurements,this instrument is appropriate for testing and calibrating other meters.展开更多
The propagation of an intense laser pulse in an under-dense plasma induces a plasma wake that is suitable for the acceleration of electrons to relativistic energies. For an ultra-intense laser pulse which has a longit...The propagation of an intense laser pulse in an under-dense plasma induces a plasma wake that is suitable for the acceleration of electrons to relativistic energies. For an ultra-intense laser pulse which has a longitudinal size shorter than the plasma wavelength, λp, instead of a periodic plasma wave, a cavity free from cold plasma electrons, called a bubble, is formed behind the laser pulse. An intense charge separation electric field inside the moving bubble can capture the electrons at the base of the bubble and accelerate them with a narrow energy spread. In the nonlinear bubble regime, due to localized depletion at the front of the pulse during its propagation through the plasma, the phase shift between carrier waves and pulse envelope plays an important role in plasma response. The carrier–envelope phase(CEP) breaks down the symmetric transverse ponderomotive force of the laser pulse that makes the bubble structure unstable. Our studies using a series of two-dimensional(2D) particle-in-cell(PIC) simulations show that the frequency-chirped laser pulses are more effective in controlling the pulse depletion rate and consequently the effect of the CEP in the bubble regime. The results indicate that the utilization of a positively chirped laser pulse leads to an increase in rate of erosion of the leading edge of the pulse that rapidly results in the formation of a steep intensity gradient at the front of the pulse. A more unstable bubble structure, the self-injections in different positions, and high dark current are the results of using a positively chirped laser pulse. For a negatively chirped laser pulse, the pulse depletion process is compensated during the propagation of the pulse in plasma in such a way that results in a more stable bubble shape and therefore, a localized electron bunch is produced during the acceleration process. As a result, by the proper choice of chirping, one can tune the number of self-injected electrons, the size of accelerated bunch and its energy spectrum to the values required for practical applications.展开更多
Using the Kane-Mele Hamiltonian, Dirac theory and self-consistent Born approximation, we investigate the effect of dilute charged impurity on the electronic heat capacity and magnetic susceptibility of two-dimensional...Using the Kane-Mele Hamiltonian, Dirac theory and self-consistent Born approximation, we investigate the effect of dilute charged impurity on the electronic heat capacity and magnetic susceptibility of two-dimensional ferromagnetic honeycomb structure of group-Ⅳ elements including silicene, germanene and stanene within the Green's function approach. We also find these quantities in the presence of applied external electric field. Our results show that the silicene(stanene) has the maximum(minimum) heat capacity and magnetic susceptibility at uniform electric fields. From the behavior of theses quantities, the band gap has been changed with impurity concentration, impurity scattering strength and electric field. The analysis on the impurity-dependent magnetic susceptibility curves shows a phase transition from ferromagnetic to paramagnetic and antiferromagnetic phases. Interestingly, electronic heat capacity increases(decreases) with impurity concentration in silicene(germanene and stanene) structure.展开更多
文摘Some kinds of break in the reactor coolant system may cause the coolant to exit rapidly from the failure site,which leads to the loss of coolant accident(LOCA).In this paper,a stress analysis of an AP1000 reactor containment is performed in an LOCA,with the passive containment cooling system(PCCS) being available and not available for cooling the wall's containment.The variations in the mechanical properties of the wall's containment,including elastic modulus,strength,and stress,are analyzed using the ABAQUS code.A general two-phase model is applied for modeling thermal-hydraulic behavior inside the containment.Obtained pressure and temperature from thermal-hydraulic models are considered as boundary conditions of the ABAQUS code to obtain distributions of temperature and stress across steel shell of the containment in the accident.The results indicate that if the PCCS fails,the peak pressure inside the containment exceeds the design value.However,the stress would still be lower than the yield stress value,and no risk would threaten the integrity of the containment.
文摘The highest thermal-hydraulic pressure in the containment occurs when reactor coolant in the first loop and steam in the secondary loop discharge simultaneously,and when the maximum amount of energy from reactor unit enters to containment volume. In this paper, we investigate temperature and pressure variations in the VVER 1000 containment compartments owing to concurrent break in the pipelines of the primary and secondary loops. A two-phase, multicellular model is applied in the presence of non-condensable gases. Convection and conduction through the main heat structures inside the containment are also considered. The predicted results agree well with available data. Maximum values of pressure and temperature in the containment are then calculated and compared to the design values. If LOCA and MSLB occur simultaneously, the maximum pressure would exceed the design value and integrity of the containment would be threatened.
文摘The present paper discusses entropy generation in fully developed turbulent flows through a subchannel,arranged in square and triangle arrays. Entropy generation is due to contribution of both heat transfer and pressure drop. Our main objective is to study the effect of key parameters such as spacer grid, fuel rod power distribution,Reynolds number Re, dimensionless heat power ω, lengthto-fuel-diameter ratio λ, and pitch-to-diameter ratio ξ on subchannel entropy generation. The analysis explicitly shows the contribution of heat transfer and pressure drop to the total entropy generation. An analytical formulation is introduced to total entropy generation for situations with uniform and sinusoidal rod power distribution. It is concluded that power distribution affects entropy generation.A smoother power profile leads to less entropy generation.The entropy generation of square rod array bundles is more efficient than that of triangular rod arrays, and spacer grids generate more entropy.
文摘Optical features of a semiconductor–dielectric photonic crystal are studied theoretically. Alternating layers of micrometer sized SiO2/In Sb slabs are considered as building blocks of the proposed ideal crystal. By inserting additional layers and disrupting the regularity, two more defective crystals are also proposed. Photonic band structure of the ideal crystal and its dependence on the structural parameters are explored at the first step. Transmittance of the defective crystals and its changes with the thicknesses of the layers are studied. After extracting the optimum values for the thicknesses of the unit cells of the crystals, the optical response of the proposed structures at different temperatures and incident angles are investigated. Changes of the defect layers’ induced mode(s) are discussed by taking into consideration of the temperature dependence of the In Sb layer permittivity. The results clearly reflect the high potential of the proposed crystals to be used at high temperature terahertz technology as a promising alternative to their electronic counterparts.
文摘The drift-flux model has a practical importance in two-phase flow analysis.In this study,a finite volume solution is developed for a transient four-equation drift-flux model through the staggered mesh,leading to the development of a fully implicit discretization method.The main advantage of the fully implicit method is its unconditional stability.Newton's scheme is a popular method of choice for the solution of a nonlinear system of equations arising from fully implicit discretization of field equations.However,the lack of convergence robustness and the construction of Jacobian matrix have created several difficulties for the researchers.In this paper,a fully implicit model is developed based on the SIMPLE algorithm for two-phase flow simulations.The drawbacks of Newton's method are avoided in the developed model.Different limiter functions are considered,and the stabilized method is developed under steady and transient conditions.The results obtained by the numerical modeling are in good agreement with the experimental data.As expected,the results prove that the developed model is not restricted by any stability limit.
文摘Measurement of dose distribution in patients during radiotherapy is impossible. The Monte Carlo simulation is an alternative method for dose calculations. In routine radiotherapy, the source-to-surface distance(SSD)method is not practical for an isocentric unit because it requires numerous values of tissue–air ratios and inverse square law. Therefore, this method is time consuming. In this paper, the curves of relative depth doses were obtained for three different SSDs using the MCNP4C Monte Carlo simulation and approximated with a single curve called calibration curve. This curve was compared to the curve obtained by published data, differing in approximately 5% in the worst case. It was also observed that the obtained results were more accurate for distances between-5 and 10 cm from source-to-axis distance.
文摘This paper deals with the investigation of the power absorption in helicon plasma excited through a half-helix antenna driven at 13.56 MHz. The simulations were carried out by means of a code,HELIC. They were carried out by taking into account different inhomogeneous radial density profiles and for a wide range of plasma densities, from 10^(11) cm^(-3) to 10^(13) cm^(-3). The magnetic field was 200, 400, 600 and 1000 G. A three-parameter function was used for generating various density profiles with different volume gradients, edge gradients and density widths. The density profile had a large effect on the efficient Trivelpiece–Gould(TG) and helicon mode excitation and antenna coupling to the plasma. The fraction of power deposition via the TG mode was extremely dependent on the plasma density near the plasma boundary. Interestingly, the obtained efficient parallel helicon wavelength was close to the anticipated value for Gaussian radial density profile.Power deposition was considerably asymmetric when the n/B_0 ratio was more than a specific value for a determined density width. The longitudinal power absorption was symmetric at approximately n_0 =10^(11) cm^(-3), irrespective of the magnetic field supposed. The asymmetry became more pronounced when the plasma density was 10^(12) cm^(-3). The ratio of density width to the magnetic field was an important parameter in the power coupling. At high magnetic fields, the maximum of the power absorption was reached at higher plasma density widths. There was at least one combination of the plasma density, magnetic field and density width for which the RF power deposition at both side of the tube reached its maximum value.
基金financially supported by Radiation Application Research School as a part of Nuclear Science and Technology Research Institute
文摘Multiphase flow meters as the potential alternatives to separation and metering techniques have been in rapid development since 1980 s.Before its field operation,the instrument should be calibrated in a standard test-facility.In spite of the known medium and large scale facilities all over the world,we developed a laboratory scale instrument for component fraction measurements.It has a two-phase flow homogenizer loop with the clamp-on potential of the meters to provide a regime independent measurement.It is capable of delivering a complete homogenization by γ-ray densitometer.With an error of±5%in component fraction measurements,this instrument is appropriate for testing and calibrating other meters.
文摘The propagation of an intense laser pulse in an under-dense plasma induces a plasma wake that is suitable for the acceleration of electrons to relativistic energies. For an ultra-intense laser pulse which has a longitudinal size shorter than the plasma wavelength, λp, instead of a periodic plasma wave, a cavity free from cold plasma electrons, called a bubble, is formed behind the laser pulse. An intense charge separation electric field inside the moving bubble can capture the electrons at the base of the bubble and accelerate them with a narrow energy spread. In the nonlinear bubble regime, due to localized depletion at the front of the pulse during its propagation through the plasma, the phase shift between carrier waves and pulse envelope plays an important role in plasma response. The carrier–envelope phase(CEP) breaks down the symmetric transverse ponderomotive force of the laser pulse that makes the bubble structure unstable. Our studies using a series of two-dimensional(2D) particle-in-cell(PIC) simulations show that the frequency-chirped laser pulses are more effective in controlling the pulse depletion rate and consequently the effect of the CEP in the bubble regime. The results indicate that the utilization of a positively chirped laser pulse leads to an increase in rate of erosion of the leading edge of the pulse that rapidly results in the formation of a steep intensity gradient at the front of the pulse. A more unstable bubble structure, the self-injections in different positions, and high dark current are the results of using a positively chirped laser pulse. For a negatively chirped laser pulse, the pulse depletion process is compensated during the propagation of the pulse in plasma in such a way that results in a more stable bubble shape and therefore, a localized electron bunch is produced during the acceleration process. As a result, by the proper choice of chirping, one can tune the number of self-injected electrons, the size of accelerated bunch and its energy spectrum to the values required for practical applications.
文摘Using the Kane-Mele Hamiltonian, Dirac theory and self-consistent Born approximation, we investigate the effect of dilute charged impurity on the electronic heat capacity and magnetic susceptibility of two-dimensional ferromagnetic honeycomb structure of group-Ⅳ elements including silicene, germanene and stanene within the Green's function approach. We also find these quantities in the presence of applied external electric field. Our results show that the silicene(stanene) has the maximum(minimum) heat capacity and magnetic susceptibility at uniform electric fields. From the behavior of theses quantities, the band gap has been changed with impurity concentration, impurity scattering strength and electric field. The analysis on the impurity-dependent magnetic susceptibility curves shows a phase transition from ferromagnetic to paramagnetic and antiferromagnetic phases. Interestingly, electronic heat capacity increases(decreases) with impurity concentration in silicene(germanene and stanene) structure.
