Uma solução híbrida para a equação de Fokker-Planck dependente da energia usando a técnica da transformada de Laplace e diferenças finitas

Neste trabalho é obtida uma solução híbrida para a equação de Fokker-Planck dependente da energia, muito utilizada em problemas de implantação iônica. A idéia consiste na aplicação da transformada de Laplace na variável de energia e aplicação de um esquema de diferenças finitas nas variáveis espacial e angular desta equação. Tal procedimento gera um problema matricial simbólico para a energia transformada. Para resolver este sistema, procede-se a inversão de Laplace da matriz (sI+A), onde s é um parâmetro complexo, I a matriz identidade e A uma matriz quadrada gerada pela discretização das variáveis espacial e angular. A matriz A não é diagonalizável, desta forma, contorna-se este problema decompondo esta matriz na soma de outras duas, onde uma delas é diagonalizável. É gerado então um método iterativo de inversão, semelhante ao método da fonte fixa associado ao método de diagonalização, do qual o resultado fornecido são os valores para o fluxo de partículas do sistema. A partir disto pode-se determinar a energia depositada no sistema eletrônico e nuclear do alvo. Para validar os resultados obtidos faz-se a simulação de implantação de íons de B em Si numa faixa energética de 1keV a 50MeV, comparam-se os resultados com simulação gerada numericamente pelo software SRIM2003.

Equação de burgers propriedades e comportamento assintótico

No presente trabalho, obtemos e analisamos diversas propriedades das soluções u(·, t) da equação de difusão linear (equação do calor em meios unidimensionais homogêneos) ut = μuxx x 2 R, t > 0 correspondentes a estados iniciais u(x, 0) = u0(x), com u0 2 Lp(R), para algum 1 p < 1; bem como da equação de Burgers ut + cuux = μuxx x 2 R, t > 0 onde c, μ são constantes dadas, sendo c 6= 0 e μ > 0 e ainda assumindo u(x, 0) = u0(x) com u0 2 Lp(R) para 1 p < 1, e limitado. Estudamos também a equação mais geral da forma ut + f(u)x = μuxx x 2 R, t > 0 discutindo várias propriedades importantes das soluções, associadas a estados iniciais u0 2 Lp(R) \ L1(R) para algum 1 p < 1. Em particular, examinamos o comportamento de ku(·, t)kLr(R), p r 1, para t >> 1, e diversas propriedades relacionadas.

Simulação numérica da transferência simultânea de energia e umidade através do solo em um sistema trocador-armazenador de calor

