ISIJ International
Online ISSN : 1347-5460
Print ISSN : 0915-1559
ISSN-L : 0915-1559
Regular Article
Reassessment of TiN(s)=Ti+N Equilibration in Liquid Iron
Jung-Mock JangSeok-Hyo SeoJoon-Suk HanDong-Sik KimYoun-Bae KangJong-Jin Pak
Author information
JOURNAL OPEN ACCESS FULL-TEXT HTML

2015 Volume 55 Issue 11 Pages 2318-2324

Details
Abstract

Thermodynamics of nitrogen solubility and TiN formation in liquid Fe–Ti alloys was reassessed over the temperature range from 1823 to 1973 K by combining new experimental data for the Ti–N–TiN relations in liquid iron at 1823–1873 K together with the authors’ previous data. Using the Wagner’s Interaction Parameter Formalism(WIPF), the first-order interaction parameters of , and , and the equilibrium constant for the TiN formation reaction in liquid iron, logKTiN were newly determined. The interaction parameters showed no temperature dependence in this temperature range. The calculated solubility product of titanium and nitrogen for TiN formation in liquid iron was in excellent agreement with the available experimental results measured at various temperatures.

1. Introduction

Titanium is often added in liquid steel during secondary refining process for various steel products to improve mechanical properties and corrosion resistance.1,2) On the other hand, an excessive formation of TiN inclusions in liquid steel can cause a nozzle clogging problem during continuous casting and various defects in final products.3) In particular, the deleterious effect of 6 μm TiN inclusion on the fatigue life of steel is equivalent to the harm of a 25 μm oxide inclusion.4)

In order to control TiN inclusion formation based on the supersaturation of [Ti] and [N] in liquid steel during cooling and casting processes, it is essential to have accurate information on thermodynamics of titanium, nitrogen and TiN formation in liquid iron over a wide range of melt temperature. In spite of such importance, there are still some uncertainties in the literature values including the recommended values of the Japan Society for the Promotion of Science (JSPS).5) Furthermore, the recommended values of JSPS for the equilibrium constant, KTiN for the formation of TiN and the interaction parameter values of e N Ti and e Ti Ti in liquid iron were taken from different sources6,7) in which the data were determined by different experimental systems. This situation will cause significant errors in predicting the reaction equilibrium of TiN formation in Fe–Ti–N melt. In the authors’ recent study,8) the thermodynamics of Fe–Ti–N–TiN system was studied using the metal-gas and the metal-gas-nitride equilibration technique in the temperature range from 1873 to 1973 K, and an extensive discussion on the accuracy of those thermodynamic parameters has been carried out. Later the similar studies were carried out for Fe–Ti–Al–N, Fe–Ti–Cr–N, Fe–Ti–Mo–N and Fe–Ti–Nb–N alloys at the same temperature range.9,10,11) The accumulation of experimental data allowed us to obtain more reliable data on the thermodynamics of Fe–Ti–N–TiN system.

In the present study, the interaction parameters between titanium and nitrogen, and the equilibrium constant for the formation of TiN in liquid iron were reassessed over the temperature range from 1823 to 1973 K by combining new experimental data for the Ti–N–TiN relations in liquid iron at 1823–1873 K together with the authors’ previous data.8,9,10,11)

2. Experimental

The metal-gas and metal-nitride-gas equilibration experiments were carried out to determine the thermodynamics of titanium, nitrogen and TiN formation in Fe–Ti–N melts. Detailed descriptions of the experimental apparatus and procedure are available in the authors’ recent studies.8,9,10,11) Five hundred grams of high purity electrolytic iron (99.99 mass% purity, 60 mass ppm O, <5 mass ppm N, 18 mass ppm C, <5 mass ppm Si, <7 mass ppm Ni, 1 mass ppm Al) contained in an Al2O3 or MgO crucible was melted in the temperature range from 1823 to 1873 K by a high frequency induction furnace. Prior to the nitrogen solubility measurement, the iron melt was deoxidized by an Ar-10%H2 gas blown onto the melt surface at a high flow rate of ~5000 ml/min for 2 h. The oxygen content in the melt decreased to a value less than 20 mass ppm. The nitrogen partial pressure in the system was controlled by a mixture of Ar-10%H2 and N2 gases.

