ISIJ International
Online ISSN : 1347-5460
Print ISSN : 0915-1559
ISSN-L : 0915-1559
Regular Article
Phase Equilibria Studies in the CaO–SiO2–Al2O3–MgO System with CaO/SiO2 Ratio of 1.10
Xiaodong Ma Dianwei ZhangZhixing ZhaoTim EvansBaojun Zhao
著者情報
ジャーナル オープンアクセス HTML

2016 年 56 巻 4 号 p. 513-519

詳細
Abstract

Phase equilibria and liquidus temperatures in the CaO–SiO2–Al2O3–MgO system with CaO/SiO2 weight ratio of 1.10 have been experimentally determined by means of high temperature equilibration and quenching technique followed by electron probe X-ray microanalysis. Isotherms in the interval of 20 K between 1613 and 1733 K were determined in the primary phase fields of dicalcium silicate, wollastonite, merwinite, periclase, spinel and melilite that are relevant to ironmaking slags. Effects of Al2O3, MgO and CaO/SiO2 ratio on the liquidus temperatures have been discussed to assist optimum operation of iron blast furnace. Compositions of the solid solutions corresponding to the liquidus have been accurately measured that will be used for development of the thermodynamic database.

1. Introduction

In modern blast furnace (BF) process, the new challenges such as low gas permeability and formation of hearth accretion will constantly arise with the increasing use of low grade iron ores and poor fuels.1,2,3) It is essential to optimise the BF slag system for smooth operation and low cost requirement. Therefore, accurate phase equilibria become more important to explain the effects of Al2O3 and MgO on operational temperature as the oxide system CaO–SiO2–Al2O3–MgO forms a base for the tapping ironmaking slag. The phase diagrams in the system CaO–SiO2–Al2O3–MgO have been summarised in Slag Atlas in the form of pseudo-ternary sections with the isotherms of 100 K interval at fixed concentrations of Al2O3 or MgO.4) In the modern industry the operating parameters can be controlled more accurately. The larger intervals of compositions and temperatures in the previous phase diagrams cannot meet the requirements of current operations. Our recent work5) showed the significant differences of the phase diagrams in the Al2O3–CaO–MgO–SiO2 system from those reported in the early studies.6,7,8) In addition, melilite solid solutions of Al2O3–CaO–MgO–SiO2 system, in which industrial slags are operated, are important for optimisation of thermodynamic database. There are lack of melilite solid solution data reported previously. Hence, accurate phase equilibria data are of importance for both industrial and scientific interests.

The binary basicity (CaO/SiO2 weight ratio) of BF slags are mostly between 1.10 and 1.30 and it does not change significantly in a particular ironmaking plant. In modern ironmaking operations, the hot metal and slag temperatures need to be controlled within a narrow range to obtain stable BF operating condition. Phase equilibria and liquidus temperatures in the CaO–SiO2–Al2O3–MgO system with CaO/SiO2 weight ratio of 1.30 have been reported recently by the present authors.5) The present study is focusing on the liquidus surfaces in the CaO–SiO2–Al2O3–MgO system with CaO/SiO2 of 1.10 to cover the composition range of iron BF final slags.

2. Experimental

The experimental procedure in this investigation has been described in details in a previous paper.5) In brief, the experiments were carried out using a vertical electric resistance furnace. High-purity reagent powders of Al2O3, SiO2, MgO and CaCO3 were mixed and pelletized according to the required compositions, then placed in the hot zone of furnace to be premelted at a temperature 20–50 K higher than the desired temperature to ensure the homogeneous slag. The temperature was then decreased to the required one to achieve equilibration. Graphite crucible was used to hold the samples and ultrahigh purity Ar atmosphere was flashed continuously to protect the graphite crucibles. SiO gas may be generated by the reaction of graphite with SiO2 in the slag. As a result, the initial bulk composition of the sample could change. However, the compositions of the liquid and solid phases of quenched sample were measured after the experiment. The change of slag bulk composition only affects the ratios of liquid to solid but does not change the phase compositions. A Pt-30 pct Rh/Pt-6 pct Rh thermocouple placed in an alumina sheath was located adjacent to the sample to monitor the temperature. The temperature of the furnace was controlled within 2 K and overall temperature uncertainty is within 5 K. The samples were directly quenched into water after the equilibration, then dried, mounted and polished for metallographic analysis. The microstructures were examined by scanning electron microscopy coupled with energy-dispersive spectroscopy analysis (SEM-EDS). Compositions of the liquid and solid phases were measured by a JEOL JXA-8200 Electron Probe X-Ray Microanalyser (EPMA) with Wavelength Dispersive Spectrometers (WDS). The EPMA measurement conditions were the same as previous work.5) The homogeneity of the phase can be confirmed by the EPMA measurements in different areas of the quenched sample. Usually 8–20 points of liquid phase and 3–5 points of solid phase were measured from different areas by EPMA. The sample with the standard deviation of the phase compositions less than 1% will be accepted and used to construct the phase diagram.