O presente trabalho pesquisa o fenômeno da transferência simultânea de calor e umidade em solos insaturados que envolvem dutos enterrados. Estes dutos e o solo envolvente compõem o chamado sistema trocador-armazenador de calor no solo, que é muitas vezes utilizado como fonte complementar de aquecimento em estufas solares agrícolas. O funcionamento do sistema trocador-armazenador de calor é baseado na energia armazenada no solo durante os horários de máxima insolação, sendo que durante a noite parte desta energia é recuperada. Durante o dia o ar quente do meio interno da estufa solar é bombeado para dentro do feixe de tubos enterrados, que devolvem o ar mais frio na outra extremidade. Por outro lado, durante a noite o ar mais frio do meio interno da estufa é bombeado para dentro dos dutos, os quais efetuam troca térmica com o solo envolvente, que neste horário possui a temperatura mais elevada do sistema. O objetivo geral deste trabalho é resolver o problema transiente periódico e tridimensional da transferência simultânea de calor e umidade em solos não-saturados, que compõem o sistema trocador-armazenador de calor, utilizando a simulação numérica. Como objetivos secundários têm-se: o melhoramento do sistema de troca térmica, a quantificação da parcela de calor transportado pela difusão da umidade no solo e a análise dos campos de temperatura e de conteúdo de umidade, sendo que a análise dos campos de umidade permite verificar se existe formação de frentes de secagem significativas na vizinhança dos dutos durante os sucessivos períodos de aquecimento e resfriamento a que o sistema é submetido. Na resolução do problema em questão é empregado o modelo clássico de Philip e De Vries para a transferência simultânea de calor e massa em meios porosos insaturados Neste modelo, as equações de conservação de energia e massa obtidas trazem explicitamente as influências combinadas dos gradientes de temperatura e de conteúdo de umidade nos processos de transporte de calor e umidade. O sistema de equações diferenciais governantes do problema em questão é resolvido numericamente utilizando o Método dos Volumes Finitos e na discretização destas equações é usada uma integração temporal totalmente implícita. Todas as propriedades difusivas e termofísicas empregadas são consideradas variáveis com a temperatura e o conteúdo de umidade. Os dutos de seção transversal circular do sistema trocadorarmazenador de calor no solo, são modelados como dutos de seção transversal quadrada de área equivalente para que coordenadas cartesianas possam ser utilizadas nos modelos analisados. Neste trabalho são simulados quatro modelos computacionais associados ao sistema trocador-armazenador de calor no solo. Estes modelos são compostos por: um duto isolado, um duto com convecção, dois dutos isolados e dois dutos com convecção. A variação da temperatura do ar na entrada do(s) escoamento(s), assim como a temperatura do meio ambiente, para os modelos com convecção, é dada por uma senóide com uma amplitude de 14 ºC. No modelo de um duto isolado, são realizadas simulações utilizando várias combinações dos parâmetros do modelo em questão e os resultados, assim obtidos, são comparados com aqueles encontrados na literatura Visando melhorar o sistema de troca térmica dos modelos computacionais investigados, são selecionados valores e intervalos de valores recomendados para os parâmetros do modelo de um duto isolado. Para este modelo, com um diâmetro de 0,1 m, são escolhidos valores (ou intervalos de valores) recomendados: de 4 m/s para a velocidade do escoamento interno dentro do duto, de 0,25 para o conteúdo de umidade do solo, de 5 até 20 metros para o comprimento do duto e de 0,20 até 0,30 m para a distância entre centros do dutos. As simulações dos quatro modelos computacionais realizadas utilizando as várias combinações dos valores recomendados para os parâmetros destes modelos, mostrou que não há diferença significativa entre os valores de calor volumétrico armazenado no solo empregando a resolução acoplada das equações de energia e de massa e a resolução da equação da temperatura. Mesmo para os modelos de um e de dois dutos com convecção a diferença percentual encontrada foi insignificante. Finalmente, são apresentados e analisados os campos de temperatura e de conteúdo de umidade para os quatro modelos computacionais avaliados. Os perfis de temperatura e de conteúdo de umidade em diferentes horários mostraram que, durante o dia, o solo absorve calor dos escoamentos internos de ar e, uma vez que, junto à superfície dos dutos tem-se regiões de maior temperatura, há, conseqüentemente, uma migração da umidade nestas regiões. Durante a noite, ocorre o contrário, o solo fornece calor aos escoamentos dentro dos dutos, e, desta forma, as regiões próximas aos dutos apresentam níveis de conteúdo de umidade superiores ao inicial. Ainda, os perfis de conteúdo de umidade para todas as situações analisadas mostraram que, não há formação de frentes de secagem significativas nas proximidades dos dutos que compõem os quatro modelos computacionais avaliados.

Solução da equação de difusão-advecção para uma CLP não-homogênea e não-estacionária pelo método GILTT

Neste trabalho é apresentada a solução da equação de difusão-advecção transiente para simular a dispersão de poluentes na Camada Limite Planetária. A solução é obtida através do método analítico GILTT (Generalized Integral Laplace Transform Technique) e da técnica de inversão numérica da quadratura de Gauss. A validação da solução é comprovada utilizando as concentraçãos obtidas a partir do modelo com as obtidas experimentalmente pelo Experimento de Copenhagen. Nesta comparação foram utilizados os perfis de vento potencial e logaritmo e os parâmetros de turbulência propostos por Degrazia et al (1997) [19] e (2002) [17]. Os melhores resultados foram obtidos utilizando o perfil de vento potencial e o coeficiente de difusão propostos por Degrazia et al (1997). A influência da velocidade vertical é mostrada através do comportamento das concentrações de poluentes na pluma. Além disso, as velocidades verticais e longitudinais geradas pelo Large Eddy Simulation (LES) foram colocadas no modelo para poder simular uma camada limite turbulenta mais realística, a qual apresentou resultados satisfatórios quando comparados com os disponíveis na literatura.

Alguns resultados para equação do calor e equações de advecção-difusão

