Monte Carlo Hamiltonian: Linear Potentials

We further study the validity of the Monte Carlo Hamiltonian method. The advantage of the method,in comparison with the standard Monte Carlo Lagrangian approach, is its capability to study the excited states. Weconsider two quantum mechanical models: a symmetric one V(x) = |x|/2; and an asymmetric one V(x) = ∞, forx ＜ 0 and V(x) = x, for x ≥ 0. The results for the spectrum, wave functions and thermodynamical observables are inagreement with the analytical or Runge-Kutta calculations.     

The Properties of Voronoi Tessellations for 2D Ising Fluid Near the Second-Order Magnetic Phase Transition

The topological and geometrical properties of Voronoi cells generated for 2D fluid of hard disks with Ising-like spins near the second-order phase transition from the paramagnetic to ferromagnetic phase are described for different disk densities. The comparison with Voronoi cells generated for the random hard disk system is given.     

The effects of flow split ratio and flow rate in manifolds

This paper reports a combined experimental and numerical investigation of three-dimensional steady turbulent flows in inlet manifolds of square cross-section. Predictions and measurements of the flows were carried out using computational fluid dynamics and laser Doppler anemometry techniques respectively. The flow structure was characterized in detail and the effects of flow split ratio and inlet flow rate were studied. These were found to cause significant variations in the size and shape of recirculation regions in the branches, and in the turbulence levels. It was then found that there is a significant difference between the flow rates through different branches. The performance of the code was assessed through a comparison between predictions and measurements. The comparison demonstrates that all important features of the flow are well represented by the predictions.     

Synthesis and characterization of novel thermoresponsive‐co‐biodegradable hydrogels composed of N‐isopropylacrylamide,poly(L‐lactic acid), and dextran

A series of novel multifunctional hydrogels that combined the merits of both thermoresponsive and biodegradable polymeric materials were designed, synthesized, and characterized. The hydrogels were copolymeric networks composed of N‐isopropylacrylamide (NIPAAM) as a thermoresponsive component, poly(L‐lactic acid) (PLLA) as a hydrolytically degradable and hydrophobic component, and dextran as an enzymatically degradable and hydrophilic component. The chemical structures of the hydrogels were characterized by an attenuated total reflection–Fourier transform infrared spectroscopy (ATR–FTIR) technique. The hydrogels were thermoresponsive, showing a lower critical solution temperature (LCST) at approximately 32 °C, and their swelling properties strongly depended on temperature changes, the balance of the hydrophilic/hydrophobic components, and the degradation of the PLLA component. The degradation of the hydrogels caused by hydrolytic cleavage of ester bonds in the PLLA component was faster at 25 °C below the LCST than at 37 °C above the LCST, determined by the ATR–FTIR technique. Due to their multifunctional properties, the designed hydrogels show great potential for biomedical applications, including drug delivery and tissue engineering. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5054–5066, 2004     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

A novel method for automated grid generation of ice shapes for local‐flow analysis

Modelling a complex geometry, such as ice roughness, plays a key role for the computational flow analysis over rough surfaces. This paper presents two enhancement ideas in modelling roughness geometry for local flow analysis over an aerodynamic surface. The first enhancement is use of the leading‐edge region of an airfoil as a perturbation to the parabola surface. The reasons for using a parabola as the base geometry are: it resembles the airfoil leading edge in the vicinity of its apex and it allows the use of a lower apparent Reynolds number. The second enhancement makes use of the Fourier analysis for modelling complex ice roughness on the leading edge of airfoils. This method of modelling provides an analytical expression, which describes the roughness geometry and the corresponding derivatives. The factors affecting the performance of the Fourier analysis were also investigated. It was shown that the number of sine–cosine terms and the number of control points are of importance. Finally, these enhancements are incorporated into an automated grid generation method over the airfoil ice accretion surface. The validations for both enhancements demonstrate that they can improve the current capability of grid generation and computational flow field analysis around airfoils with ice roughness. Copyright © 2004 John Wiley & Sons, Ltd.     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

Determination of the rheological behavior of epoxy–amine thermosets by dynamic mechanical analysis: Isothermal methods versus nonisothermal methods

The reticulation process of an epoxy resin using an amine as a cure agent was studied at different temperatures and concentrations of the cure agent with dynamic mechanical thermal analysis. The study was performed under both isothermal and nonisothermal conditions, and a temperature–time–transformation diagram was obtained. The measurements from the two modes gave similar results, although the nonisothermal mode required fewer experiments. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 1965–1977, 2003     

Basic study of the gelation of dimethacrylate‐type crosslinking agents

An investigation was made of the gelation of dimethacrylate‐type crosslinking agents in view of an application for separation media. The study mainly centered on a crosslinking agent, glycerol dimethacrylate (GDMA), which is relatively hydrophilic because of a hydroxyl group in the middle of its structure. The gelation of GDMA was compared with that of other hydrophobic crosslinking agents such as ethylene glycol dimethacrylate and 1,6‐hexanediol dimethacrylate. The diluents used in the study were toluene, toluene with methanol, and cyclohexanol. The gelation was observed in real time with a charge coupled device camera and dynamic light scattering (DLS). Also, the separated dry gels were extensively characterized with scanning electron microscopy, BET (N₂ absorption and desorption isotherm), and Fourier transform infrared. DLS analysis showed a stronger molecular interaction of GDMA gelation in toluene, whereas this interaction was much weaker in an alcoholic solvent such as toluene with methanol or cyclohexanol. This indicated that GDMA gelation might proceed through hydrogen bonding as well as a crosslinking reaction of vinyl groups. © 2005 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 44: 949–958, 2006