[pt] Mudanças das propriedades dos solos devido à ação de carregamentos dinâmicos são responsáveis por danos significativos em geo-estruturas, tais como: barragens, estruturas de concentração, fundações, taludes, etc. A ocorrência do fenômeno da liquefação, em materiais suscetíveis como areias fofas saturadas, representa um tipo de resposta desastrosa de solos. O termo liquefação tem sido empregado para descrever uma variedade de fenômenos no qual tem em comum o desenvolvimento de altas poropressões em materiais saturados sem coesão devido a carregamentos monotônicos , transientes ou ciclios. A previsão da liquefação depende de uma adequada análise do comportamento não-drenado do material, em termos do incremento de poropressões e da perda da rigidez da mistura sólido-fluido, durante e após o período de movimento. O estabelecimento das equações governantes é essencial para elaboração de um modelo matemático realista para descrever o comportamento físico deste fenômeno. As equações a srem consideradas são: equação de movimento da fase sólida, a equação do movimento da mistura sólido-fluido , a equação de continuidade da fase fluida, as equações de acoplamento das fases e as equações constitutivas desses materiais. Nesta tese a resposta dinâmica do solo foi investigada numericamente mediante a técnica dos elementos finitos. A discretização espacial das equações governantes foi feita através de método de Galerkin e a discretização temporal pelo método de Newmark Generalizado. Um modelo constitutivo elasto-plástico foi considerado para descrever o comportamento mecânico da fase sólida, desenvolvido a partir de conceitos da generealização da teoria da plasticidade, que apresenta algumas vantagens em relação aos outros modelos baseados na teoria da plasticidade clássica. A implementação computacional foi escrito em fortran 90. Exemplos numéricos analisados nesta tese comprovam tanto a eficiência do modelo constitutivo na predição do comportamento do solo sobre liquefação como a confiabilidade do programa computacional elaborado nesta pesquisa, em termos da rapidez de processamento e da boa precisão dos resultados, quando comparados com soluções analíticas e outros valores numéricos obtidos por vários autores e diferentes modelos constitutivos. / [en] Changes in soil properties due to the action of dynamic loads are responsible for significant damage of geo-structures such as dams, retaining structures,building foundations, slopes, etc. The occurrence of liquefaction phenomena in susceptible materials, such as loose saturated, represents a type of disastrous response of soil, the term liquefaction has been used to refer to a group of phenomena wich have in common the development of high pore pressures in saturated cohesionless mterial due to monotonic, transient, or cyclic loads. The prediction of soil liquefaction depends of an adequate analysis of the behavior of undrained materials, in terms of increase of pore water pressure and weakening of the solid-fluid mixture, during and after the periodic motion. The establishment of the governing equations is essential to provide a realistic mathematical model to describe the physical behavior of this phenomenon. The system of equations to be considered are: the equilibrium equation of the solid phase, the equilibrium equation of the solid-fluid mixture, the conservation mass of the fluid phase, the coupling equation of phases, and the conservation equations of materials. In this thesis the soil dynamic response was numerically investigated by the finite element method. To obtain the spatial discretization in time was the Generalized Newmark method. An elastic-plastic constitutive model was used to describe the mechanical behavior of the solid phase. This model was developed in the framework of the generalized theory of plasticity, wich has some advantages when compared with other models based on the classical plasticity theory. The computacional implementation was written in fortran 90. Numerical examples considered in this thesis demonstrate the efficiency of the constitutive model to simulated the predicted behavior of soil under liquefaction as well as the reliability of the software developed in this research, in terms of computational effort and good accuracy of the results, when compared with some analytical solutions and other numerical values obtained by various authors and different constitutive models.

[pt] Apresenta-se uma visão gerencial para o procedimento de projeto de suspensões de veículos terrestres. São descritos, em linhas gerais, os principais aspectos técnicos relativos ao projeto deste sub-sistema veicular, e tratados, com detalhes, os tópicos fundamentais para a sua administração. Discute-se, entre outras, as etapas de especificação do projeto, quesitos necessários, normas a serem aplicadas, e as etapas a cumprir, segundo a visão do gerente administrativo do projeto, responsável pela organização do grupo de técnicos que irá desenvolvê-lo. / [en] It is shown a managerial vision of the ground vehicles suspension design procedure. Are described, in general lines, the main technical aspects related to the design of this vehicular sub-system, and treats, with details, the fundamental topics for its administration. It is discussed, beside others, the design specification stages, necessary requirements, norms to be applied, and the stages to accomplish, according to the vision of the administrative design manager, responsible for the organization of the technicians' group that will develop it.

