MathSciDoc: An Archive for Mathematician ∫

We propose and study a new parallel one-shot Lagrange--Newton--Krylov--Schwarz (LNKSz) algorithm for shape optimization problems constrained by steady incompressible Navier--Stokes equations discretized by finite element methods on unstructured moving meshes. Most existing algorithms for shape optimization problems iteratively solve the three components of the optimality system: the state equations for the constraints, the adjoint equations for the Lagrange multipliers, and the design equations for the shape parameters. Such approaches are relatively easy to implement, but generally are not easy to converge as they are basically nonlinear Gauss--Seidel algorithms with three large blocks. In this paper, we introduce a fully coupled, or the so-called one-shot, approach which solves the three components simultaneously. First, we introduce a moving mesh finite element method for the shape optimization

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

We simulate blood flow in patient-specific cerebral arteries. The complicated geometry in the human brain makes the problem challenging. We use a fully unstructured three dimensional mesh to cover the artery, and Galerkin/least-squares finite element method to discretize the incompressible Navier-Stokes equations, that are employed to model the blood flow, and the resulting large sparse nonlinear system of equations is solved by a Newton-Krylov-Schwarz algorithm. From the computed flow fields, we are able to understand certain behavior of the blood flow of this particular patient before and after a stenosis is surgically removed. We also report the robustness and parallel performance of the domain decomposition based algorithm.

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.