Modelagem da dispersão de odores em ambientes urbanos

O presente trabalho investigou o problema da modelagem da dispersão de compostos odorantes em presença de obstáculos (cúbicos e com forma complexa) sob condição de estabilidade atmosférica neutra. Foi empregada modelagem numérica baseada nas equações de transporte (CFD1) bem como em modelos algébricos baseados na pluma Gausseana (AERMOD2, CALPUFF3 e FPM4). Para a validação dos resultados dos modelos e a avaliação do seu desempenho foram empregados dados de experimentos em túnel de vento e em campo. A fim de incluir os efeitos da turbulência atmosférica na dispersão, dois diferentes modelos de sub-malha associados à Simulação das Grandes Escalas (LES5) foram investigados (Smagorinsky dinâmico e WALE6) e, para a inclusão dos efeitos de obstáculos na dispersão nos modelos Gausseanos, foi empregado o modelo PRIME7. O uso do PRIME também foi proposto para o FPM como uma inovação. De forma geral, os resultados indicam que o uso de CFD/LES é uma ferramenta útil para a investigação da dispersão e o impacto de compostos odorantes em presença de obstáculos e também para desenvolvimento dos modelos Gausseanos. Os resultados também indicam que o modelo FPM proposto, com a inclusão dos efeitos do obstáculo baseado no PRIME também é uma ferramenta muito útil em modelagem da dispersão de odores devido à sua simplicidade e fácil configuração quando comparado a modelos mais complexos como CFD e mesmo os modelos regulatórios AERMOD e CALPUFF. A grande vantagem do FPM é a possibilidade de estimar-se o fator de intermitência e a relação pico-média (P/M), parâmetros úteis para a avaliação do impacto de odores. Os resultados obtidos no presente trabalho indicam que a determinação dos parâmetros de dispersão para os segmentos de pluma, bem como os parâmetros de tempo longo nas proximidades da fonte e do obstáculo no modelo FPM pode ser melhorada e simulações CFD podem ser usadas como uma ferramenta de desenvolvimento para este propósito. Palavras chave: controle de odor, dispersão, fluidodinâmica computacional, modelagem matemática, modelagem gaussiana de pluma flutuante, simulação de grandes vórtices (LES).

Estudo de falha em placas laminadas usando métodos analíticos e numéricos

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia de Manutenção

Cartografia e avaliação geotécnica de blocos rochosos em plataformas costeiras (Gaia, NW Portugal)

Mestrado em Engenharia Geotécnica e Geoambiente

Study of AC losses in medium-sized high temperature superconducting coils

The study of AC losses in superconducting pancake coils is of utmost importance for the development of superconducting devices. Due to different technical difficulties this study is usually performed considering one of two approaches: considering superconducting coils of few turns and studying AC losses in a large frequency range vs. superconducting coils with a large number of turns but measuring AC losses only in low frequencies. In this work, a study of AC losses in 128 turn superconducting coils is performed, considering frequencies ranging from 50 Hz till 1152 Hz and currents ranging from zero till the critical current of the coils. Moreover, the study of AC losses considering two different simultaneous harmonic components is also performed and results are compared to the behaviour presented by the coils when operating in a single frequency regime. Different electrical methods are used to verify the total amount of AC losses in the coil and a simple calorimetric method is presented, in order to measure AC losses in a multi-harmonic context. Different analytical and numerical methods are implemented and/or used, to design the superconducting coils and to compute the total amount of AC losses in the superconducting system and a comparison is performed to verify the advantages and drawbacks of each method.

Análisis mecánico de albañilería arqueológica de adobe bajo cargas de compresión uniaxial: El caso de Huaca de la Luna