[pt] Dentro do campo dos controladores passivos, um dos dispositivos usados pelas suas propriedades de amortecimento são os Apoios Elásticos, que constituem uma solução econômica e efetiva na supressão ou atenuação das vibrações em sistemas dinâmicos com problemas de ressonância ou instabilidade, freqüentemente pela falta de amortecimento suficiente. Este trabalho envolve o estudo de um rotor horizontal com apoios elásticos (silicone), adaptado a partir de um rotor existente, o estudo de diversos efeitos sobre a sua operação, a medição de seu movimento, a identificação dos parâmetros do problema, a medição e validação a partir de resultados simulados em um modelo numérico. Os fenômenos incluídos no estudo são o efeito giroscópio (rotor descentrado com respeito do vão), desbalanceamento do rotor e empenamento do eixo. Neste trabalho os parâmetros do sistema foram determinados usando técnicas de identificação, análise modal e otimização não linear devido à anisotropia do sistema. / [en] Viscoelastic Passive Controllers are an important field of technological research due to the development of new materials and design techniques. Damping properties allow an easy retrofit of existing machines with excessive vibration problems, developing Elastomeric Bearing Supports. They are an economic and effective solution in the suppression or attenuation of vibrations in dynamic systems suffering from instability or resonance problems, which often lack of sufficient damping. This work involves the study of an horizontal rotor with elastomeric bearing supports, adapted of another one, the study of several effects on his operation, the measurement of his motion, the identification of the problems parameters, the measurement and validation from the simulated results in a numeric model. Phenomena included in the study are the gyroscopic effect (rotor out of the middle), rotor unbalance and shaft bow. In this work the systems parameters were determined using identification, modal analysis and nonlinear optimization techniques due to the anisotropy of the system.

[pt] O controle de empuxo vetorial (TVC) é usado para o controle de atitude de foguetes aeroespaciais. No caso de propulsão usando combustível líquido, tradicionalmente o bocal é conectado ao corpo do foguete através de uma junta cardânica. Dois atuadores hidráulicos são colocados ao redor do bocal para controlar sua orientação. Nesta tese, o TVC é tratado como uma plataforma robótica de base fixa. Ao invés de usar servo-válvulas comerciais para controlar os atuadores, uma nova válvula de controle é proposta. Primeiro uma plataforma cardânica é considerada com transdutores de posição angular medindo o deslocamento da cruzeta da junta. Em seguida, uma nova configuração da plataforma é proposta substituindo o cardan por uma junta homocinética. Neste caso, a realimentação da posição da plataforma é feito usando um estimador de atitude em tempo real. Este estimador é um filtro complementar baseado em matrizes de orientação que coleta dados de uma central inercial (IMU). A modelagem do sistema começa com a cinemática. Na sequência, a modelagem dinâmica utiliza a formulação de Newton-Euler para obter a equação de movimento. A modelagem do sistema hidráulico é apresentada com o modelo da nova válvula de controle e do atuador. Inicialmente, um controlador puramente proporcional é proposto. Durante a validação experimental é mostrado que devido as características do sistema de atuação, mesmo este simples controlador é preciso e confiável. Em seguida é demonstrado um método para avaliar outras estratégias de controle. A comparação entre a plataforma cardânica e homocinética mostra que, nas condições analisadas, ambas têm um comportamento dinâmico similar. Nas duas configurações da plataforma o sistema se mostrou preciso e confiável. / [en] Thrust Vector Control (TVC) is used for the attitude control of spacecrafts. In the case of liquid-propellant fuel, the nozzle is traditionally connected to the rocket frame through a gimbal. Two hydraulic actuators are placed around the nozzle to control its orientation. In this Thesis, TVC is treated as a fixed base robotic platform. Instead of using commercial servo-valves to control the actuators, a novel control valve is proposed. First a gimbaled platform is considered with two angular position transducers to measure the angular displacement of the joint crosshead. Then, a homokinetic platform configuration is proposed replacing the gimbal by a constant velocity joint. In this case, the platform position feedback is done using a real-time attitude estimator. The estimator is a complementary filter based on orientation matrices that collects data from an inertial measurement unit (IMU). The modeling of the system begins with kinematics. Then, the dynamic modeling uses the Newton-Euler formulation to obtain the equation of motion. The modeling of the electro-hydraulic system is presented with the model of the novel control valve and the linear actuator. Initially, a full proportional controller is proposed. During the experimental validation it is shown that due to the characteristics of the actuation system, even this simple controller is accurate and reliable. Thereafter, method is demonstrated to evaluate novel control strategies. The comparison between the gimbaled and homokinetic platform shows that, under the analyzed conditions, they have a similar dynamic behavior. In both platform configurations the system is accurate and reliable.

