Direct Simulation Monte Carlo Method Applied to Aerothermodynamics

An axisymmetric supersonic flow of rarefied gas past a finite cylinder was calculated applying the direct simulation Monte Carlo method. The drag force, the coefficients of pressure, of skin friction, and of heat transfer, the fields of density, of temperature, and of velocity were calculated as function of the Reynolds number for a fixed Mach number. The variation of the Reynolds number is related to the variation of the Knudsen number, which characterizes the gas rarefaction. The present results show that all quantities in the transition regime (Knudsen number is about the unity) are significantly different from those in the hydrodynamic regime, when the Knudsen number is small.

Termodinâmica estatística de líquidos com o método de Monte Carlo. I. Metodologia

Statistical mechanics Monte Carlo simulation is reviewed as a formalism to study thermodynamic properties of liquids. Considering the importance of free energy changes in chemical processes, the thermodynamic perturbation theory implemented in the Monte Carlo method is discussed. The representation of molecular interaction by the Lennard-Jones and Coulomb potential functions is also discussed. Charges derived from quantum molecular electrostatic potential are also discussed as an useful methodology to generate an adequate set of partial charges to be used in liquid simulation.

Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Thermodynamic properties and radial distribution functions for liquid chloroform were calculated using the Monte Carlo method implemented with Metropolis algorithm in the NpT ensemble at 298 K and 1 atm. A five site model was developed to represent the chloroform molecules. A force field composed by Lennard-Jones and Coulomb potential functions was used to calculate the intermolecular energy. The partial charges needed to represent the Coulombic interactions were obtained from quantum chemical ab initio calculations. The Lennard-Jones parameters were adjusted to reproduce experimental values for density and enthalpy of vaporization for pure liquid. All thermodynamic results are in excelent agreement with experimental data. The correlation functions calculated are in good accordance with theoretical results avaliable in the literature. The free energy for solvating one chloroform molecule into its own liquid at 298 K and 1 atm was computed as an additional test of the potential model. The result obtained compares well with the experimental value. The medium effects on cis/trans convertion of a hypotetical solute in water TIP4P and chloroform solvents were also accomplished. The results obtained from this investigation are in agreement with estimates of the continuous theory of solvation.

Simulação do equilíbrio: o método de Monte Carlo

We make several simulations using the Monte Carlo method in order to obtain the chemical equilibrium for several first-order reactions and one second-order reaction. We study several direct, reverse and consecutive reactions. These simulations show the fluctuations and relaxation time and help to understand the solution of the corresponding differential equations of chemical kinetics. This work was done in an undergraduate physical chemistry course at UNIFIEO.

Simulação Monte Carlo no ensino de luminescência e cinética de decaimento de estados excitados

A software based in the Monte Carlo method have been developed aiming the teaching of important cases of mechanisms found in luminescence and in excited states decay kinetics, including: multiple decays, consecutive decays and coupled systems decays. The Monte Carlo Method allows the student to easily simulate and visualize the luminescence mechanisms, focusing on the probabilities of the related steps. The software CINESTEX was written for FreeBASIC compiler; it assumes first-order kinetics and any number of excited states, where the pathways are allowed with probabilities assigned by the user.

Simulação Monte Carlo de mecanismo de transferência de energia de excitação eletrônica: modelo de Perrin para a supressão estática da luminescência

A software based in the Monte Carlo method has been developed aiming the teaching of the Perrin´s model for static luminescence quenching. This software allows the student to easily simulate the luminescence decays of emissive molecules in the presence of quenching ones. The software named PERRIN was written for FreeBASIC compiler and it can be applied for systems where the molecules remain static during its excited state lifetime. The good agreement found between the simulations and the expected theoretical results shows that it can be used for the luminescence and excited states decay kinetic teaching.

Uma abordagem simplificada do método Monte Carlo Quântico: da solução de integrais ao problema da distribuição eletrônica

The paper presents an introductory and general discussion on the quantum Monte Carlo methods, some fundamental algorithms, concepts and applicability. In order to introduce the quantum Monte Carlo method, preliminary concepts associated with Monte Carlo techniques are discussed.

