Aujourd’hui le développement de nouveaux matériaux polymères ayant de bonnes propriétés repose de plus en plus sur des procédés de mélange ou de compoundage de polymères au lieu de recourir à la synthèse de nouvelles molécules. L’action du mélange peut fortement influer sur la morphologie des matériaux polymères multi-constituants. Les extrudeuses bi-vis (TSE) sont souvent utilisées comme mélangeurs/réacteurs pour des procédés de mélange, de compoundage et d’extrusion réactive. Cependant, l’étude sur la qualité du mélange dans les TSE demeure un grand défi en raison de la complexité géométrique et du caractère transitoire de l’écoulement. Cette thèse a pour objet de développer un nouvel instrument en line pour mesurer en temps réel la distribution des temps de séjour (DTS) qui caractérise la performance du mélange axial et la capacité de convoyage de différents types d’éléments de vis basées sur l’analyse de l’écoulement transitoire et l’évaluation systématique de la théorie de mélange dans les TSE. Le mélange distributif des polymères fondus est caractérisé par la génération de l’aire des interfaces, un paramètre difficile à mesurer expérimentalement. Alors on fait appel à des simulations numériques de type CFD / The development of new materials with improved properties seems to rely nowadays more on blending and compounding than on the synthesis of chemically new polymers. Mixing may have a great effect on the morphology and structure of multi-component polymer materials. Twin-screw extruders (TSE) are widely used as mixers/reactors for blending, compounding, and reactive processing. This work aimed at developing a new instrument to measure in real time the residence time distribution (RTD) which characterizes the axial mixing and transport abilities of different screw elements based on the analysis of the transient flow pattern and systematic evaluation of mixing theory in TSE. Distributive mixing of polymer melts is characterized by the generation of interfacial area, which is experimentally much more difficult to measure. This 3D numerical simulation based on CFD is adopted

Aujourd’hui le développement de nouveaux matériaux polymères ayant de bonnes propriétés repose de plus en plus sur des procédés de mélange ou de compoundage de polymères au lieu de recourir à la synthèse de nouvelles molécules. L’action du mélange peut fortement influer sur la morphologie des matériaux polymères multi-constituants. Les extrudeuses bi-vis (TSE) sont souvent utilisées comme mélangeurs/réacteurs pour des procédés de mélange, de compoundage et d’extrusion réactive. Cependant, l’étude sur la qualité du mélange dans les TSE demeure un grand défi en raison de la complexité géométrique et du caractère transitoire de l’écoulement. Cette thèse a pour objet de développer un nouvel instrument en line pour mesurer en temps réel la distribution des temps de séjour (DTS) qui caractérise la performance du mélange axial et la capacité de convoyage de différents types d’éléments de vis basées sur l’analyse de l’écoulement transitoire et l’évaluation systématique de la théorie de mélange dans les TSE. Le mélange distributif des polymères fondus est caractérisé par la génération de l’aire des interfaces, un paramètre difficile à mesurer expérimentalement. Alors on fait appel à des simulations numériques de type CFD / The development of new materials with improved properties seems to rely nowadays more on blending and compounding than on the synthesis of chemically new polymers. Mixing may have a great effect on the morphology and structure of multi-component polymer materials. Twin-screw extruders (TSE) are widely used as mixers/reactors for blending, compounding, and reactive processing. This work aimed at developing a new instrument to measure in real time the residence time distribution (RTD) which characterizes the axial mixing and transport abilities of different screw elements based on the analysis of the transient flow pattern and systematic evaluation of mixing theory in TSE. Distributive mixing of polymer melts is characterized by the generation of interfacial area, which is experimentally much more difficult to measure. This 3D numerical simulation based on CFD is adopted