[pt] Esta tese contém três ensaios sobre dinâmica de firmas. O primeiro ensaio avalia os efeitos de choques de oferta e demanda sobre a dinâmica de firmas e seleção no Brasil. Exploramos o fato de que o crescimento recente da China não apenas aumentou o nível de competição via importações, mas também aumentou a demanda por exportações de bens primários, um fator especialmente relevante para países em desenvolvimento. Nossos resultados mostram que firmas afetadas pelo aumento de competição proveniente de importações chinesas apresentam um aumento na probabilidade de sair do mercado, enquanto firmas em indústrias beneficiadas pela demanda por exportações para a China têm uma menor probabilidade de saída. Em ambos os casos, esses efeitos estão concentrados em firmas com um menor número de trabalhadores. O segundo ensaio descreve a relação entre a má alocação de energia e a má alocação de recursos no setor de manufaturas brasileiro, e quantifica em que medida distorções que afetam o uso eficiente de energia resultam em perdas de produto agregado. Nós encontramos que as duas medidas de má alocação são positivamente relacionadas nos setores, sugerindo que a energia é um importante componente da eficiência alocativa de recursos. Nós mostramos que a realocação de recursos entre firmas de um mesmo setor levaria a ganhos agregados significativos. Entretanto, distorções de capital são responsáveis pela maior parte dos ganhos potenciais pela realocação de recursos. O terceiro ensaio compara a dinâmica do ciclo de vida em manufaturas e serviços e encontra que o crescimento ao longo do ciclo de vida é menor para firmas do setor de serviços, mesmo controlando pelo seu tamanho inicial. Nós mostramos que esse menor crescimento ocorre devido ao padrão de seleção e à fraca relação existente entre produtividade e tamanho das firmas em serviços. Finalmente, nós investigamos o papel de duas possíveis explicações para os resultados encontrados: distorções relacionadas ao ciclo de vida e custos de monitoramento. / [en] This thesis consists of three essays on firm dynamics. The first essay evaluates the effects of supply and demand shocks on firm dynamics and selection in Brazil. We explore the fact that China’s recent growth has led not only to an increase in import competition, but also to higher export demand for commodities, which is especially relevant in developing countries. We find that firms facing greater competition from Chinese imports suffer from an increase in exit probability, while firms in industries benefiting from increased export demand have lower probability of exit. In both cases, these effects are concentrated among smaller firms. In the second article, we describe the relationship between energy misallocation and resource misallocation across manufacturing industries in Brazil, and quantify the extent to which distortions affecting energy use result in output losses at the aggregate level. We find that these two measures of misallocation are positively related across industries, which suggests that energy is an important component of resource allocation efficiency. We show that reallocating resources between firms would result in substantial aggregate output gains. However, capital distortions account for most of the potential gains in manufacturing from reallocating resources between firms. The third essay compares firm life cycle dynamics in manufacturing and services, and finds that life-cycle growth is slower for service firms, even when controlling for initial size. We show that this result arises because of the selection pattern and weaker relationship between productivity and size in service industries. Finally, we assess the role of two potential explanations for these results: age-related distortions and monitoring costs.