3. Results and Discussions

3.1. Description of the Pseudo-ternary Section

Over 150 experiments have been carried out in the CaO–SiO2–Al2O3–MgO system with the CaO/SiO2 weight ratio of 1.10. The composition range covers Al2O3 concentration up to 25 wt% and MgO up to 20 wt%. Dicalcium silicate, wollastonite, spinel, periclase, merwinite and melilite were found to be the primary phases in the composition range investigated. Typical microstructures observed in the quenched samples are presented in Fig. 1. Figure 1 shows the equilibrium of liquids with (a) merwinite, (b) dicalcium silicate (c) wollastonite, (d) melilite, (e) spinel and (f) merwinite and melilite respectively. CaO/SiO2 ratio in dicalcium silicate, merwinite and melilite are mostly higher than 1.10. Precipitation of these solid phases from the melt will result in lower CaO/SiO2 ratio in the liquid phase. The difficulty in controlling CaO/SiO2 ratio at 1.10 can be overcame by preparing initial mixtures with slightly higher CaO/SiO2 ratios. The compositions of liquid and solid phases measured by EPMA after the equilibration are presented in Table 1. EPMA measurements show that the compositions of dicalcium silicate (Ca2SiO4), wollastonite (CaSiO3), spinel (MgO·Al2O3), and merwinite (3CaO·MgO·2SiO2) are closed to their stoichiometry. Melilite is the solid solution between akermanite (2CaO·MgO·2SiO2) and gehlenite (2CaO·Al2O3·SiO2).

Fig. 1.

Typical microstructures of slags quenched from (a) merwinite, (b) dicalcium silicate (c) wollastonite, (d) melilite and (e) spinel primary phase fields and (f) merwinite and melilite phase boundary.

