Evolution of Micropores in Polymer-Polymer Composite Thin Films
Theoretical Foundations
Combined surface diffusion and evaporation-condensation problem
Define:
j = volume of matter added to per unit area of the solid surface per unit time
p = free energy reduction associated with per unit volume of matter added to per unit surface area of the solid
m = kinetic parameter in the model
Mass conservation relates the surface velocity to the fluxes of the two matter transport processes
i = the volume of matter added to per unit area of the solid surface
I = mass displacement
G = free energy
The weak statement is:
This is a weak statement including both surface diffusion and evaporation-condensation. To formulate a method for finite elements, it is convenient to use I and rn as the basic variables. Eliminating the other variables in the above weak statement, we can rewrite the weak statement as:
Finite Element Method
To approximate a solution of J and vn, the above weak statement needs to be satisfied for all virtual movements. However, to implement this process for the finite element method, this requirement can be relaxed such that the weak statement be satisfied for only a family of virtual motions.
To model the structure, we use a certain number of straight/linear elements connected by nodes. Each segment is an element, and two neighboring elements meet at a node. The motions of the nodes constitute the family of virtual motions. The length of each element can be defined arbitrarily.
In each element, we can define the two nodes based on a Cartesian coordinate system; the two end points of the element, or the two nodes, can be written as (x1, y1) and (x2, y2). The element has length, l. At the mid-point of the element, we can define that midpoint as s. Due to the surface diffusion process, the two nodal points of the element can move. As a result, we can model the incremental movement by the following equation.
Where the interpolation coefficients can be defined by the following set of equations.
Similarly, the incremental movements of the nodal points can be defined by its corresponding nodal velocities at each direction.
Due to the mass flux vector, we can model the virtual movement of the mass displacement due to diffusion at the nodes to define
The interpolation equations can be modeled by a set of quadratic equations atisfying the boundary value problem set forth for each nodal coordinate.
This same interpolation is used for the mass flux J.
From the above relations, we can formulate the element
For each element, the integration gives a bilinear form such that
and the element He matrix is
where,
Furthermore, we can specify the driving force due to surface tension acting on each element by the following equations.
With the elemental matrices formed, the finite element process assembles the global H and f matrices such that the global nodal velocities can be computed. The governing relation in terms of global coordinates can be written as the following.
The velocity at each node can then be determined and extracted.
Once the velocity has been extracted, the Runge Kutta method can be utilized to predict the movement of the nodal points. The Runge Kutta method can then be specified for a user-specified time (and time steps) to model the behavior of the surface diffusion and evaporation process.