[pt] A determinação de uma trajetória ótima não é uma tarefa simples, uma vez que ela é diretamente dependente dos limites de aceleração suportada por cada veículo. Essa pesquisa aborda um método de otimização baseado em algoritmos genéticos que identifica a trajetória que um carro deve percorrer para completar uma pista pré-definida no menor tempo. Considerando um modelo veicular de Partícula Orientada, o método otimiza os perfis de aceleração que levam o veículo a percorrer a trajetória de menor tempo. Adicionalmente, projeta-se um controlador fuzzy para emular o comportamento de um ser humano na direção do veículo ao longo da trajetória ótima. Para alimentar o controlador, foram testados dois métodos de geração de erro: o Erro Presente da Trajetória e o Erro Futuro da Trajetória (FBTE), que é a medida de posição do carro quanto a sua tendência de movimento. Resultados obtidos com controladores clássicos, como o PDD, são confrontados com os fornecidos pelo controlador fuzzy alimentado pelo procedimento de geração de Erro Futuro de Trajetória (FBTE). / [en] The definition of the minimum time trajectory in a track is not obvious, since it is directly dependent on the acceleration limits that the vehicle can withstand. This paper presents an optimization method based on Genetic Algorithms that identifies the path that a car must follow in order to complete a given circuit in minimum time. By considering an Oriented Particle model, the method optimizes the acceleration profiles that drive the vehicle along the trajectory in minimum time. In addition, a fuzzy controller is designed to emulate the behavior of a human driver controlling a high speed car along the optimized trajectory. Two different error generation procedures were tested as controller inputs: the Present Trajectory Error and the Future-based Trajectory Error (FBTE), which gives information on the car’s tendency of movement. Results obtained with other controllers in the same application, such as the PDD, are compared to those provided by the fuzzy controller fed by the FBTE procedure.

[pt] Este trabalho aplica três métodos básicos de malhas auto-adaptativas em problemas de mecânica dos fluidos e transmissão de calor. Um método unidimensional é aplicado a problemas descritos pela equação generalizada de Burgers. Um método bidimensional é aplicado ao problema da. cavidade com tampa móvel para números de Reynolds Iguais a 100, 500 e 1000 e ao problema de convecção natural em uma cavidade para números de Rayleigh iguais a 10 elevado a 3,10 elevado a 4,10 elevado a 5 e 10 elevado a 6. Um terceiro método, fundamentado em um sistema de equações elípticas é proposto e aplicado ao problema de convecção natural na cavidade, para análise de desempenho. Procura-se definir, para cada aplicação, critérios como: a) variável indicada para dirigir a adaptação; b) especiacação da intensidade desejada na adaptação e c) justificação para o uso de malhas auto-adaptativas no problema em estudo. Observa-se, nos problemas unidimensionais, que quando o problema tem uma região única de intensos gradientes, a adaptação da malha funciona convenientemente, especialmente na presença de um termo de fonte. No problema da cavidade com, tampa móvel, verifica-se que a função de corrente deve dirigir a adaptação. Para Re igual a 1000 o método empregado consegue melhorar a descrição geral do problema. O problema de convecção natural apresenta melhores resultados para Ra: 10 elevado a 6 e com a adaptação dirigida pelo gradiente da temperatura que está diretamente relacionado ao termo de fonte do problema (forças de empuxo). São consideradas as características geométricas de ortogonalidade e suavidade das malhas adaptadas. Para isso, compara-se malhas obtidas pelos métodos bidimensionais implementados. Observa-se que o método fundamentado no sistema elíptíco consegue concentrar a malha nas regiões de intensos gradientes, sem danos excessivos naquelas características. / [en] The present work is concerned to the study and implementation of self-adaptive grids to tluid mechanics and heat transfer problems. Initially, a ID method is applied to simple problems described by the generalized Burger-s equation, After that, a 2-D method is applied to the wall driven Cavity with Reynolds numbers equal to 100, 500 and 1000 and to the natural convection problem inside cavltles with Rayleigh numbers equal to 10 to the third power, 10 to the fourth power, 10 to the fifth power and 10 to the sixth power. A third method, based on a system of elliptic equations is proposed and applied to this natural convection problem, for error assessmet. For each application, one seeks to define criteria such as: a) most suitable dependent variable to drive the adaptive technique; b) intensity of the adaptive technique and c) need for the use of the technique at all. From the 1-D problems, it may be concluded that in presence of single regions of strong gradients, the technidne works line, specially shenever a source term is present. From the wall driven cavity problem, it appears that the stream function is most suitable to be the driving force for the adaptation. For instance, for Reynolds number on the order of 1000, the used method works well. The natural convectlon problem indicates best results for high Rayleigh numbers, say 10 to the sixth power, and with the use of the temperature gradient that is straight related to the source term, that is, buoyancy in the case. In the present work, geometric characteristics such as ortogonality and smooth are considered. Among the three methods, studied, it seems that the elliptic one, proposed here, is able to concentrate grid points Whenever needed, without severe penalization to those characteristics.