Les évaluations de vulnérabilité spécifiques pour les ressources d'eau souterraine sont des méthodes SIG qui établissent des indices qualitatifs spatiaux qui déterminent la sensibilité à l'infiltration de polluants superficiels. D'autre part, les fonctions de transfert, comme la Distribution de Temps de Séjour (DTS), sont utilisées pour prévoir le changement temporel de qualité de l'eau au forage, mais ils ne permettent pas de relier les concentrations aux pratiques observées sur le bassin versant. Basé sur une approche analytique, un modèle pour le transport souterrain des produits phytosanitaires liant la vulnérabilité spécifique SIG à la DTS a été développée. L’idée étant de relier les notions de vulnérabilité aux méthodes de prédiction des chroniques de la qualité des eaux. L'outil estime la qualité de l'eau à partir de l'ensemble de données des cartes de vulnérabilité. La validation de cette méthode de vulnérabilité spécifique est possible à partir des données de suivi de la qualité de l’eau au forage. Une formulation de paramètres équivalents a été proposée pour prendre en compte les caractéristiques hydrodynamiques des compartiments de sol (zone saturée et zone non saturée). Une validation théorique de l'approche est faite à l’aide de modèles souterrain à différence finie : HYDRUS et MODFLOW. Une application a été réalisée sur l'aquifère karstique du Val d'Orléans. Elle a permis de déterminer la DTS des concentrations en pesticides dans l'eau souterraine et a mis en évidence la présence non négligeable des métabolites du métazachlore dans le forage. En même temps, le modèle DTS, lié au SIG rend possible la localisation des zones contributrices dans le bassin versant. / The specific vulnerability estimations for the groundwater resources are GIS methods that establish spatial qualitative indices which determine the sensitivity of infiltration from surface contaminants. On the other hand, the transfer functions, using the Residence Time Distribution (RTD), are used to predict temporal water quality change in a borehole, but they do not integrate the spatial variability of the land use. Based on an analytic (advection / dispersion equation) approach, a simple GIS-linked RTD model for groundwater transport has been developed. The tool estimates the water quality from the vulnerability map dataset. This method enables to validate the specific vulnerability maps with the water quality monitoring at the borehole. It links the impacts of land use with the temporal evolution of the water quality. A equivalent formulation parameters is proposed to take into account the hydrodynamic characteristics of the soil compartments (unsaturated zone and Saturated Zone). A theoretical validation of the approach is made from finite-difference groundwater models: HYDRUS and MODFLOW. Also, an application of the RTD compilation was realized on the Val d’Orléans karstic aquifer. This last methodology allowed to determine the RTD of pesticides into the groundwater and highlighted the not insignificant presence of the metabolite of the metazachlor in the groundwater drilling. At the same time, the GIS-linked RTD model makes possible the localization of the contributing zones in the watershed.

Conical spouted beds have been commonly used for drying suspensions, solutions and pasty materials. They can also be utilized in many other processes, such as catalytic partial oxidation of methane to synthesis gas, coating of tablets, coal gasification and liquefaction, pyrolysis of sawdust or mixtures of wood residues. The main objectives of this work include both the experimental research and mathematical modelling of the conical spouted bed hydrodynamics. For experimental research, pressure transducers and static pressure probes were applied to investigate the evolution of the internal spout and the local static pressure distribution; optical fibre probes were utilized to measure axial particle velocity profiles and voidage profiles; the step tracer injection technique using helium as the tracer and thermal conductivity cells as detectors was used to investigate the gas mixing behaviour inside a conical spouted bed. It was found that many factors might affect calibration of the effective distance of an optical fibre probe. Therefore, a new calibration setup was designed and assembled, and a comprehensive sensitivity analysis was conducted to calibrate the optical probes used in this study. For mathematical modelling, a stream-tube model based on the bed structure inside a conical spouted bed was proposed to simulate partial spouting states. By introducing an adjustable parameter, this model is capable of predicting the total pressure drop under different operating conditions, and estimating axial superficial gas velocity profiles and gauge pressure profiles. A mathematical model based on characteristics of conical spouted beds and the commercial software FLUENT was also developed and validated using measured experimental data. The proposed new CFD model can simulate both stable spouting and partial spouting states, with an adjustable solids-phase source term. At stable spouting states, simulation results agree very well with almost all experimental data, such as static pressure profiles, axial particle velocity profiles, voidage profiles etc. A comprehensive sensitivity analysis was also conducted to investigate the effect of all possible factors on simulation results, including the fluid inlet profile, solid bulk viscosity, frictional viscosity, restitution coefficient, exchange coefficient, and solid phase source term. The proposed new CFD model was also used successfully to simulate gas mixing behaviours inside a conical spouted bed, and simulate cylindrical packed beds as well as cylindrical fluidized beds in one code package.

