A variation of the Method of Regularized Stokeslets (MRS) in three dimensions is developed for triangulated surfaces with a piecewise linear force density. The work extends the regularized Stokeslet ...segment methodology used for piecewise linear curves. By using analytic integration of the regularized Stokeslet kernel over the triangles, the regularization parameter ϵ is effectively decoupled from the spatial discretization of the surface. This is in contrast to the usual implementation of the method in which the regularization parameter is chosen for accuracy reasons to be about the same size as the spatial discretization. The validity of the method is demonstrated through several examples, including the flow around a rigidly translating/rotating sphere and in the squirmer model for ciliate self-propulsion. Notably, second order convergence in the spatial discretization for fixed ϵ is demonstrated. Considerations of mesh design and choice of regularization parameter are discussed, and the performance of the method is compared with existing quadrature-based implementations.
•A variation on the method of regularized Stokeslets in 3D is presented for triangulated surfaces with a piecewise linear force density.•The use of analytic integration effectively decouples the regularization parameter from the boundary discretization.•Second order convergence in the boundary discretization h is demonstrated numerically for fixed regularization ϵ≪h.
Summary
An integrated shape morphing and topology optimization approach based on the deformable simplicial complex methodology is developed to address Stokes and Navier‐Stokes flow problems. The ...optimized geometry is interpreted by a set of piecewise linear curves embedded in a well‐formed triangular mesh, resulting in a physically well‐defined interface between fluid and impermeable regions. The shape evolution is realized by deforming the curves while maintaining a high‐quality mesh through adaption of the mesh near the structural boundary, rather than performing global remeshing. Topological changes are allowed through hole merging or splitting of islands. The finite element discretization used provides smooth and stable optimized boundaries for simple energy dissipation objectives. However, for more advanced problems, boundary oscillations are observed due to conflicts between the objective function and the minimum length scale imposed by the meshing algorithm. A surface regularization scheme is introduced to circumvent this issue, which is specifically tailored for the deformable simplicial complex approach. In contrast to other filter‐based regularization techniques, the scheme does not introduce additional control variables, and at the same time, it is based on a rigorous sensitivity analysis. Several numerical examples are presented to demonstrate the applicability of the approach.
An integrated shape morphing and topology optimization approach based on the deformable simplicial complex (DSC) methodology is developed to address Stokes and Navier‐Stokes flow problems. The optimized geometry is interpreted by a set of piecewise linear curves embedded in a well‐formed triangular mesh, resulting in a physically well‐defined interface between fluid and impermeable regions. A surface regularization scheme is introduced to circumvent boundary oscillations and, hence, allow for stable convergence for advanced fluid design problems.
Phase field models recently gained a lot of interest in the context of tumour growth models. Typically Darcy-type flow models are coupled to Cahn–Hilliard equations. However, often Stokes or Brinkman ...flows are more appropriate flow models. We introduce and mathematically analyse a new Cahn–Hilliard–Brinkman model for tumour growth allowing for chemotaxis. Outflow boundary conditions are considered in order not to influence tumour growth by artificial boundary conditions. Existence of global-in-time weak solutions is shown in a very general setting.
A landmark of turbulence is the emergence of universal scaling laws, such as Kolmogorov’s E(q) ~ q−5∕3 scaling of the kinetic energy spectrum of inertial turbulence with the wavevector q. In recent ...years, active fluids have been shown to exhibit turbulent-like flows at low Reynolds number. However, the existence of universal scaling properties in these flows has remained unclear. To address this issue, here we propose a minimal defect-free hydrodynamic theory for two-dimensional active nematic fluids at vanishing Reynolds number. By means of large-scale simulations and analytical arguments, we show that the kinetic energy spectrum exhibits a universal scaling E(q) ~ q−1 at long wavelengths. We find that the energy injection due to activity has a peak at a characteristic length scale, which is selected by a nonlinear mechanism. In contrast to inertial turbulence, energy is entirely dissipated at the scale where it is injected, thus precluding energy cascades. Nevertheless, the non-local character of the Stokes flow establishes long-range velocity correlations, which lead to the scaling behaviour. We conclude that active nematic fluids define a distinct universality class of turbulence at low Reynolds number.Determining the properties that emerge from the equations that govern turbulent flow is a fundamental challenge in non-equilibrium physics. A hydrodynamic theory for two-dimensional active nematic fluids at vanishing Reynolds number is now put forward, revealing a universal scaling behaviour for this class of systems.
We present a new diffuse interface model for the dynamics of inextensible vesicles in a viscous fluid with inertial forces. A new feature of this work is the implementation of the local ...inextensibility condition in the diffuse interface context. Local inextensibility is enforced by using a local Lagrange multiplier, which provides the necessary tension force at the interface. We introduce a new equation for the local Lagrange multiplier whose solution essentially provides a harmonic extension of the multiplier off the interface while maintaining the local inextensibility constraint near the interface. We also develop a local relaxation scheme that dynamically corrects local stretching/compression errors thereby preventing their accumulation. Asymptotic analysis is presented that shows that our new system converges to a relaxed version of the inextensible sharp interface model. This is also verified numerically. To solve the equations, we use an adaptive finite element method with implicit coupling between the Navier–Stokes and the diffuse interface inextensibility equations. Numerical simulations of a single vesicle in a shear flow at different Reynolds numbers demonstrate that errors in enforcing local inextensibility may accumulate and lead to large differences in the dynamics in the tumbling regime and smaller differences in the inclination angle of vesicles in the tank-treading regime. The local relaxation algorithm is shown to prevent the accumulation of stretching and compression errors very effectively. Simulations of two vesicles in an extensional flow show that local inextensibility plays an important role when vesicles are in close proximity by inhibiting fluid drainage in the near contact region.
We propose and analyze an unfitted method for a dual-dual mixed formulation of a class of Stokes models with variable viscosity depending on the velocity gradient, in which the pseudoestress, the ...velocity and its gradient are the main unknowns. On a fluid domain Ω with curved boundary Γ we consider a Dirichlet boundary condition and employ an approach previously applied to the Stokes equations with constant viscosity, which consists of approximating Ω by a polyhedral computational subdomain Ωh, not necessarily fitting Ω, where a Galerkin method is applied to compute solution. Furthermore, to approximate the Dirichlet data on the computational boundary Γh, we make use of a transferring technique based on integrating the discrete velocity gradient. Then the associated Galerkin scheme can be defined by employing Raviart–Thomas of order k≥0 for the pseudostress, and discontinuous polynomials of degree k for the velocity and its gradient. For the a priori error analysis we provide suitable assumptions on the mesh near the boundary Γ ensuring that the associated Galerkin scheme is well-posed and optimally convergent with O(hk+1). Next, for the case when Γh is taken as a piecewise linear interpolation of Γ, we develop a reliable and quasi-efficient residual-based a posteriori error estimator. Numerical experiments verify our analysis and illustrate the performance of the associated a posteriori error indicator.
The thermocapillary migration of a concentric compound drop in an arbitrary viscous flow under the consideration of negligible Reynolds number is investigated. The thermocapillary effect refers to ...the migration of a drop under the influence of a temperature gradient. The thermal and hydrodynamic problems are examined. The thermal field is governed by the heat conduction equation whereas the hydrodynamic fluid velocities are governed by the linearized Navier–Stokes equations. Presence of temperature gradient results in variation of the interfacial tension which is assumed to depend on temperature linearly. Variation of interfacial gradient leads to the coupling of the hydrodynamic problem with the thermal problem through the boundary condition. A complete general solution of Stokes equations is utilized to obtain closed-form expressions for the velocity vector and pressure. The hydrodynamic forces acting on the compound drop are obtained and expressed in terms of Fax́en’s law. Some important asymptotic limiting cases of hydrodynamic drag are also derived. The hydrodynamic drag for cases of uniform flow, shear flow, and heat source with the known ambient flow are derived and it is found that in the case of shear flow, the hydrodynamic drag is contributed only by the thermal component and the shear flow has no effect on it. The obtained results for drag and torque in the limiting cases are in agreement with the existing results in the literature. Furthermore, the migration velocity of the compound drop is obtained by equating the hydrodynamic drag force to zero. The results obtained for migration velocity are explained with the aid of graphs. The migration velocity is found to be a monotonic function of the Marangoni number and the radius of the innermost drop.
Triangular surface-based 3D IIM (Immersed Interface Method) algorithms face major challenges due to the need to calculate surface derivative of jumps. This paper proposes a fast, easy-to-implement, ...surface-derivative-free of jumps, augmented IIM (AIIM) with triangulated surfaces for 3D Helmholtz interface problems for the first time, which combines the simplified AIIM with domain decomposed and embedding techniques. The computational domain is divided into sub-domains along the interface and the solutions of sub-domains are continuously extended into larger regular domains by embedding. The jumps in normal derivative of solution along the interfaces in the extended domains are introduced as unknowns to impose the original jump relations. The original problem is simplified into Helmholtz interface problems with constant coefficients by coupling them with the augmented equation, which is then solved using fast simplified AIIM. This approach eliminates the need to compute surface derivatives of jumps, making implementation of 3D IIM based on triangulated surfaces fairly simple. Numerical results demonstrate that the algorithm is efficient and can achieve the overall second-order accuracy.
•A novel, easy-to-implement and fast augmented simplified IIM is proposed for 3D Helmholtz interface problems.•It seems to be the first work on IIM with the triangular surface mesh indeed.•The method provides a fairly simple way to compute the correction terms without computing surface derivatives of jumps.•Various numerical examples verify the efficiency and accuracy of the method.