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.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Estudo de fluxos numéricos para modelos magneto-hidrodinâmicos utilizando o código Flash

A teoria magneto-hidrodinâmicos permite a estruturação de modelos computacionais, designados modelos MHDs, que são uma extensão da dinâmica dos fluidos para lidar com fluidos eletricamente carregados, tais como os plasmas, em que se precisa considerar os efeitos de forças eletromagnéticas. Tais modelos são especialmente úteis quando o movimento exato de uma partícula não é de interesse, sendo que as equações descrevem as evoluções de quantidades macroscópicas. Várias formas de modelos MHD têm sido amplamente utilizadas na Física Espacial para descrever muitos tipos diferentes de fenômenos de plasma, tais como reconexão magnética e interações de ventos estelares com diferentes objetos celestiais. Neste trabalho, o objetivo é analisar o comportamento de diversos fluxos numéricos em uma discretização de volumes finitos de um modelo numérico de MHD usando um esquema de malha entrelaçada sem separação direcional considerando alguns casos testes. Para as simulações, utiliza-se o código Flash, desenvolvido pela Universidade de Chicago, por ser um código de amplo interesse nas simulações astrofísicas e de fenômenos no espaço próximo à Terra. A metodologia consiste na inclusão de um fluxo numérico, permitindo melhoria com respeito ao esquema HLL.

Heat transfer and aerodynamics of turbine blade tips in a linear cascade

Local measurements of the heat transfer coefficient and pressure coefficient were conducted on the tip and near tip region of a generic turbine blade in a five-blade linear cascade. Two tip clearance gaps were used: 1.6% and 2.8% chord. Data was obtained at a Reynolds number of 2.3 × 10 5 based on exit velocity and chord. Three different tip geometries were investigated: a flat (plain) tip, a suction-side squealer, and a cavity squealer. The experiments reveal that the flow through the plain gap is dominated by flow separation at the pressure-side edge and that the highest levels of heat transfer are located where the flow reattaches on the tip surface. High heat transfer is also measured at locations where the tip-leakage vortex has impinged onto the suction surface of the aerofoil. The experiments are supported by flow visualisation computed using the CFX CFD code which has provided insight into the fluid dynamics within the gap. The suction-side and cavity squealers are shown to reduce the heat transfer in the gap but high levels of heat transfer are associated with locations of impingement, identified using the flow visualisation and aerodynamic data. Film cooling is introduced on the plain tip at locations near the pressure-side edge within the separated region and a net heat flux reduction analysis is used to quantify the performance of the successful cooling design. copyright © 2005 by ASME.

Swimming direction reversal of flagella through ciliary motion of mastigonemes.

Bio-inspired designs can provide an answer to engineering problems such as swimming strategies at the micron or nano-scale. Scientists are now designing artificial micro-swimmers that can mimic flagella-powered swimming of micro-organisms. In an application such as lab-on-a-chip in which micro-object manipulation in small flow geometries could be achieved by micro-swimmers, control of the swimming direction becomes an important aspect for retrieval and control of the micro-swimmer. A bio-inspired approach for swimming direction reversal (a flagellum bearing mastigonemes) can be used to design such a system and is being explored in the present work. We analyze the system using a computational framework in which the equations of solid mechanics and fluid dynamics are solved simultaneously. The fluid dynamics of Stokes flow is represented by a 2D Stokeslets approach while the solid mechanics behavior is realized using Euler-Bernoulli beam elements. The working principle of a flagellum bearing mastigonemes can be broken up into two parts: (1) the contribution of the base flagellum and (2) the contribution of mastigonemes, which act like cilia. These contributions are counteractive, and the net motion (velocity and direction) is a superposition of the two. In the present work, we also perform a dimensional analysis to understand the underlying physics associated with the system parameters such as the height of the mastigonemes, the number of mastigonemes, the flagellar wave length and amplitude, the flagellum length, and mastigonemes rigidity. Our results provide fundamental physical insight on the swimming of a flagellum with mastigonemes, and it provides guidelines for the design of artificial flagellar systems.

Numerical modelling of chirality-Induced bi-Directional swimming of artificial flagella

Biomimetic micro-swimmers can be used for various medical applications, such as targeted drug delivery and micro-object (e.g. biological cells) manipulation, in lab-on-a-chip devices. Bacteria swim using a bundle of flagella (flexible hair-like structures) that form a rotating cork-screw of chiral shape. To mimic bacterial swimming, we employ a computational approach to design a bacterial (chirality-induced) swimmer whose chiral shape and rotational velocity can be controlled by an external magnetic field. In our model, we numerically solve the coupled governing equations that describe the system dynamics (i.e. solid mechanics, fluid dynamics and magnetostatics). We explore the swimming response as a function of the characteristic dimensionless parameters and put special emphasis on controlling the swimming direction. Our results provide fundamental physical insight on the chirality-induced propulsion, and it provides guidelines for the design of magnetic bi-directional micro-swimmers. © 2013 The Author(s) Published by the Royal Society. All rights reserved.

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.

