13,003 research outputs found
Shape optimization of Stokesian peristaltic pumps using boundary integral methods
This article presents a new boundary integral approach for finding optimal
shapes of peristaltic pumps that transport a viscous fluid. Formulas for
computing the shape derivatives of the standard cost functionals and
constraints are derived. They involve evaluating physical variables (traction,
pressure, etc.) on the boundary only. By emplyoing these formulas in conjuction
with a boundary integral approach for solving forward and adjoint problems, we
completely avoid the issue of volume remeshing when updating the pump shape as
the optimization proceeds. This leads to significant cost savings and we
demonstrate the performance on several numerical examples
Development of an unstructured solution adaptive method for the quasi-three-dimensional Euler and Navier-Stokes equations
A general solution adaptive scheme based on a remeshing technique is developed for solving the two-dimensional and quasi-three-dimensional Euler and Favre-averaged Navier-Stokes equations. The numerical scheme is formulated on an unstructured triangular mesh utilizing an edge-based pointer system which defines the edge connectivity of the mesh structure. Jameson's four-stage hybrid Runge-Kutta scheme is used to march the solution in time. The convergence rate is enhanced through the use of local time stepping and implicit residual averaging. As the solution evolves, the mesh is regenerated adaptively using flow field information. Mesh adaptation parameters are evaluated such that an estimated local numerical error is equally distributed over the whole domain. For inviscid flows, the present approach generates a complete unstructured triangular mesh using the advancing front method. For turbulent flows, the approach combines a local highly stretched structured triangular mesh in the boundary layer region with an unstructured mesh in the remaining regions to efficiently resolve the important flow features. One-equation and two-equation turbulence models are incorporated into the present unstructured approach. Results are presented for a wide range of flow problems including two-dimensional multi-element airfoils, two-dimensional cascades, and quasi-three-dimensional cascades. This approach is shown to gain flow resolution in the refined regions while achieving a great reduction in the computational effort and storage requirements since solution points are not wasted in regions where they are not required
A boundary element regularised Stokeslet method applied to cilia and flagella-driven flow
A boundary element implementation of the regularised Stokeslet method of
Cortez is applied to cilia and flagella-driven flows in biology.
Previously-published approaches implicitly combine the force discretisation and
the numerical quadrature used to evaluate boundary integrals. By contrast, a
boundary element method can be implemented by discretising the force using
basis functions, and calculating integrals using accurate numerical or analytic
integration. This substantially weakens the coupling of the mesh size for the
force and the regularisation parameter, and greatly reduces the number of
degrees of freedom required. When modelling a cilium or flagellum as a
one-dimensional filament, the regularisation parameter can be considered a
proxy for the body radius, as opposed to being a parameter used to minimise
numerical errors. Modelling a patch of cilia, it is found that: (1) For a fixed
number of cilia, reducing cilia spacing reduces transport. (2) For fixed patch
dimension, increasing cilia number increases the transport, up to a plateau at
cilia. Modelling a choanoflagellate cell it is found that the
presence of a lorica structure significantly affects transport and flow outside
the lorica, but does not significantly alter the force experienced by the
flagellum.Comment: 20 pages, 7 figures, postprin
- …