Development of procedures for the simulation of atmospheric flows over complex terrain, using OpenFOAM

The main purpose of this work was the development of procedures for the simulation of atmospheric ows over complex terrain, using OpenFOAM. For this aim, tools and procedures were developed apart from this code for the preprocessing and data extraction, which were thereafter applied in the simulation of a real case. For the generation of the computational domain, a systematic method able to translate the terrain elevation model to a native OpenFOAM format (blockMeshDict) was developed. The outcome was a structured mesh, in which the user has the ability to de ne the number of control volumes and its dimensions. With this procedure, the di culties of case set up and the high computation computational e ort reported in literature associated to the use of snappyHexMesh, the OpenFOAM resource explored until then for the accomplishment of this task, were considered to be overwhelmed. Developed procedures for the generation of boundary conditions allowed for the automatic creation of idealized inlet vertical pro les, de nition of wall functions boundary conditions and the calculation of internal eld rst guesses for the iterative solution process, having as input experimental data supplied by the user. The applicability of the generated boundary conditions was limited to the simulation of turbulent, steady-state, incompressible and neutrally strati ed atmospheric ows, always recurring to RaNS (Reynolds-averaged Navier-Stokes) models. For the modelling of terrain roughness, the developed procedure allowed to the user the de nition of idealized conditions, like an uniform aerodynamic roughness length or making its value variable as a function of topography characteristic values, or the using of real site data, and it was complemented by the development of techniques for the visual inspection of generated roughness maps. The absence and the non inclusion of a forest canopy model limited the applicability of this procedure to low aerodynamic roughness lengths. The developed tools and procedures were then applied in the simulation of a neutrally strati ed atmospheric ow over the Askervein hill. In the performed simulations was evaluated the solution sensibility to di erent convection schemes, mesh dimensions, ground roughness and formulations of the k - ε and k - ω models. When compared to experimental data, calculated values showed a good agreement of speed-up in hill top and lee side, with a relative error of less than 10% at a height of 10 m above ground level. Turbulent kinetic energy was considered to be well simulated in the hill windward and hill top, and grossly predicted in the lee side, where a zone of ow separation was also identi ed. Despite the need of more work to evaluate the importance of the downstream recirculation zone in the quality of gathered results, the agreement between the calculated and experimental values and the OpenFOAM sensibility to the tested parameters were considered to be generally in line with the simulations presented in the reviewed bibliographic sources.

Controlo interno na Administração Local implantação de um novo procedimento numa Autarquia

Trabalho de Projeto Apresentado ao Instituto de Contabilidade e Administração do Porto para a obtenção do grau de Mestre em Auditoria, sob orientação do Mestre Carlos Mendes e coorientação da Engª Sónia Rodrigues

Mathematical model for HIV dynamics in HIV-specific helper cells

In this paper we study a delay mathematical model for the dynamics of HIV in HIV-specific CD4 + T helper cells. We modify the model presented by Roy and Wodarz in 2012, where the HIV dynamics is studied, considering a single CD4 + T cell population. Non-specific helper cells are included as alternative target cell population, to account for macrophages and dendritic cells. In this paper, we include two types of delay: (1) a latent period between the time target cells are contacted by the virus particles and the time the virions enter the cells and; (2) virus production period for new virions to be produced within and released from the infected cells. We compute the reproduction number of the model, R0, and the local stability of the disease free equilibrium and of the endemic equilibrium. We find that for values of R0<1, the model approaches asymptotically the disease free equilibrium. For values of R0>1, the model approximates asymptotically the endemic equilibrium. We observe numerically the phenomenon of backward bifurcation for values of R0⪅1. This statement will be proved in future work. We also vary the values of the latent period and the production period of infected cells and free virus. We conclude that increasing these values translates in a decrease of the reproduction number. Thus, a good strategy to control the HIV virus should focus on drugs to prolong the latent period and/or slow down the virus production. These results suggest that the model is mathematically and epidemiologically well-posed.

Foraging at the Edge of Chaos: Internal Clock versus External Forcing