This article presents an experimental and numerical study for the mechanical characterization under uniaxial compressive loading of the adobe masonry of one of the most emblematic archaeological complex in Peru, 'Huaca de la Luna' (100-650AD). Compression tests of prisms were carried out with original material brought to the laboratory. For measuring local deformations in the tests, displacement transducers were used which were complemented by a digital image correlation system which allowed a better understanding of the failure mechanism. The tests were then numerically simulated by modelling the masonry as a continuum media. Several approaches were considered concerning the geometrical modelling, namely 2D and 3D simplified models, and 3D refined models based on a photogrammetric reconstruction. The results showed a good approximation between the numerical prediction and the experimental response in all cases. However, the 3D models with irregular geometries seem to reproduce better the cracking pattern observed in the tests.

A FEM-based model to study the behaviour of corroded RC beams shear repaired by NSM CFRP rods technique

This paper presents the main features of finite element FE numerical model developed using the computer code FEMIX to predict the near-surface mounted NSM carbon-fiber-reinforced polymer CFRP rods shear repair contribution to corroded reinforced concrete RC beams. In the RC beams shear repaired with NSM technique, the Carbon Fibre Reinforced Polymer (CFRP) rods are placed inside pre-cut grooves onto the concrete cover of the RC beam’s lateral faces and are bonded to the concrete with high epoxy adhesive. Experimental and 3D numerical modelling results are presented in this paper in terms of load-deflection curves, and failure modes for 4 short corroded beams: two corroded beams (A1CL3-B and A1CL3-SB) and two control beams (A1T-B and A1T-SB), the beams noted with B were let repaired in bending only with NSM CFRP rods while the ones noted with SB were repaired in both bending and shear with NSM technique. The corrosion of the tensile steel bars and its effect on the shear capacity of the RC beams was discussed. Results showed that the FE model was able to capture the main aspects of the experimental load-deflection curves of the RC beams, moreover it has presented the experimental failure modes and FE numerical modelling crack patterns and both gave similar results for non-shear repaired beams which failed in diagonal tension mode of failure and for shear-repaired beams which failed due to large flexural crack at the middle of the beams along with the concrete crushing, three dimensional crack patterns were produced for shear-repaired beams in order to investigate the splitting cracks occurred at the middle of the beams and near the support.

In-plane stiffness of timber floors strengthened with CLT

Five full-scale timber floors were tested in order to analyse the in-plane behaviour of these structural systems. The main objective was an assessment of the effectiveness of in-plane strengthening using cross-laminated timber (CLT). To this end, one unstrengthened specimen (original), one specimen strengthened with a second layer of floorboards, two specimens strengthened with three CLT panels, and one specimen strengthened with two CLT panels, were tested. A numerical analysis was then performed in order to analyse the composite behaviour of the timber floors in more detail. Due to its importance as regards composite behaviour, the first phase of the experimental programme was composed of push-out tests on specimens representing the shear connection between the timber beams and the CLT panels. This paper describes the tests performed and the numerical modelling applied to evaluate the composite behaviour of the strengthened timber floors. The use of CLT panels is revealed to be an effective way to increase the in-plane stiffness of timber floors, through which the behaviour of the composite structure can be significantly changed, depending on the connection applied, or modified as required.

A FEM-based model to study the behaviour of corroded RC beams shear repaired by NSM CFRP rods technique

This paper presents the main features of finite element FE numerical model developed using the computer code FEMIX to predict the near-surface mounted NSM carbon-fiber-reinforced polymer CFRP rods shear repair contribution to corroded reinforced concrete RC beams. In the RC beams shear repaired with NSM technique, the Carbon Fibre Reinforced Polymer (CFRP) rods are placed inside pre-cut grooves onto the concrete cover of the RC beam’s lateral faces and are bonded to the concrete with high epoxy adhesive. Experimental and 3D numerical modelling results are presented in this paper in terms of load-deflection curves, failure modes and slip information of the tensile steel bars for 4 short corroded beams: two corroded beams (A1CL3-B and A1CL3-SB) and two control beams (A1T-B and A1T-SB), the beams noted with B were let repaired in bending only with NSM CFRP rods while the ones noted with SB were repaired in both bending and shear with NSM technique. The corrosion of the tensile steel bars and its effect on the shear capacity of the RC beams was discussed. Results showed that the FE model was able to capture the main aspects of the experimental load-deflection curves of the RC beams, moreover it has presented the experimental failure modes and FE numerical modelling crack patterns and both gave similar results for non-shear repaired beams which failed in diagonal tension mode of failure and for shear-repaired beams which failed due to large flexural crack at the middle of the beams along with the concrete crushing, three dimensional crack patterns were produced for shear-repaired beams in order to investigate the splitting cracks occurred at the middle of the beams and near the support.

