Derivação e solução de equações modelo da dinâmica de gases rarefeitos

Neste trabalho, duas equações modedlo na área da dinâmica de gases rarefeitos, são derivadas a partir de algumas soluções exatas da equação linearizada de Boltzmann homogênea e não homogênea. Em adição, uma versão analítca do método de ordenadas discretas é usado para resolver problemas clássicos nesta área, descritos pelo "Modelo S". Resultados numéricos são apresentados para os problemas de fluxo de Couette, fluxo de Poiseuille, "Creep" Térmico, Deslizamento Térmico e problema de Kramers.

Desenvolvimento de um método de volumes finitos de alta ordem para a simulação de escoamentos de fluidos viscoelásticos

O uso da mecânica de fluidos computacional no estudo de processos envolvendo o escoamento de fluidos poliméricos está cada vez mais presente nas indústrias de transformação de polímeros. Um código computacional voltado a esta função, para que possa ser aplicado com sucesso, deve levar a predições mais próximas possível da realidade (modelagem), de uma forma relativamente rápida e eficiente (simulação). Em relação à etapa de modelagem, o ponto chave é a seleção de uma equação constitutiva que represente bem as características reológicas do fluido, dentre as diversas opções existentes. Para a etapa de simulação, ou seja, a resolução numérica das equações do modelo, existem diversas metodologias encontradas na literatura, cada qual com suas vantagens e desvantagens. Neste tópico se enquadra o trabalho em questão, que propõe uma nova metodologia para a resolução das equações governantes do escoamento de fluidos viscoelásticos. Esta se baseia no método dos volumes finitos, usando o arranjo co-localizado para as variáveis do problema, e na utilização de aproximações de alta ordem para os fluxos médios lineares e não-lineares e para outros termos não lineares que surgem da discretização das equações constitutivas. Nesta metodologia, trabalha-se com os valores médios das variáveis nos volumes durante todo o processo de resolução, sendo que os valores pontuais são obtidos ao final do procedimento via deconvolução. A solução do sistema de equações não lineares, resultante da discretização das equações, é feita de forma simultânea, usando o método de Newton São mostrados então, resultados da aplicação da metodologia proposta em problemas envolvendo escoamentos de fluidos newtonianos e fluidos viscoelásticos. Para descrever o comportamento reológico destes últimos, são usadas duas equações constitutivas, que são o modelo de Oldroyd-B e o modelo de Phan-Thien-Tanner Simplificado. Por estes resultados pode-se ver que a metodologia é muito promissora, apresentando algumas vantagens frente às metodologias convencionais em volumes finitos. A implementação atual da metodologia desenvolvida está restrita a malhas uniformes e, consequentemente, soluções para problemas com geometrias complexas, que necessitam de refinamento localizado da malha, foram obtidas somente para baixos números de Weissenberg, devido a limitação do custo computacional. Esta restrição pode ser contornada, tornando o seu uso competitivo.

Proposta de desenvolvimento de um método de tomada de decisão usando a teoria das restrições para sistemas de produção

Este trabalho discorre sobre a aplicação de uma proposta de método de resolução de problemas utilizando a Teoria das Restrições (TOC). No desenvolvimento desta proposta procurou-se aliar a focalização dos problemas e os conceitos financeiros protagonizados pela TOC. Para que uma solução urgente de problemas possa ser conduzida de uma forma coerente, precisa e rápida, não é necessário que seja executada de maneira que os gastos ultrapassem os benefícios. Além disso, o todo que compõe o sistema analisado deve poder ser vislumbrado claramente. Esta proposta foi aplicada a uma corporação que fabrica autopeças. Ela está localizada nos Estados Unidos da América. Os resultados da aplicação do método aqui apresentado permitiram chegar às conclusões apontadas no capítulo final deste trabalho.

Caracterização das ligas zinco-níquel eletrodepositadas a partir de solução de cloreto

