Atomistic packing models for experimentally investigated swelling states induced by CO2 in glassy polysulfone and poly(ether sulfone)

Atomistic packing models have been created, which help to better understand the experimentally observed swelling behavior of glassy polysulfone and poly (ether sulfone), under CO₂ gas pressures up to 50 bar at 308 K. The experimental characterization includes the measurement of the time‐dependent volume dilation of the polymer samples after a pressure step and the determination of the corresponding gas concentrations by gravimetric gas‐sorption measurements. The models obtained by force‐field‐based molecular mechanics and molecular dynamics methods allow a detailed atomistic analysis of representative swelling states of polymer/gas systems, with respect to the dilation of the matrix. Also, changes of free volume distribution and backbone mobility are accessible. The behavior of gas molecules in unswollen and swollen polymer matrices is characterized in terms of sorption, diffusion, and plasticization. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 1874–1897, 2006     

A numerical study of an unsteady laminar flow in a doubly constricted 3D vessel

Unsteady flow dynamics in doubly constricted 3D vessels have been investigated under pulsatile flow conditions for a full cycle of period T. The coupled non‐linear partial differential equations governing the mass and momentum of a viscous incompressible fluid has been numerically analyzed by a time accurate Finite Volume Scheme in an implicit Euler time marching setting. Roe's flux difference splitting of non‐linear terms and the pseudo‐compressibility technique employed in the current numerical scheme makes it robust both in space and time. Computational experiments are carried out to assess the influence of Reynolds' number and the spacing between two mild constrictions on the pressure drop across the constrictions. The study reveals that the pressure drop across a series of mild constrictions can get physiologically critical and is also found to be sensitive both to the spacing between the constrictions and the oscillatory nature of the inflow profile. The flow separation zone on the downstream constriction is seen to detach from the diverging wall of the constriction leading to vortex shedding with 3D features earlier than that on the wall in the spacing between the two constrictions. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

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.     

Permeability,permselectivity, and penetrant‐induced plasticization in fluorinated polyimides studied by positron lifetime measurements

The ortho‐positronium (o‐Ps) lifetime τ₃ and its intensity I₃ in various fluorinated polyimides were determined by the positron annihilation technique and were studied with the spin–lattice relaxation time T₁ and the propylene permeability, solubility, diffusivity, and permselectivity for propylene/propane in them. τ₃, I₃, and the distribution of τ₃ changed when the bulky moieties in the polyimides were changed. The polyimides, having both large τ₃ and I₃ values, exhibited a short T₁ and a high permeability with a low permselectivity. The propylene permeability and diffusivity were exponentially correlated with the product of I₃ and the average free‐volume hole size estimated from τ₃. In highly plasticized states induced by the sorption of propylene, the permeability increased with the propylene pressure in excellent agreement with the change in the free‐volume hole properties probed by o‐Ps. The large and broad distribution of the free‐volume holes and increased local chain mobility for the 2,2‐bis(3,4‐decarboxyphenyl) hexafluoropropane dianhydride‐based polyimides are thought to be important physical properties for promoting penetrant‐induced plasticization. These results suggest that o‐Ps is a powerful probe of not only the free‐volume holes but also the corresponding permeation mechanism and penetrant‐induced plasticization phenomenon. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 308–318, 2003     

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.     

AN ENCOMPASSING APPROACH TO MODELING FISHERY DYNAMICS: MODELING DYNAMIC NONLINEAR SYSTEMS

Uncertainty about the nature and significance of nonlinearities and the manner in which dynamics affect future realizations makes model specification the most difficult aspect of modeling dynamic systems. By interpreting several popular fishery models as subcases of a nesting dynamic Taylor series approximation, we isolate the specification differences between these models in a way that accounts for commonalities. On the argument that the differences due to alternative nonlinear forms are likely to be small compared to more mundane considerations such as delay difference and general dynamic lag specification, we propose an alternative model that uses the terms from the first order approximation common to all models combined with a data-based determination of the appropriate lags using the methods of state space time series analysis. Finally, the success of the alternative models is judged in an application to Pacific halibut data.