Determinação da força de arrasto e da velocidade terminal de frutos de café pela técnica de elementos finitos

Determination of drag force and terminal velocity of coffee fruits using finite element technique

ARMAZENAMENTO E PROCESSAMENTO DE PRODUTOS AGRÍCOLAS

Determinação da força de arrasto e da velocidade terminal de frutos de café pela técnica de elementos finitos

Determination of drag force and terminal velocity of coffee fruits using finite element technique

Sandra M. Couto^I; Ana E. O. dos Santos^II; Stella M. J. Vieira^III; Danilo J. P. da Silva^IV

^IDEA/UFV. CEP 36571-000, Viçosa, MG. Fone: (31) 3899-1927. E-mail: [e-mail]

^IIUFV. E-mail: [e-mail]

^IIIUFV. E-mail: [e-mail]

^IVUFV. E-mail: [e-mail]

RESUMO

Este trabalho apresenta um método alternativo para a determinação de velocidades terminais e de forças de arrasto em frutos de café, ao se variar o teor de umidade do produto e sua configuração geométrica (esférica e elíptica). Especificamente, investigou-se a aplicabilidade da técnica de elementos finitos para a obtenção de: (a) forças de arrasto no produto, usando-se valores experimentais das dimensões do fruto e de sua velocidade terminal e (b) velocidades terminais do produto, empregando-se valores experimentais das dimensões do fruto e de seu peso. A técnica dos elementos finitos mostrou-se ferramenta eficiente para simular velocidades terminais e forças de arrasto em frutos de café submetidos a um fluxo de ar. Independente do teor de umidade dos frutos, os valores simulados para a velocidade terminal e para força de arrasto no produto considerado elíptico são, estatisticamente, iguais àqueles determinados experimentalmente. Os erros relativos médios, envolvidos no valores simulados para a velocidade terminal do produto, em cada um dos teores de umidade, não foram superiores a 7,5 e 13,6% para o produto considerado elíptico e esférico, respectivamente.

Palavras-chave: configuração geométrica, força de arrasto, Coffea arabica L.

ABSTRACT

This work presents an alternative method for terminal speed and drag force determinations for coffee fruits with different moisture contents and geometric configurations (spherical and elliptic). Specifically, the use of the finite element technique was investigated for obtaining: (a) product drag forces using experimental values of the fruit dimensions and its terminal velocity and (b) product terminal velocity using experimental values of fruit dimensions and weight. The use of finite element technique was shown to be an efficient tool to simulate terminal velocity and drag force on coffee fruits submitted to airflow. Independent of the fruit moisture content, simulated and experimental values for the terminal velocity and drag force did not differ statistically when the product was considered elliptic. For each product moisture content, the mean relative errors involved in the simulated values of terminal velocity were not superior to 7.5 and 13.6% for the fruits considered elliptic and spherical, respectively.

Key words: geometric configuration, drag force, Coffea arabica L.

INTRODUÇÃO

O conhecimento da velocidade terminal de um produto e da força de arrasto em partículas colocadas em fluidos em movimento é sempre necessário em muitas aplicações em engenharia, tais como transporte pneumático e hidráulico de um produto e processos de separação e classificação de materiais. Determinações experimentais da velocidade terminal de vários tipos de grãos e sementes foram realizadas no passado e, atualmente, ainda continuam sendo o objetivo de vários pesquisadores (por exemplo: Carman, 1996; Singh & Goswami, 1996; Suthar & Das, 1996; Allen & Watts, 1997; Gupta & Das, 1997; Nimkar & Chattopadhyay, 2001). Devido a importância da velocidade terminal e da força de arrasto, um grande número de pesquisas tem sido realizado objetivando o desenvolvimento de métodos de predição precisos para estimar estas grandezas; entretanto, a maioria dessas pesquisas é direcionada para materiais não agrícolas.

De acordo com Fox & McDonald (1998), a superfície de um corpo imerso em um fluido em movimento, fica submetida a uma força dinâmica, resultante de forças de pressão (normal à superfície) e de cisalhamento (tangente à superfície). Tensões cisalhantes devido à viscosidade e gradientes de velocidade na superfície do corpo causam forças tangenciais à superfície enquanto variações nas intensidades de pressão ao longo da superfície, devido aos efeitos dinâmicos, resultam em forças normais à superfície do corpo. A componente da força total exercida pelo fluido no sólido, na direção do escoamento é denominada arrasto.

