Development of the new method to investigate the phase transformations based on the computational scientific approachs.

(1)基于非线性演化方程的相分解计算机模拟（Miyazaki）。实验研究了实际合金体系（Al-Zn、Cu-Co、Fe-Mo、Fe-Al-Co、Fe-Cr和III-V化合物半导体等）相分解的时间演化过程。采用相场法对这些合金体系的微观组织变化进行了动态模拟。计算得到的微观组织形态特征与实验结果吻合较好。(2)临界现象（Mohri）的计算分析。采用聚类变分法和路径概率法对相变点附近有序-无序相变的磁化率和弛豫时间进行了数值研究。当转变类型为一阶时，磁化率向纵向无序温度发散。众所周知，临界慢化发生在二阶相变点。通过计算发现，即使在一阶相变中也会发生类似的现象。在过渡点附近的自由能的平坦性将是上述现象的起源。(3)马氏体相变（Ohta）橡胶样行为的理论研究。成功地模拟了马氏体相变的橡胶样行为。在考虑孪晶界运动的基础上，从相场型演化方程推导出应力应变场的本构方程。通过数值模拟预测了应力应变曲线及其频率依赖关系。

(1) Computer simulation of the phase decomposition based on the non-linear evolution equation (Miyazaki).The time evolution of the phase decomposition in the real alloy system (such as Al-Zn, Cu-Co, Fe-Mo, Fe-Al-Co, Fe-Cr and III-V compound semiconductor, etc.) was investigated experimentally. The dynamics of the microstructure changes for those alloy systems were simulated on the basis of the Phase-field Method developed by the present research. The morphological feature of the calculated microstructures were in good agreement with the experimental results.(2) Computational analysis of the critical phenomena (Mohri).The susceptibility and relaxation time for the order-disorder phase transition near the transition point were investigated by the Cluster Variation and Path Probability Method, numerically. The susceptibility diverged toward the spinodal disordering temperature when the transition type is first order. It is well known that the critical slowing down will occur at the second order phase transition point. We found that the similar phenomenon will take place even for the first order phase transition by our calculation. The flatness of the free energy at the vicinity of the transition point would be the origin for the phenomena above mentioned.(3) Theoretical investigation for the rubber-like behavior of Martensitic transformation (Ohta).The rubber-like behavior of Martensitic transformation was simulated, successfully. The constitutive equations for the stress and strain field were deduced from the phase field type evolution equations with considering the twin boundary motion. The stress-strain curve and its frequency dependence were predicted by the numerical simulation.