Activity rhythms in animal groups arise both from external changes in the environment, as well as from internal group dynamics. These cycles are reminiscent of physical and chemical systems with quasiperiodic and even chaotic behavior resulting from “autocatalytic” mechanisms. We use nonlinear differential equations to model how the coupling between the self-excitatory interactions of individuals and external forcing can produce four different types of activity rhythms: quasiperiodic, chaotic, phase locked, and displaying over or under shooting. At the transition between quasiperiodic and chaotic regimes, activity cycles are asymmetrical, with rapid activity increases and slower decreases and a phase shift between external forcing and activity. We find similar activity patterns in ant colonies in response to varying temperature during the day. Thus foraging ants operate in a region of quasiperiodicity close to a cascade of transitions leading to chaos. The model suggests that a wide range of temporal structures and irregularities seen in the activity of animal and human groups might be accounted for by the coupling between collectively generated internal clocks and external forcings.

Virtual Power Players Internal Negotiation and Management in MASCEM

Electricity Markets are not only a new reality but an evolving one as the involved players and rules change at a relatively high rate. Multi-agent simulation combined with Artificial Intelligence techniques may result in very helpful sophisticated tools. This paper presents a new methodology for the management of coalitions in electricity markets. This approach is tested using the multi-agent market simulator MASCEM (Multi-Agent Simulator of Competitive Electricity Markets), taking advantage of its ability to provide the means to model and simulate Virtual Power Players (VPP). VPPs are represented as coalitions of agents, with the capability of negotiating both in the market and internally, with their members in order to combine and manage their individual specific characteristics and goals, with the strategy and objectives of the VPP itself. A case study using real data from the Iberian Electricity Market is performed to validate and illustrate the proposed approach.

Cost Allocation Model for Distribution Networks Considering High Penetration of Distributed Energy Resources

The high penetration of distributed energy resources (DER) in distribution networks and the competitiveenvironment of electricity markets impose the use of new approaches in several domains. The networkcost allocation, traditionally used in transmission networks, should be adapted and used in the distribu-tion networks considering the specifications of the connected resources. The main goal is to develop afairer methodology trying to distribute the distribution network use costs to all players which are usingthe network in each period. In this paper, a model considering different type of costs (fixed, losses, andcongestion costs) is proposed comprising the use of a large set of DER, namely distributed generation(DG), demand response (DR) of direct load control type, energy storage systems (ESS), and electric vehi-cles with capability of discharging energy to the network, which is known as vehicle-to-grid (V2G). Theproposed model includes three distinct phases of operation. The first phase of the model consists in aneconomic dispatch based on an AC optimal power flow (AC-OPF); in the second phase Kirschen’s andBialek’s tracing algorithms are used and compared to evaluate the impact of each resource in the net-work. Finally, the MW-mile method is used in the third phase of the proposed model. A distributionnetwork of 33 buses with large penetration of DER is used to illustrate the application of the proposedmodel.

Multiagent System Architecture for Short-term Operation of Integrated Microgrids

The forthcoming smart grids are comprised of integrated microgrids operating in grid-connected and isolated mode with local generation, storage and demand response (DR) programs. The proposed model is based on three successive complementary steps for power transaction in the market environment. The first step is characterized as a microgrid’s internal market; the second concerns negotiations between distinct interconnected microgrids; and finally, the third refers to the actual electricity market. The proposed approach is modeled and tested using a MAS framework directed to the study of the smart grids environment, including the simulation of electricity markets. This is achieved through the integration of the proposed approach with the MASGriP (Multi-Agent Smart Grid Platform) system.

Descending pain modulation in a Kaolin-induced hydrocephalus rat model

