288 research outputs found
Applied optimal shape design
AbstractThis paper is a short survey of optimal shape design (OSD) for fluids. OSD is an interesting field both mathematically and for industrial applications. Existence, sensitivity, correct discretization are important theoretical issues. Practical implementation issues for airplane designs are critical too.The paper is also a summary of the material covered in our recent book, Applied Optimal Shape Design, Oxford University Press, 2001
Optimum Shape Design Using Automatic Differentiation in Reverse Mode
This paper shows how to use automatic differentiation in reverse mode as a powerful tool in optimization procedures. It is also shown that for aerodynamic applications the gradients have to be as accurate as possible. In particular, the effect of having the exact gradient of he first or second order spatial discretization schemes is presented. We show that the loss of precision in the gradient affects not only the convergence, but also the final shape. Both two and three dimensional configurations of transonic and supersonic flows have been investigated. These cases involve up to several thousand control parameters
An example of a global shock fitting method
AbstractBy means of an example of a single first order equation, we show how a shock can be characterized as the solution of an optimization problem. the optimization problem is solved directly by a gradient calculation and the method of steepest descent. The result is a global shock fitting method which eliminates dispersion. Extensions to the case of a system of equations are discussed
Numerical studies of the Lagrangian approach for reconstruction of the conductivity in a waveguide
We consider an inverse problem of reconstructing the conductivity function in
a hyperbolic equation using single space-time domain noisy observations of the
solution on the backscattering boundary of the computational domain. We
formulate our inverse problem as an optimization problem and use Lagrangian
approach to minimize the corresponding Tikhonov functional. We present a
theorem of a local strong convexity of our functional and derive error
estimates between computed and regularized as well as exact solutions of this
functional, correspondingly. In numerical simulations we apply domain
decomposition finite element-finite difference method for minimization of the
Lagrangian. Our computational study shows efficiency of the proposed method in
the reconstruction of the conductivity function in three dimensions
What is the optimal shape of a pipe?
We consider an incompressible fluid in a three-dimensional pipe, following
the Navier-Stokes system with classical boundary conditions. We are interested
in the following question: is there any optimal shape for the criterion "energy
dissipated by the fluid"? Moreover, is the cylinder the optimal shape? We prove
that there exists an optimal shape in a reasonable class of admissible domains,
but the cylinder is not optimal. For that purpose, we explicit the first order
optimality condition, thanks to adjoint state and we prove that it is
impossible that the adjoint state be a solution of this over-determined system
when the domain is the cylinder. At last, we show some numerical simulations
for that problem
Continuum viscoplastic simulation of a granular column collapse on large slopes: ÎĽ(I) rheology and lateral wall effects
We simulate here dry granular flows resulting from the collapse of granular columns on an inclined channel (up to 22°) and compare precisely the results with laboratory experiments. Incompressibility is assumed despite the dilatancy observed in the experiments (up to 10%). The 2-D model is based on the so-called μ(I) rheology that induces a Drucker-Prager yield stress and a variable viscosity. A nonlinear Coulomb friction term, representing the friction on the lateral walls of the channel, is added to the model. We demonstrate that this term is crucial to accurately reproduce granular collapses on slopes ≳10°, whereas it remains of little effect on the horizontal slope. Quantitative comparison between the experimental and numerical changes with time of the thickness profiles and front velocity makes it possible to strongly constrain the rheology. In particular, we show that the use of a variable or a constant viscosity does not change significantly the results provided that these viscosities are of the same order. However, only a fine tuning of the constant viscosity (η=1 Pa s) makes it possible to predict the slow propagation phase observed experimentally at large slopes. Finally, we observed that small-scale instabilities develop when refining the mesh (also called ill-posed behavior, characterized in the work of Barker et al. [“Well-posed and ill-posed behaviour of the μ(I)-rheology for granular flow,” J. Fluid Mech. 779, 794–818 (2015)] and in the present work) associated with the mechanical model. The velocity field becomes stratified and the bands of high velocity gradient appear. These model instabilities are not avoided by using variable viscosity models such as the μ(I) rheology. However we show that the velocity range, the static-flowing transition, and the thickness profiles are almost not affected by them
Strong and auxiliary forms of the semi-Lagrangian method for incompressible flows
We present a review of the semi-Lagrangian method for advection-diusion and incompressible Navier-Stokes equations discretized with high-order methods. In particular, we compare the strong form where the departure points are computed directly via backwards integration with the auxiliary form where an auxiliary advection equation is solved instead; the latter is also referred to as Operator Integration Factor Splitting (OIFS) scheme. For intermediate size of time steps the auxiliary form is preferrable but for large time steps only the strong form is stable
- …