Multiphysical modeling and meshless simulation of electric‐sensitive hydrogels

According to a multiphase mixture theory, we have mathematically developed a multiphysical model with chemoelectromechanical coupling considerations, termed the multieffect‐coupling electric‐stimulus (MECe) model, to simulate the responsive behavior of electric‐sensitive hydrogels immersed in a bath solution under an externally applied electric field. For solutions of the MECe model consisting of coupled nonlinear partial differential governing equations, a meshless Hermite–Cloud method with a hierarchical iteration technique has been used for a one‐dimensional steady‐state analysis of a hydrogel strip. The computed results are compared with the experimental data, and there is very good agreement. Simulations within the domains of both hydrogels and surrounding solutions also present distributions of the ionic concentrations and electric potential as well as the hydrogel displacement. The effects of various physical parameters on the response behavior of electric‐stimulus responsive hydrogels are discussed in detail. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1514–1531, 2004     

聚合物单晶生长的三维相场模拟研究

利用元胞自动机方法与相场模型的结合建立新型三维模拟相场模型.同时,为模拟真实的、三维的高分子结晶的过程,采用元胞自动机方法离散方程,且元胞几何形状的选取符合真实聚合物晶格扩散方式的物理规律,以及新建立的相场模型套用间规聚丙烯的实验参数.利用该模型模拟了多种三维立方体或者薄层的晶体形貌及其相互之间的演化过程,包括正方形、长方形、菱形、六边形、多层单晶等.通过模拟结果与真实形貌作对比来证明所建立的相场模型真实可靠性.     

Robust and high‐resolution simulations of nonlinear electrokinetic processes in variable cross‐section channels

We present a model and an associated numerical scheme to simulate complex electrokinetic processes in channels with nonuniform cross‐sectional area. We develop a quasi‐1D model based on local cross‐sectional area averaging of the equations describing unsteady, multispecies, electromigration‐diffusion transport. Our approach uses techniques of lubrication theory to approximate electrokinetic flows in channels with arbitrary variations in cross‐section; and we include chemical equilibrium calculations for weak electrolytes, Taylor–Aris type dispersion due of nonuniform bulk flow, and the effects of ionic strength on species mobility and on acid–base equilibrium constants. To solve the quasi‐1D governing equations, we provide a dissipative finite volume scheme that adds numerical dissipation at selective locations to ensure both unconditional stability and high accuracy. We couple the numerical scheme with a novel adaptive grid refinement algorithm that further improves the accuracy of simulations by minimizing numerical dissipation. We benchmark our numerical scheme with existing numerical schemes by simulating nonlinear electrokinetic problems, including ITP and electromigration dispersion in CZE. Simulation results show that our approach yields fast, stable, and high‐resolution solutions using an order of magnitude less grid points compared to the existing dissipative schemes. To highlight our model's capabilities, we demonstrate simulations that predict increase in detection sensitivity of ITP in converging cross‐sectional area channels. We also show that our simulations of ITP in variable cross‐sectional area channels have very good quantitative agreement with published experimental data.     

Compact adaptive-grid scheme for high numerical resolution simulations of isotachophoresis

In a previous publication we demonstrated a fast simulation tool for solution of electrophoretic focusing and separation. We here describe the novel mathematical model and numerical algorithms used to create this code. These include the representation of advection–diffusion equations on an adaptive grid, high-resolution discretization of the equations (sixth order compact), a new variational-based approach for controlling the motion of grid points, and new boundary conditions which enable solution in a moving frame of reference. We discuss the advantages of combining a high-resolution discretization with an adaptive grid in accurately resolving sharp interfaces in isotachophoresis, and provide verification against known analytical solutions and comparison with prevailing exiting numerical algorithms.     

Meshless simulation of equilibrium swelling/deswelling of pH‐sensitive hydrogels

Hydrogels have been widely used in microelectromechanical systems (MEMS) and Bio‐MEMS devices. In this article, the equilibrium swelling/deswelling of the pH‐stimulus cylindrical hydrogel in the microchannel is studied and simulated by the meshless method. The multi‐field coupling model, called multi‐effect‐coupling pH‐stimulus (MECpH) model, is presented and used to describe the chemical field, electric field, and the mechanical field involved in the problem. The partial differential equations (PDEs) describing these three fields are either nonlinear or coupled together. This multi‐field coupling and high nonlinear characteristics produce difficulties for the conventional numerical methods (e.g., the finite element method or the finite difference method), so an alternative—meshless method is developed to discretize the PDEs, and the efficient iteration technique is adopted to solve the nonlinear problem. The computational results for the swelling/deswelling diameter of the hydrogel under the different pH values are firstly compared with experimental results, and they have a good agreement. The influences of other parameters on the mechanical properties of the hydrogel are also investigated in detail. It is shown that the multi‐field coupling model and the developed meshless method are efficient, stable, and accurate for simulation of the properties of the stimuli‐sensitive hydrogel. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 326–337, 2006     

