Shape optimization of a body located in low Reynolds number flow

The purpose of this study is to perform a numerical application of the shape optimization formulation of a body located in an incompressible viscous flow field. The formulation is based on an optimal control theory in which a performance function of the fluid force is introduced. The performance function should be minimized satisfying the state equation. This problem can be transformed into the minimization problem without constraint condition by the Lagrange multiplier method and the adjoint equations using adjoint variables corresponding to the state equations. As a numerical study, the drag force minimization problem in the steady Stokes flow, which means approximated equation of the low Reynolds number Navier–Stokes equation is carried out. After that, the unsteady Navier–Stokes flow is analysed. As the minimization algorithm, the steepest descent method is successfully applied. Copyright © 2005 John Wiley & Sons, Ltd.     

Continuous and discrete adjoint approaches for aerodynamic shape optimization with low Mach number preconditioning

Discrete and continuous adjoint approaches for use in aerodynamic shape optimization problems at all flow speeds are developed and assessed. They are based on the Navier–Stokes equations with low Mach number preconditioning. By alleviating the large disparity between acoustic waves and fluid speeds, the preconditioned flow and adjoint equations are numerically solved with affordable CPU cost, even at the so‐called incompressible flow conditions. Either by employing the adjoint to the preconditioned flow equations or by preconditioning the adjoint to the ‘standard’ flow equations (under certain conditions the two formulations become equivalent, as proved in this paper), efficient optimization methods with reasonable cost per optimization cycle, even at very low Mach numbers, are derived. During the mathematical development, a couple of assumptions are made which are proved to be harmless to the accuracy in the computed gradients and the effectiveness of the optimization method. The proposed approaches are validated in inviscid and viscous flows in external aerodynamics and turbomachinery flows at various Mach numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

Optimal flow control for Navier–Stokes equations: drag minimization

Optimal control and shape optimization techniques have an increasing role in Fluid Dynamics problems governed by partial differential equations (PDEs). In this paper, we consider the problem of drag minimization for a body in relative motion in a fluid by controlling the velocity through the body boundary. With this aim, we handle with an optimal control approach applied to the steady incompressible Navier–Stokes equations. We use the Lagrangian functional approach and we consider the Lagrangian multiplier method for the treatment of the Dirichlet boundary conditions, which include the control function itself. Moreover, we express the drag coefficient, which is the functional to be minimized, through the variational form of the Navier–Stokes equations. In this way, we can derive, in a straightforward manner, the adjoint and sensitivity equations associated with the optimal control problem, even in the presence of Dirichlet control functions. The problem is solved numerically by an iterative optimization procedure applied to state and adjoint PDEs which we approximate by the finite element method. Copyright © 2007 John Wiley & Sons, Ltd.     

Vector green's function method for unsteady Navier-Stokes equations

Summary The Green's function method is reformulated in general terms to treat vector unsteady and nonlinear equations. The particular expressions of the adjoint linear operator, the Green's formula and the integral representation of the solution are derived for unsteady Navier Stokes equations. The appropriate fundamental solutions for incompressible and for certain compressible flows have been obtained in closed form. Both the positive features and the possible limits of the method are briefly outlined.

---

Sommario Si riformula in termini generali il metodo della funzione di Green per estenderlo al caso di equazioni vettoriali e non lineari. In particolare si ricavano le espressioni della formula di Green e della rappresentazione integrale della soluzione per le equazioni di Navier Stokes non stazionarie. Si ottengono le soluzioni fondamentali in forma chiusa sia per il caso di fluidi comprimibili che incomprimibili. Si discutono infine brevemente i lati positivi ed i possibili limiti della metodologia illustrata.

     

Objective function choice for control of a thermocapillary flow using an adjoint-based control strategy

The problem of suppressing flow oscillations in a thermocapillary flow is addressed using a gradient-based control strategy. The physical problem addressed is the “open boat” process of crystal growth, the flow in which is driven by thermocapillary and buoyancy effects. The problem is modeled by the two-dimensional unsteady incompressible Navier–Stokes and energy equations under the Boussinesq approximation. The goal of the control is to suppress flow oscillations which arise when the driving forces are such that the flow becomes unsteady. The control is a spatially and temporally varying temperature gradient boundary condition at the free surface. The control which minimizes the flow oscillations is found using a conjugate gradient method, where the gradient of the objective function with respect to the control variables is obtained from solving a set of adjoint equations. The issue of choosing an objective function that can be both optimized in a computationally efficient manner and optimization of which provides control that damps the flow oscillations is investigated. Almost complete suppression of the flow oscillations is obtained for certain choices of the objective function.     

Implicit application of non‐reflective boundary conditions for Navier–Stokes equations in generalized coordinates

