Modelagem do fluxo sanguíneo na aorta abdominal utilizando interação fluido-estrutura

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Desenvolvimento de uma metodologia numérica para escoamentos viscoelásticos não-isotérmicos

Pós-graduação em Matematica Aplicada e Computacional - FCT

Solução das equações de Navier-Stokes e da equação de transporte por um método de elementos finitos estabilizado pela técnica de separação baseada na característica

Pós-graduação em Engenharia Mecânica - FEIS

Analytical solution for transient one-dimensional couette flow considering constant and time-dependent pressure gradients

This paperaims to determine the velocity profile, in transient state, for a parallel incompressible flow known as Couette flow. The Navier-Stokes equations were applied upon this flow. Analytical solutions, based in Fourier series and integral transforms, were obtained for the one-dimensional transient Couette flow, taking into account constant and time-dependent pressure gradients acting on the fluid since the same instant when the plate starts it´s movement. Taking advantage of the orthogonality and superposition properties solutions were foundfor both considered cases. Considering a time-dependent pressure gradient, it was found a general solution for the Couette flow for a particular time function. It was found that the solution for a time-dependent pressure gradient includes the solutions for a zero pressure gradient and for a constant pressure gradient.

Order 10 4 speedup in global linear instability analysis using matrix formation

A unified solution framework is presented for one-, two- or three-dimensional complex non-symmetric eigenvalue problems, respectively governing linear modal instability of incompressible fluid flows in rectangular domains having two, one or no homogeneous spatial directions. The solution algorithm is based on subspace iteration in which the spatial discretization matrix is formed, stored and inverted serially. Results delivered by spectral collocation based on the Chebyshev-Gauss-Lobatto (CGL) points and a suite of high-order finite-difference methods comprising the previously employed for this type of work Dispersion-Relation-Preserving (DRP) and Padé finite-difference schemes, as well as the Summationby- parts (SBP) and the new high-order finite-difference scheme of order q (FD-q) have been compared from the point of view of accuracy and efficiency in standard validation cases of temporal local and BiGlobal linear instability. The FD-q method has been found to significantly outperform all other finite difference schemes in solving classic linear local, BiGlobal, and TriGlobal eigenvalue problems, as regards both memory and CPU time requirements. Results shown in the present study disprove the paradigm that spectral methods are superior to finite difference methods in terms of computational cost, at equal accuracy, FD-q spatial discretization delivering a speedup of ð (10 4). Consequently, accurate solutions of the three-dimensional (TriGlobal) eigenvalue problems may be solved on typical desktop computers with modest computational effort.

Modelling wake effects in large wind farms in complex terrain: the problem, the methods and the issues

Computational fluid dynamic (CFD) methods are used in this paper to predict the power production from entire wind farms in complex terrain and to shed some light into the wake flow patterns. Two full three-dimensional Navier–Stokes solvers for incompressible fluid flow, employing k − ϵ and k − ω turbulence closures, are used. The wind turbines are modeled as momentum absorbers by means of their thrust coefficient through the actuator disk approach. Alternative methods for estimating the reference wind speed in the calculation of the thrust are tested. The work presented in this paper is part of the work being undertaken within the UpWind Integrated Project that aims to develop the design tools for next generation of large wind turbines. In this part of UpWind, the performance of wind farm and wake models is being examined in complex terrain environment where there are few pre-existing relevant measurements. The focus of the work being carried out is to evaluate the performance of CFD models in large wind farm applications in complex terrain and to examine the development of the wakes in a complex terrain environment.

Influencia de un fluido en las características dinámicas de placas circulares y casicirculares