Conical spouted beds have been commonly used for drying suspensions, solutions and pasty materials. They can also be utilized in many other processes, such as catalytic partial oxidation of methane to synthesis gas, coating of tablets, coal gasification and liquefaction, pyrolysis of sawdust or mixtures of wood residues. The main objectives of this work include both the experimental research and mathematical modelling of the conical spouted bed hydrodynamics. For experimental research, pressure transducers and static pressure probes were applied to investigate the evolution of the internal spout and the local static pressure distribution; optical fibre probes were utilized to measure axial particle velocity profiles and voidage profiles; the step tracer injection technique using helium as the tracer and thermal conductivity cells as detectors was used to investigate the gas mixing behaviour inside a conical spouted bed. It was found that many factors might affect calibration of the effective distance of an optical fibre probe. Therefore, a new calibration setup was designed and assembled, and a comprehensive sensitivity analysis was conducted to calibrate the optical probes used in this study. For mathematical modelling, a stream-tube model based on the bed structure inside a conical spouted bed was proposed to simulate partial spouting states. By introducing an adjustable parameter, this model is capable of predicting the total pressure drop under different operating conditions, and estimating axial superficial gas velocity profiles and gauge pressure profiles. A mathematical model based on characteristics of conical spouted beds and the commercial software FLUENT was also developed and validated using measured experimental data. The proposed new CFD model can simulate both stable spouting and partial spouting states, with an adjustable solids-phase source term. At stable spouting states, simulation results agree very well with almost all experimental data, such as static pressure profiles, axial particle velocity profiles, voidage profiles etc. A comprehensive sensitivity analysis was also conducted to investigate the effect of all possible factors on simulation results, including the fluid inlet profile, solid bulk viscosity, frictional viscosity, restitution coefficient, exchange coefficient, and solid phase source term. The proposed new CFD model was also used successfully to simulate gas mixing behaviours inside a conical spouted bed, and simulate cylindrical packed beds as well as cylindrical fluidized beds in one code package.