Este trabalho apresenta um estudo da liga zinco-níquel obtida a partir de um banho ácido de cloretos, de formulação sugerida pela Boeing. A liga zinco níquel é uma proposta atual para a substituição de revestimentos de cádmio eletrodepositados, tanto na industria aeronáutica como na indústria automobilística. O interesse nesta substituição surge da necessidade de se encontrar um processo menos poluente ao meio ambiente e ao operador que aplica o revestimento, economicamente viável e que também atenda legislação ambiental, cada vez mais restritiva. As ligas zinco-níquel são de especial interesse, porque além de apresentarem as características descritas anteriormente, também conferem ao substrato ferroso, uma proteção do tipo catódico. O estudo aqui apresentado avaliou uma faixa de densidade de corrente de deposição, entre 0,5 e 5 A.dm-2 quanto aos teores de níquel das ligas formadas, a resistência a corrosão e outras características das ligas, através de curvas de polarização, voltametrias e ensaios de impedância eletroquímica e ensaios de névoa salina e outros ensaios . Especial atenção foi dada para a faixa de densidade de corrente considerada como ótima pela Boeing, entre 2 e 3 A.dm-2.

Um Modelo de estruturação de requisitos para o método fusion

Uma definição confiável dos requisitos de um software depende diretamente da completa e correta compreensão sobre as necessidades do sistema e sua conseqüente representação de forma adequada ao processo de desenvolvimento. Uma proposta de modelagem de requisitos deve apresentar qualidades que colaborem para a compreensão mútua das necessidades entre os envolvidos no processo e que organizem os requisitos de forma a permitir o acompanhamento no desenvolvimento do software. O presente trabalho apresenta um modelo de estruturação de requisitos fundamentado em metodologias orientadas a objetivos com utilização de cenários e preceitos da Teoria da Atividade. O modelo tem sua argumentação nas premissas que cliente e usuários normalmente expressam suas necessidades através de objetivos almejados e que a ação humana deve ser analisada dentro de um contexto para que possa fazer sentido e ser compreendida. Inserido no contexto do Projeto FILM1, cujo objetivo é expandir o Método Fusion, agregando uma etapa de modelagem de requisitos, o trabalho estabeleceu a qualidade de usabilidade como motivadora da definição de um modelo de estruturação de requisitos. A usabilidade é uma qualidade que visa facilitar a utilização do modelo como uma ferramenta de representação dos requisitos de forma inteligível, atuando tanto na especificação dos requisitos como na validação dos mesmos entre os envolvidos. Os requisitos são estruturados segundo uma abordagem voltada aos clientes e usuários do sistema. O modelo definido tem por objetivo prover a construção gradual e incremental do entendimento compartilhado entre os envolvidos sobre os domínios do problema e da solução, na concepção e no desenvolvimento do software. Metodologias orientadas a objetivos, operacionalizadas através de cenários, conjugadas a princípios da atividade oferecem um suporte adequado a estruturação de requisitos provendo usabilidade ao modelo. A avaliação da aplicabilidade do modelo é realizada com a modelagem de requisitos em três estudos de casos. Em cada caso são aplicadas técnicas de elicitação no sentido da afinar a sintonia com a estrutura do modelo de requisitos. A concepção do modelo, embasada em conceitos da Teoria da Atividade, é bastante adequado às atividades de elicitação em uma abordagem voltada aos clientes e usuários.

Solução da equação de condução de calor bidimensional, em meios multicompostos, pelos métodos nodal, com parâmetros concentrados, e a técnica da transformada de Laplace

Neste trabalho, desenvolvemos uma metodologia semi-analítica para solução de problemas de condução de calor bidimensional, não-estacionária em meios multicompostos. Esta metodologia combina os métodos nodal, com parâmetros concentrados, e a técnica da transformada de Laplace. Inicialmente, aplicamos o método nodal. Nele, a equação diferencial parcial que descreve o problema é integrada, transversalmente, em relação a uma das variáveis espaciais. Em seguida, é utilizado o método de parâmetros concentrados, onde a distribuição de temperatura nos contornos superior e inferior é substituída pelo seu valor médio. Os problemas diferenciais unidimensionais resultantes são então resolvidos com o uso da técnica da transformada de Laplace, cuja inversão é avaliada numericamente. O método proposto é usado na solução do problema de condução de calor, em paredes de edificações. A implementação computacional é feita, utilizando-se a linguagem FORTRAN e os resultados numéricos obtidos são comparados com os disponíveis na literatura.

FRIDA: um método para elicitação e modelagem de RNFs