Strong agitation of the melt by an induction furnace resulted in a fast attainment of equilibrium nitrogen solubility in liquid iron under a nitrogen partial pressure within 1 h. This was confirmed by samplings and in-situ analysis for nitrogen content in the metal samples. Pellets of titanium (99.99 mass% purity) were added into liquid iron through an 18 mm ID quartz tube after confirming the equilibrium nitrogen solubility in pure liquid iron under a nitrogen partial pressure. After each titanium addition, new nitrogen solubility equilibrium was attained within 1 h. Titanium addition and sampling were repeated until a stable TiN layer was formed on the surface of the iron melt. The formation of TiN in the iron melt could be also confirmed by a sharp decrease in nitrogen content checked by the analysis of metal samples during the experiment. In author’s previous studies,8,9,10,11) the detailed procedure for chemical analysis is available.

3. Results and Discussion

3.1. Effect of Titanium on Nitrogen Solubility in Liquid Iron

Table 1 summarizes the present experimental results of nitrogen solubility measurement with titanium additions in liquid iron. Figure 1 shows the variation of equilibrium nitrogen solubility in liquid iron with titanium additions under different nitrogen partial pressures at 1823 and 1873 K. The nitrogen solubility increases linearly as the titanium content increases in liquid iron when the melt is not saturated with TiN as shown as open symbols in the figure. When TiN is formed at higher Ti contents, the nitrogen solubility decreased significantly as shown as solid symbols in the figure. The Ti–N relation in liquid iron was not affected by the crucible material of Al2O3 or MgO.

Table 1. Experimental results of Ti and N Equilibration in liquid iron.
Temp (K)PN2 (atm)Crucible[% Ti][% N][% O][% Si][% Al][% Mg]Saturation
18230.05MgO00.01000.00170.00490.0007
0.020.01000.00080.00420.0003
0.040.01020.00110.00380.0004
0.060.01030.00050.00450.0005
0.080.01030.00060.00480.0001
0.100.01040.00060.00360.0012
0.110.01050.00050.00120.0005
0.120.01040.00080.00360.0003TiN
0.190.00680.00070.00470.0007TiN
0.240.00530.00040.00210.0002TiN
0.2Al2O300.01990.00070.00230.0008
0.020.02020.00130.00370.0001
0.050.02050.00100.00470.0011
0.060.02020.00070.00340.0007TiN
0.130.00970.00050.00140.0003TiN
0.180.00680.00080.00310.0003TiN
18730.3MgO00.02420.00180.00420.0005
0.0220.02430.00160.00210.0005
0.0450.02500.00140.00270.0007
0.0660.02570.00140.00320.0008
0.0750.02600.00180.00140.0013TiN
0.0770.02590.00110.00150.0005TiN
0.0910.02010.00090.00080.0009TiN
0.0940.01920.00090.00210.0009TiN
0.1090.01790.00070.00150.0005TiN
Al2O300.02460.00180.00320.0002
0.030.02470.00180.00440.0007
0.050.02550.00160.00380.0015
0.070.02580.00180.00450.0002
0.080.02500.00200.00480.0009TiN
0.090.02080.00160.00360.0002TiN
0.110.01660.00190.00360.0007TiN
0.140.01390.00070.00470.0018TiN
Fig. 1.

Equilibrium [%Ti]-[%N] relation in Fe–Ti–N melts at 1823–1873 K.

The dissolution of nitrogen before the saturation of TiN in liquid iron can be written as   

1 2 N 2 (g)=N (1)
  
Δ G 1 o =3   598+23.89T      J/mol 12)
  
K N = f N [%N] P N 2 1/2 (2)
where KN is the equilibrium constant for Reaction (1), [%N] is the equilibrium nitrogen content in mass% and, fN is the activity coefficient of nitrogen in the 1 mass% standard state in liquid iron. PN2 is the nitrogen partial pressure in atm over the melt surface.

In Fe–Ti–N alloy melts, the Wagner’s Interaction Parameter Formalism (WIPF) can be used to estimate the fN in Eq. (2) along with first-order interaction parameters as:   

log f N = e N Ti [%Ti]+ e N N [%N]=log K 1 -log[%N]+ 1 2 logP N 2 (3)
where the value of e N N is known to be 0 in the present experimental temperature range.13) As mentioned earlier, the oxygen content in the melt was very low, and the effect of oxygen on nitrogen was assumed to be negligible. Silicon, aluminum and magnesium contents in the melt were less than 50, 18 and 13 mass ppm, respectively, and the effects of those elements can be also neglected.

