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     

Conversion estimates for gas-solid reactions: Improved estimates for irreversible systems

Solid conversion has been of much interest to research workers and the Cumulative Gas Concentration Method (CGCM) has been developed to provide bounds for the reactant solid in irreversible systems having power-law type kinetics. This CGCM is used in conjunction with the Pseudo-Steady-State (PSS) assumption and these bounds appear to be good under certain conditions and for long times.

In this work, the PSS condition is not assumed when deriving the bounds on the cumulative gas concentration and the effect of the inclusion of the transient part is analysed. A technique, based on an iterative scheme incorporating the CGCM, is developed to provide improved estimates for the reactant solid.     

Estimation of eardrum acoustic pressure and of ear canal length from remote points in the canal

Sound pressure distributions in the human ear canal, whether unoccluded or occluded with ear molds, were studied using a probe tube technique. On average, for frequencies below 6 kHz, the measuring probe tube had to be placed within 8 mm of the vertical plane containing the top of the eardrum (TOD), determined optically, in order to obtain sound pressure magnitudes within 6 dB of "eardrum pressure." To obtain that accuracy in all of the eight subjects studied, the probe had to be within 6 mm of the TOD. Since probe location relative to the drum has to be known, a purely acoustic method was developed which can be conveniently used to localize the probe-tip position, utilizing the standing wave property of the sound pressure in the ear canal. The acoustically estimated "drum location" generally lay between the optically determined vertical planes containing the TOD and the umbo. On average, the "drum location" fell 1 mm medial to the TOD. Of the 32 estimates made acoustically in various occluded and unoccluded conditions in 14 subjects, 30 estimates lay within a +/- 2-mm range of this average.     

Optimal birth control of population dynamics

The authors studied optimal birth control policies for an age-structured population of McKendrick type which is a distributed parameter system involving 1st order partial differential equations with nonlocal bilinear boundary control. The functional analytic approach of Dubovitskii and Milyutin is adopted in the investigation. Maximum principles for problems with a free end condition and fixed final horizon are developed, and the time optimal control problems, the problem with target sets, and infinite planning horizon case are investigated.     

Resistivity exponent of two-dimensional lattice animals

We calculate the average resistanceR(L) of lattice animals spanningL×L cells on the square lattice using exact and Monte Carlo methods. The dynamical resistivity exponent, defined asR(L) L , is found to be =1.36±0.07. This contradicts the Alexander-Orbach conjecture, which predicts 0.8. Our value for differs from earlier measurements of this quantity by other methods yielding =1.17±0.05 and 1.22±0.08 by Havlin et al.On leave from the Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing, China.