A força de arrasto atuante em um corpo é, portanto, composta de duas parcelas: uma devida ao arrasto de atrito e a outra ao arrasto de pressão. A contribuição relativa de cada uma varia com o regime do escoamento e com as condições da superfície do corpo. O arrasto de pressão depende, basicamente, da forma do corpo e pode ser reduzido a um mínimo como, por exemplo, em um elipsóide de revolução exposto a um escoamento paralelo ao eixo de simetria do corpo. O arrasto de atrito não depende significativamente da forma do corpo; entretanto, é influenciado pelas características da superfície,Tabak & Wolf (1998), Dong et al.(2001).

Os arrastos de pressão e atrito são simultaneamente de grande valor apenas em poucas situações. Para objetos que não apresentam qualquer sustentação, o arrasto de pressão pode ser considerado como o arrasto total, Giles (1980). Segundo Fox & McDonald (1998), no caso de escoamento sobre uma esfera lisa, tanto o arrasto de atrito quanto o de pressão contribuem para o arrasto total. Para números de Reynolds, Re, muito baixos (Re < 1), o arrasto é predominantemente devido ao atrito. À medida que o número de Reynolds aumenta até cerca de 1000, o arrasto é uma combinação das parcelas devidas ao atrito e à pressão. A contribuição relativa do arrasto de atrito cai com o aumento do número de Reynolds; para Re » 1000, o arrasto de atrito é aproximadamente 5% do arrasto total.

A força de arrasto (F) em um corpo imerso em um fluido em movimento é comumente expressa como , em que C é o coeficiente de arrasto (adimensional); A é a área característica (m²); é a densidade do fluido (kg m^-3) e V é a velocidade uniforme do fluido (m s^-1). O coeficiente de arrasto inclui os arrastos de pressão e de atrito.

A determinação da velocidade terminal de corpos, baseando-se em expressões que relacionam o coeficiente de arrasto e o número de Reynolds, tem sido extensivamente investigada por vários pesquisadores nos últimos noventa anos. Nguyen et al. (1997) relatam que, atualmente, existem pelo menos 35 equações na literatura que expressam o coeficiente de arrasto (C) sobre uma partícula em função do número de Reynolds. A determinação da velocidade terminal através de quaisquer dessas equações propostas para "C versus Re", requer um procedimento de tentativa e erro, visto que a velocidade se encontra presente nas duas variáveis e a relação entre "C e Re" é altamente não linear.

A velocidade terminal e, conseqüentemente, o coeficiente de arrasto de uma partícula são por demais influenciados pela sua forma, tamanho e orientação, pela viscosidade do meio e pelas densidades da partícula e do fluido. Dos modelos semi-empíricos desenvolvidos para predizer forças de arrasto em corpos submersos em um fluido, grande parte considera as partículas tendo uma forma esférica, apesar das partículas envolvidas em processos de engenharia não terem, na grande maioria, esta forma. Atualmente, é possível se estimar forças de arrasto ou velocidades terminais de uma partícula esférica sob a maioria das condições de interesse prático em um meio newtoniano incompressível. Entretanto, quando se trata de partículas não esféricas, a falta de padrões para a medição da sua forma, seu tamanho e orientação durante a flutuação, de acordo com Chhabra et al. (1999), tem sido um dos grandes obstáculos para obtenção de modelos adequados para predizer forças de arrasto e no desenvolvimento de correlações aplicáveis universalmente.

Arsenijevic et al. (1999) desenvolveram um método para a determinação da velocidade terminal de partículas não esféricas, no qual usaram dados de queda de pressão em um leito recheado do produto para a determinação indireta do fator de forma para a partícula. Segundo os autores, os outros métodos existentes necessitam do conhecimento da forma da partícula, parâmetro de difícil obtenção para materiais reais.

Chhabra et al. (1999) investigaram a eficácia de correlações mais comumente usadas para estimar o coeficiente de arrasto de partículas não esféricas em fluidos viscosos incompressíveis e observaram que a influência da forma da partícula sobre o coeficiente de arrasto aumenta com o número de Reynolds e, ainda, que melhores estimativas desta grandeza poderiam ser obtidas somente se parâmetros ligados à forma da partícula fossem adicionados às correlações.

Teoricamente, a velocidade terminal de um partícula e a força de arrasto podem ser determinadas por meio de uma análise do escoamento de um fluido sobre uma superfície sólida (partícula). Para condições isotérmicas, isto envolve a resolução simultânea das equações de conservação de massa e de conservação do momentum. Na maioria dos casos, soluções analíticas para este tipo de problema não podem ser viabilizadas e se faz necessário a adoção de algum método numérico para a obtenção de soluções.

