simulation and interaction of fluid dynamics

In the fluid simulation, the fluids and their surroundings may greatly change properties such as shape and temperature simultaneously, and different surroundings would characterize different interactions, which would change the shape and motion of the fluids in different ways. On the other hand, interactions among fluid mixtures of different kinds would generate more comprehensive behavior. To investigate the interaction behavior in physically based simulation of fluids, it is of importance to build physically correct models to represent the varying interactions between fluids and the environments, as well as interactions among the mixtures. In this paper, we will make a simple review of the interactions, and focus on those most interesting to us, and model them with various physical solutions. In particular, more detail will be given on the simulation of miscible and immiscible binary mixtures. In some of the methods, it is advantageous to be taken with the graphics processing unit (GPU) to achieve real-time computation for middle-scale simulation.

Generation of internal and surface waves by seafloor movement in a two layer fluid system

Internal and surface waves generated by the deformations of the solid bed in a two layer fluid system of infinite lateral extent and uniform depth are investigated. An integral solution is developed for an arbitrary bed displacement on the basis of a linear approximation of the complete description of wave motion using a transform method (Laplace in time and Fourier in space) analogous to that used to study the generation of tsunamis by many researchers. The theoretical solutions are presented for three interesting specific deformations of the seafloor; the spatial variation of each seafloor displacement consists of a block section of the seafloor moving vertically either up or down while the time-displacement history of the block section is varied. The generation process and the profiles of the internal and surface waves for the case of the exponential bed movement are numerically illustrated, and the effects of the deformation parameters, densities and depths of the two layers on the solutions are discussed. As expected, the solutions derived from the present work include as special cases that obtained by Kervella et al. [Theor Comput Fluid Dyn 21:245-269, 2007] for tsunamis cased by an instantaneous seabed deformation and those presented by Hammack [J Fluid Mech 60:769-799, 1973] for the exponential and the half-sine bed displacements when the density of the upper fluid is taken as zero.

CFX与USAERO的水下机器人操纵性仿真计算研究

研究了借助计算流体力学软件CFX与USAERO对水下机器人的操纵性进行仿真计算的方法,并将仿真计算结果与模型试验结果对比,验证仿真计算的可行性,指出仿真计算的适用范围.该方法可以对某些现实中难以进行的操纵性试验进行仿真计算,对于研究水下机器人的操纵性尤其是在方案设计阶段的操纵性设计和降低试验成本有一定的现实意义.

A model of an isolated magnetic-flux tube in a stratified atmosphere

A two-dimensional model of a magnetic flux tube confined in a gravitational stratified atmosphere is discussed. The magnetic field in the flux tube is assumed to be force-free. By using the approximation of large scale height, the problem of a free boundary with nonlinear conditions may be reduced to one involving a fixed boundary. The two-dimensional features are obtained by applying the perturbation method and adopting the Luest-Schlueter model as the basic state. The results show that the configuration of a flux tube confined in a gravitational stratified atmosphere is divergent, and the more twisted the magnetic field, the more divergent is the flux tube.

Numerical experiments on one-dimensional model of turbulence

The initial-value problem of a forced Burgers equation is numerically solved by the Fourier expansion method. It is found that its solutions finally reach a steady state of 'laminar flow' which has no randomness and is stable to disturbances. Hence, strictly speaking, the so-called Burgers turbulence is not a turbulence. A new one-dimensional model is proposed to simulate the Navier-Stokes turbulence. A series of numerical experiments on this one-dimensional turbulence is made and is successful in obtaining Kolmogorov's (1941) k exp(-5/3) inertial-range spectrum. The (one-dimensional) Kolmogorov constant ranges from 0.5 to 0.65.

Numerical Simulation of R-M Instability

In order to capture shock waves and contact discontinuities in the field and easy to program with parallel computation a new algorithm is developed to solve the N-S equations for simulation of R-M instability problems. The method with group velocity control is used to suppress numerical oscillations, and an adaptive non-uniform mesh is used to get fine resolution. Numerical results for cylindrical shock-cylindrical interface interaction with a shock Mach number Ms=1.2 and Atwood number A=0.818, 0.961, 0.980 (the interior density of the interface/outer density p(1)/p(2) = 10, 50, 100, respectively), and for the planar shock-spherical interface interaction with Ms=1.2 and p(1)/p(2) = 14.28are presented. The effect of Atwood number and multi-mode initial perturbation on the R-M instability are studied. Multi-collisions of the reflected shock with the interface is a main reason of nonlinear development of the interface instability and formation of the spike-bubble structures In simulation with double mode perturbation vortex merging and second instability are found. After second instability the small vortex structures near the interface produced. It is important factor for turbulent mixing.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Operational Optimization of GaN Thin Film Growth Employing Numerical Simulation in a Showerhead MOCVD Reactor

