Haoxiang Luo

An efficient immersed boundary-lattice Boltzmann method for the hydrodynamic interaction of elastic filaments

We have introduced a modified penalty approach into the flow-structure interaction solver that combines an immersed boundary method (IBM) and a multi-block lattice Boltzmann method (LBM) to model an incompressible flow and elastic boundaries with finite mass. The effect of the solid structure is handled by the IBM in which the stress exerted by the structure on the fluid is spread onto the collocated grid points near the boundary. The fluid motion is obtained by solving the discrete lattice Boltzmann equation. The inertial force of the thin solid structure is incorporated by connecting this structure through virtual springs to a ghost structure with the equivalent mass. This treatment ameliorates the numerical instability issue encountered in this type of problems. Thanks to the superior efficiency of the IBM and LBM, the overall method is extremely fast for a class of flow-structure interaction problems where details of flow patterns need to be resolved. Numerical examples, including those involving multiple solid bodies, are presented to verify the method and illustrate its efficiency. As an application of the present method, an elastic filament flapping in the Kármán gait and the entrainment regions near a cylinder is studied to model fish swimming in these regions. Significant drag reduction is found for the filament, and the result is consistent with the metabolic cost measured experimentally for the live fish.

Fluid-structure interaction involving large deformations: 3D simulations and applications to biological systems

Three-dimensional fluid-structure interaction (FSI) involving large deformations of flexible bodies is common in biological systems, but accurate and efficient numerical approaches for modeling such systems are still scarce. In this work, we report a successful case of combining an existing immersed-boundary flow solver with a nonlinear finite-element solid-mechanics solver specifically for three-dimensional FSI simulations. This method represents a significant enhancement from the similar methods that are previously available. Based on the Cartesian grid, the viscous incompressible flow solver can handle boundaries of large displacements with simple mesh generation. The solid-mechanics solver has separate subroutines for analyzing general three-dimensional bodies and thin-walled structures composed of frames, membranes, and plates. Both geometric nonlinearity associated with large displacements and material nonlinearity associated with large strains are incorporated in the solver. The FSI is achieved through a strong coupling and partitioned approach. We perform several validation cases, and the results may be used to expand the currently limited database of FSI benchmark study. Finally, we demonstrate the versatility of the present method by applying it to the aerodynamics of elastic wings of insects and the flow-induced vocal fold vibration.

Effect of wing inertia on hovering performance of flexible flapping wings

Insect wings in flight typically deform under the combined aerodynamic force and wing inertia; whichever is dominant depends on the mass ratio defined as m∗=ρsh/(ρfc), where ρsh is the surface density of the wing, ρf is the density of the air, and c is the characteristic length of the wing. To study the differences that the wing inertia makes in the aerodynamic performance of the deformable wing, a two-dimensional numerical study is applied to simulate the flow-structure interaction of a flapping wing during hovering flight. The wing section is modeled as an elastic plate, which may experience nonlinear deformations while flapping. The effect of the wing inertia on lift production, drag resistance, and power consumption is studied for a range of wing rigidity. It is found that both inertia-induced deformation and flow-induced deformation can enhance lift of the wing. However, the flow-induced deformation, which corresponds to the low-mass wing, produces less drag and leads to higher aerodynamic power effi...

A Computational Study of the Effect of False Vocal Folds on Glottal Flow and Vocal Fold Vibration During Phonation

The false vocal folds are believed to be components of the acoustic filter that is responsible for shaping the voice. However, the effects of false vocal folds on the vocal fold vibration and the glottal aerodynamic during phonation remain unclear. This effect has implications for computational modeling of phonation as well as for understanding laryngeal pathologies such as glottal incompetence resulting from unilateral vocal fold paralysis. In this study, a high fidelity, two-dimensional computational model, which combines an immersed boundary method for the airflow and a continuum, finite-element method for the vocal folds, is used to examine the effect of the false vocal folds on flow-induced vibration (FIV) of the true vocal folds and the dynamics of the glottal jet. The model is notionally based on a laryngeal CT scan and employs realistic flow conditions and tissue properties. Results show that the false vocal folds potentially have a significant impact on phonation. The false vocal folds reduce the glottal flow impedance and increase the amplitude as well as the mean glottal jet velocity. The false vocal folds also enhance the intensity of the monopole acoustic sources in the glottis. A mechanism for reduction in flow impedance due to the false vocal folds is proposed.

Coupling modes of three filaments in side-by-side arrangement

A viscous flow past three filaments in side-by-side arrangement is studied by a numerical simulation and is accompanied by a previously established linear stability analysis. Other than the coupling modes reported previously, which include the in-phase mode, symmetrical mode, and out-of-phase mode, three additional modes are identified in the nonlinear regime by systematically varying the separation distance between the filaments. These modes are the half-frequency mode, irrational-frequency mode, and erratic flapping state. The dynamic characteristics of each mode at the saturated state is described, including the flapping amplitude, frequency, drag force, and mechanical energy of the filaments. Four typical vortex structures are observed in the wake of the filaments and are described as the coalesced vortices, symmetrical vortices, erratic vortices, and independent vortex streets. The vortex structures are closely related to the coupling modes and the dynamic characteristics of the filaments. As the Rey...