O método dos elementos finitos, segundo Moaveni (1999), é um procedimento numérico que pode ser usado para a obtenção de soluções para uma extensa série de problemas relacionados com engenharia, incluindo escoamento de fluidos. A vantagem deste método sobre os outros reside no fato de ser aplicável a objetos de formas irregulares, compostos de diferentes materiais e tendo condições de contorno mistas. Sendo um método numérico, as soluções, obtidas por esta técnica, aproximam as soluções exatas somente em pontos discretos do sistema, denominados nós. Para a aplicação do método o sistema é, inicialmente, discretizado em pequenas sub-regiões, denominadas elementos, conectadas através de pontos (nós). O método usa formulações integrais para criar um sistema de equações algébricas. Uma função contínua aproximada é usada para representar a solução para cada elemento. A solução completa é gerada pela conexão das soluções individuais.

Torna-se evidente, diante do exposto, a necessidade de métodos alternativos mais eficazes para a determinação de velocidades terminais e de forças de arrasto em partículas não esféricas, tais como a maioria dos produtos agrícolas. Neste sentido, o trabalho teve, por objetivo, o desenvolvimento de uma metodologia para a determinação de propriedades aerodinâmicas de frutos de café, ao se variar o teor de umidade do produto e sua configuração geométrica (esférica e elíptica). Especificamente, investigou-se a aplicabilidade da técnica de elementos finitos para a obtenção de: (a) forças de arrasto em frutos de café, usando-se valores experimentais das dimensões do fruto e de sua velocidade terminal e (b) velocidades terminais do produto através de valores experimentais das dimensões do fruto e de seu peso.

MATERIAL E MÉTODOS

Predição da força de arrasto nos frutos de café

Em geral, a superfície de um corpo imerso em um fluido, em cada ponto, faz um ângulo com a direção do escoamento e, neste caso, o arrasto total em um elemento de área da superfície, dA, é a soma das componentes, na direção do escoamento, da força cisalhante (t dA) e da força de pressão do fluido (PdA) como esquematizado na Figura 1. Considerando-se o elemento de área, dA, inclinado de um ângulo de (90º - á ) em relação à direção do escoamento, o arrasto devido ao cisalhamento (arrasto de atrito) é ô sen á dA e aquele devido à pressão (arrasto de forma ou de pressão) é igual a P cos ô dA. O arrasto total sobre um corpo é a soma das integrais dessas duas quantidades sobre toda a superfície do corpo em contato com o fluido.

A metodologia apresentada neste trabalho para a determinação de forças de arrasto em frutos de café considera:

- escoamento de ar em um duto circular (raio = 0,046 m e comprimento = 1,0 m) contendo, no seu interior, um fruto de café;

- fruto posicionado no eixo longitudinal central do duto e na metade do seu comprimento;

- velocidade do ar no duto equivalente à velocidade terminal do produto. Portanto, o fruto permanece estacionário no interior do duto, durante o escoamento de ar;

- sistema investigado corresponde à metade da seção longitudinal do duto, como esquematizado na Figura 2A, com simetria em relação ao eixo longitudinal central do duto (X), paralelo à direção do fluxo de ar;

- efeito tridimensional: obtido através da rotação do sistema em torno do eixo de simetria (X) durante as análises;

- forma geométrica do fruto de café no interior do duto: esférica ou elíptica.

- dimensões características do fruto: A e B, Figura 2A. No caso do fruto elíptico: A = a/2 e B = (c+b)/4, Figura 2B. Para fruto esférico: RE=A=B= (a+b+c)/4, em que RE é o raio equivalente esférico do fruto;

- desprezível a contribuição do arrasto de atrito para o arrasto total nos frutos, ou seja, nas análises considerou-se apenas o arrasto de pressão.

Para a determinação da força de arrasto no produto, independente da forma assumida para o fruto, na sua superfície (projeção de sua metade no plano XY) foram especificados trinta e sete pontos (defasados de cinco graus), como pode ser visualizado na Figura 3A.

Os valores de pressão total nos 37 pontos (nós) foram considerados previamente conhecidos (provenientes de análise do escoamento de ar no duto pela técnica de elementos finitos). Qualquer um dos pontos representados na Figura 3A, por exemplo: o denominado "N", pode ser localizado em relação à origem do sistema de referência, pelo ângulo "q" ou pelas suas coordenadas "X_N" e "Y_N".

