Pore network simulation of fluid imbibition into paper during coating: II. Characterization of paper's morphology and computation of its effective permeability tensor

3D high spatial resolution visualisation and quantification of interconnectivity in polymer films

A porous network acts as transport paths for drugs through films for controlled drug release. The interconnectivity of the network strongly influences the transport properties. It is therefore important to quantify the interconnectivity and correlate it to transport properties for control and design of new films. This work presents a novel method for 3D visualisation and analysis of interconnectivity. High spatial resolution 3D data on porous polymer films for controlled drug release has been acquired using a focused ion beam (FIB) combined with a scanning electron microscope (SEM). The data analysis method enables visualisation of pore paths starting at a chosen inlet pore, dividing them into groups by length, enabling a more detailed quantification and visualisation. The method also enables identification of central features of the porous network by quantification of channels where pore paths coincide. The method was applied to FIB-SEM data of three leached ethyl cellulose (EC)/hydroxypropyl cellulose (HPC) films with different weight percentages. The results from the analysis were consistent with the experimentally measured release properties of the films. The interconnectivity and porosity increase with increasing amount of HPC. The bottleneck effect was strong in the leached film with lowest porosity.

Characterization of the Interface Between Coating and Fibrous Layers of Paper

Coated paper is an example of a multi-layer porous medium, involving a coating layer along the two surfaces of the paper and a fibrous layer in the interior of the paper. The interface between these two media needs to be characterized in order to develop relevant modeling tools. After careful cutting of the paper, a cross section was imaged using focused ion beam scanning electron microscopy. The resulting image was analyzed to characterize the coating layer and its transition to the fibrous layer. Such image analysis showed that the coating layer thickness is highly variable, with a significant fraction of it being thinner than a minimum thickness required to keep ink from invading into the fibrous layer. The overall structure of the coating and fibrous layers observed in this analysis provide insights into how the system should be modeled, with the resulting conclusion pointing to a specific kind of multi-scale modeling approach.

Nonuniversality of roughness exponent of quasistatic fracture surfaces

Numerous experiments have indicated that the fracture front (in three dimensions) and crack lines (in two dimensions) in disordered solids and rocklike materials is rough. It has been argued that the roughness exponent ζ is universal. Using extensive simulations of a two-dimensional model, we provide strong evidence that if extended correlations and anisotropy-two features that are prevalent in many materials-are incorporated in the models that are used in the numerical simulation of crack propagation, then ζ will vary considerably with the extent of the correlations and anisotropy. The results are consistent with recent experiments that also indicate deviations of ζ from its supposedly universal value, as well as with the data from rock samples.

On the influence of pore structure on the free-imbibition of sessile drops into nanoporous substrates

We report here a model experimental study on the influence of pore structure on the free-imbibition of sessile drops into nanoporous substrates. The work takes advantage of the existence of distinct pore structures on the two sides of a nanoporous alumina membrane: straight parallel channels versus a denser and tortous network. We show first that the spreading which coexists with the free-imbibition predominates in the early stage well follows on both sides the power-law scaling with time predicted by the universal Tanner's law. More interestingly, we found also that the imbibition rate scales in a similar way with the time on both sides of the membrane, showing that the pore structure does not affect qualitatively the free-imbibition kinetics. On the other hand, our results clearly show that the pore structure has a quantitative impact on the imbibition rate, which increases markedly from the A side (dense network of short and tortuous pores) to the side B (straight vertical channels). This latter result shows that, as regards the free-imbibition, the topology of the pores has a preeminent impact on their volume, which is here comparable for both sides of the membrane. More unexpectedly, this quantitative impact of the pore structure on the imbibition rate seems to display a certain sensitivity to the viscosity of the liquid.

Permeability calculations in three-dimensional isotropic and oriented fiber networks

Hydraulic permeabilities of fiber networks are of interest for many applications and have been studied extensively. There is little work, however, on permeability calculations in three-dimensional random networks. Computational power is now sufficient to calculate permeabilities directly by constructing artificial fiber networks and simulating flow through them. Even with today's high-performance computers, however, such an approach would be infeasible for large simulations. It is therefore necessary to develop a correlation based on fiber volume fraction, radius, and orientation, preferably by incorporating previous studies on isotropic or structured networks. In this work, the direct calculations were performed, using the finite element method, on networks with varying degrees of orientation, and combinations of results for flows parallel and perpendicular to a single fiber or an array thereof, using a volume-averaging theory, were compared to the detailed analysis. The detailed model agreed well with existing analytical solutions for square arrays of fibers up to fiber volume fractions of 46% for parallel flow and 33% for transverse flow. Permeability calculations were then performed for isotropic and oriented fiber networks within the fiber volume fraction range of 0.3%-15%. When drag coefficients for spatially periodic arrays were used, the results of the volume-averaging method agreed well with the direct finite element calculations. On the contrary, the use of drag coefficients for isolated fibers overpredicted the permeability for the volume fraction range that was employed. We concluded that a weighted combination of drag coefficients for spatially periodic arrays of fibers could be used as a good approximation for fiber networks, which further implies that the effect of the fiber volume fraction and orientation on the permeability of fiber networks are more important than the effect of local network structure.

Construction by Molecular Dynamics Modeling and Simulations of the Porous Structures Formed by Dextran Polymer Chains Attached on the Surface of the Pores of a Base Matrix: Characterization of Porous Structures

Significant increases in the separation of bioactive molecules by using ion-exchange chromatography are realized by utilizing porous adsorbent particles in which the affinity group/ligand is linked to the base matrix of the porous particle via a polymeric extender. To study and understand the behavior of such systems, the M3B model is modified and used in molecular dynamics (MD) simulation studies to construct porous dextran layers on the surface of a base matrix, where the dextran polymer chains and the surface are covered by water. Two different porous polymer layers having 25 and 40 monomers per main polymer chain of dextran, respectively, are constructed, and their three-dimensional (3D) porous structures are characterized with respect to porosity, pore size distribution, and number of conducting pathways along the direction of net transport. It is found that the more desirable practical implications with respect to structural properties exhibited by the porous polymer layer having 40 monomers per main polymer chain, are mainly due to the higher flexibility of the polymer chains of this system, especially in the upper region of the porous structure. The characterization and analysis of the porous structures have suggested a useful definition for the physical meaning and implications of the pore connectivity of a real porous medium that is significantly different than the artificial physical meaning associated with the pore connectivity parameter employed in pore network models and whose physical limitations are discussed; furthermore, the methodology developed for the characterization of the three-dimensional structures of real porous media could be used to analyze the experimental data obtained from high-resolution noninvasive three-dimensional methods like high-resolution optical microscopy. The MD modeling and simulations methodology presented here could be used, considering that the type and size of affinity group/ligand as well as the size of the biomolecule to be adsorbed onto the affinity group/ligand are known, to construct different porous dextran layers by varying the length of the polymeric chain of dextran, the number of attachment points to the base matrix, the degree of side branching, and the number of main polymeric chains immobilized per unit surface area of base matrix. After the characterization of the porous structures of the different porous dextran layers is performed, then only a few promising structures would be selected for studying the immobilization of adsorption sites on the pore surfaces and the subsequent adsorption of the bioactive molecules onto the immobilized affinity groups/ligands.