A short note on the entrainment and exit boundary conditions

An attempt is made to find out the suitable entrainment and exit boundary conditions in laminar flow situations. Streamfunction vorticity formulation of the Navier–Stokes equations are solved by ADI method. Two‐dimensional laminar plane wall jet flow is used to test different forms of the boundary conditions. Results are compared with the experimental and similarity solution and the proper boundary condition is suggested. The Kind 1 boundary condition is recommended. It consists of zero first derivative condition for velocity variable and for streamfunction equation, mixed derivative at the entrainment and exit boundaries. Copyright © 2005 John Wiley & Sons, Ltd.     

Nanostructured thermoset epoxy resin templated by an amphiphilic poly(ethylene oxide)‐block‐poly(dimethylsiloxane) diblock copolymer

An amphiphilic poly(ethylene oxide)‐block‐poly(dimethylsiloxane) (PEO–PDMS) diblock copolymer was used to template a bisphenol A type epoxy resin (ER); nanostructured thermoset blends of ER and PEO–PDMS were prepared with 4,4′‐methylenedianiline (MDA) as the curing agent. The phase behavior, crystallization, hydrogen‐bonding interactions, and nanoscale structures were investigated with differential scanning calorimetry, Fourier transform infrared spectroscopy, transmission electron microscopy, and small‐angle X‐ray scattering. The uncured ER was miscible with the poly(ethylene oxide) block of PEO–PDMS, and the uncured blends were not macroscopically phase‐separated. Macroscopic phase separation took place in the MDA‐cured ER/PEO–PDMS blends containing 60–80 wt % PEO–PDMS diblock copolymer. However, the composition‐dependent nanostructures were formed in the cured blends with 10–50 wt % PEO–PDMS, which did not show macroscopic phase separation. The poly(dimethylsiloxane) microdomains with sizes of 10–20 nm were dispersed in a continuous ER‐rich phase; the average distance between the neighboring microdomains was in the range of 20–50 nm. The miscibility between the cured ER and the poly(ethylene oxide) block of PEO–PDMS was ascribed to the favorable hydrogen‐bonding interaction. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3042–3052, 2006     

Method of lines with boundary elements for 1‐D transient diffusion‐reaction problems

Time‐dependent differential equations can be solved using the concept of method of lines (MOL) together with the boundary element (BE) representation for the spatial linear part of the equation. The BE method alleviates the need for spatial discretization and casts the problem in an integral format. Hence errors associated with the numerical approximation of the spatial derivatives are totally eliminated. An element level local cubic approximation is used for the variable at each time step to facilitate the time marching and the nonlinear terms are represented in a semi‐implicit manner by a local linearization at each time step. The accuracy of the method has been illustrated on a number of test problems of engineering significance. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

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.     

Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

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.     

A WKB analysis of the buckling condition for a cylindrical shell of arbitrary thickness subjected to an external pressure

We consider the plane-strain buckling of a cylindrical shellof arbitrary thickness which is made of a Varga material andis subjected to an external hydrostatic pressure on its outersurface. The WKB method is used to solve the eigenvalue problemthat results from the linear bifurcation analysis. We show thatthe circular cross-section buckles into a non-circular shapeat a value of µ₁ which depends on A₁/A₂ and a mode number,where A₁ and A₂ are the undeformed inner and outer radii, andµ₁ is the ratio of the deformed inner radius to A₁. Inthe large mode number limit, we find that the dependence ofµ₁ on A₁/A₂ has a boundary layer structure: it is constantover almost the entire region of 0 < A₁/A₂ < 1 and decreasessharply from this constant value to unity as A₁/A₂ tends tounity. Our asymptotic results for A₁ – 1 = O(1) and A₁– 1 = O(1/n) are shown to agree with the numerical resultsobtained by using the compound matrix method.     

对求解非线性方程方法的探索

BaiCL在2001年提出了双曲函数法，并把解设成具体的双曲函数，即sinhw和coshw的线性组合．本文把解设成是由2个函数f(ξ)和g(ξ)组成的线性组合，它们可以有多种取法，从而使方程更直接、更有效．借助数学软件Maple，用改进后的方法和吴文俊消元法，求解Bai C L文中的一个例子，获得了包含Bai文结果在内的更为丰富、精确的行波解．     

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.     

Polymerization behaviors of racemic and chiral amphiphilic monomers in organized bilayer membranes of lamellar liquid crystalline phase

A racemic amphiphilic monomer, n‐dodecyl glyceryl itaconate (DGI), forms bilayer membranes in water in the presence of small amount of ionic cosurfactant and shows iridescent color. A chiral DGI, S‐DGI, also shows an iridescent property, but with a rather red shift in the color, which can be ascribed to the increased packing density of the monomer in the bilayer membranes. Chrial DGI has a more compact packing density than racemic one owing to closer distance between the monomer molecules; the conversion rate, however, is slower than that of racemic one when H₂O₂ is used as an initiator. When the initiator is changed to an amphiphilic one, 4‐(2‐hydroxyethoxy) phenyl‐(2‐hydroxy‐2‐propyl) ketone (Irgacure 2959), the chiral DGI shows even a little faster conversion rate than that of racemic one. The NMR chemical shift results of protons in benzene ring show that the molecules of Irgacure 2959 insert into the bilayer membranes. The molecular weights of the corresponding polymers prove that the initiation by H₂O₂ is restricted compared to that by Irgacure 2959. It is concluded that the decelerated polymerization behavior of chiral DGI initiated by H₂O₂ is a result of limited diffusion of the initiator into the lamellar bilayer structures. © 2007 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 4891–4900, 2007