Conforme mostrado na Figura 3B, em torno de cada um dos 37 nós foram construídos pontos intermediários, denominados Na (anterior ao nó N, com coordenadas X_Na e Y_Na) e Np (posterior ao nó N, com coordenadas X_Np e Y_Np), cada um defasado de 2,5 graus do nó N. O segmento de arco entre os dois pontos intermediários, vizinhos de "N", foi aproximado por um segmento de reta de comprimento dS_N, no qual a pressão correspondente ao nó "N" passou a atuar.

As equações usadas para a obtenção da força de arrasto nos frutos, considerados esferas e elipses, são apresentadas na Tabela 1, bem como as expressões e aproximações utilizadas para as grandezas envolvidas nas determinações. Os parâmetros principais presentes nas equações podem ser visualizados na Figura 3B e C.

em que:

H - número total de nós para os quais a pressão é conhecida

dS_N- magnitude do segmento de re

F - força de arrasto total sobre o fruto

F_N- força de arrasto sobre o elemento N

P_Nx - componente, na direção X, da pressão atuante no elemento N

R_E - raio equivalente da esfera

A_Ny- componente de área do elemento N, perpendicular a P_Nx

A_N- área da superfície do fruto atuada pela pressão P_N

r - dimensão que delimita a área de cada elemento (Figura x)

F_Ni - força de arrasto nos elementos sem contato com o eixo X

F_Ne- força de arrasto nos elementos em contato com o eixo X

As forças de arrasto sobre as áreas superficiais correspondentes aos nós pertencentes ao primeiro quadrante da esfera ou da elipse, foram consideradas negativas enquanto aquelas referentes aos nós do segundo quadrante, foram consideradas positivas. No caso do fruto elíptico, para nós da superfície em contato com o eixo X, os valores de X_Na e Y_Na, concernente ao nó da direita, e X_Np e Y_Np , correspondente ao nó da esquerda, foram considerados iguais a zero na equação correspondente na Tabela 1.

Escoamento de ar no duto - Análise por Elementos Finitos

Valores de pressão total nos trinta e sete nós, Figura 3(A), especificados previamente na superfície do fruto submetido a um fluxo de ar, foram obtidos através de uma análise de escoamento de fluido usando-se a técnica de elementos finitos. O problema de escoamento (isotérmico, incompressível e em regime permanente) de fluidos é definido pelas leis de conservação de massa, Eq. 1, e de momentum, eqs. 2, 3 e 4.

(1)

em que:

V_x, V_y e V_z= componentes da velocidade nas direções x, y e z, respectivamente.

Equações de momentum:

(2)

(3)

(4)

em que:

P - pressão do fluido

gx, gy e gz - componentes da aceleração da gravidade nas direções das coordenadas espaciais x, y e z, respectivamente

m_e- viscosidade efetiva

r - densidade do fluido

Para a análise do problema de escoamento usou-se um programa computacional: ANSYS (versão 5.3), especialista em soluções usando a técnica de elementos finitos. Este programa possui uma opção denominada FLOTRAN CFD (Computational Fluid Dynamics), dirigida para a resolução de problemas em escoamento de fluidos, utilizada neste trabalho. As equações diferenciais parciais que expressam as leis de conservação de massa e de momentum, são discretizadas com uma técnica baseada em elementos finitos. O processo de discretização consiste em derivar matrizes para os elementos, que serão acopladas em uma matriz global. O método de Garlekin dos resíduos ponderados é usado na obtenção das matrizes para os elementos. As matrizes dos elementos após serem formadas, são acopladas e o sistema de equações resultante é resolvido para cada grau de liberdade (velocidade e pressão), separadamente. O problema de escoamento é não linear e as equações governantes encontram-se acopladas.

FLOTRAN CFD disponibiliza dois tipos de método para a solução de sistemas de equações; o primeiro é um método rápido e aproximado: TDMA (Tri-Diagonal Matrix Algorithm) e o segundo abrange métodos semidiretos: métodos direcionais conjugados. Os dois tipos são iterativos. TDMA é usado para se obter a solução aproximada e os métodos direcionais conjugados quando soluções "exatas" são requeridas. Em certas situações, soluções exatas para as equações devem ser obtidas, enquanto em outras, soluções aproximadas são adequadas. Na solução do algoritmo seqüencial neste trabalho, o método TDMA (Tri-Diagonal Matrix Algorithm) foi usado para a solução das equações de momentum e o método direcional conjugado, gradiente conjugado precondicionado, foi usado na solução da equação de pressão.