Pain transmission at the spinal cord is modulated by descending actions that arise from supraspinal areas which collectively form the endogenous pain control system. Two key areas involved of the endogenous pain control system have a circunventricular location, namely the periaqueductal grey (PAG) and the locus coeruleus (LC). The PAG plays a crucial role in descending pain modulation as it conveys the input from higher brain centers to the spinal cord. As to the LC, it is involved in descending pain inhibition by direct noradrenergic projections to the spinal cord. In the context of neurological defects, several diseases may affect the structure and function of the brain. Hydrocephalus is a congenital or acquired disease characterized by an enlargement of the ventricles which leads to a distortion of the adjacent tissues, including the PAG and LC. Usually, patients suffering from hydrocephalus present dysfunctions in learning and memory and also motor deficits. It remains to be evaluated if lesions of the periventricular brain areas involved in pain control during hydrocephalus may affect descending pain control and, herein, affect pain responses. The studies included in the present thesis used an experimental model of hydrocephalus (the rat injected in the cisterna magna with kaolin) to study descending modulation of pain, focusing on the two circumventricular regions referred above (the PAG and the LC). In order to evaluate the effects of kaolin injection into the cisterna magna, we measured the degree of ventricular dilatation in sections encompassing the PAG by standard cytoarquitectonic stanings (thionin staining). For the LC, immunodetection of the noradrenaline-synthetizing enzyme tyrosine hydroxylase (TH) was performed, due to the noradrenergic nature of the LC neurons. In general, rats with kaolin-induced hydrocephalus presented a higher dilatation of the 4th ventricle, along with a tendency to a higher area of the PAG. Due to the validated role of detection the c-fos protooncogene as a marker of neuronal activation, we also studied neuronal activation in the several subnuclei which compose the PAG, namely the dorsomedial, dorsolateral, lateral and ventrolateral (VLPAG) parts. A decrease in the numbers of neurons immunoreactive for Fos protein (the product of activation of the c-fos protooncogene) was detected in rats injected with kaolin, whereas the remaining PAG subnuclei did not present changes in Fos-immunoreactive nuclei. Increases in the levels of TH in the LC, namely at the rostral parts of the nucleus, were detected in hydrocephalic animals. The following pain-related parameters were measured, namely 1) pain behavioural responses in a validated pain inflammatory test (the formalin test) and 2) the nociceptive activation of spinal cord neurons. A decrease in behavioral responses was detected in rats with kaolin-induced hydrocephalus was detected, namely in the second phase of the test (inflammatory phase). This is the phase of the formalin test in which the motor behaviour is less important, which is important since a semi-quantitative analysis of the motor performance of rats injected with kaolin indicates that these animals may present some motor impairments. Collectively, the results of the behavioral studies indicate that rats with kaolin-induced hydrocephalus exhibit hypoalgesia. A decrease in Fos expression was detected at the superficial dorsal layers of the spinal cord in rats with kaolin-induced hydrocephalus, further indicating that hydrocephalus decreases nociceptive responses. It remains to be ascertained if this is due to alterations in the PAG and LC in the rats with kaolin-induced hydrocephalus, which may affect descending pain modulation. It remains to be evaluated what are the mechanisms underlying the increased pain inhibition at the spinal dorsal horn in the hydrocephalus rats. Regarding the VLPAG, the decrease in neuronal activity may impair descending modulation. Since the LC has higher levels of TH in rats with kaolininduced hydrocephalus, which also appears to increase the noradrenergic innervation in the spinal dorsal horn, it is possible that an increase in the release of noradrenaline at the spinal cord accounts for pain inhibition. Our studies also determine the need to study in detail patients with hydrocephalus namely in what concerns their thresholds to pain and to perform imaging studies focused on the structure and function of pain control areas in the brain.

Optimizing potentiometric ionophore and electrode design for environmental on-site control of antibiotic drugs: Application to sulfamethoxazole

Potentiometric sensors are typically unable to carry out on-site monitoring of environmental drug contaminants because of their high limits of detection (LODs). Designing a novel ligand material for the target analyte and managing the composition of the internal reference solution have been the strategies employed here to produce for the first time a potentiometric-based direct reading method for an environmental drug contaminant. This concept has been applied to sulfamethoxazole (SMX), one of the many antibiotics used in aquaculture practices that may occur in environmental waters. The novel ligand has been produced by imprinting SMX on the surface of graphitic carbon nanostructures (CN) < 500 nm. The imprinted carbon nanostructures (ICN) were dispersed in plasticizer and entrapped in a PVC matrix that included (or not) a small amount of a lipophilic additive. The membrane composition was optimized on solid-contact electrodes, allowing near-Nernstian responses down to 5.2 μg/mL and detecting 1.6 μg/mL. The membranes offered good selectivity against most of the ionic compounds in environmental water. The best membrane cocktail was applied on the smaller end of a 1000 μL micropipette tip made of polypropylene. The tip was then filled with inner reference solution containing SMX and chlorate (as interfering compound). The corresponding concentrations were studied for 1 × 10−5 to 1 × 10−10 and 1 × 10−3 to 1 × 10−8 mol/L. The best condition allowed the detection of 5.92 ng/L (or 2.3 × 10−8 mol/L) SMX for a sub-Nernstian slope of −40.3 mV/decade from 5.0 × 10−8 to 2.4 × 10−5 mol/L.