A numerical study of automotive turbocharger mixed flow turbine inlet geometry for off design performance

Mixed flow turbines represent a potential solution to the increasing requirement for high pressure, low velocity ratio operation in turbocharger applications. While literature exists for the use of these turbines at such operating conditions, there is a lack of detailed design guidance for defining the basic geometry of the turbine, in particular, the cone angle – the angle at which the inlet of the mixed flow turbine is inclined to the axis. This investigates the effect and interaction of such mixed flow turbine design parameters.
Computational Fluids Dynamics was initially used to investigate the performance of a modern radial turbine to create a baseline for subsequent mixed flow designs. Existing experimental data was used to validate this model.
Using the CFD model, a number of mixed flow turbine designs were investigated. These included studies varying the cone angle and the associated inlet blade angle.
The results of this analysis provide insight into the performance of a mixed flow turbine with respect to cone and inlet blade angle.

Modelling and aerodynamic simulation of a vehicle and failure analysis of a laminated front fender

Master Thesis in Mechanical Engineering field of Maintenance and Production

Chaotic dynamics of electric-field domains in periodically driven superlattices

Self-sustained time-dependent current oscillations under dc voltage bias have been observed in recent experiments on n-doped semiconductor superlattices with sequential resonant tunneling. The current oscillations are caused by the motion and recycling of the domain wall separating low- and high-electric-field regions of the superlattice, as the analysis of a discrete drift model shows and experimental evidence supports. Numerical simulation shows that different nonlinear dynamical regimes of the domain wall appear when an external microwave signal is superimposed on the dc bias and its driving frequency and driving amplitude vary. On the frequency-amplitude parameter plane, there are regions of entrainment and quasiperiodicity forming Arnold tongues. Chaos is demonstrated to appear at the boundaries of the tongues and in the regions where they overlap. Coexistence of up to four electric-field domains randomly nucleated in space is detected under ac+dc driving.

Uso de la matriz extracelular como material para ingeniería de tejidos

Introducción: La evaluación de injertos vasculares de submucosa de intestino delgado para la regeneración de vasos sanguíneos ha producido una permeabilidad variable (0-100%) que ha sido concurrente con la variabilidad en las técnicas de fabricación. Metodología: Investigamos los efectos de fabricación en permeabilidad y regeneración en un diseño experimental de 22factorial que combino: 1) preservación (P) o remoción (R) de la capa estratum compactum del intestino, y 2) deshidratada (D) o hidratada (H), dentro de cuatro grupos de estudio (PD, RD, PH, RH). Los injertos fueron implantados en las Arterias Carótidas de porcinos (ID 4.5mm, N=4, 7d). Permeabilidad, trombogenicidad, reacción inflamatoria, vascularización, infiltración de fibroblastos, perfil de polarización de macrófagos y fuerza tensil biaxial fueron evaluadas. Resultados: Todos los injertos PD permanecieron permeables (4/4), pero tuvieron escasa vascularización e infiltración de fibroblastos. El grupo RD permaneció permeable (4/4), presentó una extensa vascularización e infiltración de fibroblastos, y el mayor número del fenotipo de macrófagos (M2) asociado a regeneración. El grupo RH presentó menor permeabilidad (3/4), una extensa vascularización e infiltración de fibroblastos, y un perfil dominante de M2. El grupo PH presentó el menor grado de permeabilidad, y a pesar de mayor infiltración celular que PD, exhibió un fenotipo de macrófagos dominante adverso. La elasticidad de los injertos R evolucionó de una manera similar a las Carótidas nativas (particularmente RD, mientras que los injertos P mantuvieron su rigidez inicial. Discusión: Concluimos que los parámetros de fabricación afectan drásticamente los resultados, siendo los injertos RD los que arrojaron mejores resultados.

The frequency and dynamics of stratospheric sudden warmings in the 21st century

Changes to stratospheric sudden warmings (SSWs) over the coming century, as predicted by the Geophysical Fluid Dynamics Laboratory (GFDL) chemistry climate model [Atmospheric Model With Transport and Chemistry (AMTRAC)], are investigated in detail. Two sets of integrations, each a three-member ensemble, are analyzed. The first set is driven with observed climate forcings between 1960 and 2004; the second is driven with climate forcings from a coupled model run, including trace gas concentrations representing a midrange estimate of future anthropogenic emissions between 1990 and 2099. A small positive trend in the frequency of SSWs is found. This trend, amounting to 1 event/decade over a century, is statistically significant at the 90% confidence level and is consistent over the two sets of model integrations. Comparison of the model SSW climatology between the late 20th and 21st centuries shows that the increase is largest toward the end of the winter season. In contrast, the dynamical properties are not significantly altered in the coming century, despite the increase in SSW frequency. Owing to the intrinsic complexity of our model, the direct cause of the predicted trend in SSW frequency remains an open question.