Numerical modelling of compound channel flow

Fundação para a Ciência e a Tecnologia (FCT) - SFRH/BD/64337/2009 ; projects PTDC/ECM/70652/2006, PTDC/ECM/117660/2010 and RECI/ECM-HID/0371/2012

DYNAMOD – A dynamic agent based modelling framework for digital businesses

Digital Businesses have become a major driver for economic growth and have seen an explosion of new startups. At the same time, it also includes mature enterprises that have become global giants in a relatively short period of time. Digital Businesses have unique characteristics that make the running and management of a Digital Business much different from traditional offline businesses. Digital businesses respond to online users who are highly interconnected and networked. This enables a rapid flow of word of mouth, at a pace far greater than ever envisioned when dealing with traditional products and services. The relatively low cost of incremental user addition has led to a variety of innovation in pricing of digital products, including various forms of free and freemium pricing models. This thesis explores the unique characteristics and complexities of Digital Businesses and its implications on the design of Digital Business Models and Revenue Models. The thesis proposes an Agent Based Modeling Framework that can be used to develop Simulation Models that simulate the complex dynamics of Digital Businesses and the user interactions between users of a digital product. Such Simulation models can be used for a variety of purposes such as simple forecasting, analysing the impact of market disturbances, analysing the impact of changes in pricing models and optimising the pricing for maximum revenue generation or a balance between growth in usage and revenue generation. These models can be developed for a mature enterprise with a large historical record of user growth rate as well as for early stage enterprises without much historical data. Through three case studies, the thesis demonstrates the applicability of the Framework and its potential applications.

Characterization of an OPCPA laser amplifier based on the nonlinear crystal YCOB

The work described in this thesis was performed at the Laboratory for Intense Lasers (L2I) of Instituto Superior Técnico, University of Lisbon (IST-UL). Its main contribution consists in the feasibility study of the broadband dispersive stages for an optical parametric chirped pulse amplifier based on the nonlinear crystal yttrium calcium oxi-borate (YCOB). In particular, the main goal of this work consisted in the characterization and implementation of the several optical devices involved in pulse expansion and compression of the amplified pulses to durations of the order of a few optical cycles (20 fs). This type of laser systems find application in fields such as medicine, telecommunications and machining, which require high energy, ultrashort (sub-100 fs) pulses. The main challenges consisted in the preliminary study of the performance of the broadband amplifier, which is essential for successfully handling pulses with bandwidths exceeding 100 nm when amplified from the μJ to 20 mJ per pulse. In general, the control, manipulation and characterization of optical phenomena on the scale of a few tens of fs and powers that can reach the PW level are extremely difficult and challenging due to the complexity of the phenomena of radiation-matter interaction and their nonlinearities, observed at this time scale and power level. For this purpose the main dispersive components were characterized in detail, specifically addressing the demonstration of pulse expansion and compression. The tested bandwidths are narrower than the final ones, in order to confirm the parameters of these elements and predict the performance for the broadband pulses. The work performed led to additional tasks such as a detailed characterization of laser oscillator seeding the laser chain and the detection and cancelling of additional sources of dispersion.

The thermal effects on the methanol-to-olefins reaction: A modelling and experimental approach

With the projection of an increasing world population, hand-in-hand with a journey towards a bigger number of developed countries, further demand on basic chemical building blocks, as ethylene and propylene, has to be properly addressed in the next decades. The methanol-to-olefins (MTO) is an interesting reaction to produce those alkenes using coal, gas or alternative sources, like biomass, through syngas as a source for the production of methanol. This technology has been widely applied since 1985 and most of the processes are making use of zeolites as catalysts, particularly ZSM-5. Although its selectivity is not especially biased over light olefins, it resists to a quick deactivation by coke deposition, making it quite attractive when it comes to industrial environments; nevertheless, this is a highly exothermic reaction, which is hard to control and to anticipate problems, such as temperature runaways or hot-spots, inside the catalytic bed. The main focus of this project is to study those temperature effects, by addressing both experimental, where the catalytic performance and the temperature profiles are studied, and modelling fronts, which consists in a five step strategy to predict the weight fractions and activity. The mind-set of catalytic testing is present in all the developed assays. It was verified that the selectivity towards light olefins increases with temperature, although this also leads to a much faster catalyst deactivation. To oppose this effect, experiments were carried using a diluted bed, having been able to increase the catalyst lifetime between 32% and 47%. Additionally, experiments with three thermocouples placed inside the catalytic bed were performed, analysing the deactivation wave and the peaks of temperature throughout the bed. Regeneration was done between consecutive runs and it was concluded that this action can be a powerful means to increase the catalyst lifetime, maintaining a constant selectivity towards light olefins, by losing acid strength in a steam stabilised zeolitic structure. On the other hand, developments on the other approach lead to the construction of a raw basic model, able to predict weight fractions, that should be tuned to be a tool for deactivation and temperature profiles prediction.