Computational modeling for irrigated agriculture planning. Part II: risk analysis

Techniques of evaluation of risks coming from inherent uncertainties to the agricultural activity should accompany planning studies. The risk analysis should be carried out by risk simulation using techniques as the Monte Carlo method. This study was carried out to develop a computer program so-called P-RISCO for the application of risky simulations on linear programming models, to apply to a case study, as well to test the results comparatively to the @RISK program. In the risk analysis it was observed that the average of the output variable total net present value, U, was considerably lower than the maximum U value obtained from the linear programming model. It was also verified that the enterprise will be front to expressive risk of shortage of water in the month of April, what doesn't happen for the cropping pattern obtained by the minimization of the irrigation requirement in the months of April in the four years. The scenario analysis indicated that the sale price of the passion fruit crop exercises expressive influence on the financial performance of the enterprise. In the comparative analysis it was verified the equivalence of P-RISCO and @RISK programs in the execution of the risk simulation for the considered scenario.

Neutron tomography using projection data obtained by Monte Carlo simulation for nondestructive evaluation

This work present the application of a computer package for generating of projection data for neutron computerized tomography, and in second part, discusses an application of neutron tomography, using the projection data obtained by Monte Carlo technique, for the detection and localization of light materials such as those containing hydrogen, concealed by heavy materials such as iron and lead. For tomographic reconstructions of the samples simulated use was made of only six equal projection angles distributed between 0º and 180º, with reconstruction making use of an algorithm (ARIEM), based on the principle of maximum entropy. With the neutron tomography it was possible to detect and locate polyethylene and water hidden by lead and iron (with 1cm-thick). Thus, it is demonstrated that thermal neutrons tomography is a viable test method which can provide important interior information about test components, so, extremely useful in routine industrial applications.

Monte Carlo derivation of filtered tungsten anode X-ray spectra for dose computation in digital mammography

Abstract Objective: Derive filtered tungsten X-ray spectra used in digital mammography systems by means of Monte Carlo simulations. Materials and Methods: Filtered spectra for rhodium filter were obtained for tube potentials between 26 and 32 kV. The half-value layer (HVL) of simulated filtered spectra were compared with those obtained experimentally with a solid state detector Unfors model 8202031-H Xi R/F & MAM Detector Platinum and 8201023-C Xi Base unit Platinum Plus w mAs in a Hologic Selenia Dimensions system using a direct radiography mode. Results: Calculated HVL values showed good agreement as compared with those obtained experimentally. The greatest relative difference between the Monte Carlo calculated HVL values and experimental HVL values was 4%. Conclusion: The results show that the filtered tungsten anode X-ray spectra and the EGSnrc Monte Carlo code can be used for mean glandular dose determination in mammography.

Monte Carlo thermodynamic and structural properties of the TIP4P water model: dependence on the computational conditions

Classical Monte Carlo simulations were carried out on the NPT ensemble at 25°C and 1 atm, aiming to investigate the ability of the TIP4P water model [Jorgensen, Chandrasekhar, Madura, Impey and Klein; J. Chem. Phys., 79 (1983) 926] to reproduce the newest structural picture of liquid water. The results were compared with recent neutron diffraction data [Soper; Bruni and Ricci; J. Chem. Phys., 106 (1997) 247]. The influence of the computational conditions on the thermodynamic and structural results obtained with this model was also analyzed. The findings were compared with the original ones from Jorgensen et al [above-cited reference plus Mol. Phys., 56 (1985) 1381]. It is notice that the thermodynamic results are dependent on the boundary conditions used, whereas the usual radial distribution functions g(O/O(r)) and g(O/H(r)) do not depend on them.

Estudo de reações químicas homogêneas via método de Monte Carlo

By using the Monte Carlo simulation platform with probabilistic mathematical functions of the Boltzmann type, , having activation energy and temperature as parameters, it was possible to assess important dynamic aspects of homogeneous chemical reactions of the types A → B and A B. The protocol proved a useful tool in work with the basic concepts of Kinetics and Thermodynamics allowing its application both in class activities and for assisting experimental procedures.