Os requisitos direcionam o desenvolvimento de software porque são cruciais para a sua qualidade. Como conseqüência tanto requisitos funcionais quanto não funcionais devem ser identificados o mais cedo possível e sua elicitação deve ser precisa e completa. Os requisitos funcionais exercem um papel importante uma vez que expressam os serviços esperados pela aplicação. Por outro lado, os requisitos não funcionais estão relacionados com as restrições e propriedades aplicadas ao software. Este trabalho descreve como identificar requisitos não funcionais e seu mapeamento para aspectos. O desenvolvimento de software orientado a aspectos é apontado como a solução para os problemas envolvidos na elicitação e modelagem dos requisitos não funcionais. No modelo orientado a aspectos, o aspecto é considerado o elemento de primeira ordem, onde o software pode ser modelado com classes e aspectos. As classes são comumente usadas para modelar e implementar os requisitos funcionais, já os aspectos são adotados para a modelagem e implementação dos requisitos não funcionais. Desse modo, é proposta a modelagem dos requisitos não funcionais através das fases do ciclo de vida do software, desde as primeiras etapas do processo de desenvolvimento. Este trabalho apresenta o método chamado FRIDA – From RequIrements to Design using Aspects, cujo objetivo é determinar uma forma sistemática para elicitar e modelar tanto os requisitos funcionais quanto os não funcionais, desde as fases iniciais do ciclo de desenvolvimento. Em FRIDA, a elicitação dos requisitos não funcionais é realizada usando-se checklists e léxicos, os quais auxiliam o desenvolvedor a descobrir os aspectos globais – utilizados por toda a aplicação – bem como, os aspectos parciais que podem ser empregados somente a algumas partes da aplicação. O próximo passo consiste na identificação dos possíveis conflitos gerados entre aspectos e como resolvê-los. No método FRIDA, a identificação e resolução de conflitos é tão importante quanto a elicitação de requisitos não funcionais, nas primeiras fases do ciclo de vida do software. Além disso, é descrito como usar a matriz de conflitos para automatizar esse processo sempre que possível. A extração dos aspectos e sua modelagem visual são características muito importantes, suportadas pelo método, porque elas possibilitam a criação de modelos que podem ser reutilizados no futuro. Em FRIDA, é demonstrado como transformar os requisitos em elementos da fase de projeto (classes e aspectos) e como traduzir esses elementos em código. Outra característica do método FRIDA é que a conexão entre diagramas, que pertencem a diferentes fases do processo de desenvolvimento do software, permite um alto nível de rastreabilidade. Em resumo, FRIDA requer que o desenvolvedor migre de uma visão puramente funcional para outra que contemple também os requisitos não funcionais.

Sistematização de uma abordagem da medição de uma produção diversificada e seus desempenhos num ambiente industrial pelo método das Unidades de Esforço de Produção-UEPs

As atividades ligadas à medição da produção em uma organização industrial monoprodutora são facilmente feitas em função de uma única unidade de medida: a quantidade total produzida (controle físico), com seus custos totais de fabricação associados (controle monetário). A fabricação simultânea de vários produtos diversificados aumenta a complexidade dessas atividades, principalmente pela inexistência de uma unidade de medida comum a toda a produção. Este trabalho desenvolveu-se tendo em mente procurar uma solução alternativa aos métodos usuais utilizados para medir uma produção, quais sejam: número de horas, faturamento, quantidade física, entre outras, buscando racionalizar o processo de gestão Industrial, pela definição de uma unidade de medida comum a toda a produção o qual transforma uma fábrica multiprodutora real numa fábrica de um único produto representativo e equivalente a toda a produção. A Unidade de Esforço de Produção, UEPs, responde satisfatoriamente a essa necessidade, unificando e homogeneizando as atividades desenvolvidas por uma fábrica multiprodutora. Sua definição materializa sob uma mesma base de esforços (despesas) alocados aos vários postos operativos (máquinas e equipamentos) na fabricação dos diversos produtos. A organização escolhida para que fosse desenvolvido o presente trabalho, Estudo de Caso, serve como modelo para pesquisas futuras em outros tipos de organização que produzem produtos diversificados. O objetivo principal deste estudo é a Medição de uma Produção Diversificada e de seus Desempenhos, a saber: Custos, Rendimento, Eficiência, Produtividade tanto horária quanto econômica, num ambiente industrial pelo Método das UEPs. É uma contribuição para o controle das atividades no âmbito da gestão Industrial. Os resultados apresentados no presente trabalho indicam que as organizações industriais multiprodutoras tem neste Método (UEPs) uma ferramenta gerencial para poder Medir e Controlar os seus Desempenhos fazendo comparações entre os Desempenhos Planejados (Metas) e os Desempenhos verdadeiramente ocorridos.

