A numerical investigation on the flow physics of koi carp (Cyprinus carpio koi) burst-and-coast swimming

Abstract—Burst-and-coast is the most common locomotion type in freely routine swimming of koi carps (Cyprinus carpio koi), which consists of a burst phase and a coast phase in each cycle and mostly leads to a straight-line trajectory. Combining with the tracking experiment, the flow physics of koi carp’s burst-andcoast swimming is investigated using a novel integrated CFD method solving the body-fluid interaction problem. The dynamical equations of a deforming body are formulated. Following that, the loose-coupled equations of the body dynamics and the fluid dynamics are numerically solved with the integrated method. The two burst modes, MT (Multiple Tail-beat) and HT (Half Tail-beat), which have been reported by the experiments, are investigated by numerical simulations in this paper. The body kinematics is predicted and the flow physics is visualized, which are in good agreement with the corresponding experiments. Furthermore, the optimization on the energy cost and several critical control mechanisms in burst-and-coast swimming of koi carps are explored, by varying the parameters in its selfpropelled swimming. In this paper, energetics is measured by the two mechanical quantities, total output power CP and Froude efficiency Fr. Results and discussion show that from the standpoint of mechanical energy, burst-and-coast swimming does not actually save energy comparing with steady swimming at the same average speed, in that frequently changing of speed leads to decrease of efficiency.

Experimental chemical physics of droplets and clusters

Time-of-flight measurements of energetic He atoms, field ionization of cryogenic liquid helium clusters, and time-of-flight and REMPI spectroscopy of radical salt clusters were investigated experimentally. The excited He atoms were generated in a corona discharge. Two strong neutral peaks were observed, accompanied by a prompt photon peak and a charged peak. All peaks were correlated with the pulsing of the discharge. The neutral hyperthermal and metastable atoms were formed by different mechanisms at different stages of the corona discharge. Positively charged helium droplets were produced by ionization of liquid helium in an electrostatic spraying experiment. The fluid emerging from a thin glass capillary was ionized by a high voltage applied to a needle inside the capillary. Fine droplets (less than 10 µm in diameter) were produced in showers with currents as high as 0.4 µA at 2-4 kV. The high currents resulting from field ionization in helium and the low surface tension of He I, led to charge densities that greatly exceeded the Rayleigh limit, thus resulting in coulombic explosion of the liquid. In contrast, liquid nitrogen formed a well-defined Taylor cone with droplets having diameters comparable to the jet (≈100 µm) at lower currents (10 nA) and higher voltages (8 kV). The metal-halide clusters of calcium and chlorine were generated by laser ablation of calcium metal in a Ar/CCl₄ expansion. A visible spectrum of the Ca₂Cl₃ cluster was observed from 651 to 630 nm by 1 +1' REMPI. The spectra were composed of a strong origin band at 15 350.8 cm^-1 and several weak vibronic bands. Density functional calculations predicted three minimum energy isomers. The spectrum was assigned to the ²B₂ ← X ²A₁ transition of a planar C_2V structure having a ring of two Cl and two Ca atoms and a terminal Cl atom. The ring isomer of Ca₂Cl₃ has the unpaired electron localized on one Ca²⁺ ion to form a Ca⁺ chromophore. A second electronic band of Ca₂Cl₃ was observed at 720 nm. The band is sharply different from the 650 nm band and likely due to a different isomer.

Fluid-structure Interactions of Inverted Leaves and Flags

Interactions between fluid flows and elastic bodies are ubiquitous in nature. One such phenomena that is encountered on a daily basis is the flapping and fluttering of leaves in the wind. The fluid-structure interaction that governs the physics of a leaf in the wind is poorly understood at best and has potential applications in biomechanics, vehicle design, and energy conversion. We build upon previous work on the flapping dynamics of inverted flags, which are cantilevered elastic sheets with free leading edge and fixed trailing edge that display unique large amplitude oscillatory behaviors. We model a leaf in the laboratory using modified inverted flags, experimentally probing the governing parameters behind leaf fluttering as well as shedding light on the physics behind the inverted flag phenomena. The behavior of these "inverted leaves" studied here display sensitive dependence on two biomechanically relevant parameters, stem-to-leaf rigidity and stem-to-leaf length. In addition, leaves on a tree are not often found alone. We seek to understand the complex interactions of multiple fluttering and flapping leaves by way of examining the interactions between pairs of inverted flags. Coupling through their flow fields, pairs of inverted flags exhibit striking emergent phenomena. We report these observed dynamical behaviors and the conditions upon which they arise.

Quantum Simulator for Transport Phenomena in Fluid Flows

Transport phenomena still stand as one of the most challenging problems in computational physics. By exploiting the analogies between Dirac and lattice Boltzmann equations, we develop a quantum simulator based on pseudospin-boson quantum systems, which is suitable for encoding fluid dynamics transport phenomena within a lattice kinetic formalism. It is shown that both the streaming and collision processes of lattice Boltzmann dynamics can be implemented with controlled quantum operations, using a heralded quantum protocol to encode non-unitary scattering processes. The proposed simulator is amenable to realization in controlled quantum platforms, such as ion-trap quantum computers or circuit quantum electrodynamics processors.

