Comparison of different numerical methods for the solution of the time-fractional reaction-diffusion equation with variable diffusion coefficient

In this work we perform a comparison of two different numerical schemes for the solution of the time-fractional diffusion equation with variable diffusion coefficient and a nonlinear source term. The two methods are the implicit numerical scheme presented in [M.L. Morgado, M. Rebelo, Numerical approximation of distributed order reaction- diffusion equations, Journal of Computational and Applied Mathematics 275 (2015) 216-227] that is adapted to our type of equation, and a colocation method where Chebyshev polynomials are used to reduce the fractional differential equation to a system of ordinary differential equations

High strain rate constitutive modeling forhistorical structures subjected to blast loading

Doctoral Thesis Civil Engineering

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.

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

Improvign mixing processes with the help of numerical tools

In several industrial applications, highly complex behaviour materials are used together with intricate mixing processes, which difficult the achievement of the desired properties for the produced materials. This is the case of the well-known dispersion of nano-sized fillers in a melt polymer matrix, used to improve the nanocomposite mechanical and/or electrical properties. This mixing is usually performed in twin-screw extruders, that promote complex flow patterns, and, since an in loco analysis of the material evolution and mixing is difficult to perform, numerical tools can be very useful to predict the evolution and behaviour of the material. This work presents a numerical based study to improve the understanding of mixing processes. Initial numerical studies were performed with generalized Newtonian fluids, but, due to the null relaxation time that characterize this type of fluids, the assumption of viscoelastic behavior was required. Therefore, the polymer melt was rheologically characterized, and, a six mode Phan-Thien-Tanner and Giesekus models were used to fit the rheological data. These viscoelastic rheological models were used to model the process. The conclusions obtained in this work provide additional and useful data to correlate the type and intensity of the deformation history promoted to the polymer nanocomposite and the quality of the mixing obtained.

Analysis of portuguese residential buildings’ needs and proposed solutions

The Portuguese housing sector experienced a significant growth throughout the 20th century, particularly in the last quarter, after the democratic revolution in 1974. In fact, the number of buildings built between 1970 and 1990 is more than one third of the buildings actually existing in Portugal. Therefore most of them were built before the publication of the first regulation concerning the energy efficiency in buildings. Regarding this scenario, it would be expected that rehabilitation activities would represent most of the current construction activities. However, given some remaining barriers from old social policies, this situation is not observed; actually building retrofitting is the least significant sector, accentuating the degradation level of major part of the Portuguese housing stock. Several studies show that the main problems are found in the buildings envelope elements, such as roofs and façades. Based on this context, the aim of this paper is to introduce some examples of building retrofitting systems that, adapted to the Portuguese main needs and requirements may represent sustainable solutions to overcome the identified needs of Portuguese buildings' envelope.

Global vs. local nonlinear optimization techniques for human-like movement of an anthropomorphic robot

In this paper a comparison between using global and local optimization techniques for solving the problem of generating human-like arm and hand movements for an anthropomorphic dual arm robot is made. Although the objective function involved in each optimization problem is convex, there is no evidence that the admissible regions of these problems are convex sets. For the sequence of movements for which the numerical tests were done there were no significant differences between the optimal solutions obtained using the global and the local techniques. This suggests that the optimal solution obtained using the local solver is indeed a global solution.

Anderson localization of light in disordered superlattices containing graphene layers

We theoretically investigate light propagation and Anderson localization in one-dimensional disordered superlattices composed of dielectric stacks with graphene sheets in between. Disorder is introduced either on graphene material parameters ({\it e.g.} Fermi energy) or on the widths of the dielectric stacks. We derive an analytic expression for the localization length $\xi$, and compare it to numerical simulations using transfer matrix technique; a very good agreement is found. We demonstrate that the presence of graphene may strongly attenuate the anomalously delocalised Breswter modes, and is at the origin of a periodic dependence of $\xi$ on frequency, in contrast to the usual asymptotic decay, $\xi \propto \omega^{-2}$. By unveiling the effects of graphene on Anderson localization of light, we pave the way for new applications of graphene-based, disordered photonic devices in the THz spectral range.

Electronic transport across linear defects in graphene

We investigate the low-energy electronic transport across grain boundaries in graphene ribbons and infinite flakes. Using the recursive Green’s function method, we calculate the electronic transmission across different types of grain boundaries in graphene ribbons. We show results for the charge density distribution and the current flow along the ribbon. We study linear defects at various angles with the ribbon direction, as well as overlaps of two monolayer ribbon domains forming a bilayer region. For a class of extended defect lines with periodicity 3, an analytic approach is developed to study transport in infinite flakes. This class of extended grain boundaries is particularly interesting, since the K and K0 Dirac points are superposed.

Numerical modeling and assessment of the thermal environment in an industrial premise of automotive components

Dissertação de mestrado em Engenharia Industrial

Numerical simulation of the reflow soldering process

Dissertação de mestrado integrado em Mechanical Engineering

Numerical and experimental analysis of full scale arches reinforced with GFRP materials

In this contribution, original limit analysis numerical results are presented dealing with some reinforced masonry arches tested at the University of Minho-UMinho, PT. Twelve in-scale circular masonry arches were considered, reinforced in various ways at the intrados or at the extrados. GFRP reinforcements were applied either on undamaged or on previously damaged elements, in order to assess the role of external reinforcements even in repairing interventions. The experimental results were critically discussed at the light of limit analysis predictions, based on a 3D FE heterogeneous upper bound approach. Satisfactory agreement was found between experimental evidences and the numerical results, in terms of failure mechanisms and peak load.

Bacterial vaginosis biofilms: challenges to current therapies and emerging solutions

Bacterial vaginosis (BV) is the most common genital tract infection in women during their reproductive years and it has been associated with serious health complications, such as preterm delivery and acquisition or transmission of several sexually transmitted agents. BV is characterized by a reduction of beneficial lactobacilli and a significant increase in number of anaerobic bacteria, including Gardnerella vaginalis, Atopobium vaginae, Mobiluncus spp., Bacteroides spp. and Prevotella spp.. Being polymicrobial in nature, BV etiology remains unclear. However, it is certain that BV involves the presence of a thick vaginal multi-species biofilm, where G. vaginalis is the predominant species. Similar to what happens in many other biofilm-related infections, standard antibiotics, like metronidazole, are unable to fully eradicate the vaginal biofilm, which can explain the high recurrence rates of BV. Furthermore, antibiotic therapy can also cause a negative impact on the healthy vaginal microflora. These issues sparked the interest in developing alternative therapeutic strategies. This review provides a quick synopsis of the currently approved and available antibiotics for BV treatment while presenting an overview of novel strategies that are being explored for the treatment of this disorder, with special focus on natural compounds that are able to overcome biofilm-associated antibiotic resistance.