Figure 2 shows the relation of logfN vs titanium content in mass% in liquid iron using the relation expressed by Eq. (3) from the nitrogen solubility data determined in the present study at 1823–1873 K together with the data determined in the authors’ recent studies at 1873–1973 K.8,9) The results determined at various temperatures and nitrogen partial pressures show excellent linear relationships up to 0.19 mass% Ti and the e N Ti value determined by the linear regression analysis of the data is −0.21 at all temperatures as shown in the figure. Discrepancy of experimental data obtained at low temperatures was less than those obtained at high temperatures. The standard deviations for the interaction parameter determined at each temperature are shown in Fig. 2. Temperature dependence was not evaluated for this parameter because it was found to be within experimental uncertainties. The dotted lines in Fig. 2 show the calculated relations using the e N Ti values (= −8507/T+4.18) determined in the temperature range of 1873–1973 K by Kim et al.8) The discrepancy of e N Ti value is significant at a 1823 K.

Fig. 2.

Relation of log f N Ti vs. [%Ti] in Fe–Ti–N melts at 1823–1973 K.

Figure 3 shows the same relation for the data determined at all temperatures. The e N Ti value of −0.21±0.016 can be determined by the linear regression analysis and it gives an excellent correlation of data for all temperatures. Using the WIPF’s reciprocal relationship,14) the value of e Ti N can be evaluated as −0.728±0.044. Figure 3 also compares the relations determined by previous workers at 1873 K, and there is huge scatter among their e N Ti values from −0.36 to −0.87.6,8,15,17) They all reported temperature dependence of the e N Ti values as shown in Table 2.

Fig. 3.

Relation of log f N Ti vs. [%Ti] in Fe–Ti–N melts together with previous results.

Table 2. Thermodynamic parameters of Ti and N in liquid iron at 1873 K.
e N Ti
(Temp.
Dependency)
e Ti Ti
(Temp.
Dependency)
logKTiN
(Temp.
Dependency)
Temp. (K)[% Ti] rangeMethodRef.
−0.210.048−2.742
(−12740/T+4.06)
1823–1973<0.52Sampling
(Fe–Ti–N)
Present study
−0.36
(−8507/T+4.18)
0.048−2.792
(−15780/T+5.63)
1873–1973<0.52Sampling
(Fe–Ti–N)
8
*−0.60
(−5700/T+2.45)
0.53 (1873 K)*−2.760
(−19800/T+7.78)
1873–1973<0.5Sampling
(Fe–Ti–N)
6
0.29 (1923 K)
0.26 (1973 K)
−0.53
(−6200/T+2.7)
−2.584
(−16433/T+6.19)
1873–2023<0.32Sieverts
(Fe–Ti–N)
15
−0.631973<0.3Sampling
(Fe–Ti–N)
16
−0.87
(−7240/T+3.0)
1853–1953<0.25Sampling
(Fe–Ti–N)
17
*0.0481873<8.8EMF
(Fe–Ti–O)
7
0.0481900<1Estimation from the data of Hadley et al.21) and Kelly et al.22)
(Fe–Ti–O)
18
0.0421873<12EMF
(Fe–Ti–O)
19
0.0231873<2.1Ag bath
Iso-activity
20
(*:  Recommended values of JSPS)

In the present study, a constant e N Ti value of −0.21 was obtained in the temperature range of 1823–1973 K without noticeable temperature dependence. The accuracy of the e N Ti value is important in studying thermodynamics of TiN(s) = Ti + N equilibrium in liquid iron. The validity of e N Ti values determined in the present study will be discussed more in detail in the next section.

3.2. TiN Formation in Liquid Iron

The solubility product of titanium and nitrogen in liquid iron for TiN saturation can be determined as shown as solid symbols in Fig. 1. The reaction equilibrium for the formation of pure solid TiN in liquid iron is   

TiN(s)= Ti _ + N _    ;   Δ G o (4)
  
K TiN = h Ti h N a TiN = f Ti f N [%Ti][%N] (5)

The equilibrium constant of the reaction (4) can be expressed as   