The non‐reflective boundary conditions (NRBC) for Navier–Stokes equations originally suggested by Poinsot and Lele (J. Comput. Phys. 1992; 101 :104–129) in Cartesian coordinates are extended to generalized coordinates. The characteristic form Navier–Stokes equations in conservative variables are given. In this characteristic‐based method, the NRBC is implicitly coupled with the Navier–Stokes flow solver and are solved simultaneously with the flow solver. The calculations are conducted for a subsonic vortex propagating flow and the steady and unsteady transonic inlet‐diffuser flows. The results indicate that the present method is accurate and robust, and the NRBC are essential for unsteady flow calculations. Copyright © 2005 John Wiley & Sons, Ltd.     

Direct sensitivity analysis for smooth unsteady compressible flows using complex differentiation

A method for the direct computation of the instantaneous sensitivities of unsteady compressible flows is proposed. It is based on the complex differentiation of the full compressible Navier–Stokes equations and does not require the storage of the unsteady flow solution to be differentiated. The method does not rely on any assumption on the basic Navier–Stokes solver, and can therefore be implemented in a straightforward way. The method is assessed on several cases, including a two‐dimensional subsonic mixing layer. It is observed that the sensitivity patterns can be interpreted thanks to Kovasznay's decomposition for perturbations in a compressible flow. Copyright © 2006 John Wiley & Sons, Ltd.     

Fast solvers for models of ICEO microfluidic flows

We demonstrate the performance of a fast computational algorithm for modeling the design of a microfluidic mixing device. The device uses an electrokinetic process, induced charge electroosmosis (J. Fluid Mech. 2004; 509 ), by which a flow through the device is driven by a set of polarizable obstacles in it. Its design is realized by manipulating the shape and orientation of the obstacles in order to maximize the amount of fluid mixing within the device. The computation entails the solution of a constrained optimization problem in which function evaluations require the numerical solution of a set of partial differential equations: a potential equation, the incompressible Navier–Stokes equations, and a mass‐transport equation. The most expensive component of the function evaluation (which must be performed at every step of an iteration for the optimization) is the solution of the Navier–Stokes equations. We show that by using some new robust algorithms for this task (SIAM J. Sci. Comput. 2002; 24 :237–256; J. Comput. Appl. Math. 2001; 128 :261–279), based on certain preconditioners that take advantage of the structure of the linearized problem, this computation can be done efficiently. Using this computational strategy, in conjunction with a derivative‐free pattern search algorithm for the optimization, applied to a finite element discretization of the problem, we are able to determine optimal configurations of microfluidic devices. Copyright © 2009 John Wiley & Sons, Ltd.     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

Functional and generalized separable solutions to unsteady Navier–Stokes equations

The paper studies unsteady Navier–Stokes equations with two space variables. It shows that the non-linear fourth-order equation for the stream function with three independent variables admits functional separable solutions described by a system of three partial differential equations with two independent variables. The system is found to have a number of exact solutions, which generate new classes of exact solutions to the Navier–Stokes equations. All these solutions involve two or more arbitrary functions of a single argument as well as a few free parameters. Many of the solutions are expressed in terms of elementary functions, provided that the arbitrary functions are also elementary; such solutions, having relatively simple form and presenting significant arbitrariness, can be especially useful for solving certain model problems and testing numerical and approximate analytical hydrodynamic methods. The paper uses the obtained results to describe some model unsteady flows of viscous incompressible fluids, including flows through a strip with permeable walls, flows through a strip with extrusion at the boundaries, flows onto a shrinking plane, and others. Some blow-up modes, which correspond to singular solutions, are discussed.     

Optimal shape design for the time‐dependent Navier–Stokes flow

This paper is concerned with the problem of shape optimization of two‐dimensional flows governed by the time‐dependent Navier–Stokes equations. We derive the structures of shape gradients for time‐dependent cost functionals by using the state derivative and its associated adjoint state. Finally, we apply a gradient‐type algorithm to our problem, and numerical examples show that our theory is useful for practical purposes and the proposed algorithm is feasible in low Reynolds number flows. Copyright © 2007 John Wiley & Sons, Ltd.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

Solution of the incompressible Navier–Stokes equations by the method of lines

The numerical method of lines (NUMOL) is a numerical technique used to solve efficiently partial differential equations. In this paper, the NUMOL is applied to the solution of the two‐dimensional unsteady Navier–Stokes equations for incompressible laminar flows in Cartesian coordinates. The Navier–Stokes equations are first discretized (in space) on a staggered grid as in the Marker and Cell scheme. The discretized Navier–Stokes equations form an index 2 system of differential algebraic equations, which are afterwards reduced to a system of ordinary differential equations (ODEs), using the discretized form of the continuity equation. The pressure field is computed solving a discrete pressure Poisson equation. Finally, the resulting ODEs are solved using the backward differentiation formulas. The proposed method is illustrated with Dirichlet boundary conditions through applications to the driven cavity flow and to the backward facing step flow. Copyright © 2015 John Wiley & Sons, Ltd.     

Development and assessment of an intrusive polynomial chaos expansion-based continuous adjoint method for shape optimization under uncertainties

