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.     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. 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.     

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.     

UV‐ignited frontal polymerization of an epoxy resin

By combining frontal polymerization and radical‐induced cationic polymerization, it was possible to cure thick samples of an epoxy monomer bleached by UV light. The effect of the relative amounts of cationic photoinitiator and radical initiator was thoroughly investigated and was related to the front's velocity and its maximum temperature. The materials obtained were characterized by quantitative conversion also in the deeper layers, not reached by UV light. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 2066–2072, 2004     

Photoinitiators with functional groups. VII. Covalently bonded camphorquinone—amines

Camphorquinone (CQ), a widely used photoinitiator (PI) in dental applications, was covalently bonded to aromatic amines to enhance the rate of electron and proton transfer effect due to the close vicinity of the diketone and the amine group. 10‐bromocamphorquinone and 10‐bromomethylcamphorquinone were selected as suitable precursors for esterification with the carboxyl group containing aromatic amines based on 4‐dimethylaminobenzoic acid. Properties of the new photoinitiating systems were investigated by UV spectroscopy and differential scanning photocalorimetry in lauryl acrylate. Compared to physical mixtures, in all cases similar or even better performance was obtained. Surprisingly, 10‐acetyl derivatives 7 – 9 and 18 especially, were found to be highly reactive. Compared to CQ/ethyl 4‐dimethylaminobenzoate, the rate of photopolymerization was increased by a factor of up to 2. Intramolecular reaction was confirmed by photo‐differential scanning calorimetry experiments with varying PI concentrations. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4948–4963, 2004     

Synthesis and characterization of new segmented polyurethanes with side‐chain,liquid‐crystalline chain extenders

Liquid‐crystalline, segmented polyurethanes with methoxy–biphenyl mesogens pendant on the chain extender were synthesized by the conventional prepolymer technique and esterification reaction. Two, side‐chain, liquid‐crystalline (SCLC) polyurethanes with mesogens having spacers of six and eight methylene units were prepared. The structures of the mesogenic units and SCLC polyurethanes were confirmed by Fourier transform infrared spectroscopy and ¹H NMR. Polymer properties were also examined by solubility tests, water uptakes, and inherent viscosity measurements. Differential scanning calorimetry studies indicated that the transition temperature of the isotropic to the liquid‐crystalline phase decreased with increasing spacer length. Wide‐angle X‐ray diffraction (WAXD) studies revealed the existence of liquid‐crystalline phases for both SCLC polyurethanes. Polarized optical microscopic investigations further confirmed the thermotropic liquid‐crystalline behaviors and nematic mesophases of both samples. Thermogravimetric analysis displayed better thermal stabilities for both SCLC polymers and indicated that the presence of mesogenic side chains may increase the thermal stability of segmented polyurethanes. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 290–302, 2004