A Multi-DAG Model for Real-Time Parallel Applications with Conditional Execution

The 30th ACM/SIGAPP Symposium On Applied Computing (SAC 2015). 13 to 17, Apr, 2015, Embedded Systems. Salamanca, Spain.

An automated maneuver control framework for a remotely operated vehicle

OCEANS, 2001. MTS/IEEE Conference and Exhibition (Volume:2 )

A control framework for a remotely operated vehicle

A control framework enabling the automated maneuvering of a Remotely Operate Vehicle (ROV) is presented. The control architecture is structured according to the principle of composition of vehicle motions from a minimal set of elemental maneuvers that are designed and verified independently. The principled approach is based on distributed hybrid systems techniques, and spans integrated design, simulation and implementation as the same model is used throughout. Hybrid systems control techniques are used to synthesize the elemental maneuvers and to design protocols, which coordinate the execution of elemental maneuvers within a complex maneuver. This work is part of the Inspection of Underwater Structures (IES) project whose main objective is the implementation of a ROV-based system for the inspection of underwater structures.

Simple Kinematic Design for Evading the Forced Oscillation of a Car-Wheel Suspension System

An adaptive control damping the forced vibration of a car while passing along a bumpy road is investigated. It is based on a simple kinematic description of the desired behavior of the damped system. A modified PID controller containing an approximation of Caputo’s fractional derivative suppresses the high-frequency components related to the bumps and dips, while the low frequency part of passing hills/valleys are strictly traced. Neither a complete dynamic model of the car nor ’a priori’ information on the surface of the road is needed. The adaptive control realizes this kinematic design in spite of the existence of dynamically coupled, excitable internal degrees of freedom. The method is investigated via Scicos-based simulation in the case of a paradigm. It was found that both adaptivity and fractional order derivatives are essential parts of the control that can keep the vibration of the load at bay without directly controlling its motion.

Formation Control Driven by Cooperative Object Tracking

In this paper we introduce a formation control loop that maximizes the performance of the cooperative perception of a tracked target by a team of mobile robots, while maintaining the team in formation, with a dynamically adjustable geometry which is a function of the quality of the target perception by the team. In the formation control loop, the controller module is a distributed non-linear model predictive controller and the estimator module fuses local estimates of the target state, obtained by a particle filter at each robot. The two modules and their integration are described in detail, including a real-time database associated to a wireless communication protocol that facilitates the exchange of state data while reducing collisions among team members. Simulation and real robot results for indoor and outdoor teams of different robots are presented. The results highlight how our method successfully enables a team of homogeneous robots to minimize the total uncertainty of the tracked target cooperative estimate while complying with performance criteria such as keeping a pre-set distance between the teammates and the target, avoiding collisions with teammates and/or surrounding obstacles.

Cost allocation model for distribution networks considering high penetration of distributed energy resources

The high penetration of distributed energy resources (DER) in distribution networks and the competitive environment of electricity markets impose the use of new approaches in several domains. The network cost allocation, traditionally used in transmission networks, should be adapted and used in the distribution networks considering the specifications of the connected resources. The main goal is to develop a fairer methodology trying to distribute the distribution network use costs to all players which are using the network in each period. In this paper, a model considering different type of costs (fixed, losses, and congestion costs) is proposed comprising the use of a large set of DER, namely distributed generation (DG), demand response (DR) of direct load control type, energy storage systems (ESS), and electric vehicles with capability of discharging energy to the network, which is known as vehicle-to-grid (V2G). The proposed model includes three distinct phases of operation. The first phase of the model consists in an economic dispatch based on an AC optimal power flow (AC-OPF); in the second phase Kirschen's and Bialek's tracing algorithms are used and compared to evaluate the impact of each resource in the network. Finally, the MW-mile method is used in the third phase of the proposed model. A distribution network of 33 buses with large penetration of DER is used to illustrate the application of the proposed model.