One of the most important boron minerals, colemanite is reacted with sulfuric acid to produce boric acid. During this reaction, gypsum (calcium sulfate dihydrate) is formed as a byproduct. In this study, the boric acid production was handled both in a batch and four continuously stirred slurry reactors (4-CFSSR&rsquo / s) in series system. In this reaction system there are at least three phases, one liquid and two solid phases (colemanite and gypsum). In a batch reactor all the phases have the same operating time (residence time), whereas in a continuous reactor all the phases may have different residence time distributions. The residence time of both the reactant and the product solids are very important because they affect the dissolution conversion of colemanite and the growth of gypsum crystals. The main aim of this study was to investigate the dynamic behavior of continuous flow stirred slurry reactors. By obtaining the residence time distribution of the solid and liquid components, the non-idealities in the reactors can be found. The experiments performed in the continuous flow stirred slurry reactors showed that the reactors to be used during the boric acid production experiments approached an ideal CSTR in the range of the stirring rate (500-750 rpm) studied. The steady state performance of the continuous flow stirred slurry reactors (CFSSR&rsquo / s) in series was also studied. During the studies, two colemanites having the same origin but different compositions and particle sizes were used. The boric acid production reaction consists of two simultaneous reactions, dissolution of colemanite and crystallization of gypsum. The dissolution of colemanite and the gypsum formation was followed from the boric acid and calcium ion concentrations, respectively. The effect of initial CaO/ SO42- molar ratio (1.00, 1.37 and 2.17) on the boric acid and calcium ion concentrations were searched. Also, at these initial molar ratios the colemanite feed rate was varied (5, 7.5, 10 and 15 g/min) to change the residence time of the slurry. Purity of the boric acid solution was examined in terms of the selected impurities, which were the magnesium and sulfate ion concentrations. The concentrations of them were compared at the initial molar ratios of 1.00 and 1.37 with varying colemanite feed rates. It was seen that at high initial CaO/ SO42- molar ratios the sulfate and magnesium ion concentrations decreased but the calcium ion concentration increased. The gypsum crystals formed in the reaction are in the shape of thin needles. These crystals, mixed with the insolubles coming from the mineral, are removed from the boric acid slurry by filtration. Filtration of gypsum crystals has an important role in boric acid production reaction because it affects the efficiency, purity and crystallization of boric acid. These crystals must grow to an appropriate size in the reactor. The growth process of gypsum crystals should be synchronized with the dissolution reaction. The effect of solid hold-up (0.04&ndash / 0.09), defined as the volume of solid to the total volume, on the residence time of gypsum crystals was investigated and the change of the residence time (17-60 min) on the growth of the gypsum was searched. The residence time at each reactor was kept constant in each experiment as the volumes of the reactors were equal. The growth of gypsum was examined by a laser diffraction particle size analyzer and the volume weighted mean diameters of the gypsum crystals were obtained. The views of the crystals were taken under a light microscope. It was observed that the high residence time had a positive effect on the growth of gypsum crystals. The crystals had volume weighted mean diameters of even 240 &micro / m. The gypsum crystal growth model was obtained by using the second order crystallization reaction rate equation. The residence time of the continuous reactors are used together with the gypsum growth model to simulate the continuous boric acid reactors with macrofluid and microfluid models. The selected residence times (20-240 min) were modeled for different number of CSTR&rsquo / s (1-8) and the PFR. The simulated models were, then verified with the experimental data. The experimentally found calcium ion concentrations checked with the concentrations found from the microfluid model. It was also calculated that the experimental data fitted the microfluid model with a deviation of 4-7%.