Influência das características mamárias na cintilografia miocárdica pelo método Monte Carlo

FUNDAMENTO: Pela redução da especificidade associada à perda de informação, a influência da atenuação das mamas é de fundamental importância em estudos de perfusão do miocárdio. Entretanto, apesar de vários estudos terem sido realizados ao longo dos últimos anos, pouco se tem evoluído para determinar com acurácia a influência das características das mamas sobre a qualidade da cintilografia miocárdica, evitando exposições adicionais de radiação às pacientes. OBJETIVO: O objetivo deste estudo é quantificar a atenuação de fótons pelas mamas, em estudos de perfusão do miocárdio com 99mTc, de acordo com diferentes tamanhos e composições. MÉTODOS: Cada mama foi assumida como sendo um cubo composto de tecido adiposo e fibroglandular. Os dados referentes aos fótons de 99mTc foram analisados em um modelo de Monte Carlo. Variamos a espessura e a composição das mamas e analisamos as interferências na atenuação. Foi empregado o software EGS 4 para as simulações. RESULTADOS: Fixando a espessura de uma mama, a variação da sua composição acarreta um acréscimo máximo de 2,3% no número de fótons atenuados. Em contrapartida, mantendo-se uma composição do tecido mamário fixa, a diferença na atenuação de fótons foi de 45,0%, sendo em média de 6,0% para cada acréscimo de centímetro na espessura da mama. CONCLUSÃO: A simulação por Monte Carlo demonstrou que a influência das espessuras das mamas na atenuação de fótons em cintilografias do miocárdio com 99mTc é muito maior do que a influência das suas composições.

CITRICULTURE ECONOMIC AND FINANCIAL EVALUATION UNDER CONDITIONS OF UNCERTAINTY

ABSTRACT The citriculture consists in several environmental risks, as weather changes and pests, and also consists in considerable financial risk, mainly due to the period ofreturn on the initial investment. This study was motivated by the need to assess the risks of a business activity such as citriculture. Our objective was to build a stochastic simulation model to achieve the economic and financial analysis of an orange producer in the Midwest region of the state of Sao Paulo, under conditions of uncertainty. The parameters used were the Net Present Value (NPV), the Modified Internal Rate of Return(MIRR), and the Discounted Payback. To evaluate the risk conditions we built a probabilistic model of pseudorandom numbers generated with Monte Carlo method. The results showed that the activity analyzed provides a risk of 42.8% to reach a NPV negative; however, the yield assessed by MIRR was 7.7%, higher than the yield from the reapplication of the positive cash flows. The financial investment pays itself after the fourteenth year of activity.

Dosimetria em folículos tireoidianos devido aos elétrons de baixa energia do iodo usando o método Monte Carlo

OBJETIVO: Avaliar a dose absorvida em folículos tireoidianos devido aos elétrons de baixa energia, como os elétrons Auger e os de conversão interna, além das partículas beta, para os radioisótopos de iodo (131I, 132I, 133I, 134I e 135I) usando o método Monte Carlo. MATERIAIS E MÉTODOS: O cálculo da dose foi feito ao nível folicular, simulando elétrons Auger, conversão interna e partículas beta, com o código MCNP4C. Os folículos (colóide e células foliculares) foram modelados como esferas, com diâmetros do colóide variando de 30 a 500 mm. A densidade considerada para os folículos foi a da água (1,0 g.cm-³). RESULTADOS: Considerando partículas de baixa energia, o percentual de contribuição do 131I na dose total absorvida pelo colóide é de aproximadamente 25%, enquanto os isótopos de meia-vida física curta apresentaram contribuição de 75%. Para as células foliculares, esse percentual é ainda maior, chegando a 87% para os iodos de meia-vida curta e 13% para o 131I. CONCLUSÃO: Com base nos resultados obtidos, pode-se mostrar a importância de se considerar partículas de baixa energia na contribuição para a dose total absorvida ao nível folicular (colóide e células foliculares) devido aos radioisótopos de iodo (131I, 132I, 133I, 134I e 135I).