log K TiN =- Δ G o 2.303RT =log h Ti h N =log f Ti +log f N +log[%Ti][%N] = e Ti Ti [%Ti]+ e Ti N [%N]+ e N Ti [%Ti]+log[%Ti][%N] (6)
where hi and fi are the activity and activity coefficient of i in liquid iron with respect to 1 mass% standard state, respectively. e Ti Ti is the first-order self-interaction parameter of titanium in liquid iron and, the values of e N Ti and e Ti N were determined in the preceding section. The TiN formed in the melt was identified as a pure solid stoichiometric TiN by the XRD analysis.8,9,10,11) Therefore, the activity of TiN is unity under the present experimental condition.

Then, Eq. (6) can be rearranged as   

- e Ti N [%N]- e N Ti [%Ti]-log[%Ti][%N]= e Ti Ti [%Ti]-log K TiN (7)

The left hand side of Eq. (7) can be plotted versus [%Ti], and the e Ti Ti and the logKTiN values are simultaneously obtained by taking the slope and the intercept at y-axis of the plot, provided that the experimental data could be well represented by Eq. (7). One such plot is shown in Fig. 4 using the present experimental data for the TiN solubility product in liquid iron at 1823 and 1873 K together with the previous data determined in the authors’ recent studies.8,9,10,11) As can be seen, the data points determined at different nitrogen partial pressures and temperatures show excellent linear relationships.

Fig. 4.

Relation of Eq. (7) to determine e Ti Ti and logKTiN values.

The e Ti Ti value can be obtained as 0.048±0.0024 from the slope of the lines in the temperature range from 1823 to 1973 K. Temperature dependence was not evaluated for this parameter within experimental uncertainties. The values of logKTiN are obtained as −2.935, −2.742, −2.570 and −2.402 at 1823, 1873, 1923 and 1973 K, respectively. The temperature dependence of logKTiN is shown in Fig. 5, and it can be expressed as −12740/T+4.06. The standard Gibbs free energy change of Reaction (4) is then obtained as   

Δ G TiN o =244   000-77.7T      J/mol (8)
Fig. 5.

Temperature dependence of logKTiN values.

In the present study, three important thermodynamic variables of e N Ti , e Ti Ti and logKTiN were determined simultaneously in one consistent experimental method over the temperature range from 1823 to 1973 K by combining new experimental data for the Ti–N–TiN relations in liquid iron at 1823–1873 K together with the authors’ previous data at 1873–1973 K.8,9,10,11) The first-order interaction parameters apparently showed no temperature dependence in this temperature range, and the logKTiN values showed a less temperature dependence compared to the previous results as shown in Fig. 4.

An extensive discussion on those thermodynamic parameters has been carried out in the authors’ previous study.8) Morita and Kunisada6) reported their own combination of e Ti Ti and logKTiN values in order to best represent their experimental data measured by the similar method using an induction furnace as the present study. They reported the e Ti Ti values in liquid iron as 0.53, 0.29 and 0.26 at 1873, 1923 and 1973 K, respectively, along with the e N Ti = −5700/T + 2.45 as shown in Table 2. The e Ti Ti values in liquid iron determined by various authors are compared in Table 2.6,7,8,18,19,20) In spite of different experimental techniques, the values of e Ti Ti are quite well agreed in the range of 0.042–0.048 at 1873 K except for the values obtained by Morita and Kunisada6) and Yuanchang et al.20) The recommended value for e Ti Ti in JSPS (= 0.048 at 1873 K) was taken from the data of Janke and Fischer7) who studied the titanium deoxidation equilibrium in liquid iron by the EMF method. The e Ti Ti value of 0.048 at 1823–1973 K determined in Fe–Ti–N system in the present study is well agreed with the values determined in Fe–Ti–O system.7,18,19)

Figure 6 compares the available experimental results6,8,9,10,11,15,17) and the predictions of equilibrium solubility product of titanium and nitrogen for TiN formation in liquid iron in the temperature range of 1823–1973 K. The lines were calculated using Eq. (7) from logKTiN and interaction parameters determined by the present study and other workers6,8,15) including the recommended values of JSPS5) as shown in Table 2.

Fig. 6.

Equilibrium [%Ti]-[%N] relations in Fe–Ti–N melts saturated with TiN at (a) 1823, (b) 1873, (c) 1923 and (d) 1973 K.

