一种基于粒子的牛顿流体与粘弹性流体统一模拟方法

在越来越受到人们关注的基于物理流体动画领域,目前分别模拟牛顿流体或粘弹性流体的方法很多,但很少有统一模拟两者的方法.文中基于光滑粒子流体动力学方法,通过对传统纳维-斯托克斯方程添加弹性应力项,提出了一种新的统一模拟牛顿流体和粘弹性流体的方法.通过实验说明该方法不仅有效,易于实现,而且具有良好的可控性,仅仅通过调节参数就可以模拟不同粘弹性、不同类型的流体现象.

A domain decomposition algorithm for viscous/inviscid coupling

In fluid mechanics, it is well accepted that the Euler equation is one of the reduced forms of the Navier-Stokes equation by truncating the viscous effect. There are other truncation techniques currently being used in order to truncate the Navier-Stokes equation to a reduced form. This paper describes one such technique, suitable for adaptive domain decomposition methods for the solution of viscous flow problems. The physical domain of a viscous flow problem is partitioned into viscous and inviscid subdomains without overlapping regions, and the technique is embedded into a finite volume method. Some numerical results are provided for a flat plate and the NACA0012 aerofoil. Issues related to distributed computing are discussed.

Assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Numerical modelling of fluid-structure interaction with application in hemodynamics

The thesis deals with numerical algorithms for fluid-structure interaction problems with application in blood flow modelling. It starts with a short introduction on the mathematical description of incompressible viscous flow with non-Newtonian viscosity and a moving linear viscoelastic structure. The mathematical model consists of the generalized Navier-Stokes equation used for the description of fluid flow and the generalized string model for structure movement. The arbitrary Lagrangian-Eulerian approach is used in order to take into account moving computational domain. A part of the thesis is devoted to the discussion on the non-Newtonian behaviour of shear-thinning fluids, which is in our case blood, and derivation of two non-Newtonian models frequently used in the blood flow modelling. Further we give a brief overview on recent fluid-structure interaction schemes with discussion about the difficulties arising in numerical modelling of blood flow. Our main contribution lies in numerical and experimental study of a new loosely-coupled partitioned scheme called the kinematic splitting fluid-structure interaction algorithm. We present stability analysis for a coupled problem of non-Newtonian shear-dependent fluids in moving domains with viscoelastic boundaries. Here, we assume both, the nonlinearity in convective as well is diffusive term. We analyse the convergence of proposed numerical scheme for a simplified fluid model of the Oseen type. Moreover, we present series of experiments including numerical error analysis, comparison of hemodynamic parameters for the Newtonian and non-Newtonian fluids and comparison of several physiologically relevant computational geometries in terms of wall displacement and wall shear stress. Numerical analysis and extensive experimental study for several standard geometries confirm reliability and accuracy of the proposed kinematic splitting scheme in order to approximate fluid-structure interaction problems.

Dynamic interactions between solids and viscous liquids with free surface = Interacciones dinámicas entre sólidos y líquidos viscosos con superficie libre