Table 1. Experimental results in the system CaO–SiO2–Al2O3–MgO with CaO/SiO2 ratio.
Experiment
No.
Temperature
(K)
PhasesComposition (wt%)CaO/
SiO2
CaOMgOAl2O3SiO2
Liquid
1481633Liquid45.25.96.942.01.08
1501633Liquid44.86.09.040.21.11
1791653Liquid43.96.89.939.51.11
2001653Liquid48.52.65.643.31.12
1411673Liquid43.78.39.638.41.14
1991673Liquid48.25.23.743.01.12
1541693Liquid46.87.44.141.71.12
1561693Liquid43.79.09.238.11.15
1581693Liquid40.911.711.935.61.15
2321693Liquid38.910.814.535.91.08
2271713Liquid39.09.716.634.71.12
1471733Liquid45.29.33.741.71.08
1491733Liquid42.111.27.639.11.08
1751733Liquid47.29.42.041.31.14
2331733Liquid42.10.020.937.01.14
2341733Liquid40.94.917.336.91.11
Liquid with one oxide solid
Melilite primary phase field
1741613Liquid46.45.07.940.61.14
Melilite41.310.99.038.6
1801613Liquid46.55.26.142.21.10
Melilite41.113.03.542.4
1981613Liquid45.90.013.940.21.14
Melilite41.20.036.821.9
1961613Liquid45.03.710.840.61.11
Melilite41.36.620.531.5
1441633Liquid44.43.811.540.21.11
Melilite41.05.921.731.4
1461633Liquid43.20.016.340.51.07
Melilite40.60.036.622.8
1551633Liquid46.05.88.439.71.16
Melilite42.010.110.237.7
1761653Liquid42.30.017.839.81.06
Melilite41.10.036.922.0
1771653Liquid43.33.513.639.61.09
Melilite41.73.826.927.6
2021653Liquid43.75.910.839.61.10
Melilite40.98.017.034.0
2511653Liquid44.66.49.539.51.13
Melilite41.49.911.637.1
2151673Liquid42.12.016.939.01.08
Melilite41.51.632.524.3
2161673Liquid41.56.912.938.61.07
Melilite41.56.620.231.7
2171673Liquid38.312.512.237.01.04
Melilite41.37.518.432.8
2441673Liquid43.50.017.239.31.11
Melilite40.90.036.822.2
2451673Liquid43.35.812.138.91.11
Melilite41.05.822.430.8
2491673Liquid43.33.014.539.21.11
Melilite41.22.729.926.2
2501673Liquid42.97.711.138.31.12
Melilite41.17.916.934.1
2201693Liquid42.60.019.138.31.11
Melilite41.40.036.322.3
2211693Liquid41.65.115.338.01.09
Melilite41.33.927.027.9
2301693Liquid41.61.917.938.71.08
Melilite40.91.633.124.4
2311693Liquid40.96.814.737.61.09
Melilite41.05.124.329.6
2251713Liquid42.20.020.237.61.12
Melilite41.30.036.722.1
2471713Liquid41.63.117.737.61.11
Melilite41.21.931.925.1
2481713Liquid39.77.916.336.11.10
Melilite41.23.927.027.9
2361733Liquid41.40.021.537.11.12
Melilite41.00.036.622.4
2421733Liquid40.53.919.636.01.12
Melilite41.02.031.925.1
2431733Liquid37.58.220.134.21.10
Melilite40.92.830.226.1
Merwinite primary phase field
1571633Liquid46.76.05.641.71.12
Merwinite52.112.00.135.9
1711653Liquid44.87.47.939.81.13
Merwinite51.612.30.036.0
1401673Liquid44.97.47.540.21.12
Merwinite51.812.10.136.0
1971673Liquid46.47.24.541.91.11
Merwinite51.812.20.035.9
601673Liquid37.914.213.134.71.09
Merwinite51.112.60.136.1
631673Liquid39.012.512.436.01.08
Merwinite51.312.10.136.5
421693Liquid37.914.212.335.51.07
Merwinite50.612.80.136.4
431693Liquid40.111.010.438.31.05
Merwinite50.912.70.136.3
561693Liquid38.313.912.135.61.07
Merwinite52.312.10.135.5
1731693Liquid42.011.18.738.11.10
Merwinite51.412.50.036.0
1951693Liquid46.28.72.542.61.08
Merwinite51.712.30.036.0
1531713Liquid43.110.87.838.31.13
Merwinite51.612.20.036.1
1511733Liquid39.514.69.536.41.09
Merwinite50.912.60.136.4
Merwinite52.312.10.135.5
Ca2SiO4 primary phase field
2141633Liquid48.24.93.943.01.12
Ca2SiO462.72.00.335.1
Wollastonite primary phase field
2091613Liquid46.10.011.642.31.09
Wollastonite48.10.00.051.9
2101613Liquid46.30.012.241.51.12
Wollastonite48.30.00.151.6
2131613Liquid45.32.79.142.81.06
Wollastonite48.30.00.151.6
2011633Liquid49.10.08.842.11.17
Wollastonite48.30.00.051.7
2081633Liquid48.50.08.543.01.13
Wollastonite48.20.00.151.8
2121633Liquid46.62.96.943.51.07
Wollastonite48.30.00.151.6
Spinel primary phase field
91693Liquid37.813.615.133.41.13
Spinel0.328.171.40.1
251693Liquid37.112.716.133.91.09
Spinel0.328.170.90.1
501693Liquid36.914.614.533.91.09
Spinel0.028.471.50.1
581713Liquid37.113.516.133.21.12
Spinel0.328.471.00.2
11733Liquid35.716.215.132.91.09
Spinel0.428.171.30.2
31733Liquid35.814.116.933.11.08
Spinel0.428.071.40.2
51733Liquid35.610.221.133.11.08
Spinel0.427.971.60.1
Liquid with more oxide solids
2111593Liquid45.60.013.441.01.11
Wollastonite48.30.00.051.7
Melilite41.20.036.422.3
1521633Liquid45.26.38.040.51.11
Merwinite51.811.90.136.2
Melilite41.611.86.240.4
1351673Liquid37.614.214.034.21.10
Spinel0.128.071.70.1
Merwinite51.112.40.136.4
Melilite41.23.229.026.5
2461673Liquid41.69.910.338.21.09
Merwinite50.912.50.136.6
Melilite41.09.214.135.8
1391673Liquid46.66.63.743.11.08
Ca2SiO462.82.60.134.5
Merwinite52.011.80.036.2
21693Liquid36.815.813.334.11.08
Spinel0.627.871.00.4
Merwinite50.712.70.136.5
491693Liquid36.516.513.333.61.09
Spinel0.128.571.30.1
Merwinite50.912.60.136.5
1931713Liquid46.79.12.242.11.11
Periclase61.73.40.134.8
Merwinite51.912.30.035.8
441733Liquid35.59.522.332.61.09
Spinel0.128.571.30.1
Melilite40.83.827.927.3

