Formação de nanopadrões em superfícies por sputtering iônico: Estudo numérico da equação anisotrópica amortecida de Kuramoto-Sivashinsky.

Apresenta-se uma abordagemnumérica para ummodelo que descreve a formação de padrões por sputtering iônico na superfície de ummaterial. Esse processo é responsável pela formação de padrões inesperadamente organizados, como ondulações, nanopontos e filas hexagonais de nanoburacos. Uma análise numérica de padrões preexistentes é proposta para investigar a dinâmica na superfície, baseada em ummodelo resumido em uma equação anisotrópica amortecida de Kuramoto-Sivashinsky, em uma superfície bidimensional com condições de contorno periódicas. Apesar de determinística, seu caráter altamente não-linear fornece uma rica gama de resultados, sendo possível descrever acuradamente diferentes padrões. Umesquema semi implícito de diferenças finitas com fatoração no tempo é aplicado na discretização da equação governante. Simulações foram realizadas com coeficientes realísticos relacionados aos parâmetros físicos (anisotropias, orientação do feixe, difusão). A estabilidade do esquema numérico foi analisada por testes de passo de tempo e espaçamento de malha, enquanto a verificação do mesmo foi realizada pelo Método das Soluções Manufaturadas. Ondulações e padrões hexagonais foram obtidos a partir de condições iniciais monomodais para determinados valores do coeficiente de amortecimento, enquanto caos espaço-temporal apareceu para valores inferiores. Os efeitos anisotrópicos na formação de padrões foramestudados, variando o ângulo de incidência.

Stable periodic waves in coupled Kuramoto-Sivashinsky-Korteweg-de Vries equations

Periodic waves are investigated in a system composed of a Kuramoto-Sivashinsky-Korteweg-de Vries (KS-KdV) equation linearly coupled to an extra linear dissipative one. The model describes, e.g., a two-layer liquid film flowing down an inclined plane. It has been recently shown that the system supports stable solitary pulses. We demonstrate that a perturbation analysis, based on the balance equation for the net field momentum, predicts the existence of stable cnoidal waves (CnWs) in the same system. It is found that the mean value u(0) of the wave field u in the main subsystem, but not the mean value of the extra field, affects the stability of the periodic waves. Three different areas can be distinguished inside the stability region in the parameter plane (L, u(0)), where L is the wave's period. In these areas, stable are, respectively, CnWs with positive velocity, constant solutions, and CnWs with negative velocity. Multistability, i.e., the coexistence of several attractors, including the waves with several maxima per period, appears at large value of L. The analytical predictions are completely confirmed by direct simulations. Stable waves are also found numerically in the limit of vanishing dispersion, when the KS-KdV equation goes over into the KS one.

Hydrothermal surface-wave instability and the Kuramoto-Sivashinsky equation

We consider a system formed by an infinite viscous liquid layer with a constant horizontal temperature gradient and a basic nonlinear bulk velocity profile. In the limit of long wavelength and large nondimensional surface tension we show that hydrothermal surface-wave instabilities may give rise to disturbances governed by the Kuramoto-Sivashinsky equation. A possible connection to hot-wire experiments is also discussed. © 1994.

Simulation of growth process of thin film on non-planar substrate

In (2 + 1) dimension, growth process of thin film on non-planar substrate in Kuramoto-Sivashinsky model is studied with numerical simulation approach. 15 x 15 semi-ellipsoids arranged orderly on the surface of substrate are used to represent initial rough surface. The results show that at the initial stage of growth process, the surface morphology of thin film appears to be grid-structure, and the interface width constantly decreases with the growth time, then reaches minimum. However, the grid-structure becomes ambiguous, and granules of different sizes distribute evenly on the surface of thin film with the increase of growth time. Thereafter, the average size of granules and the interface width gradually increase, and the surface morphology of thin film presents fractal properties. The numerical results of height-height correlation functions of thin film verify the surface morphology of thin film to be fractal for a longer growth time. By fitting of the height-height correlation functions of thin film with different growth times, the growth process is described quantitatively. (c) 2004 Elsevier B.V. All rights reserved.

Estimating correlated observation error statistics using an ensemble transform Kalman filter

