On accurate boundary conditions for a shape sensitivity equation method

This paper studies the application of the continuous sensitivity equation method (CSEM) for the Navier–Stokes equations in the particular case of shape parameters. Boundary conditions for shape parameters involve flow derivatives at the boundary. Thus, accurate flow gradients are critical to the success of the CSEM. A new approach is presented to extract accurate flow derivatives at the boundary. High order Taylor series expansions are used on layered patches in conjunction with a constrained least‐squares procedure to evaluate accurate first and second derivatives of the flow variables at the boundary, required for Dirichlet and Neumann sensitivity boundary conditions. The flow and sensitivity fields are solved using an adaptive finite‐element method. The proposed methodology is first verified on a problem with a closed form solution obtained by the Method of Manufactured Solutions. The ability of the proposed method to provide accurate sensitivity fields for realistic problems is then demonstrated. The flow and sensitivity fields for a NACA 0012 airfoil are used for fast evaluation of the nearby flow over an airfoil of different thickness (NACA 0015). Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of plane and axisymmetric flows of incompressible fluids with the stream tube method: Numerical simulation by trust-region optimization algorithm

New concepts for the study of incompressible plane or axisymmetric flows are analysed by the stream tube method. Flows without eddies and pure vortex flows are considered in a transformed domain where the mapped streamlines are rectilinear or circular. The transformation between the physical domain and the computational domain is an unknown of the problem. In order to solve the non-linear set of relevant equations, we present a new algorithm based on a trust region technique which is effective for non-convex optimization problems. Experimental results show that the new algorithm is more robust compared to the Newton-Raphson method.     

A generalized Adadorov theory for anisotropic thin-walled beams

This paper presents a generalized Adadorov theory for anisotropic thin—walled beams. The theory takes account of the shear strain of the middle surface, which exerts a significant influence on the anisotropic thin-walled beams. A new approach is established to solve the governing equations, which have the same form for both open and closed section beams. The numerical examples show that the effects of the shear strain cannot be neglected for this class of beams.This work was part of research project supported by the National Natural Science Foundation of China     

A note on an SOR-like method for augmented systems

Golub et al. (2001, BIT, 41, 71–85) gave a generalizedsuccessive over-relaxation method for the augmented systems.In this paper, the connection between the SOR-like method andthe preconditioned conjugate gradient (PCG) method for the augmentedsystems is investigated. It is shown that the PCG method isat least as accurate (fast) as the SOR-like method. Numericalexamples demonstrate that the PCG method is much faster thanthe SOR-like method.     

A simple four-node quadrilateral finite element for plates

A simple and accurate four-node quadrilateral finite element based on the Mindlin plate theory and Kirchhoff constraints is presented for general thin plate bending applications. The derivation of the element stiffness properties is straightforward, starting with a specified eight-node interpolation; usual discrete Kirchhoff (DK) constraints are employed to constrain out the four midside nodes of the element. The present resulting DK element passes patch tests with elements of arbitrary and even highly distorted mesh types. Numerical studies of the element convergence behaviours are undertaken for various plate bending problems so far investigated. It is indicated from comparative examples that fairly good convergence characteristics have been achieved.     

广义拟变分不等式解集的稳定性及本质连通区的存在性

在赋范线性空间下．讨论广义拟变分不等式解集的稳定性，证明了满足一定连续性和凸性条件的广义拟变分不等式问题构成的空间M中，大多数(在Baire分类意义下)广义拟变分不等式问题的解集是稳定的，并证明了M中每一个广义拟变分不等式的解集至少存在一个本质连通区。     

A discontinuous Galerkin method/HLLC solver for the Euler equations

This paper proposes a fully three‐dimensional non‐linear Euler methodology for solving aerodynamic and acoustic problems in the presence of strong shocks and rarefactions. It uses a discontinuous Galerkin method (DGM) within the element, and a Riemann solver (HLLC) at the boundaries to propagate rarefactions while preserving the entropy condition and capturing shocks with no spurious oscillations. This approach is thought to marry the best aspects of finite element and finite volume methods, achieving conservation while not requiring the solution of a large matrix. Examples in which shock and rarefaction waves are well captured are presented and the propagation of acoustic pulses is well demonstrated. Copyright © 2003 John Wiley & Sons, Ltd.     

Formation of β‐iPP in isotactic polypropylene/ethylene–propylene rubber blends: Effects of preparation method,composition, and thermal condition

The effects of preparation method, composition, and thermal condition on formation of β‐iPP in isotactic polypropylene/ethylene–propylene rubber (iPP/EPR) blends were studied using modulated differential scanning calorimeter (MDSC), wide angle X‐ray diffraction (WAXD), and phase contrast microscopy (PCM). It was found that the α‐iPP and β‐iPP can simultaneity form in the melt‐blended samples, whereas only α‐iPP exists in the solution‐blended samples. The results show that the formation of β‐iPP in the melt‐blended samples is related to the crystallization temperature and the β‐iPP generally diminishes and finally vanishes when the crystallization temperature moves far from 125 °C. The phenomena that the lower critical temperature of β‐iPP in iPP/EPR obviously increases to 114 °C and the upper critical temperature decreases to 134 °C indicate the narrowing of temperature interval, facilitating the formation of β‐iPP in iPP/EPR. Furthermore, it was found that the amount of β‐iPP in melt‐blended iPP/EPR samples is dependent on the composition and the maximum amount of β‐iPP formed when the composition of iPP/EPR blends is 85:15 in weight. The results through examining the effect of annealing for iPP/EPR samples at melt state indicate that this annealing may eliminate the susceptibility to β‐crystallization of iPP. However, only α‐iPP can be observed in solution‐blended samples subjected to annealing for different time. The PCM images demonstrate that an obvious phase‐separation happens in both melt‐blended and solution‐blended iPP/EPR samples, implying that compared with the disperse degree of EPR in iPP, the preparation method plays a dominant role in formation of β‐iPP. It is suggested that the origin of formation of β‐iPP results from the thermomechanical history of the EPR component in iPP/EPR. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1704–1712, 2007