Neste trabalho, são obtidas diversas propriedades (em especial, referentes ao comportamento ao t -+ +00) das soluções u(', t) da equação linear do calor, Ut = div(AV'u), x E JRn, t > O onde A E JRnxné uma matriz constante simétrica e positiva definida, correspondentes a estados iniciais p-somáveis, i.e., u(x, O) = uo(x), Uo E LP(JRn), onde 1 :::;p < 00. Em particular, é examinado o comportamento de Ilu(., t)IILP(lRn) ao t -+ +00, mostrando-se que Ilu(., t)IILl(lRn)-+ Ikn u(x, O)dXI quando p = 1, e Ilu(-' t)IILP(lRn)-+ O quando p > 1. São analisadas, também, as taxas de decaimento e o comportamento assintótico das soluções u(', t) de equações de advecção-difusão da forma Ut + divf(u) = div(A(u)V'u), x E JRn, t > O correspondentes a estados iniciais p-somáveis e limitados, i.e., u(x, O)= uo(x), u(', O) E LP(JRn) n LOO(JRn), onde 1 :::;p :::; 2. Novamente, é examinado o comportamento de Ilu(" t)IILP(lRn)ao t -+ +00, mostrando-se que Ilu(., t)IILl(lRn)-+ Ikn u(x, O)dxl quando p = 1, e Ilu(" t)IILP(lRn)-+ O quando p > 1. Várias outras propriedades importantes são também discutidas, seguindo principalmente [Silva, 2003], [Crandall e Tartar, 1980], [Hagstrom et al., 2003], [Zingano, 1999], [Zingano, 2004a], [Zingano, 2004b].

Avaliação das pressões respiratórias máximas em crianças e adolescentes da grande Natal: elaboração de uma equação preditiva

The human respiratory system was so designed that would allow efficient ventilation, regardless of variations in the external environment that may hinder the act of breathing, such an act involves dozens of variables, among them we find the respiratory depression, which is nothing more than respiratory muscle strength. The pressures are widely used in several cases: Neuro-muscular; evolution of pulmonary dysfunction and a predictor for discontinuation of mechanical ventilation. Therefore it was proposed to carry out evaluations of these respiratory pressures for children and adolescents aged 10 to 16 years and propose a predictive equation that involves the anthropometric variables age (A, years), body mass (BM, kilograms) and height (H, meters) with maximal respiratory pressures (maximum inspiratory and expiratory pressure). Evaluations were performed in this age group of students in public and private schools of the Grande Natal , measurements were performed using the analogue manometer, were children and adolescents and their parents gave informed consent. 517 samples were taken, and 250 for males (M), 255 for females (F) and 12 were excluded according to our exclusion criteria. The sample was subdivided into three age groups (10-11, 12-13 and 14 to 16 years old). It was found through the student s t test (p ≤ 0.05) for all variables studied, children and male adolescents had higher means than females, except for the MC. For the correlation between the variables found significant correlation (p <0.05) among all the variables when analyzed as pairs except between MIP and height for females. The development of predictive equations (for p ≤ 0.05) based on three types of strategies adopted were restricted to two association between anthropometric variables isolated, resulting in: for males: MIP = -32.29 + (-2.11*A) + (-0.52*BM), MIP = 9.99 + (-0.36*BM) + (-49.40*H); MEP = 18.54 + 3.53*A + 0, 42*BM, MEP = -33.37 + 2.78*A + 52.18* H, MEP = -17.39 + 0.33*BM + 55.04*H; and, for females we find: MEP = 24.32 + 2.59 * A + 0.24*BM

Equação de regressão linear múltipla para estimativa do erro experimental

Este trabalho teve por objetivo estimar equações de regressão linear múltipla tendo, como variáveis explicativas, as demais características avaliadas em experimento de milho e, como variáveis principais, a diferença mínima significativa em percentagem da média (DMS%) e quadrado médio do erro (QMe), para peso de grãos. Com 610 experimentos conduzidos na Rede de Ensaios Nacionais de Competição de Cultivares de Milho, realizados entre 1986 e 1996 (522 experimentos) e em 1997 (88 experimentos), estimaram-se duas equações de regressão, com os 522 experimentos, validando estas pela análise de regressão simples entre os valores reais e os estimados pelas equações, com os 88 restantes, observando que, para a DMS% a equação não estimava o mesmo valor que a fórmula original e, para o QMe, a equação poderia ser utilizada na estimação. Com o teste de Lilliefors, verificou-se que os valores do QMe aderiam à distribuição normal padrão e foi construída uma tabela de classificação dos valores do QMe, baseada nos valores observados na análise da variância dos experimentos e nos estimados pela equação de regressão.