For certain observing types, such as those that are remotely sensed, the observation errors are correlated and these correlations are state- and time-dependent. In this work, we develop a method for diagnosing and incorporating spatially correlated and time-dependent observation error in an ensemble data assimilation system. The method combines an ensemble transform Kalman filter with a method that uses statistical averages of background and analysis innovations to provide an estimate of the observation error covariance matrix. To evaluate the performance of the method, we perform identical twin experiments using the Lorenz ’96 and Kuramoto-Sivashinsky models. Using our approach, a good approximation to the true observation error covariance can be recovered in cases where the initial estimate of the error covariance is incorrect. Spatial observation error covariances where the length scale of the true covariance changes slowly in time can also be captured. We find that using the estimated correlated observation error in the assimilation improves the analysis.

Análise de erros da equação de advecção unidimensional no Método de Volumes Finitos

Uma análise utilizando a série de Taylor é apresentada para se estimar a priori os erros envolvidos na solução numérica da equação de advecção unidimensional com termo fonte, através do Método dos Volumes Finitos em uma malha do tipo uniforme e uma malha não uniforme. Também faz-se um estudo a posteriori para verificar a magnitude do erro de discretização e corroborar os resultados obtidos através da análise a priori. Por meio da técnica de solução manufaturada tem-se uma solução analítica para o problema, a qual facilita a análise dos resultados numéricos encontrados, e estuda-se ainda a influência das funções de interpolação UDS e CDS e do parâmetro u na solução numérica.

Implementação de um algoritmo numérico para solução da equação de Christoffel generalizada em acustoelasticidade.

Extensos estudos realizados nas últimas décadas sobre a propagação de ondas ultrassônicas em sólidos levaram ao desenvolvimento de técnicas não destrutivas para a avaliação da segurança e integridade de estruturas e componentes industriais. O interesse na aplicação de técnicas ultrassônicas para medição de tensões aplicadas e residuais decorre da mudança mensurável da velocidade das ondas ultrassônicas na presença de um campo de tensões, fenômeno conhecido como efeito acustoelástico. Uma teoria de acustoelasticidade fornece um meio atrativo e não destrutivo de medir a tensão média ao longo do caminho percorrido pela onda. O estudo da propagação das ondas ultrassônicas em meios homogêneos anisotrópicos sob tensão conduz a um problema não linear de autovalores dado pela equação de Christoffel generalizada. A característica não linear deste problema decorre da interdependência entre as constantes elásticas efetivas do material e as tensões atuantes. A medição experimental de tensões por técnicas ultrassônicas é um problema inverso da acustoelasticidade. Esta dissertação apresenta a implementação de um algoritmo numérico, baseado no método proposto por Degtyar e Rokhlin, para solução do problema inverso da acustoelasticidade em sólidos ortotrópicos sujeitos a um estado plano de tensões. A solução da equação de Christoffel generalizada apresenta dificuldades de natureza numérica e prática. A estabilidade e a precisão do algoritmo desenvolvido, bem como a influência das incertezas na medição experimental das velocidades das ondas ultrassônicas, foram então investigadas. Dados sintéticos para as velocidades das ondas ultrassônicas de incidência oblíqua em uma placa sujeita a um estado plano de tensões foram gerados pela solução direta da equação de Christoffel generalizada para ilustrar a aplicação do algoritmo desenvolvido. O objetivo maior desta dissertação é a disponibilização de uma nova ferramenta de cálculo para suporte às atividades experimentais de medição de tensões por ultrassom no país.

Modelagem e simulação do escoamento imiscível em meios porosos fractais descritos pela equação de Kozeny-Carman Generalizada

O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

Avaliação dos algoritmos de Picard-Krylov e Newton-Krylov na solução da equação de Richards

A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards.

Posicionamento aproximado do estado final para sistemas térmicos descritos pela equação do calor.

Neste trabalho, será considerado um problema de controle ótimo quadrático para a equação do calor em domínios retangulares com condição de fronteira do tipo Dirichlet é nos quais, a função de controle (dependente apenas no tempo) constitui um termo de fonte. Uma caracterização da solução ótima é obtida na forma de uma equação linear em um espaço de funções reais definidas no intervalo de tempo considerado. Em seguida, utiliza-se uma sequência de projeções em subespaços de dimensão finita para obter aproximações para o controle ótimo, o cada uma das quais pode ser gerada por um sistema linear de dimensão finita. A sequência de soluções aproximadas assim obtidas converge para a solução ótima do problema original. Finalmente, são apresentados resultados numéricos para domínios espaciais de dimensão 1.

Aplicação da equação de Poisson-Boltzmann modificada em sistemas biológicos: análise da partição iônica em um eritrócito