Computer‐Aided Estimation of Diffusion Coefficients in Non‐Solvent/Polymer Systems

A novel method for estimating the mutual and self‐diffusion coefficients of a non‐solvent/polymer system is proposed in this work. The idea is to study the evaporation process from non‐solvent/solvent/polymer systems as a one‐dimensional numerical experiment and to use polymer solution weight versus time data to fit the unknown parameters of the diffusion‐coefficient correlations based on free‐volume theory. For this purpose, the evaporation process is modeled as a coupled heat‐ and mass‐transfer problem with a moving boundary, and the Galerkin finite‐element method is used to solve simultaneously the non‐linear governing equations. This method is successfully applied to the estimation of water–cellulose acetate diffusion coefficients and is valid over the whole range of temperatures and concentrations for practical applications in membrane technology. Additionally, there is a detailed discussion on if water affects the morphology of the final cellulosic membrane by studying the concentration profiles of the constituents of the casting solution.     

Modeling of penetrant diffusion in glassy polymers with an integral sorption deborah number

A mathematical model was developed to explain the anomalous penetrant diffusion behavior in glassy polymers. The model equations were derived by using the linear irreversible thermodynamics theory and the kinematic relations in continuum mechanics, showing the coupling between the polymer mechanical behavior and penetrant transport. The Maxwell model was used as the stress–strain constitutive equation, from which the polymer relaxation time was defined. An integral sorption Deborah number was proposed as the ratio of the characteristic relaxation time in the glassy region to the characteristic diffusion time in the swollen region. With this definition, an integral sorption process was characterized by a single Deborah number and the controlling mechanism was identified in terms of the value of the Deborah number. The model equations were two coupled nonlinear differential equations. A finite difference method was developed for solving the model equations. Numerical simulation of integral sorption of penetrants in glassy polymers was performed. The simulation results show that (1) the present model can predict Case II transport behavior as well as the transition from Case II to Fickian diffusion and (2) the integral sorption Deborah number is a major parameter affecting the transition. © 1993 John Wiley & Sons, Inc.     

Attaining exponential convergence for the flux error with second- and fourth-order accurate finite-difference equations. II. Application to systems comprising first-order chemical reactions

This article demonstrates that exponential convergence of the flux error can be achieved for any kinetic-diffusion system comprising an arbitrary number of (pseudo) first-order chemical reactions if the underlying PDEs are discretized as outlined for the box 2 or box 4 method in the preceding part of this article. By investigating the eigenvalues and eigenvectors of the first-order kinetic coupling matrix in general form the present article demonstrates that the simulation of any multispecies first-order kinetic diffusion system can be as accurately done as the simulation of a single representative one-species system. The Fourier coefficients governing the error level of the flux are much smaller in the limiting case of kinetic control as those reported in the preceding article for the limiting case of diffusion control. The higher rate of exponential convergence predicted on the basis of the mathematical model has been fully verified by the numerical results.     

Monte Carlo simulation of diffusion effects on chain‐extension reactions

The diffusion effects on chain‐extension reactions using carboxyl‐terminated polyamide‐12 as a model reactant with bisoxazolines were investigated by the stochastic Monte Carlo method. Thus, complicated direct modeling and numerical calculations were avoided. The chain‐length dependence and detailed diffusive behavior were discussed in depth. The diffusion effects retarded the progress of chain‐extension reactions and led to lower coupling efficiency. The simulated results indicated that the diffusion effects could make the final molecular weight distributions wider. In the presence of diffusion and with the progress of the coupling efficiency, peaks in the evolution curves of the weight‐average molecular weight and valleys in the evolution curves of the polydispersity index were observed, respectively, when the coupling efficiency was low enough. These phenomena were different from those without diffusion effects and were analyzed in detail. The critical entanglement chain length had strong effects on the simulated results of the diffusion effects, especially when its value was near the average chain length. The results also showed that the effects of the reactant degradation made the molecular weight distribution of the reaction system wider and weakened the diffusion effects on the coupling reaction. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 2902–2911, 2006     

Quasi‐Spectral Method for the Solution of the Master Equation for Unimolecular Reaction Systems