The liquid compositions with CaO/SiO2 ratio close to 1.10 are used to construct liquidus isotherms of the pseudo-ternary section (CaO+SiO2)–Al2O3–MgO with a fixed CaO/SiO2 weight ratio of 1.10 as shown in Fig. 2.

Fig. 2.

Liquidus isotherms of pseudo-ternary section (CaO+SiO2)–Al2O3–MgO based on experimental results.

The previous data on the CaO–SiO2–Al2O3–MgO system6,7,8) have been used in optimisation of thermodynamic software FactSage 6.2.9) Predictions of FactSage 6.2 are shown in Fig. 3 to compare with the present experimental data. The databases selected in FactSage 6.2 are “Fact53” and “FToxide”, and the solution species selected in calculations are “FToxide-SLAGA”, “FToxide-SPINA”, “FToxide-MeO_A”, “FToxide-bC2S”, “FToxide-aC2S”, “FToxide-Mel”, and “FToxide-Merwinite”. The previous results reported by Osborn et al.,6) Cavalier and Sandrea-Deudon,8) Muan and Osborn10) are also shown in Fig. 3 for the comparison. It can be seen that the present experimental results have a good agreement with previous work reported by Osborn et al.,6) Cavalier and Sandrea-Deudon8) and Muan and Osborn10) at 1673 K. However, obvious differences can be observed between the present data and that reported by Muan and Osborn10) at 1713 and 1733 K. It can be seen that FactSage predictions show the similar overall trends as the experimental results, but the locations of the isotherms can be found significantly different. The FactSage software used the published data to optimize the databases and the accuracy of FactSage calculations rely on the accuracy of the experimental data used. The deviation of FactSage can be attributed to the limitation of the previous research techniques so that the accuracy of the data was not enough, in particular the solid solutions. Generally, the experimentally determined liquidus temperatures in the present study are approximately 20 K higher than that FactSage predictions in the melilite primary phase field. In MgO and spinel primary phase fields, the experimentally determined liquidus temperatures in the present study are usually 20 K lower than the predictions. In merwinite primary phase field, the experimentally determined liquidus temperatures agree well with the prediction values. It is noted that there is a small Ca3Si2O7 primary phase field according to the FactSage predictions. However, the Ca3Si2O7 primary phase was not observed in the present study. Instead, wollastonite (CaSiO3) primary phase field present in the composition range. Additionally, the Ca2SiO4 primary phase was observed by the present experiments at high (CaO+SiO2) corner, which is not shown by FactSage predictions.

Fig. 3.

Pseudo-ternary section (CaO+SiO2)–Al2O3–MgO with CaO/SiO2 ratio of 1.10.

To compare the effect of CaO/SiO2 ratio on primary phase fields and liquidus temperatures, the phase boundaries and isotherm at 1723 K are shown in Fig. 4 for CaO/SiO2 ratios of 1.10 and 1.30 respectively. It can be seen that the wollastonite primary phase field was not stable when the CaO/SiO2 ratio is 1.30. The dicalcium silicate primary phase area expands significantly towards Al2O3 corner when the CaO/SiO2 ratio is increased from 1.10 to 1.30. In contrast, the sizes of the melilite and merwinite primary phase fields are decreased with increasing CaO/SiO2 ratio. Full liquid area represents a composition range where the slag is fully liquid phase above a given temperature. A large full liquid area indicates that the liquidus temperature of the slag is not sensitive to the variation of the composition. The full liquid area at 1723 K is described by the isotherms from different primary phase fields. It can be seen that the full liquid area at 1723 K is narrow at CaO/SiO2 ratio of 1.30 and it is expanded significantly when the CaO/SiO2 ratio is decreased to 1.10.

