Path-following interior point algorithms for the Cartesian <Emphasis Type="Italic">P</Emphasis><Subscript>*</Subscript>(κ)-LCP over symmetric cones

In this paper, we establish a theoretical framework of path-following interior point algorithms for the linear complementarity problems over symmetric cones (SCLCP) with the Cartesian P _*(κ)-property, a weaker condition than the monotonicity. Based on the Nesterov-Todd, xy and yx directions employed as commutative search directions for semidefinite programming, we extend the variants of the short-, semilong-, and long-step path-following algorithms for symmetric conic linear programming proposed by Schmieta and Alizadeh to the Cartesian P _*(κ)-SCLCP, and particularly show the global convergence and the iteration complexities of the proposed algorithms. This work was supported by National Natural Science Foundation of China (Grant Nos. 10671010, 70841008)     

Group and phase velocities in the free-surface visco-potential flow: New kind of boundary layer induced instability

Water wave propagation can be attenuated by various physical mechanisms. One of the main sources of wave energy dissipation lies in boundary layers. The present work is entirely devoted to thorough analysis of the dispersion relation of the novel visco-potential formulation. Namely, in this study we relax all assumptions of the weak dependence of the wave frequency on time. As a result, we have to deal with complex integro-differential equations that describe transient behaviour of the phase and group velocities. Using numerical computations, we show several snapshots of these important quantities at different times as functions of the wave number. Good qualitative agreement with previous study [D. Dutykh, Eur. J. Mech. B/Fluids 28 (2009) 430] is obtained. Thus, we validate in some sense approximations made anteriorly. There is an unexpected conclusion of this study. According to our computations, the bottom boundary layer creates disintegrating modes in the group velocity. In the same time, the imaginary part of the phase velocity remains negative for all times. This result can be interpreted as a new kind of instability which is induced by the bottom boundary layer effect.     

Stress distribution and hillock formation in Au/Pd thin films as a function of aging treatment in capacitor applications

Effect of quenching in different media on hillock formation and electrical resistivity has been studied in the Au–Pd layers. Oxygen was released from substrate due to substrate relaxation process. It was suggested that hillocks appear on the triple junction grain boundaries. However, lower electrical resistivity has been seen in the sample which quenched in the air. It was concluded that grain boundary scattering decreases the conductivity of the quenched films due to higher density of dislocations.     

A ghost fluid Lattice Boltzmann method for complex geometries

A ghost fluid Lattice Boltzmann method (GF‐LBM) is developed in this study to represent complex boundaries in Lattice Boltzmann simulations of fluid flows. Velocity and density values at the ghost points are extrapolated from the fluid interior and domain boundary via obtaining image points along the boundary normal inside the fluid domain. A general bilinear interpolation algorithm is used to obtain values at image points which are then extrapolated to ghost nodes thus satisfying hydrodynamic boundary conditions. The method ensures no‐penetration and no‐slip conditions at the boundaries. Equilibrium distribution functions at the ghost points are computed using the extrapolated values of the hydrodynamic variables, while non‐equilibrium distribution functions are extrapolated from the interior nodes. The method developed is general, and is capable of prescribing Dirichlet as well as Neumann boundary conditions for pressure and velocity. Consistency and second‐order accuracy of the method are established by running three test problems including cylindrical Couette flow, flow between eccentric rotating cylinders and flow over a cylinder in a confined channel. Copyright © 2011 John Wiley & Sons, Ltd.     

Interior hybrid proximal extragradient methods for the linear monotone complementarity problem

We present new infeasible path-following methods for linear monotone complementarity problems based on Auslender, Teboulle and Ben-Tiba’s log-quadratic barrier functions. The central paths associated with these barriers are always well defined and, for those problems which have a solution, convergent to a pair of complementary solutions. Starting points in these paths are easy to compute. The theoretical iteration-complexity of these new path-following methods is derived and improved by a strategy which uses relaxed hybrid proximal-extragradient steps to control the quadratic term. Encouraging preliminary numerical experiments are presented.     