Os passos básicos envolvidos na resolução do problema em que foi usada a técnica dos elementos finitos, foram: (1) discretização do sistema - divisão do domínio em elementos e nós; (2) consideração de uma função contínua aproximada para representar a solução de um elemento; (3) desenvolvimento de equações para um elemento; (4) conexão das soluções dos elementos para se obter a solução completa do problema (construção da matriz de rigidez global); (5) aplicação das condições iniciais e de contorno no sistema, bem como as forças que nele atuam; (6) resolução simultânea de um conjunto de equações algébricas para obter as soluções nos nós. Como resultado das análises, foram obtidas velocidades e pressões em cada um dos nós dos elementos que compunham o sistema.

Para a análise, o sistema (Figura 2A) foi dividido em elementos bidimensionais do tipo FLUID141, quadrilateral com quatro nós (pertencente à biblioteca do ANSYS). Este tipo de elemento permite regiões fluidas e sólidas no sistema. Elementos com propriedades relativas ao fluido foram usados na região onde ocorria o escoamento de ar e outros, com propriedades de sólido, foram usados para a região representando o fruto de café.

Os parâmetros de entrada para a execução do programa foram: dimensões do sistema relativas ao duto de escoamento; propriedades e características físicas do produto (dimensões relevantes à forma assumida pelo fruto; peso e densidade), propriedades do ar (densidade e viscosidade); velocidade terminal do produto e condições de contorno para o sistema.

Para se estimar a força de arrasto em frutos de café submetidos a um fluxo de ar, as seguintes propriedades e características do produto foram determinadas experimentalmente: as três dimensões principais do fruto (a, b e c, Figura 2B), massa e volume. As dimensões foram medidas por meio de um paquímetro (± 0,001 cm), a massa do fruto foi obtida em uma balança semi-analítica (10^-3 g) e o volume do produto foi determinado pelo método de complementação de volume seguindo a metodologia descrita por Moreira et al. (1984/85). A densidade do produto foi calculada pela razão entre a massa do fruto e seu volume.

As dimensões principais foram também empregadas para a obtenção da esfericidade, f, do produto:

Determinou-se experimentalmente a velocidade terminal dos frutos de café utilizando-se um dispositivo constituído de um ventilador centrífugo acoplado a um tubo de acrílico transparente (0,092 m de diâmetro e 2,0 m de comprimento), posicionado na vertical. Após ser acomodado no centro de um compartimento telado posicionado a 1,0 m da extremidade inferior do tubo de acrílico, o fluxo de café era submetido a um fluxo de ar, cuja intensidade era aumentada até que ocorresse o início de flutuação do produto. A velocidade do ar no local de flutuação do fruto, medida por meio de um anemômetro de fio quente, foi considerada como a velocidade terminal do produto. Obtiveram-se as velocidades terminais de frutos cujas dimensões, massas e volumes, foram previamente determinados.

Nos experimentos para a determinação das propriedades e características físicas do produto e sua velocidade terminal, foram usados frutos de café da variedade Catuaí, no estádio de maturação correspondente à coloração "cereja", com teores de umidade iguais a 67,5; 46,9 e 18,7% b.u. Para cada teor de umidade do produto foram investigados dez frutos.

As condições impostas aos contornos do sistema (referentes às linhas 1, 2, 3 e 4, esquematizadas na Figura 2A) envolveram, respectivamente, valores nulos das componentes de velocidade do ar nas direções y e z no eixo de simetria do sistema (X); valores nulos das componentes da velocidade do ar nas paredes do duto (condição de não deslizamento); velocidade homogênea do ar, igual à velocidade terminal do produto, na entrada do duto e pressão nula (referência) na saída do duto.

Predição da velocidade terminal do produto

Obtiveram-se valores para as velocidades terminais dos frutos, através do processo de otimização. Para tal, modificações foram feitas no programa fonte (elementos finitos) usado para a determinação das forças de arrasto no fruto de café, considerado esfera e elipse. Uma rotina de otimização foi introduzida no programa e os parâmetros de entrada passaram a ser: peso real do fruto, raio equivalente esférico (fruto esférico) ou dimensões "A" e "B" (fruto elíptico), os quais envolvem menores incertezas em seus valores que as velocidades terminais do produto obtidas experimentalmente. O objetivo da otimização foi determinar um valor para a velocidade terminal do produto, tal que o erro relativo entre a força de arrasto simulada usando esta velocidade e o peso real do fruto, fosse menor ou igual 0,1%.