A detailed reaction-tran sport model was studied in a showerhead reactor for metal organic chemical vapor deposition of GaN film by using computational fluid dynamics simulation. It was found that flat flow lines without swirl are crucial to improve the uniformity of the film growth, and thin temperature gradient above the suscptor can increase the film deposition rate. By above-mentioned research, we can employ higher h (the distance from the susceptor to the inlet), P (operational pressure) and the rate of susceptor rotation to improve the film growth.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Numerical simulation of principal tidal constituents in the South China Sea, Gulf of Tonkin and Gulf of Thailand

The principal tidal constituents M-2, S-2, K-1 and O-1 in the South China Sea, Gulf of Tonkin and Gulf of Thailand are simulated simultaneously using the numerical scheme of Kwok et al. (1995 Proceedings of the 1st Asian Computational Fluid Dynamics Conference, pp. 16-19). The average differences between the computed and observed harmonic constants are mostly within 5 cm and 10 degrees for amplitudes and phase-lags, respectively. The simulated tidal regimes in the present model are believed to be more accurate than the previous numerical results. Our studies confirm that a clockwise rotating M-2 amphidromic system lies in the southeast of the Gulf of Thailand and an S-2 amphidromic system at the near-shore area of the northeast South China Sea. The linear tidal energy equation developed by Garrett (1975 Deep-Sea Research 22, 23-35) is generalized to the nonlinear case. Based on the numerical results, the energy budgets in the South China Sea and its subareas, namely the Taiwan Strait, the Gulf of Tonkin, the Gulf of Thailand and the remaining area are investigated. The tidal motion in the Taiwan Strait is maintained mainly by the energy fluxes from the East China Sea for both semidiurnal and diurnal species and partially from the Luzon Strait for semidiurnal species. For the other parts of the South China Sea, the tidal motion is mainly maintained by the energy fluxes through the Luzon Strait. The energy inputs from the tide-generating force are negative for semidiurnal species and positive for diurnal species. (C) 1999 Elsevier Science Ltd. All rights reserved.

NUMERICAL ANALYSIS OF THE FLOWFIELD IN A SUPERSONIC COIL WITH AN INTERLEAVED JET CONFIGURATION AND ITS EFFECT ON THE GAIN DISTRIBUTION

In a supersonic chemical oxygen-iodine laser (COIL) operating without primary buffer gas, the features of flowfield have significant effects on the Laser efficiency and beam quality. In this paper three-dimensional, multi-species, chemically reactive CFD technology was used to study the flowfield in mixing nozzle implemented with a supersonic interleaving jet configuration. The features of the flowfield as well as its effect on the spatial distribution of small signal gain were analyzed.

Vortex dynamics in the studies of looping in tropical cyclone tracks

This paper proposes criteria for predicting the tendency of looping in tropical cyclone tracks using the approach of vortex dynamics. We model the asymmetric structure of a cyclone by a system of vortex patches. The evolution of such system of vortices is simulated by the method of contour dynamics. A new set of exact analytic formulas for contour dynamics calculations is derived, which is shown to be more computationally effective. Based on point-vortex models, we derive analytic formulas for the criteria of looping in a cyclone track. From numerical experiments, the simulated trajectories obtained from the point-vortex system and vortex patch system agree quite well. Hence, the looping criteria obtained from the point-vortex system can be applied by forecasters to stay alert for tendency of looping in a cyclone track. To demonstrate the applicability of the proposed criteria, the trajectory of Typhoon Yancy (9012), whose field data are available from ''TCM-90'', is simulated. The case study shows that the asymmetric structure similar to the pattern of a beta gyre is responsible for its recurvature when Yancy landed Fujian Province, China on 20 August 1990.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<fluid dynamics (DFD) algorithm. Analysis shows that the basic conservative equations for discrete cells are the Euler equations, NS- and diffusion parabolized (DP) NS equations. In this paper, a new multiscale-domain decomposition method is developed for the high Reynolds number flow. First, the whole domain is decomposed to different sub-domains with the different characteristic scales. Then the different dominant equation of all sub-domains is defined according to the diffusion parabolized (DP) theory of viscous flow. Finally these different equations are solved simultaneously in whole computational region. For numerical tests of high Reynolds numerical flows, two-dimensional supersonic flows over rearward and frontward steps as well as an interaction flow between shock wave and boundary layer were solved numerically. The pressure distributions and local coefficients of skin friction on the wall are given. The numerical results obtained by the multiscale-domain decomposition algorithm are well agreement with those by NS equations. Comparing with the usual method of solving the Navier-Stokes equations in the whole flow, under the same numerical accuracy, the present multiscale domain decomposition method decreases CPU consuming about 20% and reflects the physical mechanism of practical flow more accurately.

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.

Non-relativistic conformal symmetries in fluid mechanics

The symmetries of a free incompressible fluid span the Galilei group, augmented with independent dilations of space and time. When the fluid is compressible, the symmetry is enlarged to the expanded Schrodinger group, which also involves, in addition, Schrodinger expansions. While incompressible fluid dynamics can be derived as an appropriate non-relativistic limit of a conformally invariant relativistic theory, the recently discussed conformal Galilei group, obtained by contraction from the relativistic conformal group, is not a symmetry. This is explained by the subtleties of the non-relativistic limit.