The experimental data for TiN solubility product in liquid iron measured at 1823 K are in excellent agreement with the solid line calculated in the present study. The dotted line calculated at the same temperature using the parameters determined by the author’s previous study8) is lower than the experimental data at Ti content less than 0.3 mass%. At 1873 K, the experimental data and the calculated lines are well agreed among various researchers5,6,8) except Evans and Pehlke’s data.15) Evans and Pehlke15) measured the solubility product of titanium and nitrogen for TiN formation in liquid Fe–Ti alloys by the Sieverts’ method in the temperature range of 1873–2023 K, by determining the break points in the Sieverts’ law lines of nitrogen solubility for the precipitation of TiN with increasing nitrogen partial pressure. As shown in the figure, the solubility product data measured by the Sieverts’ method shows significantly higher values than the data measured by the sampling method. As pointed out by Pehlke and Elliott in their earlier work with Fe–Ti melts,12) the formation of a stable nitride in the sieverts chamber would increase the apparent absorption of nitrogen.

The calculated lines at 1823–1973 K using the recommended values in JSPS are not in agreement with any of those experimental data.6,8,9,10,11,17) As mentioned earlier, the recommended values of e N Ti , e Ti Ti and logKTiN in JSPS were taken from the different sources6,7) where those data were determined by completely different experimental techniques. This may cause errors in predicting the reaction equilibrium in Fe–Ti–N melts.

At 1873 K, the solubility product data measured by Morita and Kunisada6) are in excellent agreement with the present author’s results.8,9,10,11) However, at 1923 and 1973 K, the solubility product data and predicted lines of Morita and Kunisada6) are higher than those determined in the present study although the experimental methods are very similar. In the Morita and Kunisada’s study, they measured the melt temperature indirectly by a thermocouple located at the bottom of crucible. In the present author’s study, the melt temperature was directly monitored by the immersion of an alumina sheathed thermocouple, and it was controlled within 2 K during experiment by the PID controller of the induction furnace. Figure 6 also shows the solubility product data measured by Ishii and Fuwa17) at 1903 and 1953 K. They are in excellent agreement with the solubility product lines at the same temperatures calculated using the thermodynamic data determined in the present study.

3.3. Temperature Dependence of Interaction Parameters

In the preceding sections, important thermodynamic interaction parameters of e N Ti and e Ti Ti were determined in the temperature range from 1823–1973 K. They apparently showed no temperature dependence in this temperature range. However, one may expect that any thermodynamic variables have temperature dependence by their inherent nature. In this section, the temperature dependence of interaction parameters is briefly discussed.

According to Darken’s quadratic formalism,21,22) in a binary system 1-2 where “1” is the solvent, Gibbs energy of the binary solution in a dilute region for the solute “2” can be expressed by the following equation.   

g= X 1 g 1 ° + X 2 ( g 2 ° + C 2 )+ RT(X 1 lnX 1 + X 2 lnX 2 )+ α 12 X 1 X 2 (9)
where C2 is introduced to change the standard state of 2 from Raoultian to Henrian, and α12 is an interaction energy between 1 and 2, and can be represented by:   
α 12 = h 12 -T s 12 (10)
where h12 and s12 are related to the excess enthalpy and the excess entropy, respectively.

The excess Gibbs energy is therefore:   

g ex = α 12 X 1 X 2 + C 2 X 2 (11)

From the definition of the excess Gibbs energy in the binary 1–2 system   

g ex =RT( X 1 ln γ 1 + X 2 ln γ 2 ) (12)

By differentiation of the Eq. (11) assuming h12 and s12 are constants, and by comparison with the Eq. (12), one obtains:   

RTln γ 1 = α 12 X 2 2 (13)
  
RTln γ 2 = C 2 + α 12 (1-2 X 2 + X 2 2 ) (14)

As lnγ2 = ln γ 2 ° + ε 2 2 X2+ ρ 2 2 X 2 2 , and neglecting the second order term,   

ε 2 2 =- 2 α 12 RT =- 2( h 12 -T s 12 ) RT (15)

The ε 2 2 is converted by utilizing the following relations derived Lupis,23) where “1” denotes the solvent.   

e 2 2 = 1 230 ( M 1 M 2 ε 2 2 - M 1 - M 2 M 2 ) (16)
  