Rating corporate groups a modelling experience

This paper aims to provide a model that allows BPI to measure the credit risk, through its rating scale, of the subsidiaries included in the corporate groups who are their clients. This model should be simple enough to be applied in practice, accurate, and must give consistent results in comparison to what have been the ratings given by the bank. The model proposed includes operational, strategic, and financial factors and ends up giving one of three results: no support, partial support, or full support from the holding to the subsidiary, and each of them translates in adjustments in each subsidiary’s credit rating. As it would be expectable, most of the subsidiaries should have the same credit rating of its parent company.

The potential field method and the nonlinear attractor dynamics approach: what are the differences?

One of the most popular approaches to path planning and control is the potential field method. This method is particularly attractive because it is suitable for on-line feedback control. In this approach the gradient of a potential field is used to generate the robot's trajectory. Thus, the path is generated by the transient solutions of a dynamical system. On the other hand, in the nonlinear attractor dynamic approach the path is generated by a sequence of attractor solutions. This way the transient solutions of the potential field method are replaced by a sequence of attractor solutions (i.e., asymptotically stable states) of a dynamical system. We discuss at a theoretical level some of the main differences of these two approaches.

Plastic-damage smeared crack model to simulate the behaviour of structures made by cement based materials

This work proposes a constitutive model to simulate nonlinear behaviour of cement based materials subjected to different loading paths. The model incorporates a multidirectional fixed smeared crack approach to simulate crack initiation and propagation, whereas the inelastic behaviour of material between cracks is treated by a numerical strategy that combines plasticity and damage theories. For capturing more realistically the shear stress transfer between the crack surfaces, a softening diagram is assumed for modelling the crack shear stress versus crack shear strain. The plastic damage model is based on the yield function, flow rule and evolution law for hardening variable, and includes an explicit isotropic damage law to simulate the stiffness degradation and the softening behaviour of cement based materials in compression. This model was implemented into the FEMIX computer program, and experimental tests at material scale were simulated to appraise the predictive performance of this constitutive model. The applicability of the model for simulating the behaviour of reinforced concrete shear wall panels submitted to biaxial loading conditions, and RC beams failing in shear is investigated.

Numerical Assessment of the out-of-plane response of a brick masonry structure without box behaviour

This paper presents the assessment of the out-of-plane response due to seismic loading of a masonry structure without rigid diaphragm. This structure corresponds to real scale brick masonry specimen with a main façade connected to two return walls. Two modelling approaches were defined for this evaluation. The first one consisted on macro modelling, whereas the second one on simplified micro modelling. As a first step of this study, static nonlinear analyses were conducted to the macro model aiming at evaluating the out-of-plane response and failure mechanism of the masonry structure. A sensibility analyses was performed in order to assess the mesh size and material model dependency. In addition, the macro models were subjected to dynamic nonlinear analyses with time integration in order to assess the collapse mechanism. Finally, these analyses were also applied to a simplified micro model of the masonry structure. Furthermore, these results were compared to experimental response from shaking table tests. It was observed that these numerical techniques simulate correctly the in-plane behaviour of masonry structures. However, the

Seismic vulnerability of existing masonry buildings: Nonlinear parametric analysis

Existing masonry structures are usually associated to a high seismic vulnerability, mainly due to the properties of the materials, weak connections between floors and load-bearing walls, high mass of the masonry walls and flexibility of the floors. For these reasons, the seismic performance of existing masonry structures has received much attention in the last decades. This study presents the parametric analysis taking into account the deviations on features of the gaioleiro buildings - Portuguese building typology. The main objective of the parametric analysis is to compare the seismic performance of the structure as a function of the variations of its properties with respect to the response of a reference model. The parametric analysis was carried out for two types of structural analysis, namely for the non-linear dynamic analysis with time integration and for the pushover analysis with distribution of forces proportional to the inertial forces of the structure. The Young's modulus of the masonry walls, Young's modulus of the timber floors, the compressive and tensile non-linear properties (strength and fracture energy) were the properties considered in both type of analysis. Additionally, in the dynamic analysis, the influences of the vis-cous damping and of the vertical component of the earthquake were evaluated. A pushover analysis proportional to the modal displacement of the first mode in each direction was also carried out. The results shows that the Young's modulus of the masonry walls, the Young's modulus of the timber floors and the compressive non-linear properties are the pa-rameters that most influence the seismic performance of this type of tall and weak existing masonry structures. Furthermore, it is concluded that that the stiffness of the floors influences significantly the strength capacity and the collapse mecha-nism of the numerical model. Thus, a study on the strengthening of the floors was also carried out. The increase of the thickness of the timber floors was the strengthening technique that presented the best seismic performance, in which the reduction of the out-of-plane displacements of the masonry walls is highlighted.

