Detecting instabilities in flows of viscoelastic fluids

There is a growing interest in developing numerical tools to investigate the onset of physical instabilities observed in experiments involving viscoelastic flows, which is a difficult and challenging task as the simulations are very sensitive to numerical instabilities. Following a recent linear stability analysis carried out in order to better understand qualitatively the origin of numerical instabilities occurring in the simulation of flows viscoelastic fluids, the present paper considers a possible extension for more complex flows. This promising method could be applied to track instabilities in complex (i.e. essentially non‐parallel) flows. In addition, results related to transient growth mechanism indicate that it might be responsible for the development of numerical instabilities in the simulation of viscoelastic fluids. Copyright © 2003 John Wiley & Sons, Ltd.     

A parallel cell‐based DSMC method on unstructured adaptive meshes

A parallel DSMC method based on a cell‐based data structure is developed for the efficient simulation of rarefied gas flows on PC‐clusters. Parallel computation is made by decomposing the computational domain into several subdomains. Dynamic load balancing between processors is achieved based on the number of simulation particles and the number of cells allocated in each subdomain. Adjustment of cell size is also made through mesh adaptation for the improvement of solution accuracy and the efficient usage of meshes. Applications were made for a two‐dimensional supersonic leading‐edge flow, the axi‐symmetric Rothe's nozzle, and the open hollow cylinder flare flow for validation. It was found that the present method is an efficient tool for the simulation of rarefied gas flows on PC‐based parallel machines. Copyright © 2004 John Wiley & Sons, Ltd.     

A numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. 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.     

A virtual boundary method with improved computational efficiency using a multi‐grid method

The flow around spherical, solid objects is considered. The boundary conditions on the solid boundaries have been applied by replacing the boundary with a surface force distribution on the surface, such that the required boundary conditions are satisfied. The velocity on the boundary is determined by extrapolation from the flow field. The source terms are determined iteratively, as part of the solution. They are then averaged and are smoothed out to nearby computational grid points. A multi‐grid scheme has been used to enhance the computational efficiency of the solution of the force equations. The method has been evaluated for flow around both moving and stationary spherical objects at very low and intermediate Reynolds numbers. The results shows a second order accuracy of the method both at creeping flow and at Re=100. The multi‐grid scheme is shown to enhance the convergence rate up to a factor 10 as compared to single grid approach. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. Copyright © 2004 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

Synthesis and structural characterization of novel,fluorinated poly(phthalazinone ether)s containing perfluorophenylene moieties

The synthesis and structural characterization of a series of novel, fluorinated poly(phthalazinone ether)s containing perfluorophenylene moieties are described. The monomers, 4‐(4′‐hydroxyaryl)phthalazin‐1(2H)‐ones ( 2a – 2d ), were conveniently and efficiently synthesized from phenols and phthalic anhydride in two steps via 2‐(4′‐hydroxybenzoyl)benzoic acids, which were first obtained by the Friedel–Crafts reaction in good yields and with high stereoselectivity and were then converted into 2a – 2d by fusion with hydrazine. All the polymers were prepared by nucleophilic aromatic substitution (S_NAr) polycondensation between the compounds perfluorobiphenyl and 4‐(4′‐hydroxyaryl)phthalazin‐1(2H)‐ones ( 2a ‐ 2d ). The resulting fluorinated polymers were readily soluble in common organic solvents (e.g., CHCl₃, tetrahydrofuran, dimethylformamide, dimethyl sulfoxide, N‐methylpyrrolidone, etc.) at room temperature. Their weight‐average molecular weights and the polydispersities ranged from (7.96–18.25) × 10³ to 1.31–2.71, respectively. Their glass‐transition temperatures varied from 213 to 263 °C. They were all stable up to 390 °C both in air and in argon. The 5% weight‐loss temperatures of these polymers in air and argon ranged from 393–487 to 437–509 °C, respectively. Wide‐angle X‐ray diffraction studies indicated they were all amorphous and could be attributed to the presence of kink nonplanar moiety, phenyl phthalazinone along the polymer backbone. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 925–932, 2004     

Intramolecular energy transferable (2,2‐diphenylvinyl)phenyl substituted poly(p‐phenylenevinylene) derivative with efficient photoluminescence and electroluminescence

A poly(p‐phenylenevinylene) (PPV) derivative containing a bulky (2,2‐diphenylvinyl)phenyl group in the side chain, EHDVP‐PPV, was synthesized by Gilch route. The reduced tolane‐bisbenzyl (TBB) defects, as well as the structure of the polymer, was confirmed by various spectroscopic methods. The intramolecular energy transfer from the (2,2‐diphenylvinyl)phenyl side group to the PPV backbone was studied by UV‐vis and photoluminescence (PL) of the obtained polymer and model compound. The polymer film showed maximum absorption and emission peaks at 454 and 546 nm, respectively, and high PL efficiency of 57%. A yellow electroluminescence (λ_max = 548 nm) was obtained with intensities of 6479 cd/m² when the light‐emitting diodes of ITO/PEDOT/EHDVP‐PPV/LiF/Al were fabricated. The maximum power efficiency of the devices was 0.729 lm/W with a turn‐on voltage of 3.6 V. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5636–5646, 2004