Neste trabalho, a partição iônica e o potencial de membrana em um eritrócito são analisados via equação de Poisson-Boltzmann modificada, considerando as interações não eletrostáticas presentes entre os íons e macromoléculas, assim como, o potencial β. Este potencial é atribuído à diferença de potencial químico de referência entre os meios intracelular e extracelular e ao transporte ativo de íons. O potencial de Gibbs-Donnan via equação de Poisson-Boltzmann na presença de carga fixa em um sistema contendo uma membrana semipermeável também é estudado. O método de aproximação paraboloide em elementos finitos em um sistema estacionário e unidimensionalé aplicado para resolver a equação de Poisson-Boltzmann em coordenadas cartesianas e esféricas. O parâmetro de dispersão relativo às interações não eletrostáticas écalculado via teoria de Lifshitz. Os resultados em relação ao potencial de Gibbs-Donnan mostram-se adequados, podendo ser calculado pela equação de Poisson-Boltzmann. No sistema contendo um eritrócito, quando o potencial β é considerado igual a zero, não se verifica a diferença iônica observada experimentalmente entre os meios intracelular e extracelular. Dessa forma, os potenciais não eletrostáticos calculados via teoria de Lifshitz têm apenas uma pequena influência no que se refere à alta concentração de íon K+ no meio intracelular em relação ao íon Na+

Equação de referência do Teste do Degrau medida em distância para indivíduos saudáveis e sedentários

O teste de caminhada em seis minutos (TC6M) avalia a capacidade respiratória durante o exercício. Recentemente, o teste do degrau em seis minutos (TD6M) está sendo estudado como uma proposta para essa mesma avaliação. Diante do exposto, o desenvolvimento de uma equação de referência se torna importante. O objetivo desse estudo foi desenvolver uma equação de referência padrão para o Teste do Degrau em Seis minutos. Esse estudo foi do tipo transversal, em que foram selecionados 452 indivíduos. Após a aplicação dos critérios de inclusão/exclusão, foram selecionados 326 sujeitos saudáveis e sedentários com idade entre 20 e 80 anos. Para serem considerados saudáveis, os participantes não podiam ter história de doenças (exceto hipertensão arterial sistêmica ou diabetes mellitus em tratamento) e foram submetidos à radiografia de tórax, espirometria e eletrocardiograma, que deveriam ser normais. O nível de sedentarismo foi avaliado através do International Physical Activity Questionnaire (IPAQ). Foram coletados os seguintes dados demográficos: idade, peso e altura. Todos os indivíduos realizaram TD6M na sua própria cadência (autocadenciado) em um degrau de 16,5 cm de altura, 65 de largura e 30 cm de comprimento. O número de subidas e descidas foi contado por um pedômetro digital. Foram mensuradas a pressão arterial, a frequência cardíaca e a saturação de oxigênio, antes e depois do TD6M. A análise estatística foi realizada pelo software STATA 12.0, e as equações foram desenvolvidas pelo modelo estatístico de regressão linear múltipla. Como resultado, observou-se que os participantes apresentaram exames normais, sendo 135 homens e 191 mulheres. O IPAQ demonstrou 157 ativos, 114 irregularmente ativos e 40 inativos, sendo que 14 indivíduos não responderam. A análise da distância alcançada em relação à idade e à diferença da frequência cardíaca tanto para homens quanto para mulheres mostrou significância estatística, demonstrando a sua importância para o desenvolvimento das equações para cada gênero. Para mulheres: Distância (m) = 88,83 - [(FC final- FC inicial)* 0,23] (Idade *0,37); para homens: Distância (m)= 110,20 - [(FC final- FC inicial)*0,18] (Idade * 0,59). O estudo concluiu que as equações de referência desenvolvidas nesse estudo foram realizadas em uma população de indivíduos saudáveis e sedentários e pode ser usada como padrão de referência para o TD6M.

Equação de volume para araucárias centenárias da Reserva Florestal Embrapa/Epagri.

2010

PASSADO E FUTURO NO PRESENTE DOS PERCURSOS DE QUALIFICAÇÃO: A EQUAÇÃO DA APRENDIZAGEM

É ponto assente que o desenvolvimento social e económico de um país dependedo conhecimento e do nível educativo dos seus cidadãos. Num momento de crise económica profunda em que a sustentabilidade do ensino superior se encontra ameaçada e a sua organização com rumo incerto, importa, mais do que nunca, criar um espaço de reflexão em que a ciência se pronuncie, aproximando e colocando em evidência os contributos que em vários domínios se vão produzindo.