La influencia de un fluido en las características dinámicas de estructuras se ha estudiado desde hace tiempo. Sin embargo muchos estudios se refieren a aplicaciones bajo el agua, como es el caso del sonar de un submarino por lo que el fluido circundante se considera líquido (sin efectos de compresibilidad). Más recientemente en aplicaciones acústicas y espaciales tales como antenas o paneles muy ligeros, ha sido estudiada la influencia en las características dinámicas de una estructura rodeada por un fluido de baja densidad. Por ejemplo se ha mostrado que el efecto del aire en el transmisor-reflector del Intelsat VI C-B con un diámetro de 3,2 metros y con un peso de sólo 34,7 kg disminuye la primera frecuencia en torno a un 20% con respecto a su valor en vacío. Por tanto es importante en el desarrollo de estas grandes y ligeras estructuras disponer de un método con el que estimar el efecto del fluido circundante sobre las frecuencias naturales de éstas. De esta manera se puede evitar el ensayo de la estructura en una cámara de vacío que para el caso de una gran antena o panel puede ser difícil y costoso. Se ha desarrollado un método de elementos de contorno (BEM) para la determinación del efecto del fluido en las características dinámicas de una placa circular. Una vez calculados analíticamente los modos de vibración de la placa en vacío, la matriz de masa añadida debido a la carga del fluido se determina por el método de elementos de contorno. Este método utiliza anillos circulares de manera que el número de elementos para obtener unos resultados precisos es muy bajo. Se utiliza un procedimiento de iteración para el cálculo de las frecuencias naturales del acoplamiento fluido-estructura para el caso de fluido compresible. Los resultados del método se comparan con datos experimentales y otros modelos teóricos mostrando la precisión y exactitud para distintas condiciones de contorno de la placa. Por otro lado, a veces la geometría de la placa no es circular sino casi-circular y se ha desarrollado un método de perturbaciones para determinar la influencia de un fluido incompresible en las características dinámicas de placas casi-circulares. El método se aplica a placas con forma elíptica y pequeña excentricidad. Por una parte se obtienen las frecuencias naturales y los modos de deformación de la placa vibrando en vacío. A continuación, se calculan los coeficientes adimensionales de masa virtual añadida (factores NAVMI). Se presentan los resultados de estos factores y el efecto del fluido en las frecuencias naturales. ABSTRACT The influence of the surrounding fluid on the dynamic characteristics of structures has been well known for many years. However most of these works were more concerned with underwater applications, such as the sonar of a submarine and therefore the surrounding fluid was considered a liquid (negligible compressibility effects). Recently for acoustical and spatial applications such as antennas or very light panels the influence on the dynamic characteristics of a structure surrounded by a fluid of low density has been studied. Thus it has been shown that the air effect for the Intelsat VI C-B transmit reflector with a diameter of 3,2 meters and weighting only 34,7 kg decreases the first modal frequency by 20% with respect to the value in vacuum. It is important then, in the development of these light and large structures to have a method that estimates the effect that the surrounding fluid will have on the natural frequencies of the structure. In this way it can be avoided to test the structure in a vacuum chamber which for a large antenna or panel can be difficult and expensive A BEM method for the determination of the effect of the surrounding fluid on the dynamic characteristics of a circular plate has been developed. After the modes of the plate in vacuum are calculated in an analytical form, the added mass matrix due to the fluid loading is determined by a boundary element method. This method uses circular rings so the number of elements to obtain an accurate result is very low. An iteration procedure for the computation of the natural frequencies of the couple fluid-structure system is presented for the case of the compressibility effect of air. Comparisons of the present method with various experimental data and other theories show the efficiency and accuracy of the method for any support condition of the plate. On the other hand, sometimes the geometry of the plate is not circular but almost-circular, so a perturbation method is developed to determine the influence of an incompressible fluid on the dynamic characteristics of almost-circular plates. The method is applied to plates of elliptical shape with low eccentricity. First, the natural frequencies and the mode shapes of the plate vibrating in vacuum are obtained. Next, the nondimensional added virtual mass coefficients (NAVMI factors) are calculated. Results of this factors and the effect of the fluid on the natural frequencies are presented.

The forced oscillations of submerged bodies