[pt] Correntes de turbidez ocorrem tanto na natureza como em situações criadas pelo homem. Segundo alguns pesquisadores, grande parte das reservas de petróleo conhecidas no mundo estão armazenadas nos reservatórios de hidrocarbonetos formados a partir dos sistemas turbidíticos. Tendo em vista a importância dessas correntes, este trabalho pretende propor um modelo numérico consistente e com baixo custo computacional, capaz de auxiliar na previsão de sedimentação de bacias nos processos de modelagem estratigráfica no contexto do simulador numérico de sedimentação tridimensional com ênfase nos processos deposicionais em ambientes de talude plataforma e bacia, chamado Steno, e validar o modelo numérico proposto através de simulações físicas conduzidas no Instituto de Pesquisas Hidráulicas da Universidade Federal do Rio Grande do Sul. O modelo numérico proposto é baseado nas equações de Navier-Stokes, que são resolvidas tirando-se a média na profundidade do perfil de velocidades característico. Além disso, uma aplicação gráfica tridimensional para simulação numérica de correntes de turbidez, chamada Turb3D, foi desenvolvida com base nas equações propostas. O programa possui uma interface gráfica amigável para o usuário no que diz respeito à entrada de dados, solução e visualização dos resultados. Os resultados obtidos mostraram que apesar do modelo numérico possuir um bom desempenho computacional, se comparado com os complexos modelos numéricos existentes, não foi possível reproduzir com precisão a evolução das correntes de turbidez e seus depósitos. Contudo, acredita-se que a abordagem dada ao problema de simulação numérica de correntes de turbidez apresentada neste trabalho pode fornecer bons resultados. Entretanto, é fundamental que seja dada continuidade a pesquisa, de modo que ajustes e melhorias ao modelo proposto possam ser realizadas. / [en] Turbidity currents occur in both natural and man-made situations. In agreement with some researchers, most of the world’s oil reserves are stored in hydrocarbon reservoir built by turbidity systems. Because of the importance of these currents, this work proposes a consistent and efficient numerical method for simulations of turbidity currents for basin sedimentations predictions in the stratigraphic modelling process, which will be incorporated to the numerical simulator of sedimentation 3D called Steno. In order to validate the proposed numerical model experiments were conducted at the Hydraulics Research Institute of Federal University of Rio Grande do Sul. The proposed numerical model is based on Navier-Stokes equations that are solved in the depth-average. Moreover, a three-dimensional graphic application for numerical simulations of turbidity currents called Turb3D was developed. The application user interface provides a common, user-friendly, graphical environment for pre-processing, solution and post-processing. Despite the good computational performance achieved by using this approach, the method presented did not reproduce accurately the evolution of turbidity current and their deposits. However, it is believed that the approach to the numeric simulation of the turbidity current problem given in this work can provide better results, although this research should continue and improvements should be made.