Waste stabilisation ponds (WSP) are a popular form of wastewater treatment worldwide, especially for rural-based manufacturing plants and small community sewage treatment. Ponds offer a robust and operationally simple technology, which are inexpensive where land is available, and have the potential to provide a considerable degree of treatment. However the continued use of WSP is being undermined by their inconsistent performance relative to current discharge requirements, particularly with respect to suspended solids, pathogen and nutrient removal. In a climate of increased public awareness of pollution, and the ever more stringent environmental protection regulations, novel pond designs need to be developed, and existing ponds retro-fitted, to improve their performance. This dissertation investigated the hydraulic modelling of non-mechanically mixed ponds, and produced a modelling framework from which improved pond designs could be evaluated. Computational fluid dynamics (CFD) simulations were used to develop models which were able to predict the hydraulics of arbitrarily shaped, non-mechanically mixed ponds under controlled conditions. The models represent an important departure from traditional pond modelling techniques, which are based on either historical experience or simple hydraulic and reaction models. The CFD approach overcomes the main limitation of these models, as it accounts for spatial variations of parameters within a pond such as fluid velocity, or pollutant concentration. This allows for the prediction of pond hydraulics based on the pond geometry (such as inlet configuration, pond shape or baffle placement), pond inlet boundary conditions and the fluid properties. Thus CFD models allow the rapid investigation of the effect of design modifications on pond performance. The WSP models were designed using a two stage process. The first stage, a steady state simulation, calculated the velocity and turbulence fields for the pond; the second stage, a transient numerical tracer, utilised the underlying steady state results to calculate the advection and diffusion of a tracer species. The species concentration at the outlet was then integrated to produce residence time distributions (RTD) and other quantities which were used to characterise the pond hydraulics, and quantitatively compare the models with experimental results to assess the pond¡¦s performance. These techniques could be applied to any numerical pond flow model, and are a discerning test of the model¡¦s consistency. RTD generated from two-dimensional (2-D) CFD simulations were compared to experimental RTD derived by Mangelson and Watters (1972). In one of the three geometries simulated, the 2-D CFD model successfully predicted the experimental RTD. However, the flow patterns in the other two geometries were not well described, due to the difficulty of representing a three dimensional (3-D) inlet in the 2-D CFD model. As no general relationship could be found for approximating a 2-D inlet in 3-D, full 3-D simulations were used to model the unsuccessful cases. The 3 D simulations provided much improved results, predicting all the major features of the RTD over the first residence time, and matching exponential decay of the RTD after this period. Due to the uncertainty in the exact experimental inlet dimensions, a range of inlet depths were simulated. This showed that the CFD model was sensitive to changes in the inlet configuration, and using the appropriate inlet depth, the simulated RTD matched the experimental results well. A sensitivity analysis of the effect of the inlet turbulent boundary conditions and tracer molecular diffusivity for the k-ƒÕ turbulent model, showed the RTD was insensitive to these properties, thereby confirming similar results in related systems (Benelmouffok, 1989; De Vantier and Larock, 1987). This is significant for future pond modelling, as these properties are difficult to measure experimentally or predict reliably. Tracer studies were performed in this dissertation on five full-scale pond systems. In Tasmania three identical sewage ponds with different inlet and baffle configurations were investigated. However wind conditions in this locality masked any effect of these modifications. Tracer studies were also performed on sugar mill ponds near Mackay. While the models predicted qualitative consistent RTD results, they did not match the experimentally measured RTD due to uncontrolled environmental mixing factors and the long residence times of these ponds. A preliminary investigation of the effect of wind mixing was undertaken by imposing a velocity to the top surface of the model. These results confirmed the strong influence of even small wind velocities due to the large surface area of the ponds. Practical experience has indicated that the pond hydraulics are often the limiting factor in pond performance. Both experimental and simulation results have confirmed this through the presence of short circuiting and dead zones within the pond. Three baffle designs were assessed, all of which improved the pond hydraulics by either dispersing the inlet jet, or utilising the jet to generate specific pond mixing. Finally the work in this thesis has highlighted a number of other areas for future investigation. These include reservations over the use of RTD to characterise full-scale pond hydraulics, and considerations regarding the most efficient use of the inlet mixing in ponds. The hydraulic models developed in this dissertation can be extended to include solids, stratification and reaction models, which would enable the direct validation of the model based on physical or chemical parameters. In addition, a coupled flow and reaction model would provide a tool that could be used to truly optimise pond performance. This offers the possibility of tailoring the design of ponds for specific reactions, such as improved biological nutrient removal.

La torréfaction est un traitement thermique à basse température (250 à 300 °C) en atmosphère inerte qui permet de modifier les propriétés de la biomasse. La biomasse torréfiée est alors plus dense énergétiquement, plus hydrophobe et plus fragile. Dans cette étude, un modèle numérique de torréfaction en four tournant à une dimension a été développé. Le transport des plaquettes de bois, les transferts thermiques, le séchage ainsi que les cinétiques de torréfaction ont été modélisés séparément. Après confrontation aux résultats expérimentaux, ces différents sous-modèles ont été assemblés dans un modèle global. Le modèle prédit alors l’évolution de la température et de la perte de masse des plaquettes le long du four. Les résultats numériques montrent une adéquation satisfaisante avec les valeurs obtenues lors d’expériences de torréfaction sur un four tournant pilote. Les solides torréfiés ont été analysés et leurs propriétés ont été corrélées à la perte de masse. Il a en particulier été démontré que l'énergie de broyage de la biomasse torréfiée était fortement réduite. / Torrefaction is a thermal treatment at low temperature (250-300°C) used to improve biomass properties. Torrefied biomass has a higher energy density, it is more hydrophobic and more brittle. In this study, a one-dimensional numerical model of torrefaction in a rotary kiln has been developed. The wood chips flow, the thermal transfers, the drying step and the torrefaction kinetics have been modelled separately. These submodels have been experimentally validated before being implemented together. The model can thus predict the temperature and the mass loss of wood chips along the kiln. These results are in good agreement with values obtained during torrefaction experiments in the pilot-scale rotary kiln. In parallel, torrefied biomass has been analysed in terms of composition, heating value and structural properties with emphasis on the decrease of grinding energy consumption.

