NetLogo Models Library:
Sample Models/Chemistry & Physics/Materials Science

(back to the library)

Solid Diffusion

If you download the NetLogo application, this model is included. You can also Try running it in NetLogo Web

WHAT IS IT?

This model describes how diffusion occurs between two adjacent solids.

Diffusion is one of the most important phenomena in fields such as biology, chemistry, geology, chemistry, engineering and physics. Interestingly, before becoming a famous for the Relativity Laws, Albert Einstein wrote extensively about diffusion, and was one of the first to connect diffusion to the Brownian motion of atoms.

Diffusion can take place in gases, liquids, or solids. In solids, particularly, diffusion occurs due to thermally-activated random motion of atoms - unless the material is at absolute zero temperature (zero Kelvin), individual atoms keep vibrating and eventually move within the material. One of the possible net effects of diffusion is that atoms move from regions of high concentration of one element to regions with low concentration, until the concentration is equal throughout the sample.

This model demonstrates a solid diffusion couple, such as copper and nickel. In a real laboratory, such experiment would take place at very high temperatures, for the process to take place in a reasonable amount of time (note that the diffusion coefficient varies exponentially with the inverse of the temperature). There are many mechanisms for diffusion in solids. In this model we demonstrate one of them, which is caused by missing atoms in the metal crystal. The locations, of the missing atoms are often called vacancies. Therefore, this type of diffusion mechanism is referred to as "vacancy diffusion". The extent to which the diffusion can happen depends on the temperature and the number of vacancies in the crystal.

In addition, there are various other conditions that are needed for solid diffusion to occur. Some examples of these are similar atomic size, similar crystal structure, and similar electronegativity. This model assumes all of these conditions are present.

HOW IT WORKS

There are two types of atoms, green and blue. At the beginning, all green atoms are on the left and the blue atoms are on the right. All the vacancies start out between the two metals. As atoms move into vacancies, the vacancies disperse. In most real-world scenarios, vacancies are scattered in the material to begin with. In this model, for simplification purposes, we assume that the materials have no vacancies in the beginning, and that all the vacancies start off in between the two materials.

In this model we also assume that the heat is evenly distributed throughout the metals. Therefore, each atom has an equal chance of breaking bonds with its neighbors and moving to a vacancy.

HOW TO USE IT

To run the model, first press the SETUP button, then press the GO button.

"Atoms by Column" is a distribution diagram of the two atom types. The other graph is a maximum diffusion distance, squared, versus time. If the model runs long enough, this plot will show an approximately linear relationship between the squared distance and time, following the known equation (for one-dimensional diffusion):

> x² = 2 * D * t

where x is the maximum diffusion distance, D is the diffusion coefficient, and t is elapsed time.

THINGS TO NOTICE

If you run the model for a few hundred ticks, the distribution graph should look like two interleaving curves. The far edges remain purely one color, while the middle is about 50-50.

The other graph should be generally linear. The "diffusion coefficient" of the system is proportional to the slope, and can be easily calculated using the above equation.

THINGS TO TRY

Let the model run for a long time. (You can use the speed slider to make the model run faster.) Do you think the metal will ever become completely diffused?

Try increasing the dimensions of the world. Does the behavior change at all?

EXTENDING THE MODEL

The model uses a very simple initial state in which there is always exactly one column of vacancies and they are all located in the middle. Try adding settings that dictate how many vacancies there are and where they start out.

Give the two metals, or the two sides of the world, different characteristics. For example, a temperature difference could be simulated by making atomic movements on one side happen less often than on the other.

Try changing the crystal structure of the atoms. In close-packed atoms in two dimensions, atoms actually have six neighbors (hexagonal) instead of four (square).

NETLOGO FEATURES

This model uses a non-wrapping world.

RELATED MODELS

MaterialSim Grain Growth GasLab Two Gas

CREDITS AND REFERENCES

Thanks to James Newell for his work on this model.

For additional information:

Porter, D.A., and Easterling, K.E., Phase Transformations in Metals and Alloys, 2nd ed., Chapman & Hall, 1992

Shewmon, P.G., Diffusion in solids, 2nd ed., TMS, 1989

HOW TO CITE

If you mention this model or the NetLogo software in a publication, we ask that you include the citations below.

For the model itself:

• Wilensky, U. (2007). NetLogo Solid Diffusion model. http://ccl.northwestern.edu/netlogo/models/SolidDiffusion. Center for Connected Learning and Computer-Based Modeling, Northwestern University, Evanston, IL.

Please cite the NetLogo software as:

• Wilensky, U. (1999). NetLogo. http://ccl.northwestern.edu/netlogo/. Center for Connected Learning and Computer-Based Modeling, Northwestern University, Evanston, IL.

COPYRIGHT AND LICENSE

Copyright 2007 Uri Wilensky.

This work is licensed under the Creative Commons Attribution-NonCommercial-ShareAlike 3.0 License. To view a copy of this license, visit https://creativecommons.org/licenses/by-nc-sa/3.0/ or send a letter to Creative Commons, 559 Nathan Abbott Way, Stanford, California 94305, USA.

Commercial licenses are also available. To inquire about commercial licenses, please contact Uri Wilensky at uri@northwestern.edu.

(back to the NetLogo Models Library)