• Energy Research
• Closed Access close
• JP close
• IN close
• University of Tokyo close

The concept of a high temperature fast reactor cooled by supercritical water (SCFR-H) was developed for achieving high thermal efficiency and a compact reactor system. The core characteristics were obtained from single channel thermal-hydraulic analysis. Thus, it is necessary to carry out subchannel analysis to estimate the effect of local power peaking and cross flows. For this purpose, a subchannel analysis code is developed. It is verified by comparing the results with experimental data of High Conversion Pressurized Water Reactor (HCPWR). Sensitivities of the outlet coolant and cladding temperature to the subchannel flow area and local power peaking are high. One of the reasons is that the ratio of the coolant flow rate of SCFR-H to the power is smaller than that of LWR. Another reason is that, temperature of supercritical water is more sensitive to the enthalpy change above 450°C. The outlet coolant temperature distribution can be flattened by reducing the area of the peripheral subchannels and by enhancing the mixing between the subchannels.

The concept of a high temperature fast reactor cooled by supercritical water (SCFR-H) was developed for achieving high thermal efficiency and a compact reactor system. The core characteristics were obtained from single channel thermal-hydraulic analysis. Thus, it is necessary to carry out subchannel analysis to estimate the effect of local power peaking and cross flows. For this purpose, a subchannel analysis code is developed. It is verified by comparing the results with experimental data of High Conversion Pressurized Water Reactor (HCPWR). Sensitivities of the outlet coolant and cladding temperature to the subchannel flow area and local power peaking are high. One of the reasons is that the ratio of the coolant flow rate of SCFR-H to the power is smaller than that of LWR. Another reason is that, temperature of supercritical water is more sensitive to the enthalpy change above 450°C. The outlet coolant temperature distribution can be flattened by reducing the area of the peripheral subchannels and by enhancing the mixing between the subchannels.

Abstract An experimental study has been conducted on the mechanisms of flame propagation through combustible solid particle clouds of 1-octadecanol. The combustible particle cloud is ignited in its centre by an electric spark, and the growth of flame kernel is observed with a CCD video camera. The direct light emission and schlieren images of propagating flame and the laser light scattering images of particles have been simultaneously recorded. After ignition, a flame kernel is observed to grow with a yellow luminous zone whose outline is of an irregular shape. At the same time, a smooth shaped schlieren front is observed to propagate at 4–8 mm ahead of the outline of the yellow luminous zone. Inside the schlieren front, dispersed blue flame spots appear but no smaller particles can be seen, and only bigger particles are observed in the border region near the schlieren front. Across the schlieren front, smaller particles (most of them are about 10–20 μm in diameter) rapidly gasify just behind the schlieren front, while the gasification of particles with a diameter larger than 80 μm is delayed and the vapour lumps formed behind the schlieren front ignite to form circular dispersed blue flames. It has also been revealed that the average propagation velocity of the schlieren front increases with the number density of smaller particles, while it is scarcely affected by the mean diameter of combustible particle clouds. This fact implies that flame propagation is mainly supported by the combustion of smaller particles gasifying across the schlieren front.

Abstract An experimental study has been conducted on the mechanisms of flame propagation through combustible solid particle clouds of 1-octadecanol. The combustible particle cloud is ignited in its centre by an electric spark, and the growth of flame kernel is observed with a CCD video camera. The direct light emission and schlieren images of propagating flame and the laser light scattering images of particles have been simultaneously recorded. After ignition, a flame kernel is observed to grow with a yellow luminous zone whose outline is of an irregular shape. At the same time, a smooth shaped schlieren front is observed to propagate at 4–8 mm ahead of the outline of the yellow luminous zone. Inside the schlieren front, dispersed blue flame spots appear but no smaller particles can be seen, and only bigger particles are observed in the border region near the schlieren front. Across the schlieren front, smaller particles (most of them are about 10–20 μm in diameter) rapidly gasify just behind the schlieren front, while the gasification of particles with a diameter larger than 80 μm is delayed and the vapour lumps formed behind the schlieren front ignite to form circular dispersed blue flames. It has also been revealed that the average propagation velocity of the schlieren front increases with the number density of smaller particles, while it is scarcely affected by the mean diameter of combustible particle clouds. This fact implies that flame propagation is mainly supported by the combustion of smaller particles gasifying across the schlieren front.