Three-dimensional flow and lift characteristics of a hovering ruby-throated hummingbird.

A three-dimensional computational fluid dynamics simulation is performed for a ruby-throated hummingbird (Archilochus colubris) in hovering flight. Realistic wing kinematics are adopted in the numerical model by reconstructing the wing motion from high-speed imaging data of the bird. Lift history and the three-dimensional flow pattern around the wing in full stroke cycles are captured in the simulation. Significant asymmetry is observed for lift production within a stroke cycle. In particular, the downstroke generates about 2.5 times as much vertical force as the upstroke, a result that confirms the estimate based on the measurement of the circulation in a previous experimental study. Associated with lift production is the similar power imbalance between the two half strokes. Further analysis shows that in addition to the angle of attack, wing velocity and surface area, drag-based force and wing–wake interaction also contribute significantly to the lift asymmetry. Though the wing–wake interaction could be beneficial for lift enhancement, the isolated stroke simulation shows that this benefit is buried by other opposing effects, e.g. presence of downwash. The leading-edge vortex is stable during the downstroke but may shed during the upstroke. Finally, the full-body simulation result shows that the effects of wing–wing interaction and wing–body interaction are small.

Analysis of flow-structure interaction in the larynx during phonation using an immersed-boundary method

A recently developed immersed-boundary method is used to model the flow-structure interaction associated with the human phonation. The glottal airflow is modeled as a two-dimensional incompressible flow driven by a constant subglottal pressure, and the vocal folds are modeled as a pair of three-layered, two-dimensional, viscoelastic structures. Both the fluid dynamics and viscoelasticity are solved on fixed Cartesian grids using a sharp-interface immersed boundary method. It is found that the vibration mode and frequency of the vocal fold model are associated with the eigenmodes of the structures, and that the transition of the vibration mode takes place during onset of the sustained vibration. The computed glottal waveforms of the volume flux, velocity, and pressure are reasonably realistic. The glottal flow features an unsteady jet whose direction is deflected by the large-scale vortices in the supraglottal region. A detailed analysis of the flow and vocal fold vibrations is conducted in order to gain insights into the biomechanics of phonation.

Stability of axisymmetric core–annular flow in the presence of an insoluble surfactant

The effect of an insoluble surfactant on the stability of the core-annular flow of two immiscible fluids is investigated by a normal-mode linear analysis and by numerical simulations based on the immersed-interface method for axisymmetric perturbations. The results reveal that, although the Marangoni stress due to surfactant concentration variations is unable to initiate a new type of instability as in the case of two-dimensional two-layer channel flow, it does destabilize the interface by broadening the range of growing wavenumbers and by raising the growth rate of unstable perturbations. Numerical simulations for large-amplitude disturbances reveal that the surfactant plays an important role in determining the morphology of the interfacial structures developing in the nonlinear stages of the motion.

Numerical study on dielectrophoretic chaining of two ellipsoidal particles

Electric field-induced assembly of biological and synthetic particles has proven useful in two- and three-dimensional fabrication of composite materials, microwires, photonic crystals, artificial tissues, and more. Biological particles are typically irregularly shaped, and using non-spherical synthetic particles has the ability to expand current applications. However, there is much to be understood about the dielectrophoretic (DEP) interaction that takes place between particles of general shape. In this work, we numerically study the DEP interaction between two prolate spheroid particles suspended in an unbounded fluid. The boundary-element method (BEM) is applied to solve the coupled electric field, Stokes flow, and particle motion, and the DEP forces are obtained by integrating the Maxwell stress tensor over the particle surfaces. Effects of the initial configuration and aspect ratio are investigated. Results show that the particles go through a self-rotation process, that is, electro-orientation, while translating slowly to form a chain pair. The final formation resembles the chaining pattern observed previously in experiments using densely distributed ellipsoidal particles. Thus, the transient behavior and particle-particle interaction exhibited in the current study could be used as the fundamental mechanism to explain the phenomenon in the experiment.

Onset of instability of a flag in uniform flow

This paper numerically and analytically studies the onset of instability of a flag in uniform flow. The three-dimensional (3D) simulation is performed by using an immersed-boundary method coupled with a nonlinear finite element method. The global stability, bistability and instability are identified in the 3D simulations. The Squires theorem is extended to analyze the stability of the fluid-flag system with 3D initial perturbations. It is found that if a parallel flow around the flag admits an unstable 3D disturbance for a certain value of the flutter speed, then a two-dimensional (2D) disturbance at a lower flutter speed is also admitted. In addition, the growth rate of 2D disturbance is larger than that of the 3D disturbance.