The increasing use of very light structures in aerospace applications are given rise to the need of taking into account the effects of the surrounding media in the motion of a structure (as for instance, in modal testing of solar panels or antennae) as it is usually performed in the motion of bodies submerged in water in marine applications. New methods are in development aiming at to determine rigid-body properties (the center of mass position and inertia properties) from the results of oscillations tests (at low frequencies during modal testing, by exciting the rigid-body modes only) by using the equations of the rigid-body dynamics. As it is shown in this paper, the effect of the surrounding media significantly modifies the oscillation dynamics in the case of light structures and therefore this effect should be taken into account in the development of the above-mentioned methods. The aim of the paper is to show that, if a central point exists for the aerodynamic forces acting on the body, the motion equations for the small amplitude rotational and translational oscillations can be expressed in a form which is a generalization of the motion equations for a body in vacuum, thus allowing to obtain a physical idea of the motion and aerodynamic effects and also significantly simplifying the calculation of the solutions and the interpretation of the results. In the formulation developed here the translational oscillations and the rotational motion around the center of mass are decoupled, as is the case for the rigid-body motion in vacuum, whereas in the classical added mass formulation the six motion equations are coupled. Also in this paper the nonsteady motion of small amplitude of a rigid body submerged in an ideal, incompressible fluid is considered in order to define the conditions for the existence of the central point in the case of a three-dimensional body. The results here presented are also of interest in marine applications.

Reconstruction of an unsteady flow from incomplete boundary data:Part I. Theory

Purpose – To propose and investigate a stable numerical procedure for the reconstruction of the velocity of a viscous incompressible fluid flow in linear hydrodynamics from knowledge of the velocity and fluid stress force given on a part of the boundary of a bounded domain. Design/methodology/approach – Earlier works have involved the similar problem but for stationary case (time-independent fluid flow). Extending these ideas a procedure is proposed and investigated also for the time-dependent case. Findings – The paper finds a novel variation method for the Cauchy problem. It proves convergence and also proposes a new boundary element method. Research limitations/implications – The fluid flow domain is limited to annular domains; this restriction can be removed undertaking analyses in appropriate weighted spaces to incorporate singularities that can occur on general bounded domains. Future work involves numerical investigations and also to consider Oseen type flow. A challenging problem is to consider non-linear Navier-Stokes equation. Practical implications – Fluid flow problems where data are known only on a part of the boundary occur in a range of engineering situations such as colloidal suspension and swimming of microorganisms. For example, the solution domain can be the region between to spheres where only the outer sphere is accessible for measurements. Originality/value – A novel variational method for the Cauchy problem is proposed which preserves the unsteady Stokes operator, convergence is proved and using recent for the fundamental solution for unsteady Stokes system, a new boundary element method for this system is also proposed.

Modélisation d'éléments finis pour les interactions fluide-structure

Cette thèse concerne la modélisation des interactions fluide-structure et les méthodes numériques qui s’y rattachent. De ce fait, la thèse est divisée en deux parties. La première partie concerne l’étude des interactions fluide-structure par la méthode des domaines fictifs. Dans cette contribution, le fluide est incompressible et laminaire et la structure est considérée rigide, qu’elle soit immobile ou en mouvement. Les outils que nous avons développés comportent la mise en oeuvre d’un algorithme fiable de résolution qui intégrera les deux domaines (fluide et solide) dans une formulation mixte. L’algorithme est basé sur des techniques de raffinement local adaptatif des maillages utilisés permettant de mieux séparer les éléments du milieu fluide de ceux du solide que ce soit en 2D ou en 3D. La seconde partie est l’étude des interactions mécaniques entre une structure flexible et un fluide incompressible. Dans cette contribution, nous proposons et analysons des méthodes numériques partitionnées pour la simulation de phénomènes d’interaction fluide-structure (IFS). Nous avons adopté à cet effet, la méthode dite «arbitrary Lagrangian-Eulerian» (ALE). La résolution fluide est effectuée itérativement à l’aide d’un schéma de type projection et la structure est modélisée par des modèles hyper élastiques en grandes déformations. Nous avons développé de nouvelles méthodes de mouvement de maillages pour aboutir à de grandes déformations de la structure. Enfin, une stratégie de complexification du problème d’IFS a été définie. La modélisation de la turbulence et des écoulements à surfaces libres ont été introduites et couplées à la résolution des équations de Navier-Stokes. Différentes simulations numériques sont présentées pour illustrer l’efficacité et la robustesse de l’algorithme. Les résultats numériques présentés attestent de la validité et l’efficacité des méthodes numériques développées.