RESULTADOS E DISCUSSÃO

Os valores experimentais das dimensões principais volume, densidade real e esfericidade de cada um dos dez frutos de café, usados nas simulações nos teores de umidade (TU) de 67,8; 46,9 e 18,7% b.u., são apresentados na Tabela 2.

Os valores de densidade do produto, apresentados na Tabela 2, de peso e de velocidade terminal (VT) expostos na Tabela 3, foram usados como entrada no programa computacional para estimar a força de arrasto no produto, usando-se a técnica de elementos finitos. Nas simulações, os valores de densidade e de viscosidade do ar escoando no duto foram considerados iguais a 1,2050 kg m^-3 e 1,8135x10^-5 Pa.s, respectivamente.

Para a análise em que se empregou a técnica de elementos finitos, o sistema (Figura 2) foi dividido em 634 elementos, resultando em 693 nós. O sistema discretizado em elementos pode ser visto na Figura 4: (A) visão global e (B) visão ampliada da região do corpo sólido (fruto de café).

Sendo o número de elementos disponíveis para a discretização do sistema (versão 5.3 do ANSYS) limitado, construíram-se duas regiões contornando o fruto, como mostrado na Figura 4C, para se obter elementos de tamanhos menores nos locais de maiores gradientes de velocidade e de pressão. Os valores das dimensões A, A1, A2, B, B1 e B2 usados para cada fruto de café, são apresentados na Tabela 4. No caso de fruto esférico, as dimensões A e B foram iguais ao raio da esfera equivalente (RE) e as duas regiões foram círculos de raios Ra1 = A1 = B1 e Ra2 = A2 = B2.

Neste trabalho, a introdução de uma viscosidade artificial (VA) para a execução do programa foi necessária para se obter convergência na solução. Os valores desta grandeza diferem para os frutos simulados e são mostrados na Tabela 3.

É interessante notar, na Tabela 3 que, em alguns casos, a consideração esférica para os frutos (células sombreadas), resultou em melhores resultados que a elíptica; entretanto, na maioria dos casos a consideração de uma forma elíptica para o produto produziu erros menores. Observa-se que os erros envolvidos nos valores das forças de arrasto para frutos com teores de umidade de 67,8; 46,9 e 18,7% b.u. variaram, respectivamente, de: (a) 9,0 a 90,2%; 6,6 a 42,8% e 8,3 a 55,5%, para o fruto com forma esférica e (b) 0,14 a 40,3%; 0,54 a 27,4% e 3,5 a 25,6%, para o fruto com forma elíptica. Em média, os erros envolvidos na estimativa da força de arrasto pela técnica de elementos finitos, para frutos considerados esfera e elipse foram, respectivamente, iguais a 34,0 ± 31,4% e 10,711,5%, para frutos a 67,8% b.u.; 20,014,8% e 9,510,2%, para frutos a 46,9% b.u.; e 25,416,4% e 13,18,2%, para frutos a 18,7% b.u.

Forças de arrasto simuladas

Os valores simulados das forças de arrasto usando-se a técnica de elementos finitos para frutos de café considerados esferas e elipses, são apresentados na Tabela 3. Os erros relativos envolvidos nos valores dessas forças encontram-se, também, na referida tabela. No cálculo desses erros considerou-se que o peso do fruto, determinado experimentalmente (Tabela 3), é o valor mais correto da força de arrasto.

Um dos problemas mais comuns encontrados na solução de problemas envolvendo escoamento de um fluido, quando se usa o ANSYS, é a divergência ou oscilação da pressão e/ou campo de velocidade. Uma das maneiras de se reduzir o impasse é a introdução de uma viscosidade artificial visando suavizar a solução para a velocidade em regiões com grandes gradientes. No caso de fluidos incompressíveis, a introdução desta viscosidade não altera a solução final. Na prática, o valor da viscosidade artificial deve ser menor que 1000 vezes a viscosidade efetiva para problemas incompressíveis turbulentos (ANSYS Theory Manual, 1996).

Em termos de média em relação aos dez frutos, os valores da força de arrasto obtidos experimentalmente e por simulação, considerando-se o fruto esférico e elíptico foram, respectivamente, iguais a 0,015 ± 0,003; 0,0190,004 e 0,0150,005 N, para frutos a 67,8% b.u.; 0,0090,001; 0,0090,002 e 0,0090,001 N, para frutos a 46,9% b.u. e 0,00560,0006; 0,0070,001 e 0,0060,001 N, para frutos a 18,7% b.u.