A systematic method to calculate anharmonic force constants of crystals is presented. The method employs the direct-method approach, where anharmonic force constants are extracted from the trajectory of first-principles molecular dynamics simulations at high temperature. The method is applied to Si where accurate cubic and quartic force constants are obtained. We observe that higher-order correction is crucial to obtain accurate force constants from the trajectory with large atomic displacements. The calculated harmonic and anharmonic force constants are, then, combined with the Boltzmann transport equation (BTE) and non-equilibrium molecular dynamics (NEMD) methods in calculating the thermal conductivity. The BTE approach successfully predicts the lattice thermal conductivity of bulk Si, whereas NEMD shows considerable underestimates. To evaluate the linear extrapolation method employed in NEMD to estimate bulk values, we analyze the size dependence in NEMD based on BTE calculations. We observe strong nonlinearity in the size dependence of NEMD in Si, which can be ascribed to acoustic phonons having long mean-free-paths and carrying considerable heat. Subsequently, we also apply the whole method to a thermoelectric material Mg2Si and demonstrate the reliability of the NEMD method for systems with low thermal conductivities.

A systematic method to calculate anharmonic force constants of crystals is presented. The method employs the direct-method approach, where anharmonic force constants are extracted from the trajectory of first-principles molecular dynamics simulations at high temperature. The method is applied to Si where accurate cubic and quartic force constants are obtained. We observe that higher-order correction is crucial to obtain accurate force constants from the trajectory with large atomic displacements. The calculated harmonic and anharmonic force constants are, then, combined with the Boltzmann transport equation (BTE) and non-equilibrium molecular dynamics (NEMD) methods in calculating the thermal conductivity. The BTE approach successfully predicts the lattice thermal conductivity of bulk Si, whereas NEMD shows considerable underestimates. To evaluate the linear extrapolation method employed in NEMD to estimate bulk values, we analyze the size dependence in NEMD based on BTE calculations. We observe strong nonlinearity in the size dependence of NEMD in Si, which can be ascribed to acoustic phonons having long mean-free-paths and carrying considerable heat. Subsequently, we also apply the whole method to a thermoelectric material Mg2Si and demonstrate the reliability of the NEMD method for systems with low thermal conductivities.

The concept of resilience is a crucial part in crafting visions of desirable futures designed to withstand the widest variety of external shocks to the system. Backcasting scenarios are widely used to envision desirable futures with a discontinuous change from the present in mind. However, less effort has been devoted to developing theoretical frameworks and methods for building backcasting scenarios with a particular focus on resilience, although resilience has been explored in related sustainability fields. This paper proposes a method that helps design backcasting scenarios for resilient futures. A characteristic of the method is to delineate “collapse” futures, based upon which resilient futures are described to avoid the various collapsed states. In the process of designing backcasting scenarios, fault tree analysis (FTA) is used to support the generation of various risk factors and countermeasures to improve resilience. In order to test the effectiveness of the proposed method, we provide a case study to describe resilient energy systems for a Japanese community to 2030. Four expert workshops involving researchers from different disciplines were organized to generate diversified ideas on resilient energy systems. The results show that three scenarios of collapsed energy systems were described, in which policy options to be taken toward achieving resilient energy systems were derived.

The concept of resilience is a crucial part in crafting visions of desirable futures designed to withstand the widest variety of external shocks to the system. Backcasting scenarios are widely used to envision desirable futures with a discontinuous change from the present in mind. However, less effort has been devoted to developing theoretical frameworks and methods for building backcasting scenarios with a particular focus on resilience, although resilience has been explored in related sustainability fields. This paper proposes a method that helps design backcasting scenarios for resilient futures. A characteristic of the method is to delineate “collapse” futures, based upon which resilient futures are described to avoid the various collapsed states. In the process of designing backcasting scenarios, fault tree analysis (FTA) is used to support the generation of various risk factors and countermeasures to improve resilience. In order to test the effectiveness of the proposed method, we provide a case study to describe resilient energy systems for a Japanese community to 2030. Four expert workshops involving researchers from different disciplines were organized to generate diversified ideas on resilient energy systems. The results show that three scenarios of collapsed energy systems were described, in which policy options to be taken toward achieving resilient energy systems were derived.

Lightning surge characteristics of a transmission line comprising a tower, a shield wire and phase conductors are studied with the help of NEC-2. A tower struck by a lightning stroke behaves as an antenna until a traveling wave makes several roundtrips in the tower. During this interval, the tower footing impedance appears to be higher than the footing resistance, and the coupling coefficient between a shield wire and a phase conductor is much lower than that in the TEM mode. These affect the estimation of insulator voltages of transmission towers subject to lightning currents having short rise times. On the basis of these findings, the parameters in the multistory tower model, which has been used widely in EMTP multi-conductor analyzes, are newly proposed.

Lightning surge characteristics of a transmission line comprising a tower, a shield wire and phase conductors are studied with the help of NEC-2. A tower struck by a lightning stroke behaves as an antenna until a traveling wave makes several roundtrips in the tower. During this interval, the tower footing impedance appears to be higher than the footing resistance, and the coupling coefficient between a shield wire and a phase conductor is much lower than that in the TEM mode. These affect the estimation of insulator voltages of transmission towers subject to lightning currents having short rise times. On the basis of these findings, the parameters in the multistory tower model, which has been used widely in EMTP multi-conductor analyzes, are newly proposed.