Fig. 4.

A comparison in pseudo-ternary sections (CaO+SiO2)–Al2O3–MgO between CaO/SiO2 ratio of 1.10 and 1.305) at 1723 K.

Melilite is the solid solution between akermanite (2CaO·MgO·2SiO2) and gehlenite (2CaO·Al2O3·SiO2). The distributions of MgO and Al2O3 between melilite and liquid are shown in Figs. 5 and 6 respectively. The solid lines in Figs. 5 and 6 are fitted from the experimental data that may have slightly different CaO/SiO2 ratios. The predictions from FactSage 6.2 are also shown in the figures for comparison. As seen in Fig. 5, the MgO concentrations in melilite solid solution increase with increasing that in liquid up to 15 mol% MgO. At 1673 K, the MgO concentrations in the melilite reach a maximum at 15 mol% MgO and then decrease with increasing MgO in the corresponding liquid. It can be seen that the MgO concentrations in melilite decrease with increasing temperature at a given MgO concentration in liquid. At 1613 K, the MgO concentration in the melilite is approximately as twice as that in the corresponding liquid. At 1673 K, the MgO concentration in the melilite is almost the same as that in the corresponding liquid. In contrast, the MgO concentration in the melilite is only half as that in the corresponding liquid at 1733 K. The predicted MgO concentrations in melilite by FactSage 6.2 show a good agreement with experimental results at 1613, 1673 and 1733 K up to 15 mol% MgO in the liquid. However, FactSage predictions show that merwinite is a stable phase at 1673 K over 15 mol% MgO in the liquid which is different from the present results.

Fig. 5.

Distribution of MgO between melilite and liquid for temperatures from 1613 to 1733 K.

Fig. 6.

Distribution of Al2O3 between melilite and liquid for temperatures from 1613 to 1733 K.

It can be seen from Fig. 6 that the Al2O3 concentrations in melilite solid solution increase with increasing that in the corresponding liquid at 1613 and 1673 K. The bending orientation of curves gradually change from 1613 to 1733 K which reflects the same shape of the isotherms in the melilite primary phase field. Similar retrograde curves were also reported at CaO/SiO2 ratio of 1.30 from 1723 to 1773 K.5) The Al2O3 concentration in melilite is much higher than that in the corresponding liquid phase at a given temperature. It is interesting to see that FactSage can predict the general trend of the distributions but the values are significantly different which may explain the difference in liquidus temperatures predicted.

3.2. Application of Phase Diagram in Ironmaking Process

3.2.1. Effect of MgO on Liquidus Temperature

Recently BF operation is facing the challenges to lower the cost by increased utilization of pulverized coal injection and low-grade ores, and reduced flux addition such as MgO. These changes would result in higher Al2O3 or lower MgO in the slag.11,12,13) The Figs. 7 and 8 were redrawed in the form of pseudo-binary section for better understanding from pseudo-ternary phase diagram of Fig. 3. The data was collected from Fig. 3 when fixing the Al2O3 concentration of 15% and 20%, or MgO concentration of 5% and 10%. Figure 7 shows the relationship between the liquidus temperatures and MgO concentration at fixed Al2O3 of 15 and 20 wt% and CaO/SiO2 ratio of 1.10 and 1.30. 15 wt% Al2O3 in slag represents an average in the current BF operations and 20 wt% Al2O3 in slag is a result of low grade iron ore utilization. It can be seen from Fig. 7 that the stable primary phase field is melilite at low MgO concentrations and spinel at high MgO concentrations. The liquidus temperatures always increase with increasing MgO concentration in the spinel primary phase field. In the melilite primary phase field the liquidus temperatures have different behaviours depending on the Al2O3 concentration and CaO/SiO2 ratio. It can be seen that at 15 wt% Al2O3 and CaO/SiO2 ratio of 1.10, the liquidus temperatures first increase rapidly and decrease with increasing MgO. At 20 wt% Al2O3 and CaO/SiO2 ratio of 1.10, the trend of the liquidus temperature changes with MgO is similar but it is less sensitive than 15 wt% Al2O3. At 20 wt% Al2O3 and CaO/SiO2 ratio of 1.30, the liquid temperatures do not change with MgO at MgO below 5 wt%. However, the liquid temperatures decrease with increasing MgO above 5 wt%. The stability of the spinel phase seems not sensitive to CaO/SiO2 ratio but to the Al2O3 concentration. For example, at 20 wt% Al2O3, the boundary between the melilite and spinel primary phase fields is the same at 10 wt% MgO for the CaO/SiO2 ratios of 1.10 and 1.30. However, at the same CaO/SiO2 ratio of 1.10, the spinel phase will be formed above 10 and 13 wt% MgO for the Al2O3 concentrations 20 and 15 wt% respectively. The liquidus temperatures in the melilite primary phase field increase with increasing both Al2O3 concentration and CaO/SiO2 ratio. It seems that the increment of the melilite liquidus temperature with CaO/SiO2 ratio is more sensitive at lower MgO concentrations. In contrast, the liquidus temperatures in the spinel primary phase field are more sensitive to the Al2O3 concentration and less sensitive to the CaO/SiO2 ratio. The accurate information given in the present study will provide useful guide to BF operators to work at temperatures high enough to avoid the precipitation of solid phases when complex feeds are used in ironmaking.