Solução LTSn da equação adjunta de transporte de nêutrons com fonte arbitrária para elevada ordem de quadratura numa placa homogênea

O objetivo deste trabalho consiste em estender o método LTSN à solução do problema adjunto de transporte de nêutrons. A solução adjunta é interpretada fisicamente como uma função importância que designa a capacidade de contribuição de cada cela do espaço de fase para um funcional resposta. A derivação desta interpretação, através do princípio variacional, está sucintamente apresentada. Surgida da necessidade de generalização da fonte adjunta, também propõe-se uma nova formulação LTSN capaz de resolver problemas de transporte, tanto direto quanto adjunto, com fonte arbitrária, para elevada ordem de quadratura em geometria de placa. Esta nova formulção inspira-se na propriedade de invariância de projeção dos meios isotrópicos mas também é válida para os meios anisotrópicos. Todos os resultados apresentados pelas simulações numéricas de problemas adjuntos são calculados pela nova formulação LTSN e são comparados ou com a definição de função importância ou pelas relações de reciprocidade ou pelo código ANISN.

Ensaio clínico randomizado da preservação de enxertos hepáticos com solução de perfusão de UW e HTK no transplante ortotópico de fígado

Introdução: A solução de UW (University of Wisconsin) tem sido utilizada como padrão para preservação de enxertos hepáticos para transplante ortotópico de fígado (TOF) desde 1990. Seu custo é alto, e não existem, em nosso meio, estudos clínicos que comparem a sua efetividade com a de outras soluções. A solução de HTK (histidina-triptofano-cetoglutarato) foi desenvolvida inicialmente para cardioplegia, porém estudos experimentais e alguns estudos clínicos retrospectivos demonstraram sua eficácia na preservação hepática, mesmo com tempos de isquemia prolongados. Objetivos: Comparar a efetividade das soluções de preservação de órgãos HTK e UW com relação as variáveis: disfunção primária do enxerto (DPE), tempo de isquemia fria (TIF), complicações de via biliar (CVB), alterações de provas funcionais hepáticas (PFH) e sobrevida do enxerto e dos pacientes. Método: Foram estudados os fígados de doadores de múltiplos órgãos, implantados nos receptores pela técnica de piggyback, segundo ordem cronológica de ingresso em lista de espera no RS no período janeiro de 2003 a agosto de 2004. As soluções de preservação HTK e UW foram utilizadas de forma randomizada em blocos. A perfusão na aorta foi feita com 4 litros de HTK ou 2 litros de UW e a perfusão venosa portal com 1 litro de ambas as soluções, utilizando-se 500 ml para perfusão venosa e arterial adicional ex-situ e armazenagem do enxerto. Realizou-se biópsia hepática em cunha do lobo esquerdo quando a esteatose macroscópica estava presente. A análise bioquímica sérica foi diária na primeira semana e, em 15 e 30 dias pós-operatórios. Resultados: Foram estudados 102 pacientes submetidos ao TOF, sendo 65 no grupo UW (63,7%) e 37 no grupo HTK (36,3%). As frequências de sexo, raça, estado hemodinâmico, o uso de vasopressores e a presença de esteatose nos doadores foram igual nos dois grupos (pα>,05). A idade média dos doadores foi de 38,1 anos (DP +-14,4) no grupo UW e de 44,6 anos (DP +-14,2) no HTK (pα=,036). A distribuição de sexo, raça, idade, etiologia da cirrose, re-transplante, hepatite fulminante, trombose portal e escore de Child-Pugh dos receptores foi igual nos dois grupos (pα>,05). O grupo HTK teve 8 casos (25,8%) de CVB (4 estenose, 2 fístulas e 2 lesões do tipo isquêmica) contra 5 casos (8,6%) do grupo UW (pα=,033) em 89 pacientes que completaram 4 meses de seguimento (OR=2,0; IC 95%=1,2 a 3,5). A média do TIF nos dois grupos foi semelhante (UW= 579,2 min.; HTK= 527,9 min. pα>,05) e não houve diferenças nas incidências de CVB, DPE e óbito com relação a TIF estratificados entre os grupos. Não houve variação nas medianas das PFH (pα>,05 para BT, AST, ALT, FA, GGT, LDH, Fator 5 e TP). A incidência de óbito foi similar em ambos os grupos: UW= 6 (9,4%) e HTK= 4 (11,1%). A incidências de DPE foi de 2,8% no grupo HTK (1 caso) e 9,4% no grupo UW (6 casos) (pα=0,15), dos quais 5 (71,4%) evoluíram para o óbito por não funcionamento primário do enxerto, com ou sem outras morbidades. Conclusão: As soluções de UW e HTK foram igualmente efetivas na preservação dos enxertos hepáticos de doadores cadavéricos na amostra analisada, considerando-se aspectos clínicos, laboratoriais e sobrevida dos pacientes e dos enxertos. A utilização rotineira da solução de HTK poderá diminuir os custos do TOF. A média de idade maior dos doadores e a utilização de um volume reduzido de solução podem ter contribuído para uma incidência maior de CVB no grupo HTK.