Evaluation of seismic vulnerability assessment parameters for portuguese vernacular constructions with nonlinear numerical analysis

Considering that vernacular architecture may bear important lessons on hazard mitigation and that well-constructed examples showing traditional seismic resistant features can present far less vulnerability than expected, this study aims at understanding the resisting mechanisms and seismic behavior of vernacular buildings through detailed finite element modeling and nonlinear static (pushover) analysis. This paper focuses specifically on a type of vernacular rammed earth constructions found in the Portuguese region of Alentejo. Several rammed earth constructions found in the region were selected and studied in terms of dimensions, architectural layout, structural solutions, construction materials and detailing and, as a result, a reference model was built, which intends to be a simplified representative example of these constructions, gathering the most common characteristics. Different parameters that may affect the seismic response of this type of vernacular constructions have been identified and a numerical parametric study was defined aiming at evaluating and quantifying their influence in the seismic behavior of this type of vernacular buildings. This paper is part of an ongoing research which includes the development of a simplified methodology for assessing the seismic vulnerability of vernacular buildings, based on vulnerability index evaluation methods.

Thermal conditions affecting heat transfer in FDM/FFE: a contribution towards the numerical modelling of the process

The performance of parts produced by Free Form Extrusion (FFE), an increasingly popular additive manufacturing technique, depends mainly on their dimensional accuracy, surface quality and mechanical performance. These attributes are strongly influenced by the evolution of the filament temperature and deformation during deposition and solidification. Consequently, the availability of adequate process modelling software would offer a powerful tool to support efficient process set-up and optimisation. This work examines the contribution to the overall heat transfer of various thermal phenomena developing during the manufacturing sequence, including convection and radiation with the environment, conduction with support and between adjacent filaments, radiation between adjacent filaments and convection with entrapped air. The magnitude of the mechanical deformation is also studied. Once this exercise is completed, it is possible to select the material properties, process variables and thermal phenomena that should be taken in for effective numerical modelling of FFE.

Modelling Amazonian forest eddy covariance data: a comparison of big leaf versus sun/shade models for the C-14 tower at Manaus I. Canopy photosynthesis

In this study, we concentrate on modelling gross primary productivity using two simple approaches to simulate canopy photosynthesis: "big leaf" and "sun/shade" models. Two approaches for calibration are used: scaling up of canopy photosynthetic parameters from the leaf to the canopy level and fitting canopy biochemistry to eddy covariance fluxes. Validation of the models is achieved by using eddy covariance data from the LBA site C14. Comparing the performance of both models we conclude that numerically (in terms of goodness of fit) and qualitatively, (in terms of residual response to different environmental variables) sun/shade does a better job. Compared to the sun/shade model, the big leaf model shows a lower goodness of fit and fails to respond to variations in the diffuse fraction, also having skewed responses to temperature and VPD. The separate treatment of sun and shade leaves in combination with the separation of the incoming light into direct beam and diffuse make sun/shade a strong modelling tool that catches more of the observed variability in canopy fluxes as measured by eddy covariance. In conclusion, the sun/shade approach is a relatively simple and effective tool for modelling photosynthetic carbon uptake that could be easily included in many terrestrial carbon models.

Modelling integral viscoelastic flows in OpenFOAM

[Excerpt] A large number of constitutive equations were developed for viscoelastic fluids, some empirical and other with strong physical foundations. The currently available macroscopic constitutive equations can be divided in two main types: differential and integral. Some of the constitutive equations, e.g. Maxwell are available both in differential and integral types. However, relevant in tegral models, like K - BKZ, just possesses the integral form. (...)

Spatio-temporal modelling of environmental data

Programa Doutoral em Matemática e Aplicações.