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.

Instabilities of a liquid film flowing down an inclined porous plane

The problem of a film flowing down an inclined porous layer is considered. The fully developed basic flow is driven by gravitation. A careful linear instability analysis is carried out. We use Darcy's law to describe the porous layer and solve the coupling equations of the fluid and the porous medium rather than the decoupled equations of the one-sided model used in previous works. The eigenvalue problem is solved by means of a Chebyshev collocation method. We compare the instability of the two-sided model with the results of the one-sided model. The result reveals a porous mode instability which is completely neglected in previous works. For a falling film on an inclined porous plane there are three instability modes, i.e., the surface mode, the shear mode, and the porous mode. We also study the influences of the depth ratio d, the Darcy number delta, and the Beavers-Joseph coefficient alpha(BJ) on the instability of the system.

Design support for turbomachinery

From the steam turbines which provide most of our electricity to the jet engines which have shrunk our World, turbomachines undoubtedly play a major role in life today. Competition in the turbomachinery industry is fiercely strong [Wisler, 1998], hence good aerodynamic design is vital. However, with efficiency levels already close to their theoretical maxima, companies are increasingly looking to reduce costs and increase reliability through improved design practice. Computational Fluid Dynamics (CFD) can make a strong contribution to assisting this process as it has the potential to increase performance while reducing cost. The situation is, however, complicated by an ever decreasing number of engineers with sufficient design experience to reap the full benefits offered by CFD. With the large risks involved, novice designers of today are increasingly confined to refining old designs rather than gaining experience, like their forebears, through 'clean sheet' exercises. Hence it is desirable to capture the knowledge and experience of older designers, before it is lost, to assist the engineers of tomorrow. It is therefore the aim of this project to produce a design support tool which will not only store the appropriate CFD codes, but also provide a dynamic signpost (based on elicited knowledge and experience) to advise the engineer in their use. The signposting methodology developed for the aerospace industry [Clarkson and Hamilton, 1997] will provide the basic framework for the tool. This paper reviews current turbomachinery design practice (including an examination of the relevant CFD) in order to establish the important issues which a support tool must address. Current design support methodologies and their propriety are then reviewed, followed by a detailed description of the signposting concept. It then sets out a clear statement of the objectives for the research and the methods proposed to meet them. The paper concludes with a timetable of the work.

Estimativa do erro de discretização analítico na solução de equações diferenciais utilizando o Método de Volumes Finitos

As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado

Effects of density differences on lateral mixing in open-channel flows

This study investigates lateral mixing of tracer fluids in turbulent open-channel flows when the tracer and ambient fluids have different densities. Longitudinal dispersion in flows with longitudinal density gradients is investigated also.

Lateral mixing was studied in a laboratory flume by introducing fluid tracers at the ambient flow velocity continuously and uniformly across a fraction of the flume width and over the entire depth of the ambient flow. Fluid samples were taken to obtain concentration distributions in cross-sections at various distances, x, downstream from the tracer source. The data were used to calculate variances of the lateral distributions of the depth-averaged concentration. When there was a difference in density between the tracer and the ambient fluids, lateral mixing close to the source was enhanced by density-induced secondary flows; however, far downstream where the density gradients were small, lateral mixing rates were independent of the initial density difference. A dimensional analysis of the problem and the data show that the normalized variance is a function of only three dimensionless numbers, which represent: (1) the x-coordinate, (2) the source width, and (3) the buoyancy flux from the source.

A simplified set of equations of motion for a fluid with a horizontal density gradient was integrated to give an expression for the density-induced velocity distribution. The dispersion coefficient due to this velocity distribution was also obtained. Using this dispersion coefficient in an analysis for predicting lateral mixing rates in the experiments of this investigation gave only qualitative agreement with the data. However, predicted longitudinal salinity distributions in an idealized laboratory estuary agree well with published data.