Bookmark and Share FiPy: A Finite Volume PDE Solver Using Python
Version 3.0.1-dev139-ge5d2233

This Page

Contact

FiPy developers
Jonathan Guyer
Daniel Wheeler
James Warren

Join our mailing list

100 Bureau Drive, M/S 6555
Gaithersburg, MD 20899

301-975-5329 Telephone
301-975-4553 Facsimile

examples.levelSet.electroChem.goldΒΆ

Model electrochemical superfill of gold using the CEAC mechanism.

This input file is a demonstration of the use of FiPy for modeling gold superfill. The material properties and experimental parameters used are roughly those that have been previously published [NIST:damascene:2005].

To run this example from the base fipy directory type:

$ python examples/levelSet/electroChem/gold.py

at the command line. The results of the simulation will be displayed and the word finished in the terminal at the end of the simulation. The simulation will only run for 10 time steps. To run with a different number of time steps change the numberOfSteps argument as follows,

>>> runGold(numberOfSteps=10, displayViewers=False) 
1

Change the displayViewers argument to True if you wish to see the results displayed on the screen. This example has a more realistic default boundary layer depth and thus requires gmsh to construct a more complex mesh.

There are a few differences between the gold superfill model presented in this example and in examples.levelSet.electroChem.simpleTrenchSystem. Most default values have changed to account for a different metal ion (gold) and catalyst (lead). In this system the catalyst is not present in the electrolyte but instead has a non-zero initial coverage. Thus quantities associated with bulk catalyst and catalyst accumulation are not defined. The current density is given by,

i = \frac{c_m}{c_m^{\infty}} \left( b_0 + b_1 \theta \right).

The more common representation of the current density includes an exponential part. Here it is buried in b_0 and b_1. The governing equation for catalyst evolution includes a term for catalyst consumption on the interface and is given by

\dot{\theta} = J v \theta - k_c v \theta

where k_c is the consumption coefficient (consumptionRateConstant). The trench geometry is also given a slight taper, given by taperAngle.

If the MayaVi plotting software is installed (see Installation) then a plot should appear that is updated every 10 time steps and will eventually resemble the image below.

catalyst coverage as a function of time during gold superfill