Rate constants of elementary reactions involving unimolecular steps can be calculated from molecular data in a most general way by solving appropriate master equations. The conventional numerical solution requires rather a fine discretization applied over a sufficiently large energy range to achieve a reasonable accuracy. This leads to linear but very high‐dimensional systems of differential equations. We propose a quasi‐spectral method that uses Gaussian radial basis functions to establish a low‐dimensional linear model to speed up the numerical integration. The combination with an iterative adaptation provides a further improvement of computational efficiency. The suggested approach is illustrated and exemplified by means of the unimolecular decomposition of 2,3‐dihydro‐2,5‐dimethylfuran‐3‐yl, an intermediate radical occurring in the pyrolysis and oxidation of 2,5‐dimethylfuran. A comparison of the conventional and the proposed method is presented to validate the novel approach and to demonstrate its performance.     

Range of Validity of a Simplified Model for Diffuse Charge Dynamics

The modelling of electrochemical processes often requires the solution of the Poisson‐Nernst‐Planck (PNP) equations. In complex geometries, such as porous electrodes, that is challenging due to the presence of disparate length scales, ranging from the Debye screening length (～nm) to the device length scale (～cm). To overcome this difficulty, one often assumes that the electric double layer (EDL) is at quasi‐equilibrium to construct a simplified model that accounts for ion diffusion in the electro‐neutral bulk of the electrolyte while replacing the EDLs with appropriate boundary conditions. Various researchers have demonstrated that such an approach is valid in the asymptotic limit of a thin EDL and moderate electrode potentials. In this note, we explore the range of validity of this approximation by considering a one‐dimensional electrolytic cell with blocking electrodes subjected to a step change and time‐periodic alternations in the electrodes’ potentials by calculating the errors associated with the approximate approach as functions of the EDL thickness and electric field frequency and intensity. Additionally, we delineate numerical instabilities associated with the numerical solutions of the bulk equations with the nonlinear boundary condition peculiar to this problem.     

Soliton‐like regimes,echo and concave spiral waves in mathematical models of biological excitable media

On the basis of numerical simulations with partial reaction‐diffusion equations describing the dynamics of electrical processes in biological excitable tissues, such as a nerve axons and cardiac Purkinje fibers, we have revealed unusual – soliton‐like – regimes of interaction of propagating nonlinear excitation waves: reflection of the colliding waves instead of their annihilation. The specific nature of wave dynamics due to the reflection effects in two‐dimensional excitable media is considered.     

Improving the accuracy of the spatial discretization in finite-difference electrochemical kinetic simulations, by means of the extended Numerov method

The fourth-order accurate, three-point finite-difference Numerov spatial discretization provides accurate and efficient solutions to the time-dependent governing differential equations of electrochemical kinetics in one-dimensional space geometry, when the equations contain first time derivatives of the solution, second spatial derivatives, and homogeneous reaction terms only. However, the original Numerov discretization is not applicable when the governing equations involve first spatial derivative terms. To overcome this limitation, an appropriately extended Numerov discretization is required. We examine the utility of one of such extensions, first described by Chawla. Relevant discrete formulae are outlined for systems of linear governing equations involving first derivative terms, and applied to five representative example models of electrochemical transient experiments. The extended Numerov discretization proves to have an accuracy and efficiency comparable to the original Numerov scheme, and its accuracy is typically up to four orders of magnitude higher, compared to the conventional, second-order accurate spatial discretization, commonly used in electrochemistry. This results in a considerable improvement of efficiency. Therefore, the application of the extended Numerov discretization to the electrochemical kinetic simulations can be fully recommended.     

A mathematical model for stress-driven diffusion in polymers

A model for case II diffusion into polymers is presented. The addition of stress terms to the Fickian flux is used to produce the characteristics progressive front. The stress in turn obeys a concentration-dependent evolution equation. The model equations are analyzed in the limit of small diffusivity for the problem of penetration into a semiinfinite medium. Provided that the coefficient functions obey two monotonicity conditions, the solvent concentration profile is shown to have a steep front that progresses into the medium. The formulas governing the progression of the front are developed. After the front decays away, the long time behavior of the solution is shown to be a similarity solution as in Fickian diffusion. Two techniques for approximating the solvent concentration and the front position are presented. The first approximation method is a series expansion; formulas are given for the initial speed and deceleration of the front. The second approximation method uses a portion of the long time similarity solution to represent the short time solution behind the front.     

Dynamic evolution of gases in the γ‐ and helium‐ion radiolyses of solid polymers

