• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Energetics comparison between zebrafish C-shaped turning and escape: self-propelled simulation with novel curvature models*

    2019-07-26 11:27:02WANGZhongweiYUYongliang

    WANG Zhongwei, YU Yongliang

    (School of Engineering Science, University of Chinese Academy of Sciences, Beijing 100049, China)

    Abstract The zebrafish C-start occurs not only in the escape response (C-escape), but also in the non-escape response (C-turn). They have different deforming modes and durations, and put in distinct kinematics performances. With the same bending amplitude, C-turns achieve larger turning angles, while C-escapes have higher velocities. However, little is known about their hydrodynamic mechanisms and energetics features. We proposed two novel curvature models based on the deformation characteristics of C-turn and C-escape, respectively. An optimization algorithm was used to determine the model parameters according to experimental data. Through self-propelled numerical simulation, we found that the positive moment at stage 1 of C-turn was large and the negative moment at stage 2 was small. So C-turn achieved a large turning angle. However, C-escape had a large thrust at stage 2, which led to the negative moment increasing and the turning angle reducing. We also found that the duration had little effect on the turning angle, but had a significant influence on escape velocity and energy consumption. Therefore, in order to save energy the duration of C-turn is usually over 100 ms, and the duration of C-escape is near 50 ms in order to pursue high speed.

    Keywords C-start; turning; curvature model; self-propelled swimming; energetics

    C-start is a type of burst maneuver pattern, in which the fish starts from rest and bends its body into a C shape, after which the tail flaps back and moves in a new direction[1]. It is widely believed that C-start occurs in the fish’s escape response[2]. In order to avoid predatory threats or external stimuli, fish can rapidly turn and escape by using C-start[3-4]. However, experimental studies have shown that, in the absence of predatory threats and external stimuli, the zebrafish also performs a C-start to make a turn starting from rest[5]. According to the purpose of swimming, this kind of C-start is called a C-turn, whereas the C-start in escape response is called a C-escape. Zebrafish C-turns and C-escapes have different deforming modes and durations, which result in the difference in the locomotion performance: with the same bending amplitude, the speed of C-turn is not as fast as that of C-escape, but it achieves a larger turning angle. In order to explain the difference in their swimming performance, it is necessary to compare the hydrodynamics and the energetics of zebrafish’s C-turn and C-escape.

    Researchers previously conducted extensive studies on the hydrodynamics of C-escape. Hu et. al. numerically studied the tail-flapping model of C-start[6]. Jing et al. conducted in-vivo experimental observations and theoretical analyses of C-escapes of crucian carp and yellow catfish[7-8]. Tytell and Lauder conducted the first in-vivo PIV experiment on C-escape and found that there were three jets in the flow field[9]. Witt et al.[10]studied the effects of flexibility on C-escape with simple physical models. Liu et al.[11]performed the self-propelled simulation of C-escape for the first time, and revealed the important role of the traveling curvature wave in the C-escape. Subsequently, Gazzola et al.[12]used evolutionary optimization to prove that C-escape is the best way to achieve maximum swimming distance. Li et al.[13-14]studied the effect of the head passing through a wake vortex on the trajectory of the locomotio. Borazjani et al.[15-16]studied the role of caudal and dorsal fins. Song et al.[17]studied the rapid turning problem in larval fish’s escape response. However, to the best of our knowledge, research on the hydrodynamics and energetics of zebrafish C-turn has not been reported.

    Self-propelled numerical simulation is an important method to study on hydrodynamics of C-starts, which has been used for researching C-escape[11-14]. During the self-propelled simulation, the deformation action of the fish body is needed as input conditions. For steady swimming, people commonly use sine waves to describe the undulation of the fish body[18-20]. For C-escape, due to the complexity of the deformation, many previous studies used the discrete experimental data directly[13-17]. However, this requires spatio-temporal interpolation and fairness, in which the calculation is easy to divergence and repeatability is poor. Moreover it is difficult to grasp the principal characters of deformation. Liu et al.[11]and Gazzola et al.[12]established a simplified curvature model for C-escape by superimposing sinusoidal traveling waves and standing waves. However this model may not be applicable to C-turns because of the difference in the deformation pattern between C-turns and C-escapes. How does one select the model parameters so that the model is in agreement with the experiment? This problem has not been resolved either.

    In this paper, we propose novel curvature models to describe the deformations of C-turn and C-escape of zebrafish. Parameters of curvature models are determined based on experimental data with an optimization algorithm. Then we perform self-propelled numerical simulation in order to explore the hydrodynamic mechanism and energetics features of C-turns and C-escapes. The influences of both the deforming mode and duration on the C-start swimming performance are discussed.

    1 Materials and methods

    1.1 The fish body model

    The contour of zebrafish’s body is symmetrical up and down, and it can be described by the function of the half width of the fish body:

    (1)

    wheresis the curve coordinate of the fish midline,L

    is the length of the fish body, andw1,w2,w3,w4,w5are coefficients whose values are fitted according to the real fish body contour and listed in Table 1. The cross section of the real fish body is approximately elliptical. According to the geometric properties of the ellipse, the average half width of the fish body isw(s)=wm(s)π/4. In this paper, a two-dimensional fish model is adopted, whose half width is equal to the average half width of the real fish.

    Table 1 Coefficients in the half width function

    1.2 Establishment of the curvature model

    Figure 1 gives experimental images of typical C-turn and C-escape of zebrafish, which are provided by Zhang Bingbing[5]. We can find that both the C-turn and C-escape have two stages, bending of the body and backward flapping of the tail, which is the main deformation characteristic of the C-start. The difference is that C-turns only flap back to the straight-line state at most, and the body is always bent to only one side and there is no reverse bending. However, C-escape has reversed bending after 43.7 ms. Domenici and Blake[3]classify C-start into single-bend and double-bend types based on the absence or presence of reverse bending. The C-turn is of the single-bend type, and the C-escape is of the double-bend type. In addition, the durations of C-turn and C-escape are significantly different, and the former is nearly three times the latter.

    Fig.1 Experimental images of zebrafish C-turn and C-escape

    The unilateral bending of C-turns means that the curvature of the fish midline is always positive. Therefore,we establish a curvature space modelκ(s) similar to a solitary wave as shown in Fig. 2(a). It consists of five parameters.Ais the amplitude of curvature, ands1,s2,s3, ands4are the coordinates along the fish midline. This model can be described by the function

    κ(s)=f(s,s1,s2,s3,s4,0,A),

    (2)

    wherefis a piecewise function defined as follows:

    f(ζ,ζ1,ζ2,ζ3,ζ4,f1,f2)=

    f1+(f2-f1)[ξ-sin(2πξ)/2π],

    (3a)

    (3b)

    whereζis an independent variable, andζ1,ζ2,ζ3,ζ4,f1, andf2are parameters. The five parameters of the curvature space model are all related to time. Their tendency to change over time is approximately illustrated in Fig. 2(b), i.e., the curvature time model, which reflects the growth and propagation of curvature. They can also be described by the functionfas follows:

    A(t)=f(t,0,ta,T,T,T,Am),

    (4a)

    s1(t)=f(t,t1,10T,10T,10T,s11,s12),

    (4b)

    s2(t)=f(t,t2,T,T,T,s21,L),

    (4c)

    s3(t)=f(t,0,t3,T,T,s31,L),

    (4d)

    s4(t)=f(t,0,t4,T,T,s41,L),

    (4e)

    whereTis the duration of C-turns;ta,Am,t1,s11,t2,s21,t3,s31,t4, ands41are parameters of the curvature time model. Each curve has two parameters, and their meanings are shown in Fig. 2(b) in detail. In order to makes1(t)=L,s12should be:

    Fig.2 Curvature models and flow chart of the optimization

    s12=s11+(L-s11)/

    For the C-escape, it is seen in Fig. 1 that the deformation in the bendding phase is very similar to the C-turn. The curvature is established gradually and propagates backward. Liu et al.[11]and Gazzola et al.[12]established a C-escape curvature model by superimposing standing waves and sinusoidal traveling waves. However, we find that this model can not reflect the deformation of the zebrafish C-escape accurately, because the standing wave does not meet the characteristics of backward propagation of curvature. Even so, the idea of superimposing sinusoidal traveling waves is worth learning. We utilise the curvature model of C-turn (i.e., Eqs. (2)-(4)), and add a sinusoidal traveling wave to model the curvature of C-escape. The sine curvature traveling wave is defined as follows:

    (5a)

    (5b)

    (5c)

    whereAsinrepresents the amplitude of the sine wave,λis the wavelength,Tsinis the period,φ0is the initial phase,a0anda1are coefficients of the amplitude distribution functionσ(s),shindicates the head length of the fish, andtsandteare parameters in the time functionτ(t).

    1.3 Determination of model parameters

    There are 10 and 19 parameters in the C-turn and C-escape curvature models, respectively. The curvature models have to be consistent with the experiment. In order to identify their parameters, we used an optimization algorithm to fit the experimental data. First, we extract several discrete points on the midline of the fish body from the experimental image and convert them to the fish head coordinate system. Then we calculate the model midline with the curvature model in the fish head system. The optimal objective function is

    (6)

    The optimization process is shown in Fig. 2(c), wherepiis theithparameter andδiexpresses the iterative step length ofpi, whose value is 0.01 time the unit ofpi. For the C-turn, in each iteration, all the parameters are traversed in three cases,pi-δi,pi, andpi+δi. There are 310(59 049) cases, and the parameter combination of the least error is selected. If the deviation does not converge, then the iterations continue; If the opposite is true, then the iterations terminate. For the C-escape, the parameters of the ‘C-turn’ model are firstly determined by the above method, and then the parameters of sine wave. The parameters should meet certain physical constraints, such as 0

    Table 2 displays the parameters of the optimized curvature model. The midline of the fish body is calculated by the curvature model and compared with the experimental points under the fish head coordinate system, as shown in Fig. 3(a), 3(b). It is observed that there is some error in the model curve compared to the experiment. At the end of the C-turn, the fish body is not completely straightened, but the model curve is fully straightened. The error for C-escape is largest att/T=1/11, as shown in Fig. 3(b). However, as a whole, the curvature model describes the deformation of zebrafish C-turn and C-escape very well. Their errors are only 2.18% and 1.44%, respectively.

    Table 2 Values of the curvature model parameters

    Fig.3 Deformations of C-turn and C-escape in the fish head and center-of-mass frames

    Although the curvature model established in this paper has 19 parameters at most, it is much lower than the 204 discrete experimental points (12 maps, 17 points for each map). Moreover, the curvature model grasps the main features of deformation, and the curvature is smooth, which facilitates numerical simulation. The curvature model has practical guiding significance for robotic fish[21].

    1.4 Governing equations of the self-propelled swimming

    The self-propelled swimming of fish is jointly controlled by the equations of fluid mechanics anddeforming body dynamics. The fluid control equations are incompressible Navier-Stokes equations

    ·u=0,

    (7a)

    (7b)

    whereu,p,ρandνare the velocity, pressure, density, and kinematic viscosity of the fluid, respectively.

    In the fish center-of-mass coordinate system, the swimming velocity of the fishu, can be decomposed into the translational velocityuc, the rotation velocityΩ, and the deformation velocityu′:

    u=uc+Ω×r′+u′,

    (8)

    wherer′ is the coordinate of the fish body in the center-of-mass reference system. The deformation satisfies the linear and angular momentum conservation laws

    What Fig.3(a), 3(b) shows is under the head system. To satisfy Eq. (9), it should be transformed into the center-of-mass system, as shown in Fig. 3(c). The rotational angle of the center-of-mass reference system relative to the head system is computed using Newton’s iteration method according to Eq. (9b). The deformation speed is determined by the deformation under the center-of-mass system. The translation and rotation of the fish body are governed by Newton’s second law of motion

    The finite-volume method,artificial compression, and overlapping grids were used to solve the hydrodynamic equations. The Newton iteration method was used to couple the hydrodynamics and the deforming body dynamics. Those methods have been verified and validated, and has been used in the study of fishlike C-escape[11].

    1.5 Dimensionless parameters

    In the experiment shown in Fig.1, the fish body length isL=3.31 cm, the duration of the C-turn isTturn=145.6 ms, and the duration of the C-escape isTescape=53.4 ms. TakeLandTescapeas the characteristic length and time, then the characteristic velocity isU=L/Tescape=0.62 m/s, the Reynolds number isRe=UL/ν=2.05×104, and the dimensionless duration of C-turn is 2.73. Whereν=1.003 ×10-6m2/s is the kinematic viscosity coefficient of water. According to the two-dimensional fish body model defined by Eq. (1), dimensionless fish mass and initial moment of inertia are 6.36×10-2and 4.01×10-3, respectively.

    To measure hydrodynamic and energetics properties, we define the thrust, torque, and power coefficients as follows:

    Where thrust is the component of the fluid force along the direction of the fish movement, the sign of the torque indicates whether it is in the same direction with the angular speed, and Power is the mechanical power consumed during the fish’s self-propelled swimming, which is defined as

    (11)

    2 Results and Discussions

    Both thedeformation and duration of C-turns and C-escapes are different. In order to explore their effects separately, we investigated four cases, as shown in Table 3. Case1 and case3 are the C-turn and C-escape of zebrafish in the experiment shown in Fig. 1. Case2 and case4 have the same deformation with case1 and case2, respectively, but the duration is different.

    Table 3 The setups and numerical simulation results of the four cases

    2.1 Numerical simulation results and comparisons with experiment

    Figure 4 shows the trajectory of the fish midline of C-turn (case 1) and C-escape (case 3) obtained from experiments and numerical simulations. We use the position and angle of the fish head of the end moment to represent errors between numerical simulations and experiments. In Fig. 4(a), the dashed line is the fish midline of numerical simulation, and the black line is that of experiment. It can be seen that the head position error of C-turn is 9.3%, the turning angle error is 7.1% (1-78.5/84.5), and the position error of C-escape is 9.6%, the angle error is 4.8% (55.1/52.6-1). All the error was lower than 10%. Because what we used in this paper was a two-dimensional and uniform density fish model, and the pectoral and caudal fins, three-dimensional fish body and flow field, and the density distribution was not taken into consideration, the error within 10% is acceptable. Overall, numerical simulations were basically consistent with experiments, which reflected the kinematics characteristics of C-turns and C-escapes.

    Table 3 demonstrates that, case1 has the lowest energy consumption and the highest relative propulsion efficiency and relative turning efficiency, while case 3 achieves the fastest swimming velocity. They have different deformation and duration, whose influence will be discussed separately in the succeeding sections.

    2.2 Effects of deformation patterns

    C-turn (case 2) and C-escape (case 3) are numerically simulated with the same duration (53.4 ms). As shown in Table 3, C-turn still have greater turning angles and lower energy consumption than C-escape, even with the same duration. Figure 5 displays the time curves of the velocity, angular velocity, turning angle, thrust coefficient, torque coefficient and power coefficient of C-turn (case 2) and C-escape (case 3). According to the definition in Fig. 3,t/T=0-0.4 is the bending phase (stage 1), andt/T=0.4-1 is the backward flapping phase (stage 2).

    Figure 5 (a) shows that velocity of C-turn and C-escape is very close at stage 1. Aftert/T=0.4, velocity of C-escape increase rapidly, peak at aroundt/T=0.7 (17.5L/s), and then decline slightly. However, the velocity enhancement of C-turn is not as significant as C-escape, whose peak value is 11.2L/s at aroundt/T=0.9. At the end time, the speed of C-escape (13.9L/s) is 40.4% higher than that of C-turn (9.9L/s). The difference in speed is due to the thrust of the fluid. It can be seen from Fig. 5 (d), both the thrust coefficients of C-turn and C-escape peak at aboutt/T=0.57 in stage 2. The peak of the C-escape (0.46) is 2.4 times that of the C-turn (0.19).

    Fig.4 Comparison of the fish midline trajectory between the experiment and numerical simulation

    Fig.5 Time curves of C-turn (case 2) and C-escape (case 3)

    As shown in Fig.5 (e), both the torque coefficients of C-turn and C-escape present positive torque at first and then negative torque. Positive torque has a large number and last a short time, however negtive torque has a relative little number and last a long time. This leads the angular velocity in Fig. 5 (b) to raise rapidly first and then decrease slowly. Time integration of the angular velocity results in the change of turning angle in Fig. 5 (c). Although the positive torque peak of C-turn (0.12) is slightly lower than C-escape (0.14), the duration of positive torque of C-turn (t/T=0-0.3) is longer than C-escape (t/T=0-0.22), which result that C-turn has a greater peak angular velocity (3.6(°)/ms) than C-escape (3.1(°)/ms). On the other hand, C-escape has a larger negative torque peak (-0.1) than C-turn (-0.048 ), which results that the angular velocity of C-turn decline slowly and almost stays positive all the time, however the angular velocity of C-escape fall quickly and reduce to negative att/T=0.7. So the turning angle of C-turn increase constantly, while the turning angle of C-escape begins to decrease att/T=0.7. Figure 5 (c) shows that the gap of turning angle between C-turn and C-escape expands aftert/T=0.7. At last, the turning angle of C-turn (85.7°) is 46% larger than C-escape (58.7°).

    Figure 5 (f) shows that the peak power coefficient at stage 1 is much larger than stage 2. Integrate the power coefficient, we could see that the energy consumed by C-turn and C-escape at stage 1 account for 77% and 63% of the total energy, respectively. This indicates that energy is mainly consumed at stage 1 in C-start, especially in C-turn. On the other hand, both the peak power coefficient of C-escape at stage 1 and stage 2 are larger than C-turn. At stage 1 the peak power coefficient of C-escape (3.36) is 1.78 times that of C-turn (1.89), and at stage 2 the peak power coefficient of C-escape (0.94) is 2.76 times that of C-turn (0.34).

    In addition, we found that there were two peaks in each curve of C-turn at stage 1. This is due to the two mechanismsof C-bend: curvature growth and propagation, as shown in the curvature-time model in Fig. 2 (b). The twin peaks in the time curves of C-turn reflect the curvature growth and curvature propagation, respectively. Parametertaindicates the speed of curvature growth.ta=0.15T, 0.3Tfor C-turn and C-escape, respectively, as shown in Table 2. Due totaof C-turn is small, before the curvature of the C-turn begin to propagate, it has reached a peak. However,taof C-escape is large, so when the curvature starts to propagate, it’s still increasing, and the growth and spread of curvature overlap with each other.

    In order to further understand the hydrodynamics and energetics mechanism of C-turns and C-escapes, we ploted the vorticity of the flow field, the center-of-mass reference system and the fluid stress distribution in Fig.6. Since the Reynolds number in this paper is 2.05×104, so the viscous stress is much smaller than the pressure stress and we only show the pressure distribution. The pressure on the two sides of the fish was concentrated to the midline, and the length of the arrow indicates the relative magnitude of the pressure.

    Fig.6 The vorticity field, the fluid stress distribution, and the center-of-mass frame

    As was stated above, we decomposed the fish swimming into translation, rotation and deformation. The translation and rotation is the motion of the center-of-mass coordinate system in Fig. 6, and the deformation is under the center-of-mass system, which is exactly the same with that in Fig. 3(c). The so-called self-propelled numerical simulation is to determine the acceleration and angular acceleration of the fish body, and then the velocity, angular velocity, position and angle of the center-of-mass system, according to the fluid force and torque, at the given predetermined deformation under the center-of-mass system.

    Figure 6 demonstrates that C-turn and C-escape have similar flow field. At stage 1, a positive vortex sheds from the head and a negative vortex sheds from the tail. At stage 2, another positive vortex sheds from the tail. However, due to the reversed bending at stage 2 of C-escape, a negative vortex generates att/T=0.9 near to the tail. Moreover, the positive vortex of C-escape near to the tail is stronger than that of C-turn (for example,t/T=0.5, 0.6), and the time of shedding is earlier. The positive vortex sheds from the tail completely att/T=0.8 in C-escape, while at the end time in C-turn. The vortex will result in the generation of jet in the flow field, which will induce the fluid force and torque on the fish.

    In Fig.5(f), power coefficient peaks att/T=0.2 in C-turn, and att/T=0.1 in C-escape. Form Fig. 6 we found that the fluid stress also peaks att/T=0.2 and 0.1 in C-turn and C-escape, respectively. This indicates that the power at stage 1 is mainly used to overcome the fluid stress generated by the C-shaped bending. At stage 1, there are three alternating reversed fluid stress on the head, the middle of the fish body, and the tail. Att/T=0-0.4 in C-turn andt/T=0-0.2 in C-escape, the stress on the head causes a negative torque and the stress on the tail causes a positive torque. Although the stress on the head is greater, the range of stress on the tail is larger, and the tail is further away from the center-of-mass so the arm of force is longer. The positive torque on the tail will overcome the negative torque on the head, and the fish will subject to a positive torque overall. However, att/T=0.5-0.8 in C-turn andt/T=0.3-0.7 in C-escape, the direction of stress will reverse. The stress on the head will generate a positive torque and the stress on the tail will generate a negative torque, and the negative torque at the tail will dominate for the same reason as above. The thrust coefficient of C-escape peaks att/T=0.6 in Fig. 5(d), and the stress on the tail also reaches the peak of stage 2 and along the direction of body movement. The stress peak at stage 2 of C-escape is greater than C-turn, which leads to a large enhancement of the velocity in C-escape. However, the larger stress will generate a larger negative torque. Att/T=0.7, the stress on the tail in C-escape is almost perpendicular to the body, so the arm of force is large, and the negative torque reaches a peak. However, the stress on the tail in C-turn is smaller and the arm is short, so the negative torque of C-turn is smaller than C-escape.

    The Strouhal numbers (St=A/uT) in fish steady swimming is in the range of 0.25-0.35, where A is the width of the wake. In zebrafish C-start, we define A as the amplitude of the tail’s deformation in the center-of-mass frame, as shown in Fig. 3(c). Then the Strouhal number of C-turn and C-escape is 0.59 and 0.57, respectively, as listed in Table 4. Such Strouhal numbers are very high, which will induce a sharp thrust and a low propulsive efficiency. Strouhal numbers of C-turn and C-escape have little difference, and the velocity is proportional to the amplitude of the tail. This indicates that due to the reverse bending in C-escape the deformation amplitude of tail increases, and the velocity will increase accordingly.

    Table 4 Strouhal numbers of C-turn and C-escape

    2.3 Effects of duration

    The effect of duration on the velocity, the turning angle and the energy consumption of C-turn and C-escape were studied in this section.

    Table 3 shows that, if the duration of C-turn increase from 53.4 ms (case 2) to 145.6 ms (case 1), the swimming velocity will reduce by 69% (1-3.1/9.9), but the energy consumption can reduce by 86% (1-4.47/32.9), and the turning angle only reduce by 8% (1-78.5/85.7). However, if the duration of C-escape decrease from 145.6 ms (case 4) to 53.4 ms (case 3), the energy consumption will increase by 6.3 times (46.9/6.44-1), swimming velocity will increase by 2.2 times (13.9/4.3-1), and the turning angle only increase by 10% (58.7/53.3-1). This demonstrates that the duration has little effect on the turning angle, but has a significant effect on the velocity and energy consumption. Increasing duration of C-turn can significantly reduce energy consumption, and decreasing duration of C-escape can greatly increase the velocity. The reason may be that the turning angle is a dimensionless quantity, but the dimension of velocity (L/T) is proportional to the reciprocal of time and the dimension of energy consumption (ρL4/T2) is proportional to the square of time’s reciprocal.

    By comparing the four cases in Table 3, we find that the C-turn in reality (case 1) consumes the least amount of energy, and its relative propulsion and turning efficiency are the highest, but it has the slowest swimming velocity. However, the C-escape in reality (case 3) has the fastest swimming velocity, but meanwhile it consumes the most energy, and the relative propulsion and turning efficiency are the lowest. Moreover, the highest turning angle is in case 2, and the second is case1. This indicates that zebrafish pursues relatively high turning angles and to save energy in C-turn, while pursues the fastest swimming velocity at the cost of high energy consumption in C-escape.

    3 Conclusion

    In conclusion, we established models of zebrafish’s C-turn and C-escape. Parameters were derived with an optimization algorithm according to the experimental data. The deformation obtained by the curvature model accorded with the experiment very well. We performed self-propelled numerical simulation, and the patterns of the midline of fish was also consistent with experiment(the error was within 10%).

    By comparing C-turn and C-escape models, there are three points worth noting. Firstly,the difference in swimming velocity and turning angle between the C-turn and C-escape mainly occurred at stage 2. Secondly, the mechanism of the turning was that a large positive torque was generated at stage 1 and a relative small reverse torque was generated at stage 2, which led to the positive angular velocity to maintain for a long period of time and to obtain turning angle. Furthermore, the reverse bending of C-escape at stage 2 increased the tail’s deformation amplitude, resulting in a strong fluid stress and a high velocity. Meanwhile, the reverse bending also led to the negative torque to increase, which resulted in that the angular velocity decreased rapidly and the turning angle reduced.

    We conclude that the duration had little effect on the turning angle, but had a significant influence on the swimming speed and energy consumption. Improving the duration of C-turn significantly reduced the energy consumption and maintained a relatively high turning angle at the same time. Decreasing the duration of C-escape dramatically increased the swimming velocity, but it required much more energy.

    In fact, fishes need to adjust their turning angles and the swimming velocity to perform various swimming behavior,and fishes have different shapes in general. In the future, we will study the impact of different curvature model parameters and fish shapes on the performance of C-start.

    Supplementarymaterial: Original experimental images and animations of numerical simulation are available at https:∥pan.baidu.com/s/1TlOU3_-ycWdqdtXIUeBAVw.

    The authors thank Zhang Bingbing for providing zebrafish experimental images.

    久久97久久精品| 永久免费av网站大全| 国产亚洲精品第一综合不卡| 国产精品久久久久久精品电影小说| www.熟女人妻精品国产| 街头女战士在线观看网站| 中文字幕人妻丝袜制服| 深夜精品福利| 亚洲精华国产精华液的使用体验| 精品视频人人做人人爽| 欧美日韩亚洲高清精品| 精品卡一卡二卡四卡免费| 国产精品秋霞免费鲁丝片| 亚洲av成人精品一二三区| 蜜桃在线观看..| 九色亚洲精品在线播放| 欧美日韩国产mv在线观看视频| 中文精品一卡2卡3卡4更新| 99re6热这里在线精品视频| 欧美av亚洲av综合av国产av | 99国产精品免费福利视频| 日韩 欧美 亚洲 中文字幕| 天天躁夜夜躁狠狠躁躁| 亚洲三区欧美一区| 免费日韩欧美在线观看| 久久韩国三级中文字幕| 国产精品免费大片| 日本vs欧美在线观看视频| 只有这里有精品99| 亚洲第一区二区三区不卡| 熟女av电影| 国产免费又黄又爽又色| 国产精品蜜桃在线观看| 国产精品久久久久久精品电影小说| 99久久99久久久精品蜜桃| 青春草亚洲视频在线观看| 亚洲精品乱久久久久久| 婷婷色综合大香蕉| 午夜日本视频在线| 最近中文字幕2019免费版| 尾随美女入室| 无遮挡黄片免费观看| 美女视频免费永久观看网站| 国产片特级美女逼逼视频| 日日爽夜夜爽网站| 欧美久久黑人一区二区| 18禁裸乳无遮挡动漫免费视频| 日本色播在线视频| 男女下面插进去视频免费观看| 国产精品.久久久| 久久天堂一区二区三区四区| 日本91视频免费播放| 韩国av在线不卡| 午夜91福利影院| 久久青草综合色| 成人影院久久| 99久久精品国产亚洲精品| 哪个播放器可以免费观看大片| 满18在线观看网站| 成人亚洲欧美一区二区av| 国产又色又爽无遮挡免| 午夜福利视频在线观看免费| 飞空精品影院首页| 亚洲熟女毛片儿| 久久久久网色| 91精品国产国语对白视频| 青青草视频在线视频观看| 一区二区三区精品91| 国产免费现黄频在线看| 青春草视频在线免费观看| 久久韩国三级中文字幕| 国产男女内射视频| 视频在线观看一区二区三区| 天天添夜夜摸| 精品少妇黑人巨大在线播放| 国产淫语在线视频| 丁香六月天网| 中文字幕亚洲精品专区| 妹子高潮喷水视频| 一级a爱视频在线免费观看| 亚洲成人手机| 观看美女的网站| 男女边吃奶边做爰视频| 狂野欧美激情性xxxx| 亚洲人成电影观看| 午夜福利网站1000一区二区三区| 日本黄色日本黄色录像| 久久婷婷青草| 午夜福利视频在线观看免费| 天天躁夜夜躁狠狠躁躁| 热re99久久精品国产66热6| 国产老妇伦熟女老妇高清| 精品国产国语对白av| 天美传媒精品一区二区| 丰满乱子伦码专区| 免费久久久久久久精品成人欧美视频| 日韩制服骚丝袜av| 免费观看a级毛片全部| 一级,二级,三级黄色视频| 国产免费福利视频在线观看| 久久精品人人爽人人爽视色| 男女床上黄色一级片免费看| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产成人一精品久久久| 日韩免费高清中文字幕av| 日本欧美国产在线视频| 久久久久精品性色| 亚洲国产欧美一区二区综合| 国产精品国产三级国产专区5o| 国产精品免费视频内射| 午夜激情av网站| 国产精品熟女久久久久浪| 人妻人人澡人人爽人人| 王馨瑶露胸无遮挡在线观看| 亚洲,一卡二卡三卡| 日日撸夜夜添| 日韩精品有码人妻一区| 青草久久国产| 久久精品亚洲av国产电影网| 大话2 男鬼变身卡| 制服丝袜香蕉在线| 亚洲成人av在线免费| 深夜精品福利| 亚洲视频免费观看视频| 国产不卡av网站在线观看| 亚洲人成网站在线观看播放| 搡老岳熟女国产| 美女午夜性视频免费| 亚洲第一av免费看| www.自偷自拍.com| 乱人伦中国视频| 人人妻,人人澡人人爽秒播 | 国产熟女欧美一区二区| av.在线天堂| 高清黄色对白视频在线免费看| av.在线天堂| 国产午夜精品一二区理论片| 一个人免费看片子| 欧美激情极品国产一区二区三区| 波野结衣二区三区在线| 免费观看a级毛片全部| 成人漫画全彩无遮挡| 又黄又粗又硬又大视频| 欧美成人午夜精品| 日韩 亚洲 欧美在线| 成年人午夜在线观看视频| 狂野欧美激情性bbbbbb| 97人妻天天添夜夜摸| 亚洲一区中文字幕在线| 欧美精品亚洲一区二区| 又黄又粗又硬又大视频| 色婷婷久久久亚洲欧美| 咕卡用的链子| 亚洲国产最新在线播放| 亚洲欧美清纯卡通| 老司机影院成人| 狂野欧美激情性bbbbbb| 热re99久久精品国产66热6| 女的被弄到高潮叫床怎么办| www.自偷自拍.com| 巨乳人妻的诱惑在线观看| 国产精品香港三级国产av潘金莲 | 日本av手机在线免费观看| 观看美女的网站| 一级片'在线观看视频| 精品国产乱码久久久久久小说| 亚洲一码二码三码区别大吗| 亚洲国产日韩一区二区| 亚洲第一青青草原| 国产精品二区激情视频| 国产一区二区激情短视频 | 欧美国产精品一级二级三级| svipshipincom国产片| 国产精品女同一区二区软件| av不卡在线播放| 人人妻,人人澡人人爽秒播 | 欧美久久黑人一区二区| 亚洲精品久久午夜乱码| 亚洲精品av麻豆狂野| 午夜福利视频精品| 一级毛片我不卡| 日韩 亚洲 欧美在线| 妹子高潮喷水视频| 国产欧美日韩一区二区三区在线| 国产成人系列免费观看| 欧美在线黄色| 亚洲精华国产精华液的使用体验| 亚洲久久久国产精品| 九草在线视频观看| 亚洲成人av在线免费| 91国产中文字幕| 国产精品偷伦视频观看了| 日本色播在线视频| 久久久精品国产亚洲av高清涩受| 亚洲av成人精品一二三区| 日韩大片免费观看网站| 黄色怎么调成土黄色| 国产精品久久久久久精品电影小说| 超碰97精品在线观看| 亚洲三区欧美一区| 久久久久久人人人人人| 亚洲图色成人| 久久国产精品男人的天堂亚洲| 国产一区亚洲一区在线观看| 亚洲成国产人片在线观看| 国产精品 欧美亚洲| 两个人看的免费小视频| 黄色毛片三级朝国网站| 咕卡用的链子| 亚洲欧美中文字幕日韩二区| 少妇 在线观看| 亚洲精品国产区一区二| 乱人伦中国视频| 国产色婷婷99| 亚洲男人天堂网一区| 黄色视频不卡| 免费观看av网站的网址| 中文天堂在线官网| 久久精品亚洲av国产电影网| 观看av在线不卡| 国产精品亚洲av一区麻豆 | 国产一级毛片在线| 九九爱精品视频在线观看| 国产精品久久久av美女十八| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av涩爱| 国产熟女欧美一区二区| 国产精品久久久久久人妻精品电影 | 色视频在线一区二区三区| av卡一久久| 热99国产精品久久久久久7| 天堂俺去俺来也www色官网| 亚洲av在线观看美女高潮| 天天影视国产精品| 午夜福利,免费看| 免费高清在线观看日韩| 交换朋友夫妻互换小说| 日本午夜av视频| 十分钟在线观看高清视频www| 亚洲欧洲日产国产| 欧美日韩视频高清一区二区三区二| 亚洲五月色婷婷综合| 免费黄频网站在线观看国产| 成人免费观看视频高清| 超色免费av| 国产男人的电影天堂91| 最近中文字幕高清免费大全6| 欧美国产精品va在线观看不卡| 亚洲精品一区蜜桃| 亚洲成人国产一区在线观看 | 搡老岳熟女国产| 国产又色又爽无遮挡免| 欧美日韩一区二区视频在线观看视频在线| 婷婷色综合大香蕉| 久久久久人妻精品一区果冻| 国产精品熟女久久久久浪| 乱人伦中国视频| 十八禁网站网址无遮挡| 精品福利永久在线观看| 又大又黄又爽视频免费| 国产一区二区在线观看av| 欧美日韩成人在线一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产男女内射视频| 一区二区三区乱码不卡18| 亚洲欧美激情在线| 一区福利在线观看| 久久人妻熟女aⅴ| 精品久久久久久电影网| 啦啦啦中文免费视频观看日本| 欧美成人午夜精品| 久久国产精品男人的天堂亚洲| 女人久久www免费人成看片| 一级片免费观看大全| 国产麻豆69| 久久久久人妻精品一区果冻| 超碰成人久久| 人人妻,人人澡人人爽秒播 | 97精品久久久久久久久久精品| 午夜影院在线不卡| 免费高清在线观看日韩| 在线观看免费午夜福利视频| 我的亚洲天堂| 国产日韩欧美视频二区| 久久久精品免费免费高清| 久久久久久久久久久久大奶| 2018国产大陆天天弄谢| 久久久久久久久免费视频了| 少妇人妻精品综合一区二区| 少妇被粗大的猛进出69影院| 久久久亚洲精品成人影院| av网站在线播放免费| 欧美日韩精品网址| 日韩一区二区视频免费看| 国产精品一区二区在线观看99| 日本wwww免费看| 国精品久久久久久国模美| 侵犯人妻中文字幕一二三四区| 新久久久久国产一级毛片| 国产精品二区激情视频| 天天躁夜夜躁狠狠久久av| 中国三级夫妇交换| 亚洲综合精品二区| 国产精品久久久av美女十八| 亚洲av日韩精品久久久久久密 | 亚洲国产毛片av蜜桃av| 久久久久久久久久久免费av| 亚洲av中文av极速乱| 成人18禁高潮啪啪吃奶动态图| 如何舔出高潮| 一本一本久久a久久精品综合妖精| 亚洲精品,欧美精品| av在线播放精品| 欧美日韩福利视频一区二区| 久久青草综合色| 国产精品久久久久久精品古装| 制服丝袜香蕉在线| 亚洲国产av影院在线观看| 亚洲国产最新在线播放| 91成人精品电影| 国产免费福利视频在线观看| www.精华液| 在线观看www视频免费| 九草在线视频观看| 国产麻豆69| 中文字幕精品免费在线观看视频| 天天影视国产精品| 日韩av在线免费看完整版不卡| 黑人欧美特级aaaaaa片| 亚洲国产精品成人久久小说| tube8黄色片| 丝袜美腿诱惑在线| 亚洲国产av影院在线观看| 一二三四中文在线观看免费高清| 亚洲av综合色区一区| 伦理电影大哥的女人| av免费观看日本| 日本爱情动作片www.在线观看| 亚洲五月色婷婷综合| 国产午夜精品一二区理论片| 交换朋友夫妻互换小说| 中文字幕最新亚洲高清| 国产成人欧美| 国产免费现黄频在线看| 免费在线观看黄色视频的| 中文字幕最新亚洲高清| 久久久久国产精品人妻一区二区| 国产淫语在线视频| 国产成人精品久久二区二区91 | 制服人妻中文乱码| 国产亚洲av高清不卡| 成人午夜精彩视频在线观看| 国产女主播在线喷水免费视频网站| 操美女的视频在线观看| 99九九在线精品视频| 日韩 欧美 亚洲 中文字幕| 亚洲av日韩精品久久久久久密 | 久热这里只有精品99| 亚洲七黄色美女视频| 桃花免费在线播放| 午夜福利影视在线免费观看| 新久久久久国产一级毛片| 美女脱内裤让男人舔精品视频| 欧美人与性动交α欧美精品济南到| 天美传媒精品一区二区| 婷婷成人精品国产| 亚洲,欧美,日韩| 深夜精品福利| 肉色欧美久久久久久久蜜桃| 最近最新中文字幕免费大全7| 大片免费播放器 马上看| 亚洲图色成人| 国产高清不卡午夜福利| 日日撸夜夜添| 建设人人有责人人尽责人人享有的| 免费观看av网站的网址| 精品少妇久久久久久888优播| 丰满饥渴人妻一区二区三| 亚洲精品成人av观看孕妇| 国产一区二区激情短视频 | 亚洲少妇的诱惑av| 国产黄频视频在线观看| 国产精品免费视频内射| 久久人妻熟女aⅴ| 美女福利国产在线| 赤兔流量卡办理| 大码成人一级视频| 中文字幕制服av| 母亲3免费完整高清在线观看| av卡一久久| 欧美黄色片欧美黄色片| 啦啦啦在线观看免费高清www| 狠狠精品人妻久久久久久综合| 久久天躁狠狠躁夜夜2o2o | 在线免费观看不下载黄p国产| 国产精品av久久久久免费| 久久热在线av| 男女无遮挡免费网站观看| 又大又黄又爽视频免费| 国产日韩欧美视频二区| 最近最新中文字幕大全免费视频 | a 毛片基地| 精品国产一区二区三区久久久樱花| 亚洲情色 制服丝袜| 久久精品熟女亚洲av麻豆精品| 日韩一卡2卡3卡4卡2021年| 中文字幕亚洲精品专区| 久久99精品国语久久久| 9色porny在线观看| 一级爰片在线观看| 成人黄色视频免费在线看| 国产精品 国内视频| 久久久久国产一级毛片高清牌| 黄色视频在线播放观看不卡| 午夜福利,免费看| bbb黄色大片| 欧美 亚洲 国产 日韩一| 久久av网站| 一级毛片 在线播放| 国产精品一区二区精品视频观看| 久久av网站| 国产欧美日韩综合在线一区二区| netflix在线观看网站| 巨乳人妻的诱惑在线观看| 母亲3免费完整高清在线观看| 成人漫画全彩无遮挡| 亚洲精品国产av蜜桃| 狠狠婷婷综合久久久久久88av| 欧美国产精品一级二级三级| 天天躁夜夜躁狠狠躁躁| 欧美精品人与动牲交sv欧美| 久久韩国三级中文字幕| h视频一区二区三区| 国产精品偷伦视频观看了| 久久精品亚洲熟妇少妇任你| 在线精品无人区一区二区三| 成人国语在线视频| 亚洲精品中文字幕在线视频| 黑丝袜美女国产一区| 精品第一国产精品| 蜜桃国产av成人99| 肉色欧美久久久久久久蜜桃| 亚洲国产欧美在线一区| 国产精品免费视频内射| 国产成人精品在线电影| 黑人猛操日本美女一级片| avwww免费| 哪个播放器可以免费观看大片| 亚洲欧美精品综合一区二区三区| 永久免费av网站大全| 婷婷色综合www| 国产有黄有色有爽视频| 咕卡用的链子| 黄色怎么调成土黄色| 久久人人爽人人片av| 看十八女毛片水多多多| 中文精品一卡2卡3卡4更新| 亚洲成人一二三区av| 51午夜福利影视在线观看| 观看美女的网站| 国产深夜福利视频在线观看| 欧美日韩成人在线一区二区| 热99国产精品久久久久久7| 女性生殖器流出的白浆| 免费av中文字幕在线| 老汉色∧v一级毛片| 永久免费av网站大全| 国精品久久久久久国模美| h视频一区二区三区| 日韩制服骚丝袜av| 妹子高潮喷水视频| 老司机影院毛片| 在线观看国产h片| 欧美97在线视频| 日本wwww免费看| 天堂俺去俺来也www色官网| 精品免费久久久久久久清纯 | 久久亚洲国产成人精品v| 成人三级做爰电影| 999久久久国产精品视频| 亚洲一区二区三区欧美精品| 99久久99久久久精品蜜桃| 国产精品久久久人人做人人爽| 少妇人妻 视频| 国产精品一区二区精品视频观看| 啦啦啦 在线观看视频| 啦啦啦在线观看免费高清www| 18禁动态无遮挡网站| 一级片免费观看大全| 午夜日韩欧美国产| 欧美激情 高清一区二区三区| 亚洲精品乱久久久久久| 中文字幕制服av| 母亲3免费完整高清在线观看| 亚洲精品国产av成人精品| 亚洲精品美女久久久久99蜜臀 | 性少妇av在线| 久热爱精品视频在线9| 精品国产露脸久久av麻豆| 天堂8中文在线网| avwww免费| 操美女的视频在线观看| 纯流量卡能插随身wifi吗| 久久免费观看电影| netflix在线观看网站| 1024香蕉在线观看| 狂野欧美激情性xxxx| 女的被弄到高潮叫床怎么办| 电影成人av| 一级毛片我不卡| 亚洲,欧美,日韩| 国产精品三级大全| 亚洲精品一区蜜桃| 别揉我奶头~嗯~啊~动态视频 | 欧美精品一区二区大全| av卡一久久| 看十八女毛片水多多多| 国精品久久久久久国模美| 国产毛片在线视频| 制服人妻中文乱码| 久久久欧美国产精品| 操出白浆在线播放| 国产精品嫩草影院av在线观看| 精品国产露脸久久av麻豆| 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 亚洲av电影在线观看一区二区三区| 国产欧美日韩一区二区三区在线| 看十八女毛片水多多多| 国产片内射在线| 妹子高潮喷水视频| 成人亚洲欧美一区二区av| 一区二区三区四区激情视频| 国产一卡二卡三卡精品 | 亚洲四区av| 国产高清不卡午夜福利| 日日摸夜夜添夜夜爱| 99国产精品免费福利视频| 久久ye,这里只有精品| 精品视频人人做人人爽| 成年女人毛片免费观看观看9 | 国产成人精品久久二区二区91 | 性色av一级| 欧美精品一区二区大全| 老司机深夜福利视频在线观看 | 91精品国产国语对白视频| 国产毛片在线视频| 久久ye,这里只有精品| 日韩欧美一区视频在线观看| 18禁裸乳无遮挡动漫免费视频| 免费黄网站久久成人精品| 69精品国产乱码久久久| 亚洲精品中文字幕在线视频| av在线观看视频网站免费| 青春草国产在线视频| 亚洲精品自拍成人| 日韩 亚洲 欧美在线| 国产女主播在线喷水免费视频网站| 美国免费a级毛片| 大陆偷拍与自拍| 王馨瑶露胸无遮挡在线观看| www.av在线官网国产| 日本av手机在线免费观看| 亚洲,欧美精品.| svipshipincom国产片| avwww免费| 国产亚洲欧美精品永久| av在线老鸭窝| 色吧在线观看| 国产亚洲av片在线观看秒播厂| 亚洲伊人久久精品综合| 亚洲图色成人| 亚洲美女视频黄频| 日韩制服丝袜自拍偷拍| 国产国语露脸激情在线看| 欧美 日韩 精品 国产| 日本午夜av视频| 国产成人欧美| 视频区图区小说| 一个人免费看片子| 自拍欧美九色日韩亚洲蝌蚪91| h视频一区二区三区| 免费av中文字幕在线| 午夜福利乱码中文字幕| 国产精品免费大片| h视频一区二区三区| √禁漫天堂资源中文www| 两个人免费观看高清视频| 日韩av免费高清视频| 欧美在线黄色| 成年人免费黄色播放视频| 欧美日韩视频高清一区二区三区二| 久久av网站| 精品少妇久久久久久888优播| 日本黄色日本黄色录像| 亚洲,一卡二卡三卡| av有码第一页| 亚洲婷婷狠狠爱综合网| 波野结衣二区三区在线| av国产久精品久网站免费入址| 成人亚洲欧美一区二区av| 久久精品国产综合久久久| 国产av码专区亚洲av| 亚洲四区av| 日韩人妻精品一区2区三区| 亚洲欧美日韩另类电影网站| 国产极品天堂在线| 丰满饥渴人妻一区二区三| 制服诱惑二区| 中文字幕亚洲精品专区| 久久久久精品久久久久真实原创| 嫩草影视91久久| e午夜精品久久久久久久| 丁香六月天网| 久久久国产一区二区|