Os processos oxidativos avançados são promissores para a degradação de compostos orgânicos resistentes aos tratamentos convencionais, como o fenol. A fluidodinâmica computacional (CFD) tornou-se uma poderosa ferramenta para analisar processos fotoquímicos por resolver os balanços acoplados de quantidade de movimento, de massa e de radiação. O objetivo deste trabalho é investigar o processo UV/H2O2 num reator fotoquímico anular usando CFD e um modelo cinético mais realista. O modelo em CFD foi criado de forma progressiva. Inicialmente, foram determinados os campos de velocidade para três vazões (30, 60 e 100 L/h). Considerou-se dois diâmetros de lâmpada para reproduzir a configuração experimental do sistema. A discretização foi feita com malhas tetraédricas variando entre 390 000 e 1 200 000 elementos. Quatro modelos de turbulência RANS foram analisados: k-e, k-w, o shear stress transport (SST) e o modelo de tensões de Reynolds (RSM). O campo de velocidades foi validado comparando a DTR com seu levantamento experimental. A próxima etapa foi incluir o mecanismo de degradação de fenol proposto por Edalatmanesh, Dhib e Mehrvar (2008) no modelo em CFD. Trata-se de um modelo cinético baseado em equações dinâmicas para todas as espécies. O campo de radiação foi calculado pelo modelo radial e pela solução da equação de transporte de radiação através do método discrete transfer. As simulações reproduziram dados experimentais abrangendo uma larga gama de concentrações iniciais de fenol, razões molares H2O2/fenol e três potências de emissão das lâmpadas. O campo de velocidades obtido era dependente da vazão: o fluido pode manter movimento helicoidal sobre toda a extensão do reator ou se desenvolver como um escoamento pistonado. O modelo k-e não reproduziu bem o escoamento por não ser adequado para escoamentos rotativos. Os outros modelos geraram curvas de DTR com bom ajuste aos dados experimentais, especialmente o modelo k-w. O desvio médio entre as simulações de degradação de fenol e os dados experimentais é inferior a 8%. Verificou-se que, devido ao escoamento rotativo, os reagentes ficavam concentrados próximos à parede externa e migravam para a região da lâmpada ao longo do reator. A elevada intensidade de radiação na superfície da lâmpada criou uma camada ao seu redor na qual a fotólise do H2O2 ocorreu com grande taxa. Os radicais OH gerados nessa camada eram transportados para a região das paredes por convecção. Isso fez com que a maior parte do fenol fosse atacada na segunda metade do reator e gerou acúmulo do radical próximo à lâmpada na seção de saída do reator, já que o poluente já fora oxidado nessa área. O método discrete transfer previu intensidades de radiação maiores que o modelo radial, e, consequentemente, maior concentração de radicais OH. Os resultados satisfatórios indicam que CFD foi uma ferramenta adequada para analisar este escoamento reativo. / Advanced oxidation processes are a promising technology for degradation of organic compounds resistant to conventional treatments such as phenol. Computational fluid dynamics (CFD) has recently emerged as a powerful tool that allows a deeper understanding of photochemical processes in reactor engineering by solving the coupled momentum, mass and radiation balances. This work aimed to investigate the UV/H2O2 process in an annular photoreactor using CFD and a more realistic kinetic model. A progressive approach was used to develop the CFD reactor model. First, the velocity fields were determined for three volumetric flow rates (30, 60 and 100 L/h). Two lamp diameters were considered to reflect the experimental configuration of the system. Tetrahedral meshes varying form 390,000 to 1,200,000 elements were analyzed to achieve grid independence. For accounting turbulence effects, four RANS models were tested: k-e, k-w, the Shear Stress Transport (SST) and the Reynolds Stress models (RSM). The velocity field was validated through comparison to RTD experimental data. Next step was introducing the mechanism of phenol degradation proposed by Edalatmanesh, Dhib and Mehrvar (2008) into the CFD model. This kinetic model is based on dynamic equations for all species. The fluence rate field was calculated by the radial model and by solving the radiation transport equation with the discrete transfer method. Simulations reproduced experimental data spanning a wide range of initial phenol concentrations, H2O2/phenol molar ratios and three values for lamp power. It was found that the velocity field depends on the volumetric flow rate: either it maintains a swirling motion through the whole reactor or might develop like a plug flow. The k-e model did not represent the RTD data accurately, and the velocity field therefore, since it is not appropriate for swirling flows. The other turbulence models showed good match of RTD, especially the k-w model. Simulations of phenol degradation deviated less than 8% from experimental data. It was possible verified that, due to the swirling inlet effects, reactants got concentrated close to the outer wall and migrated on the lamp direction along the reactor path. High radiation intensities close to the lamp surface created a layer around it where photolysis of H2O2 took place with higher rates. OH radicals were generated in that layer and transported towards the outer wall by convection. This caused most of phenol to be consumed in the second half of the reactor and accumulation of the radical near the lamp and the reactor outlet, since the pollutant in this area was already oxidized. The discrete transfer method predicted higher incident radiation intensity than the radial model, and higher concentrations of OH radicals as a consequence. Satisfactory results indicated that CFD was an appropriate tool for analyzing this reactive flow.

