Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

Absorbing boundary conditions for reaction-diffusion equations

We treat here of the question of absorbing boundary conditionsfor nonlinear diffusion equations. We use the conditions designedfor the linear equation, we prove them to be well posed forthe nonlinear problem, and through numerical experiments thatthey are well suited for reaction–diffusion equations.     

Synthesis and refractive‐index control of fluoropolymers by the polyaddition of bis(epoxide)s and triazine di(aryl ether)s

The polyaddition of fluorine‐containing bis(epoxide)s and fluorine‐containing triazine di(aryl ether)s were examined to give the corresponding fluorine‐containing poly(cyanurate)s. It was observed that the synthesized fluoropolymers had good thermal stabilities and good film‐forming properties. The glass transition temperatures (T_g's) and refractive‐indices (n_D's) of synthesized polymers were determined by differential scanning calorimetry and ellipsometry, respectively, and it was found that the values of T_g's and n_D's were supported by their fluorine containing ratios and skeletons. © 2007 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 4421–4429, 2007     

A test problem for molecular dynamics integrators

We derive a test problem for evaluating the ability of time-steppingmethods to preserve the statistical properties of systems inmolecular dynamics. We consider a family of deterministic systemsconsisting of a finite number of particles interacting on acompact interval. The particles are given random initial conditionsand interact through instantaneous energy- and momentum-conservingcollisions. As the number of particles, the particle density,and the mean particle speed go to infinity, the trajectory ofa tracer particle is shown to converge to a stationary Gaussianstochastic process. We approximate this system by one describedby a system of ordinary differential equations and provide numericalevidence that it converges to the same stochastic process. Wesimulate the latter system with a variety of numerical integrators,including the symplectic Euler method, a fourth-order Runge-Kuttamethod, and an energyconserving step-and-project method. Weassess the methods' ability to recapture the system's limitingstatistics and observe that symplectic Euler performs significantlybetter than the others for comparable computational expense.     

Reference values for drag and lift of a two‐dimensional time‐dependent flow around a cylinder

This paper presents a numerical study of a two‐dimensional time‐dependent flow around a cylinder. Its main objective is to provide accurate reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time. In addition, the accuracy of these values obtained with different time stepping schemes and different finite element methods is studied. 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     

Thermoplastic poly(ester‐imide)s derived from anhydride‐terminated polyester prepolymer and diisocyanate

The phase‐separation behavior of thermoplastic poly(ester‐imide) [P(E‐I)] multiblock copolymers, (A‐B)_n, was investigated by a stepwise variation of the imide content. All the multiblock copolymers were synthesized by solution polycondensation with dimethylformamide as a solvent. P(E‐I)s were prepared with anhydride‐terminated polyester prepolymer and diisocyanates. Polyester prepolymers were prepared by the reaction of pyromellitic dianhydride and two different polyols [poly(tetramethylene oxide glycol) (PTMG) and polycaprolactone diol (PCL)]. Structural determination was done with Fourier transform infrared spectroscopy and Fourier transform NMR, and the molecular weight was determined by gel permeation chromatography. The effect of the imide content on the thermal properties of the synthesized P(E‐I)s was investigated by thermogravimetric analysis and differential scanning calorimetry. The polymers were also characterized for static and dynamic mechanical properties. Thermal analysis data indicated that the polymers based on PTMG were stable up to 330 °C in nitrogen atmosphere and exhibited phase‐separated morphology. Polymers based on PCL showed multistage decomposition, and the films derived from them were too fragile to be characterized for static and dynamic mechanical properties. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 341–350, 2004     

New liquid‐crystalline poly(ester amide)s: The role of nitro groups in the phase behavior

New hydrogen‐bonded liquid‐crystalline poly(ester amide)s (PEA)s were obtained from 1,4‐terephthaloyl[bis‐(3‐nitro‐N‐anthranilic acid)] (5) or 1,4‐terephthaloyl[bis‐(N‐anthranilic acid)] (6), with or without nitro groups, respectively, through the separate condensation of each with hydroquinone or dihydroxynaphthalene. The dicarboxylic monomers were synthesized from 2‐aminobenzoic acid. The phase behavior of the monomers and polymers were studied with differential scanning calorimetry, polarized light microscopy, and wide‐angle X‐ray diffraction methods. Monomer 5, containing nitro groups, exhibited a smectic liquid‐crystalline phase, whereas the texture of monomer 6 without nitro groups appeared to be nematic. The PEAs containing nitro groups exhibited polymorphism (smectic and nematic), whereas those without nitro groups exhibited only one phase transition (a nematic threaded texture). The changes occurring in the phase behavior of the polymers were explained by the introduction of nitro groups. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1289–1298, 2004