Fig. 7.

The effect of MgO on the liquidus temperature in which Al2O3=15 and 20 wt%, CaO/SiO2=1.10 and 1.30.

Fig. 8.

The effect of Al2O3 on the liquidus temperature in which MgO=5 and 10 wt%, CaO/SiO2=1.10 and 1.30.

3.2.2. Effect of Al2O3 on Liquidus Temperature

Figure 8 shows the relationship between liquidus temperature and Al2O3 concentration at fixed MgO of 5 and 10 wt% and CaO/SiO2 ratio of 1.10 and 1.30. 10 wt% MgO in slag represents an average in the current operations and 5 wt% MgO in slag is a result of low slag rate operation. It can be seen that Al2O3 has significant influence in primary phase field and liquidus temperature in the composition range investigated. Dicalcium silicate and merwinite are the major stable phases at low Al2O3 concentrations and melilite is the major stable phase at high Al2O3 concentrations. Wollastonite (CaSiO3) also appears at 5 wt% MgO and CaO/SiO2 ratio of 1.10, but disappears when the MgO concentration is 10 wt%. The spinel also appears at 10 wt% MgO and CaO/SiO2 ratio of 1.10. In the melilite and spinel primary phase fields the liquidus temperatures increase significantly with increasing Al2O3 concentration. In contrast, the liquidus temperatures in the primary phase fields of dicalcium silicate, merwinite and wollastonite decrease significantly with increasing Al2O3 concentration. At CaO/SiO2 ratio of 1.10, the liquidus temperature with 5 wt% MgO is lower than that with 10 wt% MgO below 16 wt% Al2O3. At 5 wt% MgO concentration for low slag volume operation, the liquidus temperatures are increased by approximate 40 K in the melilite primary phase field above 12 wt% Al2O3 when the CaO/SiO2 ratio increases from 1.10 to 1.30. Clearly the operating range of Al2O3 where the slag is fully liquid is much wider at CaO/SiO2 ratio of 1.10. The BF operation at high CaO/SiO2 ratio and low MgO concentration will be very difficult.

4. Conclusions

The phase equilibria and liquidus temperatures in the CaO–SiO2–Al2O3–MgO system with CaO/SiO2 ratio of 1.10 have been experimentally determined in the composition range related to blast furnace slags. Providing a fixed MgO or Al2O3 concentration closed to the composition of industrial slag, the liquidus temperatures increases with increasing Al2O3 concentration, but not sensitive with increasing MgO concentration. It is confirmed that low-slag volume BF operation can be achieved by decreasing MgO concentration at certain conditions. Extensive solid solutions of melilite have been accurately determined that will provide invaluable data for optimisation of thermodynamic database.

Acknowledgements

The authors would like to thank Ms. Jie Yu for the lab assistance in the high temperature experiments and financial support from University of Queensland, Shougang and Rio Tinto. The authors also would like to thank Mr Ron Rasch and Ms Ying Yu in Centre for Microscopy and Microanalysis (CMM) for technical support of EPMA and SEM.

References
 
© 2016 by The Iron and Steel Institute of Japan
feedback
Top