Análise espacial dos fatores da equação universal de perda de solo em área de nascentes

O objetivo deste trabalho foi avaliar a perda de solo de área de nascentes da Microbacia do Córrego do Tijuco, SP. Foi utilizada a análise espacial dos fatores da equação universal da perda de solo (EUPS), em integração com análise de componentes principais e geoestatística. A perda de solo média, estimada para a área, foi de 118,5 Mg ha‑1 por ano, considerada alta. Próximo à zona urbana, houve alta interação dos fatores erosividade da chuva e práticas conservacionistas, o que evidencia grande perda de solo, em razão da concentração da água proveniente da camada impermeabilizada urbana, com alta velocidade de escoamento. Nos divisores de águas, a atuação da erodibilidade foi proeminente, em contraste com o fator topográfico. Foram observadas áreas com atuação conjunta destes fatores, inclusive em locais de inclinação suave, porém com alto potencial natural de erosão. A interação das análises multivariadas e geoestatística permite a estratificação da área, identifica locais com propriedades específicas quanto à perda de solo, e espacializa os fatores do processo erosivo e suas interações ao longo do relevo.

Estimativas de mortalidade para a região nordeste do Brasil em 2010: uma associação do método demográfico equação geral de balanceamento, com o estimador bayesiano empírico

One of the greatest challenges of demography, nowadays, is to obtain estimates of mortality, in a consistent manner, mainly in small areas. The lack of this information, hinders public health actions and leads to impairment of quality of classification of deaths, generating concern on the part of demographers and epidemiologists in obtaining reliable statistics of mortality in the country. In this context, the objective of this work is to obtain estimates of deaths adjustment factors for correction of adult mortality, by States, meso-regions and age groups in the northeastern region, in 2010. The proposal is based on two lines of observation: a demographic one and a statistical one, considering also two areas of coverage in the States of the Northeast region, the meso-regions, as larger areas and counties, as small areas. The methodological principle is to use the General Equation and Balancing demographic method or General Growth Balance to correct the observed deaths, in larger areas (meso-regions) of the states, since they are less prone to breakage of methodological assumptions. In the sequence, it will be applied the statistical empirical Bayesian estimator method, considering as sum of deaths in the meso-regions, the death value corrected by the demographic method, and as reference of observation of smaller area, the observed deaths in small areas (counties). As results of this combination, a smoothing effect on the degree of coverage of deaths is obtained, due to the association with the empirical Bayesian Estimator, and the possibility of evaluating the degree of coverage of deaths by age groups at counties, meso-regions and states levels, with the advantage of estimete adjustment factors, according to the desired level of aggregation. The results grouped by State, point to a significant improvement of the degree of coverage of deaths, according to the combination of the methods with values above 80%. Alagoas (0.88), Bahia (0.90), Ceará (0.90), Maranhão (0.84), Paraíba (0.88), Pernambuco (0.93), Piauí (0.85), Rio Grande do Norte (0.89) and Sergipe (0.92). Advances in the control of the registry information in the health system, linked to improvements in socioeconomic conditions and urbanization of the counties, in the last decade, provided a better quality of information registry of deaths in small areas

Rotina computacional e equação simplificada para modelar o transporte de sedimentos num Latossolo Vermelho Distrófico

Os objetivos do presente trabalho foram desenvolver rotina computacional para a solução da equação de Yalin e do diagrama de Shields e avaliar uma equação simplificada para modelar a capacidade de transporte de sedimento num Latossolo Vermelho Distrófico que possa ser utilizada no Water Erosion Prediction Project - WEPP, assim como em outros modelos de predição da erosão do solo. A capacidade de transporte de sedimento para o fluxo superficial foi representada como função-potência da tensão cisalhante, a qual revelou ser aproximação da equação de Yalin. Essa equação simplificada pôde ser aplicada em resultados experimentais oriundos de topografia complexa. A equação simplificada demonstrou acuracidade em relação à equação de Yalin, quando calibrada utilizando-se da tensão média cisalhante. Testes de validação com dados independentes demonstraram que a equação simplificada foi eficiente para estimar a capacidade de transporte de sedimento.