Unsteady incompressible boundary layer flow of a micropolar fluid at a stagnation point

The flow, heat and mass transfer on the unsteady laminar incompressible boundary layer in micropolar fluid at the stagnation point of a 2-dimensional and an axisymmetric body have been studied when the free stream velocity and the wall temperature vary arbitrarily with time. The partial defferential equations governing the flow have been solved numerically using a quasilinear finite-difference scheme. The skin friction, microrotation gradient and heat transfer parameters are found to be strongly dependent on the coupling parameter, mass transfer and time, whereas the effect of the microrotation parameter on the skin friction and heat transfer is rather weak, but microrotation gradient is strongly affected by it. The Prandtl number and the variation of the wall temperature with time affect the heat-transfer very significantly but the skin friction and micrortation gradient are unaffected by them.

Incompressible micropolar fluid flow over a semi-infinite plate

Micropolar fluid flow over a semi-infinite flat plate has been described by using the parabolic co-ordinates and the method of series truncation in order to study the flow for low to large Reynolds numbers. These co-ordinates permit to study the flow regime at the leading edge. Numerical results have been presented for different Reynolds numbers. Results show a reduction in skin friction.

Intrinsic equations of steady rotating, incompressible viscous fluid flow

The equations governing the flow of a steady rotating incompressible viscous fluid are expressed in intrinsic form along the vortex lines and their normals. Using these equations the effects of rotation on the geometric properties of viscous fluid flows are studied. A particular flow in which the vortex lines are right circular helices is discussed.

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

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

Two-dimensional bow and stern flows in a finite-depth fluid with constant vorticity

This thesis is concerned with two-dimensional free surface flows past semi-infinite surface-piercing bodies in a fluid of finite-depth. Throughout the study, it is assumed that the fluid in question is incompressible, and that the effects of viscosity and surface tension are negligible. The problems considered are physically important, since they can be used to model the flow of water near the bow or stern of a wide, blunt ship. Alternatively, the solutions can be interpreted as describing the flow into, or out of, a horizontal slot. In the past, all research conducted on this topic has been dedicated to the situation where the flow is irrotational. The results from such studies are extended here, by allowing the fluid to have constant vorticity throughout the flow domain. In addition, new results for irrotational flow are also presented. When studying the flow of a fluid past a surface-piercing body, it is important to stipulate in advance the nature of the free surface as it intersects the body. Three different possibilities are considered in this thesis. In the first of these possibilities, it is assumed that the free surface rises up and meets the body at a stagnation point. For this configuration, the nonlinear problem is solved numerically with the use of a boundary integral method in the physical plane. Here the semi-infinite body is assumed to be rectangular in shape, with a rounded corner. Supercritical solutions which satisfy the radiation condition are found for various values of the Froude number and the dimensionless vorticity. Subcritical solutions are also found; however these solutions violate the radiation condition and are characterised by a train of waves upstream. In the limit that the height of the body above the horizontal bottom vanishes, the flow approaches that due to a submerged line sink in a $90^\circ$ corner. This limiting problem is also examined as a special case. The second configuration considered in this thesis involves the free surface attaching smoothly to the front face of the rectangular shaped body. For this configuration, nonlinear solutions are computed using a similar numerical scheme to that used in the stagnant attachment case. It is found that these solution exist for all supercritical Froude numbers. The related problem of the cusp-like flow due to a submerged sink in a corner is also considered. Finally, the flow of a fluid emerging from beneath a semi-infinite flat plate is examined. Here the free surface is assumed to detach from the trailing edge of the plate horizontally. A linear problem is formulated under the assumption that the elevation of the plate is close to the undisturbed free surface level. This problem is solved exactly using the Wiener-Hopf technique, and subcritical solutions are found which are characterised by a train of sinusoidal waves in the far field. The nonlinear problem is also considered. Exact relations between certain parameters for supercritical flow are derived using conservation of mass and momentum arguments, and these are confirmed numerically. Nonlinear subcritical solutions are computed, and the results are compared to those predicted by the linear theory.