Ajuste de parâmetros de transporte de solutos no solo utilizando Matlab 6.5

O sucesso na utilização de modelos matemáticos no estudo do transporte de íons no solo está intimamente ligado à precisão com que os parâmetros de transporte envolvidos neste processo são estabelecidos. De maneira geral, tais parâmetros são determinados mediante a resolução de um problema de otimização não linear em que os dados experimentais, obtidos em ensaios de deslocamento miscível, são ajustados a um modelo teórico. Neste sentido, a utilização de softwares de alta performance no ajuste destes parâmetros mostra-se vantajosa, uma vez que, além da consistência e da disponibilidade de ferramentas numéricas preexistentes, possibilita a incorporação de novas rotinas de acordo com o fenômeno que se queira simular. Sendo assim, este trabalho teve como objetivo desenvolver, em ambiente MATLAB 6.5, uma rotina computacional para a otimização dos seguintes parâmetros de transporte: fator de retardamento (R) e coeficiente de dispersão (D). A rotina desenvolvida foi aplicada a dados experimentais de três ensaios de deslocamento miscível do íon potássio em colunas preenchidas com um Latossolo Vermelho- -Amarelo, fase arenosa. A qualidade dos ajustes obtidos foi avaliada utilizando-se do coeficiente de exatidão. Concluiu-se que a rotina proposta apresentou ótimo desempenho, o que, além de reforçar a consistência do método numérico utilizado, indica que a rotina proposta neste trabalho pode contribuir com o avanço dos estudos teóricos da dinâmica da água e de solutos em meios porosos não saturados.

BEM/FEM non-linear model applied to transient analysis with viscous damping

In this article a two-dimensional transient boundary element formulation based on the mass matrix approach is discussed. The implicit formulation of the method to deal with elastoplastic analysis is considered, as well as the way to deal with viscous damping effects. The time integration processes are based on the Newmark rhoand Houbolt methods, while the domain integrals for mass, elastoplastic and damping effects are carried out by the well known cell approximation technique. The boundary element algebraic relations are also coupled with finite element frame relations to solve stiffened domains. Some examples to illustrate the accuracy and efficiency of the proposed formulation are also presented.