Method for simulating impurity diffusion in semiconductor with oxidation

A method for simulating a concentration of impurities within a semiconductor device while said semiconductor device is being oxidized, comprising the steps of:

• setting a first triangular mesh configuration in a semiconductor device, said first triangular mesh configuration having first grids arranged in first triangles, one of said first triangles being divided into three first areas each of which subtends a respective first control volume, each of said first control volumes including a circumcenter and one grid of a respective one of said first triangles;
• deforming said first triangular mesh configuration to a second triangular mesh configuration after said semiconductor device is assumed to be oxidized for a time step period, said second triangular mesh configuration having second grids arranged in second triangles being divided into three second areas each of which subtends a respective second control volume, each of said second control volumes including a circumcenter and one grid of a respective one of said second triangles;
• changing impurity concentrations of said second control volumes in accordance with ratios of said second control volumes to respective ones of said first control volumes;
• regenerating a third triangular mesh configuration in said semiconductor device after said impurity concentration changing step is carried out, said third triangular mesh configuration having third grids arranged in third triangles, one of said third triangles being divided into three third areas each of which subtends a respective third control volume, each of said third control volumes including a circumcenter and one grid of a respective one of said third triangles;
• calculating impunity concentrations in said third grids by solving an impurity diffusion equation in accordance with product areas between said third control volumes and said second control volumes and impurity concentrations of said second control volumes.