Recently there have two different effective methods proposed by Kanzow et al. in (Kanzow, 2001 [8]) and (Kanzow and Petra, 2004 [9]), respectively, which commonly use the FischerBurmeister (FB) function to recast the mixed complementarity problem (MCP) as a constrained minimization problem and a nonlinear system of equations, respectively. They all remark that their algorithms may be improved if the FB function is replaced by other NCP functions. Accordingly, in this paper, we employ the generalized FischerBurmeister (GFB) where the 2-norm in the FB function is relaxed to a general p-norm (p> 1) for the two methods and investigate how much the improvement is by changing the parameter p as well as which method is influenced more when we do so, by the performance profiles of iterations and function evaluations for the two methods with different p on MCPLIB collection.

In the paper [4], the authors proposed a derivative-free descent algorithm for the nonlinear complementarity problems (NCPs) by the generalized Fischer-Burmeister merit function: p (a, b)= 1

We show that the gradient mapping of the squared norm of FischerBurmeister function is globally Lipschitz continuous and semismooth, which provides a theoretical basis for solving nonlinear second-order cone complementarity problems via the conjugate gradient method and the semismooth Newtons method.

Shape optimization aims to optimize an objective function by changing the shape of the computational domain. In recent years, shape optimization has received considerable attentions. On the theoretical side there are several publications dealing with the existence of solution and the sensitivity analysis of the problem; see e.g., [6] and references therein.

Wind power is an increasingly popular renewable energy. In the design process of the wind turbine blade, the accurate aerodynamic simulation is important. In the past, most of the wind turbine simulations were carried out with some low fidelity methods, such as the blade element momentum method [9]. Recently, with the rapid development of the supercomputers, high fidelity simulations based on 3D unsteady Navier-Stokes (N-S) equations become more popular. For example, Sorensen et al. studied the 3D wind turbine rotor using the Reynolds-Averaged Navier-Stokes (RANS) framework where a finite volume method and a semi-implicit method are used for the spatial and temporal discretization, respectively [17]. Bazilevs et al. investigated the aerodynamic of the NREL 5MW offshore baseline wind turbine rotor using large eddy simulation built with a deforming-spatial-domain/stabilized space-time