En esta tesis se investiga la interacción entre un fluido viscoso y un cuerpo sólido en presencia de una superficie libre. El problema se expresa teóricamente poniendo especial atención a los aspectos de conservación de energía y de la interacción del fluido con el cuerpo. El problema se considera 2D y monofásico, y un desarrollo matemático permite una descomposición de los términos disipativos en términos relacionados con la superficie libre y términos relacionados con la enstrofía. El modelo numérico utilizado en la tesis se basa en el método sin malla Smoothed Particle Hydrodynamics (SPH). De manera análoga a lo que se hace a nivel continuo, las propiedades de conservación se estudian en la tesis con el sistema discreto de partículas. Se tratan también las condiciones de contorno de un cuerpo que se mueve en un flujo viscoso, implementadas con el método ghost-fluid. Se ha desarrollado un algoritmo explícito de interacción fluido / cuerpo. Se han documentado algunos casos de modo detallado con el objetivo de comprobar la capacidad del modelo para reproducir correctamente la disipación de energía y el movimiento del cuerpo. En particular se ha investigado la atenuación de una onda estacionaria, comparando la simulación numérica con predicciones teóricas. Se han realizado otras pruebas para monitorizar la disipación de energía para flujos más violentos que implican la fragmentación de la superficie libre. La cantidad de energía disipada con los diferentes términos se ha evaluado en los casos estudiados con el modelo numérico. Se han realizado otras pruebas numéricas para verificar la técnica de modelización de la interacción fluido / cuerpo, concretamente las fuerzas ejercidas por las olas en cuerpos con formas simples, y el equilibrio de un cuerpo flotante con una forma compleja. Una vez que el modelo numérico ha sido validado, se han realizado simulaciones numéricas para obtener una comprensión más completa de la física implicada en casos (casi) realistas sobre los había aspectos que no se conocían suficientemente. En primer lugar se ha estudiado el el flujo alrededor de un cilindro bajo la superficie libre. El estudio se ha realizado con un número de Reynolds moderado, para un rango de inmersiones del cilindro y números de Froude. La solución numérica permite una investigación de los patrones complejos que se producen. La estela del cilindro interactúa con la superficie libre. Se han identificado algunos inestabilidades características. El segundo estudio se ha realizado sobre el problema de sloshing, tanto experimentalmente como numéricamente. El análisis se restringe a aguas poco profundas y con oscilación horizontal, pero se ha estudiado un gran número de condiciones, lo que lleva a una comprensión bastante completa de los sistemas de onda involucradas. La última parte de la tesis trata también sobre un problema de sloshing pero esta vez el tanque está oscilando con rotación y hay acoplamiento con un sistema mecánico. El sistema se llama pendulum-TLD (Tuned Liquid Damper - con líquido amortiguador). Este tipo de sistema se utiliza normalmente para la amortiguación de las estructuras civiles. El análisis se ha realizado analíticamente, numéricamente y experimentalmente utilizando líquidos con viscosidades diferentes, centrándose en características no lineales y mecanismos de disipación. ABSTRA C T The subject of the present thesis is the interaction between a viscous fluid and a solid body in the presence of a free surface. The problem is expressed first theoretically with a particular focus on the energy conservation and the fluid-body interaction. The problem is considered 2D and monophasic, and some mathematical development allows for a decomposition of the energy dissipation into terms related to the Free Surface and others related to the enstrophy. The numerical model used on the thesis is based on Smoothed Particle Hydrodynamics (SPH): a computational method that works by dividing the fluid into particles. Analogously to what is done at continuum level, the conservation properties are studied on the discrete system of particles. Additionally the boundary conditions for a moving body in a viscous flow are treated and discussed using the ghost-fluid method. An explicit algorithm for handling fluid-body coupling is also developed. Following these theoretical developments on the numerical model, some test cases are devised in order to test the ability of the model to correctly reproduce the energy dissipation and the motion of the body. The attenuation of a standing wave is used to compare what is numerically simulated to what is theoretically predicted. Further tests are done in order to monitor the energy dissipation in case of more violent flows involving the fragmentation of the free-surface. The amount of energy dissipated with the different terms is assessed with the numerical model. Other numerical tests are performed in order to test the fluid/body interaction method: forces exerted by waves on simple shapes, and equilibrium of a floating body with a complex shape. Once the numerical model has been validated, numerical tests are performed in order to get a more complete understanding of the physics involved in (almost) realistic cases. First a study is performed on the flow passing a cylinder under the free surface. The study is performed at moderate Reynolds numbers, for various cylinder submergences, and various Froude numbers. The capacity of the numerical solver allows for an investigation of the complex patterns which occur. The wake from the cylinder interacts with the free surface, and some characteristical flow mechanisms are identified. The second study is done on the sloshing problem, both experimentally and numerically. The analysis is restrained to shallow water and horizontal excitation, but a large number of conditions are studied, leading to quite a complete understanding of the wave systems involved. The last part of the thesis still involves a sloshing problem but this time the tank is rolling and there is coupling with a mechanical system. The system is named pendulum-TLD (Tuned Liquid Damper). This kind of system is normally used for damping of civil structures. The analysis is then performed analytically, numerically and experimentally for using liquids with different viscosities, focusing on non-linear features and dissipation mechanisms.

A method for computing nonsteady, incompressible, viscous fluid flows /

"Contract W-7405-ENG. 36 with the U.S. Atomic Energy Commission."