The dynamic evolution of gaseous hydrogen, methane, and carbon dioxide in the γ‐ and ⁴He‐ion radiolyses of solid polymers was investigated. The polymers used include low‐density and high‐density polyethylene, polypropylene, polystyrene, poly(methyl methacrylate), Nylon 11, Nylon 6, and poly(dimer acid‐co‐alkyl polyamine). An inline quadrupole mass spectrometer was utilized to monitor the dynamic profiles of the gases produced in the radiolysis. One‐ and two‐dimensional numerical diffusion models were developed to simulate and extract optimum diffusion coefficients and gas yields from the experimental dynamic gas profiles. It was found that the dynamic evolution of molecular hydrogen from the bulk polymer is controlled by its diffusion in most cases, such as CO₂ in poly(methyl methacrylate). In the γ radiolysis of some polymers such as low‐density polyethylene and polypropylene, the dynamic evolution of methane is only partially controlled by the diffusion process, and some other postirradiation process is a factor. It is concluded that the simulation method developed in this article is helpful in understanding and predicting the mechanisms of gas evolution in the radiolysis of solid polymers. © 2001 John Wiley & Sons, Inc. J Polym Sci B Part B: Polym Phys 39: 1449–1459, 2001     

A New Numerical Method for Accurate Simulation of Fast Cyclic Adsorption Processes

In the simulation of fast cyclic adsorption processes, to apply the Fickian diffusion model it is necessary to include an increasing number of numerical discretization points as the cycle time is reduced in comparison to the characteristic diffusional time constant. We propose a new numerical method based on the definition of two distinct regions within an adsorbent particle: an outer layer where the concentration varies significantly with large internal gradients leading to enhanced mass fluxes, and an internal region where the concentration profile is virtually flat. The proposed method leads to the automated generation of a numerical grid that has a constant number of elements independent of the process cycle time. The procedure is demonstrated on a model for the simulation of a heatless dryer pressure swing adsorption process.     

Dynamics in the crystalline polymorphic forms I and II and form III of isotactic poly‐1‐butene

Following an earlier study of the ¹H relaxation and NMR line shapes, we have carried out selective one‐dimensional and two‐dimensional ¹³C solid‐state NMR studies that yield to detailed interpretation of the dynamics in form I, II, and III polymorphs of isotactic poly‐1‐butene. A specific defect diffusion along the side group is proposed to account for the temperature dependence of the ¹³C spectra in form I. The backbone of the helix in forms II and III is shown to undergo large angle motions above the glass‐transition temperature. High‐resolution solid‐state ¹³C two‐dimensional exchange NMR under magic‐angle spinning with cross‐polarization techniques demonstrates the existence of slow rotational jumps of the helices in form III with typical jump rates of about 10 s⁻¹. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 2611–2624, 2000     

Diffusion behavior in polymer–clay nanocomposites

This paper reviews our previous studies on the diffusion behavior in polymers clay nanocomposites. A geometric model for predicting the effective diffusivity through this type of systems as a function of clay sheets orientation, volume fraction, polymer clay interaction, and aspect ratio is proposed. Model predictions are compared to the effective diffusivity generated using random walk simulations as well as with predictions obtained from already existing theoretical models. Fair agreement is found between the model prediction and the results of numerical simulations. With respect to the already existing theoretical models, the present mathematical derivation seems more adequate to describe diffusion behavior in conventional nanocomposites systems (i.e. when fillers present very low values of volume to surface ratio). Experimental diffusion tests are discussed and interpreted with the aid of the proposed model. In addition to the aspect ratio and clay concentration, the polymer clay interactions as well as the sheets orientation are the factors controlling the barrier properties of polymer‐layered silicate nanocomposites. Good agreement was found in the case of samples containing exfoliated clay, whereas the model fails in the case of micro‐composites, in which the inorganic lamellae are agglomerated in clusters. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 265–274, 2006     

Model for anomalous moisture diffusion through a polymer–clay nanocomposite

Experimental data are reported on moisture diffusion and the elastoplastic response of an intercalated nanocomposite with vinyl ester resin matrix and montmorillonite clay filler at room temperature. Observations in diffusion tests showed that water transport in the neat resin is Fickian, whereas it becomes anomalous (non‐Fickian) with the growth of the clay content. This transition is attributed to immobilization of penetrant molecules on the surfaces of hydrophilic clay layers. Observations in uniaxial tensile tests demonstrate that the response of vinyl ester resin is strongly elastoplastic, whereas an increase in the clay content results in a severe decrease of plastic strains observed as a noticeable reduction in the curvatures of the stress‐strain diagrams. This is explained by slowing down the molecular mobility in the host matrix driven by confinement of chains in galleries between platelets. Constitutive equations are developed for moisture diffusion through and the elastoplastic behavior of a nanocomposite. Adjustable parameters in these relations are found by fitting the experimental data. Fair agreement is demonstrated between the observations and the results of numerical simulation. A striking similarity is revealed among changes in diffusivity, ultimate water uptake, and the rate of plastic flow with an increased clay content. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 476–492, 2003