=- 2 230 ( M 1 M 2 h 12 R ) 1 T + 2 230 ( M 1 M 2 s 12 R ) - M 1 - M 2 230 M 2 (17)

The temperature dependence of the e 2 2 is   

d e 2 2 d(1/T) =- 2 230 ( M 1 M 2 ) ( h 12 R ) (18)

Similarly, in a 1–2–3 ternary solution, the excess Gibbs energy can be expressed by the following equation.21,22)   

g ex = C 2 X 2 + C 3 X 3 + α 12 X 1 X 2 + α 23 X 2 X 3 + α 31 X 3 X 1 (19)
and by following similar steps shown for the binary solution,   
RTln γ 1 = α 12 X 2 2 + α 31 X 3 2 +( α 12 + α 31 - α 23 ) X 2 X 3 (20)
  
RTln γ 2 =( α 12 + C 2 )-2 α 12 X 2 -( α 12 + α 31 - α 23 ) X 3 +[ α 12 X 2 2 + α 31 X 3 2 +( α 12 + α 31 - α 23 ) X 2 X 3 ] (21)
  
RTln γ 3 =( α 31 + C 3 )-2 α 31 X 3 -( α 12 + α 31 - α 23 ) X 2 +[ α 12 X 2 2 + α 31 X 3 2 +( α 12 + α 31 - α 23 ) X 2 X 3 ] (22)
where - ( α 12 + α 31 - α 23 ) RT corresponds to ε 2 3 which can be expressed by the following equation:   
ε 2 3 =- ( α 12 + α 31 - α 23 ) RT (23)
  
=- ( h 12 -T s 12 )+( h 31 -T s 31 )-( h 23 -T s 23 ) RT (24)
  
=- h 12 + h 31 - h 23 R 1 T + s 12 + s 31 - s 23 R (25)

When, Fe, Ti and N are substituted into 1, 2 and 3, respectively, ε Ti N becomes:   

ε Ti N =- h FeTi + h FeN - h TiN R 1 T + s FeTi + s FeN - s TiN R = ε N Ti (26)

The e N Ti is expressed as:   

e N Ti = 1 230 ( M Fe M Ti ε N Ti - M Fe - M Ti M Ti ) (27)
  
=- 1 230 ( M Fe M Ti ) ( h FeTi + h FeN - h TiN R ) 1 T + 1 230 ( M Fe M Ti ) ( s FeTi + s FeN - s TiN R ) - M Fe - M Ti 230 M Ti (28)

Then, the temperature dependence of the e N Ti (= ATiN/T + BTiN) is   

A TiN = d e N Ti d(1/T) =- 1 230R ( M Fe M Ti ) ( h FeTi + h FeN - h TiN ) (29)

Therefore, if hFeTi + hFeNhTiN ≈ 0, then the temperature dependence of e N Ti becomes almost zero.

The values of hFeTi, hFeN and hTiN have been reported in literature, but those are not just constant, but were expressed as functions on compositions.24,25,26) By differentiating the Eq. (19) taking into account the composition dependence of hFeTi, hFeN and hTiN, the values corresponding to the term in the last parenthesis in the Eq. (29) becomes 772 kJ/mol. Inserting these values into the Eq. (29) gives the ATiN ≈ −471 J/mol. This value is not zero, but is much smaller than the reported values (−8507,8) −5700,6) −620015) and −7240 J/mol17) as shown in Table 2). The above estimation assumes no ternary interaction in the Fe–Ti–N system. The value estimated in this section shows the validity of e N Ti and e Ti Ti values with negligible temperature dependence determined in the present study, compared to those reported in the literature.6,8,15,17)

4. Conclusions

The values of Wagner’s interaction parameters of e N Ti , e Ti N and e Ti Ti and the equilibrium constant for the TiN formation reaction in liquid iron, logKTiN were newly determined over the temperature range from 1823 to 1973 K as follows:   

e N Ti =-0.21±0.016
  
e Ti N =-0.728±0.044
  
e Ti Ti =0.048±0.0024
  
log K TiN =- 12   740 T +4.06

Acknowledgments

This work was supported by the research fund of Hanyang University (HY-2013-P).

References
 
© 2015 by The Iron and Steel Institute of Japan

This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs license.
https://creativecommons.org/licenses/by-nc-nd/4.0/
feedback
Top