This article contributes to the development of methods for shape optimization under uncertainties, associated with the flow conditions, based on intrusive Polynomial Chaos Expansion (iPCE) and continuous adjoint. The iPCE to the Navier–Stokes equations for laminar flows of incompressible fluids is developed to compute statistical moments of the Quantity of Interest which are, then, compared with those obtained through the Monte Carlo method. The optimization is carried out using a continuous adjoint-enabled, gradient-based loop. Two different formulations for the continuous adjoint to the iPCE PDEs are derived, programmed, and verified. Intrusive PCE methods for the computation of the statistical moments require mathematical development, derivation of a new system of governing equations and their numerical solution. The development is presented for a chaos order of two and two uncertain variables and can be used as a guide to those willing to extend this development to a different set of uncertain variables or chaos order. The developed method and software, programmed in OpenFOAM, is applied to two optimization problems pertaining to the flow around isolated airfoils with uncertain farfield conditions.     

A general methodology for investigating flow instabilities in complex geometries: application to natural convection in enclosures

This paper presents a general methodology for studying instabilities of natural convection flows enclosed in cavities of complex geometry. Different tools have been developed, consisting of time integration of the unsteady equations, steady state solving, and computation of the most unstable eigenmodes of the Jacobian and its adjoint. The methodology is validated in the classical differentially heated cavity, where the steady solution branch is followed for vary large values of the Rayleigh number and most unstable eigenmodes are computed at selected Rayleigh values. Its effectiveness for complex geometries is illustrated on a configuration consisting of a cavity with internal heated partitions. We finally propose to reduce the Navier–Stokes equations to a differential system by expanding the unsteady solution as the sum of the steady state solution and of a linear combination of the leading eigenmodes. The principle of the method is exposed and preliminary results are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Evaluation of an energy relaxation method for the simulation of unsteady,viscous, real gas flows

This study investigates a new energy relaxation method designed to capture the dynamics of unsteady, viscous, real gas flows governed by the compressible Navier–Stokes equations. We focus on real gas models accounting for inelastic molecular collisions and yielding temperature‐dependent heat capacities. The relaxed Navier–Stokes equations are discretized using a mixed finite volume/finite element method and a high‐order time integration scheme. The accuracy of the energy relaxation method is investigated on three test problems of increasing complexity: the advection of a periodic set of vortices, the interaction of a temperature spot with a weak shock, and finally, the interaction of a reflected shock with its trailing boundary layer in a shock tube. In all cases, the method is validated against benchmark solutions and the numerical errors resulting from both discretization and energy relaxation are assessed independently. Copyright © 2004 John Wiley & Sons, Ltd.     

An implicit meshless method for application in computational fluid dynamics

An implicit meshless scheme is developed for solving the Euler equations, as well as the laminar and Reynolds‐averaged Navier–Stokes equations. Spatial derivatives are approximated using a least squares method on clouds of points. The system of equations is linearised, and solved implicitly using approximate, analytical Jacobian matrices and a preconditioned Krylov subspace iterative method. The details of the spatial discretisation, linear solver and construction of the Jacobian matrix are discussed; and results that demonstrate the performance of the scheme are presented for steady and unsteady two dimensional fluid flows. Copyright © 2012 John Wiley & Sons, Ltd.     

Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.     

Development of new finite volume schemes on unstructured triangular grid for simulating the gas–liquid two‐phase flow

This paper presents a coupled finite volume inner doubly iterative efficient algorithm for linked equations (IDEAL) with level set method to simulate the incompressible gas–liquid two‐phase flows with moving interfaces on unstructured triangular grid. The finite volume IDEAL method on a collocated grid is employed to solve the incompressible two‐phase Navier–Stokes equations, and the level set method is used to capture the moving interfaces. For the sake of mass conservation, an effective second‐order accurate finite volume scheme is developed to solve the level set equation on triangular grid, which can be implemented much easier than the classical high‐order level set solvers. In this scheme, the value of level set function on the boundary of control volume is approximated using a linear combination of a high‐order Larangian interpolation and a second‐order upwind interpolation. By the rotating slotted disk and stretching and shrinking of a circular fluid element benchmark cases, the mass conservation and accuracy of the new scheme is verified. Then the coupled method is applied to two‐phase flows, including a 2D bubble rising problem and a 2D dam breaking problem. The computational results agree well with those reported in literatures and experimental data. Copyright © 2015 John Wiley & Sons, Ltd.     

A Chebyshev collocation method for the Navier–Stokes equations with application to double-diffusive convection

A Chebyshev collocation method for solving the unsteady two-dimensional Navier–Stokes equations in vorticity–streamfunction variables is presented and discussed. The discretization in time is obtained through a class of semi-implicit finite difference schemes. Thus at each time cycle the problem reduces to a Stokes-type problem which is solved by means of the influence matrix technique leading to the solution of Helmholtz-type equations with Dirichlet boundary conditions. Theoretical results on the stability of the method are given. Then a matrix diagonalization procedure for solving the algebraic system resulting from the Chebyshev collocation approximation of the Helmholtz equation is developed and its accuracy is tested. Numerical results are given for the Stokes and the Navier–Stokes equations. Finally the method is applied to a double-diffusive convection problem concerning the stability of a fluid stratified by salinity and heated from below.