Em termos de média em relação aos dez frutos, os valores da força de arrasto obtidos experimentalmente e por simulação, considerando-se o fruto esférico e elíptico foram, respectivamente, iguais a 0,015 ± 0,003; 0,0190,004 e 0,0150,005 N, para frutos a 67,8% b.u.; 0,0090,001; 0,0090,002 e 0,0090,001 N, para frutos a 46,9% b.u. e 0,00560,0006; 0,0070,001 e 0,0060,001 N, para frutos a 18,7% b.u.

Dos parâmetros que podem ter influenciado nos erros envolvidos nos valores simulados das forças de arrasto, citam-se: (a) utilização de formas geométricas inadequadas para modelar o fruto de café; (b) considerações inadequadas na formulação do problema (forças cisalhantes e rugosidade na superfície do fruto desprezadas) e (c) incertezas nos valores experimentais da velocidade terminal (VT). Em média, as esfericidades dos frutos são iguais ou maiores que 85; 77 e 74% para frutos com teores de umidade de 67,8; 46,9 e 18,7% b.u., respectivamente, e seus valores parecem não ter uma relação direta com os erros nos valores simulados das forças de arrasto de frutos considerados esféricos.

Acredita-se que as incertezas nos valores experimentais da velocidade terminal foram um dos parâmetros que mais influenciaram nos valores da força de arrasto. Essas incertezas se devem, principalmente, à rotação do fruto durante sua flutuação e, também, a desvios no posicionamento do sensor de velocidade na seção transversal do duto.

Testes de comparação de médias (Newman-Keuls, a um nível de significância de 5%) mostraram que os valores simulados para força de arrasto dos frutos a 18,7 e 46,9% b.u., considerados elípticos e esféricos, não diferem dos pesos do produto determinados experimentalmente; entretanto, para o produto a 67,8% b.u., apenas a consideração elíptica para o fruto proporcionou forças de arrasto que não diferiram estatisticamente dos pesos experimentais do produto.

Velocidades terminais simuladas

Os valores das velocidades terminais provenientes das otimizações com os respectivos erros relativos, considerando-se o valor experimental da velocidade terminal do produto como o mais correto, são apresentados na Tabela 5, para frutos de café, tidos como esferas e elipses, a diferentes teores de umidade.

É interessante notar, na Tabela 5, que os valores das velocidades terminais otimizadas dos frutos considerados esferas são, na maioria das vezes, menores que os obtidos experimentalmente para as velocidades terminais do produto. Os valores simulados para as velocidades terminais de 100, 50 e 70% dos frutos a 18,7; 46,9 e 67,8% b.u., respectivamente, tendem a subestimar os valores das velocidades obtidas experimentalmente.

Em termos de média para os dez frutos, as velocidades determinadas experimentalmente e por simulação para frutos elípticos e esféricos foram, respectivamente iguais a: (a) 15,2 ± 1,3; 15,51,8 e 13,72,1 m s^-1 para frutos a 67,8% b.u.; (b) 13,90,9; 14,30,81 e 13,71,4 m s^-1 para o produto a 46,9% b.u. e (c) 11,70,8; 11,61,1 e 10,20,66 m s^-1 para frutos a 18,7% b.u. Os menores erros relativos médios (para os dez frutos) aconteceram quando o fruto foi considerado elíptico. Os maiores erros relativos médios envolvidos nas simulações das velocidades terminais de frutos, nos três teores de umidade, considerados elípticos e esféricos, foram iguais a 7,5 e 13,6%, respectivamente.

Testes de comparação de médias (Newman-Keuls, a um nível de significância de 5%), aplicados aos dados experimentais, mostraram que os valores simulados de velocidade terminal para o produto elíptico a 18,7 e 67,8% b.u., não diferem daqueles determinados de forma experimental; entretanto, estes valores são diferentes daqueles simulados ao se considerar o fruto como uma esfera. Quando os frutos se encontravam a 46,9% b.u., não se detectou diferença significativa entre os valores das velocidades terminais obtidas experimentalmente nem por simulações, ao se considerar os frutos elípticos e esféricos.

Portanto, valores de velocidade terminal de frutos de café a diferentes teores de umidade podem ser obtidos, por meio de simulação, conhecendo-se apenas o peso e as dimensões principais do produto.

CONCLUSÕES