[pt] Motivado pela riqueza de fenômenos produzidos pelos seres vivos, este trabalho busca estudar a formação de padrões espaciais de populações biológicas. De um ponto de vista mesoscópico, definimos os processos básicos que podem ocorrer na dinâmica, construindo uma equação diferencial parcial para a evolução da distribuição da população. Essa equação incorpora duas generalizações de um modelo pre-existente para a dinâmica de um espécie, que leva em conta interações de longo alcance (não locais). A primeira generalização consiste em considerar que a difusão é não linear, isto é, é afetada pela densidade local de tal modo que o coeficiente de difusão segue uma lei de potência. Por outro lado, visto a alta complexidade envolvida na natureza dos parâmetros do modelo, introduzimos como segunda generalização parâmetros que flutuam no tempo. Idealizamos estas flutuações como um ruído descorrelacionado temporalmente e que obedece uma distribuição gaussiana (ruído branco). Para estudar o modelo resultante, utilizamos uma abordagem analítica e numérica. As ferramentas analíticas se baseiam na linearização da equação de evolução e portanto são aproximadas. Todavia, complementadas com resultados numéricos, conseguimos extrair conclusões relevantes. A não localidade das interações induz a formação de padrões. O alcance dessas interações é o que determina o modo dominante presente nos padrões. Assim, para valores dos parâmetros acima de um limiar crítico, emergem padrões. Analiticamente, mostramos que, mesmo abaixo desse limiar, as flutuações nos parâmetros podem induzir a aparição de ordem espacial. Os efeitos da difusão não-linear são captados superficialmente pela análise linear. Numericamente, mostraremos que sua presença modifica a forma dos padrões. Observamos, especialmente, a existência de uma transição quando alternamos entre o caso em que a difusão é facilitada por altas densidades e o caso oposto. Para o primeiro caso, verificamos que os padrões se tornam fragmentados, ou seja, a população é agora composta de sub-grupos desconectados. / [en] Motivated by the richness of phenomena produced by living beings, this work aims to study the formation of spatial patterns in biological populations. From the mesoscopic point of view, we define the basic processes that may occur in the dynamics, building a partial differential equation for the evolution of the population distribution. This equation incorporates two generalizations of a pre-existing model for the dynamics of one species, which takes into account long-range (nonlocal) interactions. The first generalization is to consider that diffusion is nonlinear, i.e., it is affected by the local density such that the diffusion coeficient follows a power law. On the other hand, because of the high complexity involved in the nature of model parameters, we introduced as a second generalization time-fluctuating parameters. We idealize these fluctuations as Gaussian temporally uncorrelated (white) noises. To study the resulting model, we use an analytical and numerical approach. Analytical tools are based on the linearization of the evolution equation and are therefore approximate. However, as evidenced by numerical results, we draw important conclusions. The nonlocal feature of the interaction is the main mechanism which induces pattern formation. We show that the extent of these interactions is what characterizes the dominant mode. Thus, for parameter values above a critical threshold patterns emerge. Analytically, we also show that even below this threshold, fluctuations in the parameters can induce the appearance of spatial order. The effects of nonlinear diffusion are only superficially captured by the linear analysis. Numerically, we show that their presence modifies the patterns shape. We mainly observed the existence of a qualitative difference between the cases when diffusion is facilitated or not by high densities. In the first case, we note that the patterns become fragmented, that is, population becomes composed of disconnected clusters.

[pt] O estudo da interação de dibenzo[a,h]antraceno (DBahA) e de seu metabólito com modelos de membrana celular e surfactante pulmonar foi realizado através de dinâmica molecular. Os modelos de membrana celular e de surfactante pulmonar são geralmente misturas de dipalmitoil fosfatidilcolina (DPPC), dipalmitoil fosfatidilglicerol (DPPG), e colesterol. No caso do modelo de surfactante pulmonar pode ser incluido as proteínas surfactantes (SP-A, SP-B, SP-C e SP-D). Neste projeto, o dibenzo[a,h]antraceno (DBahA) foi simulado com o DPPC sozinho e com uma mistura 32/32/1 de DPPC/DPPG/Colesterol. DBahA é encontrado nos gases de exaustão de veículos automotores (especialmente os movidos a diesel), na fumaça do cigarro e da madeira, além de alimentos grelhados na brasa. Ele é capaz de ser metabolizado pelo citocromo P450 e seu metabólito interage com o DNA, sendo então mutagênico e altamente carcinogênico. Os principais resultados mostram que o DBahA se difunde para o interior dos modelos e forma aglomerados. Quando o DBahA está em concentração elevada na parte exterior dos modelos, este não consegue se difundir facilmente para o interior dos modelos na escala de tempo simulado e forma aglomerados na interface água/modelo. O metabólito age similarmente, no entanto prefere ficar mais próximo da cabeça polar dos modelos. / [en] The study of the interaction of dibenz[a,h]anthracene (DBahA) and its metabolite with cell membrane and pulmonary surfactant models was performed by molecular dynamics. The cell membrane and pulmonary surfactant models usually are mixtures of dipalmitoyl phosphatidylcholine (DPPC), dipalmitoyl phosphatidylglycerol (DPPG), and cholesterol. In the case of pulmonary surfactant, the models may include surfactant proteins (SP-A, SP-B, SP-C and SP-D). In this project, the DBahA was simulated with DPPC and with a 32/32/1 mixture of DPPC/DPPG/Cholesterol. DBahA is found in automotive vehicles (especially diesel vehicles), in cigarette and wood smoke, and grilled food. The DBahA molecule is metabolized by cytochrome P450 and its metabolite interacts with DNA, being mutagenic and highly carcinogenic. The results show that the DBahA diffuses into the interior of the models forming clusters. In the simulated time scale, when the DBahA is in high concentration in the outer part of the models, it may not spread easily to the inner side of the models because it forms clusters in the water/model interface. The metabolite acts similarly, but prefers to stay closer to the polar head of the models.