Shaking table tests of a two-storey log house

This paper presents the findings of an experimental campaign that was conducted to investigate the seismic behaviour of log houses. A two-storey log house designed by the Portuguese company Rusticasa® was subjected to a series of shaking table tests at LNEC, Lisbon, Portugal. The paper contains the description of the geometry and construction of the house and all the aspects related to the testing procedure, namely the pre-design, the setup, instrumentation and the testing process itself. The shaking table tests were carried out with a scaled spectrum of the Montenegro (1979) earthquake, at increasing levels of PGA, starting from 0.07g, moving on to 0.28g and finally 0.5g. The log house did not suffer any major damage and remained in working condition throughout the entire process. The preliminary analysis of the overall behaviour of the log house is also discussed.

Tunnel engineering – influence of the type and the quantity of measurements in the back analysis of geomechanical parameters

The monitoring data collected during tunnel excavation can be used in inverse analysis procedures in order to identify more realistic geomechanical parameters that can increase the knowledge about the interested formations. These more realistic parameters can be used in real time to adapt the project to the real structure in situ behaviour. However, monitoring plans are normally designed for safety assessment and not especially for the purpose of inverse analysis. In fact, there is a lack of knowledge about what types and quantity of measurements are needed to succeed in identifying the parameters of interest. Also, the optimisation algorithm chosen for the identification procedure may be important for this matter. In this work, this problem is addressed using a theoretical case with which a thorough parametric study was carried out using two optimisation algorithms based on different calculation paradigms, namely a conventional gradient-based algorithm and an evolution strategy algorithm. Calculations were carried for different sets of parameters to identify several combinations of types and amount of monitoring data. The results clearly show the high importance of the available monitoring data and the chosen algorithm for the success rate of the inverse analysis process.

Development of an integral viscoelastic flow solver in OpenFOAM®

The usual high cost of commercial codes, and some technical limitations, clearly limits the employment of numerical modelling tools in both industry and academia. Consequently, the number of companies that use numerical code is limited and there a lot of effort put on the development and maintenance of in-house academic based codes. Having in mind the potential of using numerical modelling tools as a design aid, of both products and processes, different research teams have been contributing to the development of open source codes/libraries. In this framework, any individual can take advantage of the available code capabilities and/or implement additional features based on his specific needs. These type of codes are usually developed by large communities, which provide improvements and new features in their specific fields of research, thus increasing significantly the code development process. Among others, OpenFOAM® multi-physics computational library, developed by a very large and dynamic community, nowadays comprises several features usually only available in their commercial counterparts; e.g. dynamic meshes, large diversity of complex physical models, parallelization, multiphase models, to name just a few. This computational library is developed in C++ and makes use of most of all language capabilities to facilitate the implementation of new functionalities. Concerning the field of computational rheology, OpenFOAM® solvers were recently developed to deal with the most relevant differential viscoelastic rheological models, and stabilization techniques are currently being verified. This work describes the implementation of a new solver in OpenFOAM® library, able to cope with integral viscoelastic models based on the deformation field method. The implemented solver is verified through the comparison of the predicted results with analytical solutions, results published in the literature and by using the Method of Manufactured Solutions.

Development of a new 3D OpenFOAM® solver to model the cooling stage in profile extrusion