1. A técnica dos elementos finitos mostrou-se ferramenta eficiente para simular forças de arrasto em frutos de café submetidos a um fluxo de ar.

2. Os valores de força de arrasto para o produto considerado com forma elíptica não diferiram estatisticamente dos pesos dos frutos determinados de maneira experimental.

3. A velocidade terminal de frutos de café pode ser obtida através de simulações usando-se a técnica dos elementos finitos; para tal, são necessários apenas dois parâmetros: o peso do produto e suas dimensões principais.

4. Independente do teor de umidade dos frutos, os valores simulados para a velocidade terminal e para força de arrasto no produto considerado elíptico são estatisticamente iguais àqueles determinados experimentalmente.

5. Os erros relativos médios (para os dez frutos) envolvidos no valores simulados para a velocidade terminal do produto em cada um dos teores de umidade, não foram superiores a 7,5 e 13,6% para o produto considerado elíptico e esférico, respectivamente.

LITERATURA CITADA

Protocolo 129 - 2/9/2002 - Aprovado em 16/2/2004

• Allen, C.A.W.; Watts, K.C. Properties of cowpeas (var Minica Beans). Journal of Agricultural Engineering Research, Silsoe, v.68, n.2, p.159-167. 1997.
• ANSYS Theory Reference. ANSYS Inc., Houston, PA, 1996.
• Arsenijevic, Z.Lj.; Grbavcic, Z.B.; Garic-Grulovic, R.V.; Zdanski, F.K. Determination of non-spherical particle terminal velocity using particulate expansion data. Powder Technology, Amsterdam, v.103, n.3, p.265-273. 1999.
• Carman, K. Some physical properties of lentil seeds. Journal of Agricultural Engineering Research, Silsoe, v.63, n.2, p. 87-92. 1996.
• Chhabra , R.P.; Agarwal, L.; Sinha, N.K. Drag on non-spherical particles: an evaluation of available methods. Powder Technology, Elsevier, v.101, n.3, p. 288-295. 1999.
• Dong, Z.; Gau, S.; Fryrear, D.W. Drag coefficients, roughness length and zero-plane displacement height as disturbed by artificial standing vegetation. Journal of Arid Environments, San Diego, v.49, n.3, p.485-505. 2001.
• Fox, R.W.; McDonald, A.T. Introdução à mecânica dos fluidos. 4.ed. Rio de Janeiro: Livros Técnicos e Científicos Editora S.A., 1998. 662p.
• Giles, R.V. Mecânica dos fluidos e hidráulica. 2.ed. São Paulo: McGraw-Hill do Brasil Ltda, 1980. 401p.
• Gupta, R.K.; Das, S.K. Physical properties of sunflower seeds. Journal of Agricultural Engineering Research, Silsoe, v.66, n.1, p.1-8. 1997.
• Moaveni, S. Finite element analysis. Theory and Application with ANSYS. 1.ed. Upper Saddle River: Prentice-Hall, 1999. 527p.
• Moreira, S.M.C.; Chaves, M.A.; Oliveira, L.M. Comparação da eficiência de líquidos na determinação da massa específica aparente de grãos agrícolas. Revista Brasileira de Arma-zenamento, Viçosa, v.9, n.1e2, p.22-24, 1984/85.
• Nguyen, A.V.; Stechemesser, H.; Zobel, G.; Schulze, H.J. An improved formula for terminal velocity of rigid spheres. International Journal of Mineral Processing, Amsterdam, v.50, n.1-2, p.53-61. 1997.
• Nimkar, P.M.; Chattopadhyay, P.K. Some physical properties of green gram. Journal of Agricultural Engineering Research, Silsoe, v.80, n.2, p.183-189. 2001.
• Singh, K.K.; Goswami, T.K. Physical properties of cumin seed. Journal of Agricultural Engineering Research, Silsoe, v.64, n.2, p.93-98. 1996.
• Suthar, S.H.; Das, S.K. Some physical properties of karingda [Citrullus lanatus (Thumb)] seeds. Journal of Agricultural Engineering Research, Silsoe, v.65, n.1, p.15-22. 1996.
• Tabak, S.; Wolf, D. Aerodynamics properties of cottonseeds. Journal of Agricultural Engineering Research, Silsoe, v.70, n.3, p.257-265. 1998.

Datas de Publicação

• Publicação nesta coleção
  26 Jan 2005
• Data do Fascículo
  Dez 2004

Histórico

• Aceito
  16 Fev 2004
• Recebido
  02 Set 2002