High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Nano-tribological study on a super-hydrophobic film formed on rough aluminium substrates

A novel super-hydrophobic stearic acid (STA) film with a water contact angle of 166^o was prepared by chemical adsorption on aluminum wafer coated with polyethyleneimine (PEI) film. The micro-tribological behavior of the super-hydrophobic STA monolayer was compared with that of the polished and PEI-coated Al surfaces. The effect of relative humidity on the adhesion and friction was investigated as well. It was found that the STA monolayer showed decreased friction, while the adhesive force was greatly decreased by increasing the surface roughness of the Al wafer to reduce the contact area between the atomic force microscope (AFM) tip and the sample surface to be tested. Thus the friction and adhesion of the Al wafer was effectively decreased by generating the STA monolayer, which indicated that it could be feasible and rational to prepare a surface with good adhesion resistance and lubricity by properly controlling the surface morphology and the chemical composition. Both the adhesion and friction decreased as the relative humidity was lowered from 65% to 10%, though the decrease extent became insignificant for the STA monolayer. The project supported by the National Natural Science Foundation of China (50375151, 50323007, 10225209) and the Chinese Academy of Sciences (KJCX-SW-L2)     

Fracture toughening of epoxy matrices with blends of resins of different molecular weights and other modifiers

The microstructure and fracture behavior of epoxy mixtures containing two monomers of different molecular weights were studied. The variation of the fracture toughness by the addition of other modifiers was also investigated. Several amounts of high‐molecular‐weight diglycidyl ether of bisphenol A (DGEBA) oligomer were added to a nearly pure DGEBA monomer. The mixtures were cured with an aromatic amine, showing phase separation after curing. The curing behavior of the epoxy mixtures was investigated with thermal measurements. A significant enhancement of the fracture toughness was accompanied by slight increases in both the rigidity and strength of the mixtures that corresponded to the content of the high‐molecular‐weight epoxy resin. Dynamic mechanical and atomic force microscopy measurements indicated that the generated two‐phase morphology was a function of the content of the epoxy resin added. The influence of the addition of an oligomer or a thermoplastic on the morphologies and mechanical properties of both epoxy‐containing mixtures was also investigated. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 3920–3933, 2004     

Control of self‐organized low‐dimensional morphology in poly(styrene‐b‐4vinylpyridine)/polystyrene blend thin films

The surface morphologies of poly(styrene‐b‐4vinylpyridine) (PS‐b‐P4VP) diblock copolymer and homopolystyrene (hPS) binary blend thin films were investigated by atomic force microscopy as a function of total volume fraction of PS (?_PS) in the mixture. It was found that when hPS was added into symmetric PS‐b‐P4VP diblock copolymers, the surface morphology of this diblock copolymer was changed to a certain degree. With ?_PS increasing at first, hPS was solubilized into the corresponding domains of block copolymer and formed cylinders. Moreover, the more solubilized the hPS, the more cylinders exist. However, when the limit was reached, excessive hPS tended to separate from the domains independently instead of solubilizing into the corresponding domains any longer, that is, a macrophase separation occurred. A model describing transitions of these morphologies with an increase in ?_PS is proposed. The effect of composition on the phase morphology of blend films when graphite is used as a substrate is also investigated. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 3496–3504, 2004     

On General Mixed Quasivariational Inequalities

In this paper, we suggest and analyze several iterative methods for solving general mixed quasivariational inequalities by using the technique of updating the solution and the auxiliary principle. It is shown that the convergence of these methods requires either the pseudomonotonicity or the partially relaxed strong monotonicity of the operator. Proofs of convergence is very simple. Our new methods differ from the existing methods for solving various classes of variational inequalities and related optimization problems. Various special cases are also discussed.     

Augmented Lagrangian Algorithms Based on the Spectral Projected Gradient Method for Solving Nonlinear Programming Problems

The spectral projected gradient method SPG is an algorithm for large-scale bound-constrained optimization introduced recently by Birgin, Martínez, and Raydan. It is based on the Raydan unconstrained generalization of the Barzilai-Borwein method for quadratics. The SPG algorithm turned out to be surprisingly effective for solving many large-scale minimization problems with box constraints. Therefore, it is natural to test its perfomance for solving the sub-problems that appear in nonlinear programming methods based on augmented Lagrangians. In this work, augmented Lagrangian methods which use SPG as the underlying convex-constraint solver are introduced (ALSPG) and the methods are tested in two sets of problems. First, a meaningful subset of large-scale nonlinearly constrained problems of the CUTE collection is solved and compared with the perfomance of LANCELOT. Second, a family of location problems in the minimax formulation is solved against the package FFSQP.     

Diffusion in poro‐plastic media

A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial‐boundary‐value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto‐visco‐plastic type. The variational form of this problem in Hilbert space is a non‐linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi‐static momentum equation. The essential sufficient conditions for the well‐posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd.     

Stability Conditions of the MMAP [ K ]/ G [ K ]/1/ LCFS Preemptive Repeat Queue

In this paper, we study the stability conditions of the MMAP[K]/G[K]/1/LCFS preemptive repeat queue. We introduce an embedded Markov chain of matrix M/G/1 type with a tree structure and identify conditions for the Markov chain to be ergodic. First, we present three conventional methods for the stability problem of the queueing system of interest. These methods are either computationally demanding or do not provide accurate information for system stability. Then we introduce a novel approach that develops two linear programs whose solutions provide sufficient conditions for stability or instability of the queueing system. The new approach is numerically efficient. The advantages and disadvantages of the methods introduced in this paper are analyzed both theoretically and numerically.     

晶体微观结构的数值计算

晶体微观结构是晶体材料在特定物理条件下其多个能量极小平衔态在空间形成的某种微尺度的规则分布．几何非线性的连续介质力学理论可以用能量极小化原理来解释晶体微观结构的形成，并用Young测度来刻画平衡态各变体在空间的概率分布．定性的理解与定量地分析和计算晶体材料的微观结构对于发展和改进高级晶体功能材料，如形状记忆合金、铁电体、磁至伸缩材料等，有重要的意义．本文回顾了近年来晶体微观结构数值计算方面的最新进展．介绍了计算晶体微观结构的几种数值方法及有关的数值分析结果。     

美式期权定价中非局部问题的有限元方法

在本文中 ,我们关心的是美式期权的有限元方法 .首先 ,根据 [4 ]我们对所讨论的问题引进一个新奇的实用的方法 ,它涉及到对原问题重新形成准确的数学公式 ,使得数值解的计算可以在非常小的区域上进行 ,从而该算法计算速度快精度高 .进而 ,我们利用超逼近分析技术得到了有限元解关于 L2 -模的最优估计 .