Two solutions for a singular elliptic equation with critical growth at infinity

We look for positive solutions for the singular equation $-\Delta u-\frac{1}{2}\left(x\cdot \nabla u\right)=\mu h\left(x\right)u^{q-1}+\lambda u+u^{(N+2)/(N-2)},\text{in }{\mathbb{R}}^N,$ where $N\ge 3$, $\lambda >0$, $\mu >0$ is a parameter, $0<q<1$ and $h$ has some summability properties. By using a perturbation method and critical point theory, we obtain two solutions when $max\{1,N/4\}<\lambda <N/2$ and the parameter $\mu >0$ is small.     

Prediction of binary interaction coefficients in grain boundary segregation

A simple semiempirical method is proposed for prediction of the values of the binary interaction coefficients in grain boundary segregation (Fowler interaction). The method is based on the volume solid solubility, the maximum possible concentration of the solute at the grain boundary, and an additional parameter that represents an activity exponent of both the solute and the host metal. Values of the interaction coefficients for numerous solutes in α‐iron are predicted. The comparison with available literature data shows reasonable agreement.     

Free (open) boundary condition: some experiences with viscous flow simulations

The free (or open) boundary condition (FBC, OBC) was proposed by Papanastasiou et al. (A new outflow boundary condition, International Journal for Numerical Methods in Fluids, 1992; 14:587–608) to handle truncated domains with synthetic boundaries where the outflow conditions are unknown. In the present work, implementation of the FBC has been tested in several benchmark problems of viscous flow in fluid mechanics. The FEM is used to provide numerical results for both cases of planar and axisymmetric domains under laminar, isothermal or non‐isothermal, steady‐state conditions, for Newtonian fluids. The effects of inertia, gravity, compressibility, pressure dependence of the viscosity, slip at the wall, and surface tension are all considered individually in the extrudate‐swell benchmark problem for a wide range of the relevant parameters. The present results extend previous ones regarding the applicability of the FBC and show cases where the FBC is inappropriate, namely in the extrudate‐swell problem with gravity or surface‐tension effects. Particular emphasis has been given to the pressure at the outflow, which is the most sensitive quantity of the computations. In all cases where FBC is appropriate, excellent agreement has been found in comparisons with results from very long domains. The formulation for Picard‐type iterations is given in some detail, and the differences with the Newton–Raphson formulation are highlighted regarding some computational aspects. Copyright © 2011 John Wiley & Sons, Ltd.     

Simulation of viscoelastic fluids in a 2D abrupt contraction by spectral element method

This study presents the vortex structure and numerical instability increase occurring when the level of elasticity is enhanced in inertial flows in planar contraction configuration for finitely extensible nonlinear elastic model by Peterlin (FENE‐P) fluid 1 . The re‐entrant corner effect on corner vortices is also considered. The calculations are performed using extended matrix logarithm formulation described in a previous paper: A. Jafari et al. A new extended matrix logarithm formulation for the simulation of viscoelastic fluids by spectral elements. Computer & Fluids 2010; 39 (9):1425–1438. In that reference, the proposed algorithm has been tested for simple geometry such as Poiseuille flow. In this study, we are interested in the capability of this algorithm for more complex geometry. This formulation helps to reach higher values of the Weissenberg number when compared with the classical one. Copyright © 2015 John Wiley & Sons, Ltd.     

The parameter uniform numerical method for singularly perturbed parabolic reaction–diffusion problems on equidistributed grids

This article presents the study of singularly perturbed parabolic reaction–diffusion problems with boundary layers. To solve these problems, we use a modified backward Euler finite difference scheme on layer adapted nonuniform meshes at each time level. The nonuniform meshes are obtained by equidistribution of a positive monitor function, which involves the second-order spatial derivative of the singular component of the solution. The equidistributing monitor function at each time level allows us to use this technique to non-linear parabolic problems. The truncation error and the stability analysis are obtained. Parameter–uniform error estimates are derived for the numerical solution. To support the theoretical results, numerical experiments are carried out.