Capillary break-up rheometry of low-viscosity elastic fluids

We investigate the dynamics of the capillary thinning and break-up process for low viscosity elastic fluids such as dilute polymer solutions. Standard measurements of the evolution of the midpoint diameter of the necking fluid filament are augmented by high speed digital video images of the break up dynamics. We show that the successful operation of a capillary thinning device is governed by three important time scales (which characterize the relative importance of inertial, viscous and elastic processes), and also by two important length scales (which specify the initial sample size and the total stretch imposed on the sample). By optimizing the ranges of these geometric parameters, we are able to measure characteristic time scales for tensile stress growth as small as 1 millisecond for a number of model dilute and semi-dilute solutions of polyethylene oxide (PEO) in water and glycerol. If the final aspect ratio of the sample is too small, or the total axial stretch is too great, measurements are limited, respectively, by inertial oscillations of the liquid bridge or by the development of the well-known beads-on-a-string morphology which disrupt the formation of a uniform necking filament. By considering the magnitudes of the natural time scales associated with viscous flow, elastic stress growth and inertial oscillations it is possible to construct an operability diagram characterizing successful operation of a capillary break-up extensional rheometer. For Newtonian fluids, viscosities greater than approximately 70 mPas are required; however for dilute solutions of high molecular weight polymer, the minimum Viscosity is substantially lower due to the additional elastic stresses arising from molecular extension. For PEO of molecular weight 2.10(6) g/mol, it is possible to measure relaxation times of order 1 ms in dilute polymer solutions with zero-shear-rate viscosities on the order of 2-10 mPas.

Two-dimensional Stokes flow driven by elliptical paddles

A fast and accurate numerical technique is developed for solving the biharmonic equation in a multiply connected domain, in two dimensions. We apply the technique to the computation of slow viscous flow (Stokes flow) driven by multiple stirring rods. Previously, the technique has been restricted to stirring rods of circular cross section; we show here how the prior method fails for noncircular rods and how it may be adapted to accommodate general rod cross sections, provided only that for each there exists a conformal mapping to a circle. Corresponding simulations of the flow are described, and their stirring properties and energy requirements are discussed briefly. In particular the method allows an accurate calculation of the flow when flat paddles are used to stir a fluid chaotically.

The hemodynamic effects of in-tandem carotid artery stenosis: Implications for carotid endarterectomy

Objectives: It remains controversial whether patients with severe disease of the internal carotid artery and a coexisting stenotic lesion downstream would benefit from a carotid endarterectomy (CEA) of the proximal lesion. The aim of this study was to simulate the hemodynamic and wall shear effects of in-tandem internal carotid artery stenosis using a computational fluid dynamic (CFD) idealized model to give insight into the possible consequences of CEA on these lesions. Methods: A CFD model of steady viscous flow in a rigid tube with two asymmetric stenoses was introduced to simulate blood flow in arteries with multiple constrictions. The effect of varying the distance between the two stenoses, and the severity of the upstream stenosis on the pressure and wall shear stress (WSS) distributions on the second plaque, was investigated. The influence of the relative positions of the two stenoses was also assessed. Results: The distance between the plaques was found to have minimal influence on the overall hemodynamic effect except for the presence of a zone of low WSS (range -20 to 30 dyne/cm2) adjacent to both lesions when the two stenoses were sufficiently close (<4 times the arterial diameter). The upstream stenosis was protective if it was larger than the downstream stenosis. The relative positions of the stenoses were found to influence the WSS but not the pressure distribution. Conclusions: The geometry and positions of the lesions need to be considered when considering the hemodynamic effects of an in-tandem stenosis. Low WSS is thought to cause endothelial dysfunction and initiate atheroma formation. The fact that there was a flow recirculation zone with low WSS in between the two stenoses may demonstrate how two closely positioned plaques may merge into one larger lesion. Decision making for CEA may need to take into account the hemodynamic situation when an in-tandem stenosis is found. CFD may aid in the risk stratification of patients with this problem.

Theoretical study of a 16-µm CO2 downstream-mixing gasdynamic laser: A two-dimensional approach