Este estudo apresenta a modelagem fluidodinâmica de um reator solar utilizado para Processos Oxidativos Avançados (POA). Desenvolveu-se um modelo que considera a fluidodinâmica, o campo de radiação e cinética da reação de actinometria química (ferrioxalato), em regime transiente. Essa modelagem foi feita utilizando-se o código de fluidodinâmica computacional PHOENICS. Para análise dos resultados de simulações com o modelo, consideraram-se os dados experimentais de actinometria química para um trecho do reator, constituído de dois tubos (hairpin), e de distribuição de tempos de residência (DTR), para o reator completo, constituído de dez tubos. Os dados experimentais foram obtidos por RIBEIRO (2009). O resultado da análise da distribuição do tempo de residência do reator completo mostrou que o modelo baseado em escoamento laminar apresentou uma maior aderência aos dados experimentais de DTR. Como os experimentos de actinometria foram realizados em trecho com dois tubos do reator, construiu-se a geometria do hairpin que apresentou uma DTR mais aderente aos modelos teóricos. Outra simplificação foi necessária para a modelagem do campo de radiação de forma mais precisa, adaptando-se o hairpin para um tubo reto simples. A partir dos resultados de actinometria química foi possível estimar, pelo modelo, a taxa de fótons incidentes na parede do reator. / This study consists of the fluid dynamic modeling of a solar reactor used in Advanced Oxidation Processes (AOP). The model was developed by considering fluid dynamics, radiation field and the kinetics of the chemical actinometry reaction (ferrioxalate) in transient regime. This modeling was developed using computational fluid dynamics (CFD) in PHOENICS. Simulation results based on the model were analyzed by comparing them with a set of chemical actinometry experimental data obtained by RIBEIRO (2009). This considered a reactor section constituted by two pipes (hairpin), and the residence time distribution (RTD) of the complete reactor, composed of ten pipes. Residence time distribution results showed that the laminar flow model presented a better fitting to experimental data. Since the actinometry experiments were carried out in a reactor section with two pipes, a new geometry was designed, which resulted in a better fitting of RTD results with theoretical models. In order to obtain a more precise radiation field model, another simplification was necessary, which consisted of assuming a straight cylindrical pipe geometry. The use of the chemical actinometry and the radiation field model enabled the estimation of the incident photons rate at reactor wall.