An improved numerical simulator for multiphase flow in porous media

In basin modelling the thermodynamics of a multicomponent multiphase fluid flux are computationally too expensive when derived from an equation of state and the Gibbs equality constraints. In this article we present a novel implicit molar mass formulation technique using binary mixture thermodynamics. The two proposed solution methods, with and without cross derivative terms between components, are based on a preconditioned Newton‐GMRES scheme for each time‐step with analytical computation of the derivatives. These new algorithms reduce significantly the numerical effort for the computation of the molar masses, and we illustrate the behavior of these methods with numerical computations. Copyright © 2004 John Wiley & Sons Ltd.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

Assessment of two‐equation turbulence modelling for high Reynolds number hydrofoil flows

This paper presents an evaluation of the capability of turbulence models available in the commercial CFD code FLUENT 6.0 for their application to hydrofoil turbulent boundary layer separation flow at high Reynolds numbers. Four widely applied two‐equation RANS turbulence models were assessed through comparison with experimental data at Reynolds numbers of 8.284×106 and 1.657×107. They were the standard k–εmodel, the realizable k–εmodel, the standard k–ωmodel and the shear‐stress‐transport (SST) k–ωmodel. It has found that the realizable k–εturbulence model used with enhanced wall functions and near‐wall modelling techniques, consistently provides superior performance in predicting the flow characteristics around the hydrofoil. Copyright © 2004 John Wiley & Sons, Ltd.     

Robust handling of non‐linear constraints for GA optimization of aerodynamic shapes

A new approach to the robust handling of non‐linear constraints for GAs (genetic algorithms) optimization is proposed. A specific feature of the approach consists of the change in the conventional search strategy by employing search paths which pass through both feasible and infeasible points (contrary to the traditional approach where only feasible points may be included in a path). The method (driven by full Navier–Stokes computations) was applied to the problem of multiobjective optimization of aerodynamic shapes subject to various geometrical and aerodynamic constraints. The results demonstrated that the method retains high robustness of conventional GAs while keeping CFD computational volume to an acceptable level, which allowed the algorithm to be used in a demanding engineering environment. Copyright © 2004 John Wiley & Sons, Ltd.     

Synthesis and characterization of oxidized cellulose

Samples of oxidized cellulose (OC) with various carboxyl contents and degrees of crystallinity were obtained by the oxidation of native and mercerized cellulose with a solution of nitrogen(IV) oxide in CCl₄. A detailed characterization of these OC samples was performed. The effect of oxidation conditions (concentration of N₂O₄ in the solution and oxidation time) and starting cellulose material on OC characteristics (carboxyl, carbonyl and nitrogen content, degree of crystallinity and polymerization, surface area and swelling, and acidic properties) was investigated. Reactivity in the oxidation process was higher in mercerized cellulose than in native cellulose. The action of dilute solutions (10–15%) of N₂O₄ did not affect the degree of crystallinity of cellulose samples. Under these conditions, the oxidation took place mainly in amorphous regions and on the surface of crystallites. Oxidation in a concentrated (40%) N₂O₄ solution led to the destruction of crystallites, which increased the surface area and swelling of cellulose in water. The surface area and the swelling of OC samples increased with a decrease in the index of crystallinity. The acidic properties of OC were shown to increase with an increase of swelling in water. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4785–4791, 2004     

Photochemical attachment of polymers on planar surfaces with a covalently anchored monolayer of a novel naphthyl ketone photochemical radical generator

A monolayer of covalently anchored, novel, binaphthyl ketone is used as a surface‐confined photochemical radical generator (PRG) for anchoring a variety of polymers to silicon surfaces. The precursor PRG is synthesized by the application of a facile and novel method for the oxidation of sterically hindered benzylic hydrocarbons to carbonyl compounds. Oxidation was carried out with a stoichiometric amount of potassium peroxydisulfate, in the presence of a catalytic amount of copper sulfate in an acetonitrile/water mixture. The PRG synthesized is characterized by ¹H NMR, UV, and Fourier transform infrared (FTIR). The covalently attached monolayers are characterized by X‐ray photoelectron spectroscopy, ellipsometry, and water contact angle measurements. The method developed is applicable to the preparation of a monolayer of a variety of polymers on a wide range of substrates carrying surface hydroxyl groups. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5413–5423, 2004     

Antimicrobial activity on the copolymers of 2,4‐dichlorophenyl methacrylate with methyl methacrylate: Synthesis and characterization

Copolymers of monomers 2,4‐dichlorophenyl methacrylate (2,4‐DMA) and methyl methacrylate (MMA) were synthesized with different monomer feed ratios using toluene as a solvent and 2,2′‐azobisisobutyronitrile (AIBN) as an initiator at 70 °C. The copolymers were characterized by IR‐spectroscopy, and copolymer composition was determined with UV‐spectroscopy. The linearization method of Fineman–Ross was employed to obtain the monomer reactivity ratios. The molecular weights and polydispersity indexes were determined by gel permeation chromatography (GPC). Thermogravimetric analyses of polymers were carried out in nitrogen atmosphere. The homo‐ and copolymers were tested for their antimicrobial properties against selected microorganisms. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5227–5234, 2004     

Dark cure studies of cationic photopolymerizations of epoxides: Characterization of kinetic rate constants at high conversions

The effective propagation rate constant (k_p; averaged over all the propagating active centers) was characterized for solvent‐free cationic photopolymerizations of phenyl glycidyl ether over the entire range of conversions, including the high conversion regime in which mass transfer limitations become important. The profile for the k_p as a function of conversion was found to exhibit a constant plateau value at low to intermediate conversions, followed by a monotonic increase above a threshold value of conversion. To explain this trend, it is proposed that at high conversion the diffusional mobility of the photoinitiator counterion is reduced whereas the mobility of the cationic active center remains high because of reactive diffusion. Therefore, with increasing conversion, the average distance between the active centers and counterions may increase, resulting in an increase in the propagation rate constant. The profiles for the k_p values were investigated as a function of the temperature, photoinitiator anion, and photoinitiator concentration. As the photoinitiator concentration was increased, the plateau value of the effective propagation rate constant decreased whereas the threshold conversion increased. All of the experimental trends are consistent with the proposed increase in ion separation at high conversions. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4409–4416, 2004     

Photopolymerization of 1,6‐hexanedioldiacrylate initiated by three‐component systems based on N‐arylphthalimides

Three‐component photoinitiators comprised of an N‐arylphthalimide, a diarylketone, and a tertiary amine were investigated for their initiation efficiency of acrylate polymerization. The use of an electron‐deficient N‐arylphthalimide resulted in a greater acrylate polymerization rate than an electron‐rich N‐arylphthalimide. Triplet energies of each N‐arylphthalimide, determined from their phosphorescence spectra, and the respective rate constants for triplet quenching by the N‐arylphthalimide derivatives (acquired via laser flash photolysis) indicated that an electron–proton transfer from an intermediate radical species to the N‐arylphthalimide (not energy transfer from triplet sensitization) is responsible for generating the initiating radicals under the conditions and species concentrations used for polymerization. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4009–4015, 2004