Solução LTSn para problemas de transferência radiativa com polarização em geometria plana

O método LTSN tem sido utilizado na resolução de uma classe abrangente de problemas de transporte de partículas neutras que são reduzidos a um sistema linear algébrico depois da aplicação da transformada de Laplace. Na maioria dos casos estudados os autovalores associados são reais e simétricos. Para o problema de criticalidade os autovalores associados são reais ou imaginários puros e simétricos, e para o o problema de multigrupo podem aparecer autovalores complexos. O objetivo deste trabalho consiste na generalização da formulação LTSN para problemas de transporte com autovalores complexos. Por esse motivo é focada a solução de um problema radiativo de transporte com polarização em uma placa plana. A solução apresentada fundamenta-se na aplicação da transformada de Laplace ao conjunto de equações SN dos problemas resultantes da decomposição da equação de transferência radiativa com polarização em série de Fourier, seguindo o procedimento de Chandrasekhar. Esse procedimento gera 2L + 2 sistemas lineares de ordem 4N dependentes do parâmetro complexo "s". Aqui, L é o grau de anisotropia e N a ordem de quadratura. A solução desse sistema simbólico é obtida através da aplicação da transformada inversa de Laplace depois da inversão da matriz simbólica pelo método da diagonalização. Para a obtenção das constantes de integração é assumido que os componentes do vetor de Stokes são reais e as matrizes dos autovalores e autovetores são separadas em suas partes real e imaginária. A solução LTSN para autovalores complexos é validada através da comparação da solução para uma placa com espessura unitária, grau de anisotropia L = 13, albedo de espalhamento simples $ = 0:99, coe ciente de re exão de Lambert ¸0 = 0:1 e N = 150, segundo dados da literatura consultada.

Síntese de pós de alumina nanocristalina por combustão em solução

Esta dissertação teve como objetivo o estudo e desenvolvimento da técnica de síntese por combustão em solução para obtenção de nanocristais de alumina, e a caracterização desse material por diferentes técnicas. O estudo da obtenção incluiu o modelamento termodinâmico das reações investigadas e a influência da razão combustíveloxidante no tamanho de partícula. A síntese da alumina por combustão foi realizada a partir do par precursor químico glicina-nitrato de alumínio em água, para diferentes razões combustível-oxidante, temperaturas de chama e número total de moles gasosos, relacionando esses parâmetros com tamanho de cristalito no pó, deformação de rede, área superficial e porosidade total. As reações químicas de combustão em solução resultaram na formação in situ de fases metaestáveis de Al2O3 para o pó como-sintetizado. A fase Al2O3-α, cuja formação não ocorreu nem na presença de um extra-oxidante, somente foi obtida após um tratamento térmico a 1100°C do pó como-sintetizado. O modelamento termodinâmico da reação de combustão em solução mostrou que quando a razão combustível-oxidante aumenta, obtém-se uma elevação da temperatura de chama adiabática e da quantidade de gás produzida, definindo características do particulado como morfologia, tamanho de cristalito, área superficial e nível de aglomeração. A caracterização do pó utilizou técnicas como análise termodiferencial (ATD) e termogravimétrica (ATG), granulometria por dispersão de laser (GDL), BET para análise de área superficial e porosidade total, difração de raios X (DRX), microscopia eletrônica de varredura (MEV), e microscopia eletrônica de transmissão (MET). Os pós de alumina alfa obtidos apresentaram-se como constituídos de cristalitos nanométricos, dispostos na forma de agregados de tamanho micrométrico. O tamanho de cristalito médio, calculado pelo método de Williamson-Hall, a partir de dados de análises por DRX, foi de 97 nm, e deformação média calculada de 5x10-3 % A análise granulométrica dos pós como-sintetizados verificou um tamanho de aglomerado de 40 µm, com distribuição monomodal. As análises por MEV e MET confirmaram a aglomeração das partículas do pó como-sintetizado, a partir de partículas com morfologia irregular.

