The algebraic Bethe ansatz for the Izergin–Korepin model with open boundary conditions

We present the procedure of exactly solving the Izergin–Korepin model with open boundary conditions by using the algebraic Bethe ansatz, which include constructing the multi-particle state and achieving the eigenvalue of the transfer matrix and corresponding Bethe equations. We give a proof about our conclusions on the multi-particle state based on an assumption. When the model is U_q(su(2)) quantum invariant, our results agree with that obtained by analytic Bethe ansatz method.     

广义C-R组的性质及两个边值问题

本文考虑广义C-R组^[1](H)的解f=u iv jw∈C^2的一些性质，提出与之相关的两个边值问题，用积分方程方法和调和函数性质证明了边值问题解的存在唯一性，并写出解的积分表达式．     

Methods for implementing the stream boundary condition in DSMC computations

In the direct simulation Monte‐Carlo (DSMC) method for simulating rarefied gas flows, the velocities of simulator particles that cross a simulation boundary and enter the simulation space are typically generated using the acceptance–rejection procedure that samples the velocities from a truncated theoretical velocity distribution that excludes low and high velocities. This paper analyses an alternative technique, where the velocities of entering particles are obtained by extending the simulation procedures to a region adjacent to the simulation space, and considering the movement of particles generated within that region during the simulation time step. The alternative method may be considered as a form of acceptance–rejection procedure, and permits the generation of all possible velocities, although the population of high velocities is depleted with respect to the theoretical distribution. Nevertheless, this is an improvement over the standard acceptance–rejection method. Previous implementations of the alternative method gave a number flux lower than the theoretical number required. Two methods for obtaining the correct number flux are presented. For upstream boundaries in high‐speed flows, the alternative method is more computationally efficient than the acceptance–rejection method. However, for downstream boundaries, the alternative method is extremely inefficient. The alternative method, with the correct theoretical number flux, should therefore be used in DSMC computations in favour of the acceptance–rejection method for upstream boundaries in high‐speed flows. Copyright © 2003 John Wiley & Sons, Ltd.     

Critical Point Theorems and Applications to Differential Equations

This paper contains a generalization of the well-known Palais-Smale and Cerami compactness conditions. The compactness condition introduced is used to prove some general existence theorems for critical points. Some applications are given to differential equations.     

On the application of slip boundary condition on curved boundaries

Hydrodynamic simulations of sloshing phenomena often involve the application of slip boundary condition at the wetted surfaces. If these surfaces are curved, the ambiguous nature of the normal vector in the discretized problem can interfere with the application of such a boundary condition. Even the use of consistent normal vectors, preferred from the point of view of conservation, does not assure good approximation of the continuum slip condition in the discrete problem, and non‐physical recirculating flow fields may be observed. As a remedy, we consider the Navier slip condition, and more successfully, the so‐called BC‐free boundary condition. 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.     

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.     

Free boundary viscous flows at micro and mesolevel during liquid composites moulding process

Numerical simulation aspects, related to low Reynolds number free boundary viscous flows at micro and mesolevel during the resin impregnation stage of the liquid composite moulding process (LCM), are presented in this article. A free boundary program (FBP), developed by the authors, is used to track the movement of the resin front accurately by accounting for the surface tension effects at the boundary. Issues related to the global and local mass conservation (GMC and LMC) are identified and discussed. Unsuitable conditions for LMC and consequently GMC are uncovered at low capillary numbers, and hence a strategy for the numerical simulation of such flows is suggested. FBP encompasses a set of subroutines that are linked to modules in ANSYS. FBP can capture the void formation dynamics based on the analysis developed. We present resin impregnation dynamics in two dimensions. Extension to three dimensions is a subject for further research. Several examples are shown and efficiency of different stabilization techniques are compared. 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.     

Desorption behavior of a surfactant in a linear low‐density polyethylene blend at elevated temperatures

The desorption behavior of a surfactant in a linear low‐density polyethylene (LLDPE) blend at elevated temperatures of 50, 70, and 80 °C was studied with Fourier transform infrared spectroscopy. The composition of the LLDPE blend was 70:30 LLDPE/low‐density polyethylene. Three different specimens (II, III, and IV) were prepared with various compositions of a small molecular penetrant, sorbitan palmitate (SPAN‐40), and a migration controller, poly(ethylene acrylic acid) (EAA), in the LLDPE blend. The calculated diffusion coefficient (D) of SPAN‐40 in specimens II, III, and IV, between 50 and 80 °C, varied from 1.74 × 10^?11 to 6.79 × 10^?11 cm²/s, from 1.10 × 10^?11 to 5.75 × 10^?11 cm²/s, and from 0.58 × 10^?11 to 4.75 × 10^?11 cm²/s, respectively. In addition, the calculated activation energies (E_D) of specimens II, III, and IV, from the plotting of ln D versus 1/T between 50 and 80 °C, were 42.9, 52.7, and 65.6 kJ/mol, respectively. These values were different from those obtained between 25 and 50 °C and were believed to have been influenced by the interference of Tinuvin (a UV stabilizer) at elevated temperatures higher than 50 °C. Although the desorption rate of SPAN‐40 increased with the temperature and decreased with the EAA content, the observed spectral behavior did not depend on the temperature and time. For all specimens stored over 50 °C, the peak at 1739 cm^?1 decreased in a few days and subsequently increased with a peak shift toward 1730 cm^?1. This arose from the carbonyl stretching vibration of Tinuvin, possibly because of oxidation or degradation at elevated temperatures. In addition, the incorporation of EAA into the LLDPE blend suppressed the desorption rate of SPAN‐40 and retarded the appearance of the 1730 cm^?1 peak. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1114–1126, 2004