Linear stability of miscible two-fluid flow down an incline

We show that miscible two-layer free-surface flows of varying viscosity down an inclined substrate are different in their stability characteristics from both immiscible two-layer flows, and flows with viscosity gradients spanning the entire flow. New instability modes arise when the critical layer of the viscosity transport equation overlaps the viscosity gradient. A lubricating configuration with a less viscous wall layer is identified to be the most stabilizing at moderate miscibility (moderate Peclet numbers). This also is in contrast with the immiscible case, where the lubrication configuration is always destabilizing. The co-existence that we find under certain circumstances, of several growing overlap modes, the usual surface mode, and a Tollmien-Schlichting mode, presents interesting new possibilities for nonlinear breakdown. © 2013 AIP Publishing LLC.

Novel fluid grid and voidage calculation techniques for a discrete element model of a 3D cylindrical fluidized bed

A discrete element model (DEM) combined with computational fluid dynamics (CFD) was developed to model particle and fluid behaviour in 3D cylindrical fluidized beds. Novel techniques were developed to (1) keep fluid cells, defined in cylindrical coordinates, at a constant volume in order to ensure the conditions for validity of the volume-averaged fluid equations were satisfied and (2) smoothly and accurately measure voidage in arbitrarily shaped fluid cells. The new technique for calculating voidage was more stable than traditional techniques, also examined in the paper, whilst remaining computationally-effective. The model was validated by quantitative comparison with experimental results from the magnetic resonance imaging of a fluidised bed analysed to give time-averaged particle velocities. Comparisons were also made between theoretical determinations of slug rise velocity in a tall bed. It was concluded that the DEM-CFD model is able to investigate aspects of the underlying physics of fluidisation not readily investigated by experiment. © 2014 The Authors.

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.

Molecular weight dependence of phase Behavior of PEO/P(EO-b-DMS) blends: Application of Sanchez-Lacombe lattice fluid theory

Molecular weight dependence of phase separation behavior of the Poly (ethylene oxide) (PEO)/Poly(ethylene oxide-block-dimethylsiloxane) (P(EO-b-DMS)) blends was investigated by both experimental and theoretical methods. The cloud point curves of PEO/P(EO-b-DMS) blends were obtained by turbidity method. Based on Sanchez-Lacombe lattice fluid theory (SLLFT), the adjustable parameter, epsilon*(12)/k (quantifying the interaction energy between different components), was evaluated by fitting the experimental data in phase diagrams. To calculate the spinodals, binodals, and the volume changes of mixing for these blends, three modified combining rules of the scaling parameters for the block copolymer were introduced.

Pressure effects on the thermodynamics of trans-decahydronaphthalene/polystyrene polymer solutions: Application of the Sanchez-Lacombe lattice fluid theory

The cloud-point temperatures (T-cl's) of trans-decahydronaphthalene(TD)/polystyrene (PS, (M) over bar (w) = 270 000) solutions were determined by light scattering measurements over a range of temperatures (1-16degreesC), pressures (100-900 bar), and compositions (4.2-21.6 vol.-% polymer). The system phase separates upon cooling and T-cl was found to increase with rising pressure for constant composition. In the absence of special effects, this finding indicates positive excess volume for the mixing. Special attention was paid to the demixing temperatures as a function of pressure for different polymer solutions and the plots in the T-phi plane (where phi signifies volume fractions). The cloud-point curves of polymer solutions under different pressures were observed for different compositions, which demonstrated that pressure has a greater effect on the TD/PS solutions when far from the critical point as opposed to near the critical point. The Sanchez-Lacombe lattice fluid theory (SLLFT) was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalpy of mixing, and the volume changes of mixing. The calculated results show that modified PS scaling parameters can describe the thermodynamics of the TD/PS system well. Moreover the SLLFT describes the experimental results well.

Statistical thermodynamics of polydisperse polymer systems in the framework of lattice fluid model: Effect of molecular weight and its distribution on the spinodal in polymer solution

With the aid of thermodynamics of Gibbs, the expression of the spinodal was derived for the polydisperse polymer-solvent system in the framework of Sanchez-Lacombe Lattice Fluid Theory (SLLFT). For convenience, we considered that a model polydisperse polymer contains three sub-components. According to our calculation, the spinodal depends on both weight-average ((M) over bar (w)) and number-average ((M) over bar (n)) molecular weights of the polydisperse polymer, but the z-average molecular weight ((M) over bar (z)) dependence on the spinodal is invisible. The dependence of free volume on composition, temperature, molecular weight, and its distribution results in the effect of (M) over bar (n) on the spinodal. Moreover, it has been found that the effect of changing (M) over bar (w) on the spinodal is much bigger than that of changing (M) over bar (n) and the extrema of the spinodal increases with the rise of the weight-average molecular weight of the polymer in the solutions with upper critical solution temperature (UCST). However, the effect of polydispersity on the spinodal can be neglected for the polymer with a considerably high weight-average molecular weight. A more simple expression of the spinodal for the polydisperse polymer solution in the framework of SLLFT was also derived under the assumption of upsilon(*)=upsilon(1)(*)=upsilon(2)(*) and (1/r(1)(0))-(1/r(2i)(0))-->(1/r(1)(0)).