Desenvolvimento de método analítico para determinação de Mg, Ti e Cl em catalisadores ziegler-natta por WDXRF

O desenvolvimento de um método analítico para determinação de Mg, Ti e Cl emcatalisadores suportados Ziegler-Natta (ZN) por espectroscopia de fluorescência de raios Xcom dispersão em comprimento de onda (WDXRF) foi avaliado. Três famílias decatalisadores com distintos teores de Mg, Ti e Cl foram utilizados. Para fins de comparação,os analitos foram também determinados por métodos clássicos, a saber, Mg porcomplexometria com EDTA, Ti por espectrofotometria e Cl por argentimetria. A naturezahigroscópica desses catalisadores impossibilitou a análise direta por WDXRF. Melhoresresultados na preparação da amostra foram obtidos através da calcinação (1000 0C), seguida de mistura com aglomerante (H3BO3), na proporção 1:4, para a determinação da razãoMg/Ti. O estudo do efeito da pressão sobre a intensidade do sinal analítico de Mg e de Ti evidencioua necessidade de prensagem das amostras em pressões superiores a 200 MPa. No caso de Mg e Ti, a determinação por Parâmetros Fundamentais (FP), a partir da medida da linha espectral do Ti, forneceu resultados equivalentes àqueles obtidos porcalibração univariada e por métodos clássicos. No entanto, a utilização de FP a partir damedida da linha espectral do Mg mostrou-se inadequada, provavelmente, devido ao baixo valor do número atômico do Mg. Nas três famílias de catalisadores ZN a razão Mg/Tiobservada foi 0,44 ± 0,01; 1,8 ± 0,05 e 4,9 ± 0,02. A determinação de Cl foi estudada através de extração em meio aquoso, seguida pordeposição de um alíquota da solução sobre diferentes suportes (filmes finos) e posteriorevaporação do solvente. Foram avaliados os suportes papel filtro quantitativo, membranasMillipore FHLC Teflon®, de difluoreto de polivinilideno, celulose, Nylon®, fita Kapton®(poliamida) e Mylar® (poliéster). Devido à natureza do Cl, houve perdas do elemento químicopor volatilização, durante as medidas, possivelmente devido ao vácuo e aquecimento daamostra dentro da câmara de análise do equipamento. A tentativa de fixação do Cl comoAgCl, seguida por deposição sobre a membrana Millipore FHLC foi considerada adequada.No entanto, a média dos resultados obtidos por WDXRF mostrou-se inferior quando comparada ao método argentimétrico utilizado como referência, para um nível de certeza de95%.

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.

Aerodinâmica de automóveis baseado no método dos contornos virtuais

Este trabalho desenvolve um método numérico para a solução de escoamentos bidimensionais em torno de geometrias automobilísticas utilizando o método de diferenças finitas. O código computacional resolve as equações de Navier-Stokes e de Euler para uma distribuição adequada dos pontos discretos na malha. O método de integração empregado baseia-se no esquema explícito de Runge-Kutta de 3 estágios para as equações da quantidade de movimento e no de sub-relaxações sucessivas para a pressão na base Gauss-Seidel. Utilizou-se a técnica dos contornos virtuais em coordenadas cartesianas para resolver o escoamento sobre uma geometria simplificada, com a superfície coincidente com a malha computacional, e uma geometria automobilística mais complexa (BMW). Para a certificação da técnica empregada, optou-se pela utilização da teoria do escoamento potencial e pela comparação com dados experimentais encontrados na literatura e outros coletados em túnel de vento em escala reduzida. Houve dificuldade nesta comparação devido à falta de artigos relativos às simulações numéricas de escoamentos sobre automóveis e na aplicação da técnica dos contornos virtuais em geometrias complexas. Os resultados foram satisfatórios, com boas perspectivas para trabalhos futuros, contribuindo assim para o desenvolvimento da área.