A 16-µm CO2-N2 downstream-mixing gasdynamic laser, where a cold CO2 stream is mixed with a vibrationally excited N2 stream at the exit of the nozzle, is studied theoretically. The flow field is analyzed using a two-dimensional, unsteady, laminar and viscous flow model including appropriate finite-rate vibrational kinetic equations. The analysis showed that local small-signal gain up to 21.75 m−1 can be obtained for a N2 reservoir temperature of 2000 K and a velocity ratio of 1:1 between the CO2 and N2 mixing streams. Applied Physics Letters is copyrighted by The American Institute of Physics.

Transition from radiatively inefficient to cooling dominated phase in two temperature accretion disks around black holes

We investigate the transition of a radiatively inefficient phase of a viscous two temperature accreting flow to a cooling dominated phase and vice versa around black holes. Based on a global sub-Keplerian accretion disk model in steady state, including explicit cooling processes self-consistently, we show that general advective accretion flow passes through various phases during its infall towards a black hole. Bremsstrahlung, synchrotron and inverse Comptonization of soft photons are considered as possible cooling mechanisms. Hence the flow governs a much lower electron temperature similar to 10(8) - 10(9.5) K compared to the hot protons of temperature similar to 10(10.2) - 10(11.8) K in the range of the accretion rate in Eddington units 0.01 less than or simiar to (M) over dot less than or similar to 100. Therefore, the solutions may potentially explain the hard X-rays and the gamma-rays emitted from AGNs and X-ray binaries. We finally compare the solutions for two different regimes of viscosity and conclude that a weakly viscous flow is expected to be cooling dominated compared to its highly viscous counterpart which is radiatively inefficient. The flow is successfully able to reproduce the observed minosities of the under-fed AGNs and quasars (e.g. Sgr A*), ultra-luminous X-ray sources (e.g. SS433), as well as the highly luminous AGNs and ultra-luminous quasars (e.g. PKS 0743-67) at different combinations of the mass accretion rate and ratio of specific heats.

On the ignition and extinction problems in forced convection systems

In this paper the problem of ignition and extinction has been formulated for the flow of a compressible fluid with Prandtl and Schmidt numbers taken as unity. In particular, the problems of (i) a jet impinging on a wall of combustible material and (ii) the opposed jet diffusion flame have been studied. In the wall jet case, three approximations in the momentum equation namely, (i) potential flow, (ii) viscous flow, (ii) viscous incompressible with k = 1 and (iii) Lees' approximation (taking pressure gradient terms zero) are studied. It is shown that the predictions of the mass flow rates at extinction are not very sensitive to the approximations made in the momentum equation. The effects of varying the wall temperature in the case (i) and the jet temperature in the case (ii) on the extinction speeds have been studied. The effects of varying the activation energy and the free stream oxidant concentration in case (ii), have also been investigated.

Effect of surfactants on drop breakage in turbulent liquid dispersions

The existingm odels of drop breakage in stirred dispersions grossly overpredict the maximum drop size when surface active agents are present inspite of using the lowered value of interfacial tension. It is shown that the difference in the values of dynamic and static interfacial tension, aids the turbulent stresses in drop breakage. When the difference is zero, e.g. for pure liquids and for high concentration of surfactants, the influence of the addition of surfactant is merely to reduce the interfacial tension and can be accounted for by existingm odels. A modified model has been developed, where the drop breakage is assumed to be represented by a Voigt element. The deforming stresses are due to turbulence and the difference between dynamic and static interfacial tensions. The resisting stresses arise due to interfacial tension and the viscous flow inside the drop. The model yields the existing expressions for dmax as special cases. The model has been found to be satisfactory when tested against experimental results using the styrene-water-teepol system.

Power extraction studies in a gasdynamic laser

A generalized two‐dimensional flow‐radiation coupled model to extract power from a gasdynamic laser is proposed. The model is used for the study of power extraction from a 9.4‐μm CO2 downstream‐mixing gasdynamic laser, where a cold CO2+H2 stream is mixed with a vibrationally excited N2 stream at the nozzle exits. This model is developed by coupling radiation with the two‐dimensional, unsteady, laminar and viscous flow modeling needed for such systems. The analysis showed that the steady‐state value of 9.4‐μm intensity as high as 5×107 W/m2 can be obtained from the system studied. The role of H2 relaxant in the power extraction process has also been investigated.