[Excerpt] The advantages resulting from the use of numerical modelling tools to support the design of processing equipment are almost consensual. The design of calibration systems in profile extrusion is not an exception . H owever , the complex geome tries and heat exchange phenomena involved in this process require the use of numerical solvers able to model the heat exchange in more than one domain ( calibrator and polymer), the compatibilization of the heat transfer at the profile - calibrator interface and with the ability to deal with complex geometries. The combination of all these features is usually hard to find in commercial software. Moreover , the dimension of the meshes required to ob tain accurate results, result in computational times prohibitive for industrial application. (...)

Using OpenFOAM to aid the design of complex profile extrusion dies

[Extrat] Thermoplastic profiles are very attractive due to their inherent design freedom. However, the usual methodologies employed to design extrusion forming tools, based on experimental based trial–and–error procedures, are highly dependent on the designer’s experience and lead to high resources consumption. Despite of the relatively low cost of the raw materials employed on the production of this type of profiles, the resources involved in the die design process significantly increase their cost. These difficulties are even more evident when a complex geometry profile has to be produced and there is no previous experience with similar geometries. Therefore, novel design approaches are required, in order to reduce the required resources and guarantee a good performance for the produced profile. (...)

Understanding the dispersion of fillers in polymer matrices

Understanding the mixing process of complex composite materials is fundamental in several industrial processes. For instance, the dispersion of fillers in polymer melt matrices is commonly employed to manufacture polymer composites, using a twin-screw extruder. The effectiveness of the filler dispersion depends not only on the complex flow patterns generated, but also on the polymer melt rheological behavior. Therefore, the availability of a numerical tool able to predict mixing, taking into account both fluid and particles phases would be very useful to increase the process insight, and thus provide useful guidelines for its optimization. In this work, a new Eulerian-Lagrangian numerical solver is developed OpenFOAM® computational library, and used to better understand the mechanisms determining the dispersion of fillers in polymer matrices. Particular attention will be given to the effect of the rheological model used to represent the fluid behavior, on the level of dispersion obtained. For the Eulerian phase the averaged volume fraction governing equations (conservation of mass and linear momentum) are used to describe the fluid behavior. In the case of the Lagrangian phase, Newton’s second law of motion is used to compute the particles trajectories and velocity. To study the effect of fluid behavior on the filler dispersion, several systems are modeled considering different fluid types (generalized Newtonian or viscoelastic) and particles volume fraction and size. The results obtained are used to correlate the fluid and particle characteristics on the effectiveness of mixing and morphology obtained.

A new integral viscoelastic flow solver in OpenFOAM®

The usual high cost of commercial codes, and some technical limitations, clearly limits the employment of numerical modelling tools in both industry and academia. Consequently, the number of companies that use numerical code is limited and there a lot of effort put on the development and maintenance of in-house academic based codes . Having in mind the potential of using numerical modelling tools as a design aid, of both products and processes, different research teams have been contributing to the development of open source codes/libraries. In this framework, any individual can take advantage of the available code capabilities and/or implement additional features based on his specific needs. These type of codes are usually developed by large communities, which provide improvements and new features in their specific fields of research, thus increasing significantly the code development process. Among others, OpenFOAM® multi-physics computational library, developed by a very large and dynamic community, nowadays comprises several features usually only available in their commercial counterparts; e.g. dynamic meshes, large diversity of complex physical models, parallelization, multiphase models, to name just a few. This computational library is developed in C++ and makes use of most of all language capabilities to facilitate the implementation of new functionalities. Concerning the field of computational rheology, OpenFOAM® solvers were recently developed to deal with the most relevant differential viscoelastic rheological models, and stabilization techniques are currently being verified. This work describes the implementation of a new solver in OpenFOAM® library, able to cope with integral viscoelastic models based on the deformation field method. The implemented solver is verified through the comparison of the predicted results with analytical solutions, results published in the literature and by using the Method of Manufactured Solutions