Numerical solution of a PDE system with non-linear steady state conditions that translates the air stripping pollutants removal

This work deals with the numerical simulation of air stripping process for the pre-treatment of groundwater used in human consumption. The model established in steady state presents an exponential solution that is used, together with the Tau Method, to get a spectral approach of the solution of the system of partial differential equations associated to the model in transient state.

The use of the boundary element method in the analysis of single lap joints

The most common techniques for stress analysis/strength prediction of adhesive joints involve analytical or numerical methods such as the Finite Element Method (FEM). However, the Boundary Element Method (BEM) is an alternative numerical technique that has been successfully applied for the solution of a wide variety of engineering problems. This work evaluates the applicability of the boundary elem ent code BEASY as a design tool to analyze adhesive joints. The linearity of peak shear and peel stresses with the applied displacement is studied and compared between BEASY and the analytical model of Frostig et al., considering a bonded single-lap joint under tensile loading. The BEM results are also compared with FEM in terms of stress distributions. To evaluate the mesh convergence of BEASY, the influence of the mesh refinement on peak shear and peel stress distributions is assessed. Joint stress predictions are carried out numerically in BEASY and ABAQUS®, and analytically by the models of Volkersen, Goland, and Reissner and Frostig et al. The failure loads for each model are compared with experimental results. The preparation, processing, and mesh creation times are compared for all models. BEASY results presented a good agreement with the conventional methods.

Strap repairs using embedded patches: numerical analysis and experimental results

Adhesively bonded repairs offer an attractive option for repair of aluminium structures, compared to more traditional methods such as fastening or welding. The single-strap (SS) and double-strap (DS) repairs are very straightforward to execute but stresses in the adhesive layer peak at the overlap ends. The DS repair requires both sides of the damaged structures to be reachable for repair, which is often not possible. In strap repairs, with the patches bonded at the outer surfaces, some limitations emerge such as the weight, aerodynamics and aesthetics. To minimize these effects, SS and DS repairs with embedded patches were evaluated in this work, such that the patches are flush with the adherends. For this purpose, in this work standard SS and DS repairs, and also with the patches embedded in the adherends, were tested under tension to allow the optimization of some repair variables such as the overlap length (LO) and type of adhesive, thus allowing the maximization of the repair strength. The effect of embedding the patch/patches on the fracture modes and failure loads was compared with finite elements (FE) analysis. The FE analysis was performed in ABAQUS® and cohesive zone modelling was used for the simulation of damage onset and growth in the adhesive layer. The comparison with the test data revealed an accurate prediction for all kinds of joints and provided some principles regarding this technique.

On the analytical solutions of the Hindmarsh-Rose neuronal model

In this article we analytically solve the Hindmarsh-Rose model (Proc R Soc Lond B221:87-102, 1984) by means of a technique developed for strongly nonlinear problems-the step homotopy analysis method. This analytical algorithm, based on a modification of the standard homotopy analysis method, allows us to obtain a one-parameter family of explicit series solutions for the studied neuronal model. The Hindmarsh-Rose system represents a paradigmatic example of models developed to qualitatively reproduce the electrical activity of cell membranes. By using the homotopy solutions, we investigate the dynamical effect of two chosen biologically meaningful bifurcation parameters: the injected current I and the parameter r, representing the ratio of time scales between spiking (fast dynamics) and resting (slow dynamics). The auxiliary parameter involved in the analytical method provides us with an elegant way to ensure convergent series solutions of the neuronal model. Our analytical results are found to be in excellent agreement with the numerical simulations.

On the numerical computation of the Mittag-Leffler function

Recently simple limiting functions establishing upper and lower bounds on the Mittag-Leffler function were found. This paper follows those expressions to design an efficient algorithm for the approximate calculation of expressions usual in fractional-order control systems. The numerical experiments demonstrate the superior efficiency of the proposed method.

Desenvolvimento de um procedimento de Numerical Site Calibration para um Parque Eólico usando Dinâmica de Fluidos Computacional

Esta dissertação descreve o desenvolvimento e avaliação de um procedimento de \Numerical Site Calibration" (NSC) para um Parque Eólico, situado a sul de Portugal, usando Dinâmica de Fluídos Computacional (CFD). O NSC encontra-se baseado no \Site Calibration" (SC), sendo este um método de medição padronizado pela Comissão Electrónica Internacional através da norma IEC 61400. Este método tem a finalidade de quantificar e reduzir os efeitos provocados pelo terreno e por possíveis obstáculos, na medição do desempenho energético das turbinas eólicas. Assim, no SC são realizadas medições em dois pontos, no mastro referência e no local da turbina (mastro temporário). No entanto, em Parques Eólicos já construídos, este método não é aplicável visto ser necessária a instalação de um mastro de medição no local da turbina e, por conseguinte, o procedimento adequado para estas circunstâncias é o NSC. O desenvolvimento deste método é feito por um código CFD, desenvolvido por uma equipa de investigação do Instituto Superior de Engenharia do Porto, designado de WINDIETM, usado extensivamente pela empresa Megajoule Inovação, Lda em aplicações de energia eólica em todo mundo. Este código é uma ferramenta para simulação de escoamentos tridimensionais em terrenos complexos. As simulações do escoamento são realizadas no regime transiente utilizando as equações de Navier-Stokes médias de Reynolds com aproximação de Bussinesq e o modelo de turbulência TKE 1.5. As condições fronteira são provenientes dos resultados de uma simulação realizada com Weather Research and Forecasting, WRF. Estas simulações dividem-se em dois grupos, um dos conjuntos de simulações utiliza o esquema convectivo Upwind e o outro utiliza o esquema convectivo de 4aordem. A análise deste método é realizada a partir da comparação dos dados obtidos nas simulações realizadas no código WINDIETM e a coleta de dados medidos durante o processo SC. Em suma, conclui-se que o WINDIETM e as suas configurações reproduzem bons resultados de calibração, ja que produzem erros globais na ordem de dois pontos percentuais em relação ao SC realizado para o mesmo local em estudo.

On the numerical computation of the Mittag-Leffler function

Recently simple limiting functions establishing upper and lower bounds on the Mittag-Leffler function were found. This paper follows those expressions to design an efficient algorithm for the approximate calculation of expressions usual in fractional-order control systems. The numerical experiments demonstrate the superior efficiency of the proposed method.

An efficient numerical scheme for solving multi-dimensional fractional optimal control problems with a quadratic performance index

The shifted Legendre orthogonal polynomials are used for the numerical solution of a new formulation for the multi-dimensional fractional optimal control problem (M-DFOCP) with a quadratic performance index. The fractional derivatives are described in the Caputo sense. The Lagrange multiplier method for the constrained extremum and the operational matrix of fractional integrals are used together with the help of the properties of the shifted Legendre orthonormal polynomials. The method reduces the M-DFOCP to a simpler problem that consists of solving a system of algebraic equations. For confirming the efficiency and accuracy of the proposed scheme, some test problems are implemented with their approximate solutions.

A filter inexact-restoration method for nonlinear programming

A new iterative algorithm based on the inexact-restoration (IR) approach combined with the filter strategy to solve nonlinear constrained optimization problems is presented. The high level algorithm is suggested by Gonzaga et al. (SIAM J. Optim. 14:646–669, 2003) but not yet implement—the internal algorithms are not proposed. The filter, a new concept introduced by Fletcher and Leyffer (Math. Program. Ser. A 91:239–269, 2002), replaces the merit function avoiding the penalty parameter estimation and the difficulties related to the nondifferentiability. In the IR approach two independent phases are performed in each iteration, the feasibility and the optimality phases. The line search filter is combined with the first one phase to generate a “more feasible” point, and then it is used in the optimality phase to reach an “optimal” point. Numerical experiences with a collection of AMPL problems and a performance comparison with IPOPT are provided.

Heat recovery from wastewater: numerical modelling of sewer systems

Dissertação para obtenção do Grau de Mestre em Engenharia do Ambiente, perfil Engenharia Sanitária

Numerical Evaluation of Ventilation Performance in Children Day Care Centres

Modelling of ventilation is strongly dependent on the physical characteristics of the building of which precise evaluation is a complex and time consuming task. In the frame of a research project, two children day care centres (CDCC) have been selected in order to measure the envelope air permeability, the flow rate of mechanical ventilation systems and indoor and outdoor temperature. The data obtained was used as input to the computer code CONTAM for ventilation simulations. The results obtained were compared with direct measurements of ventilation flow from short term measurements with CO2 tracer gas and medium term measurements with perfluorocarbon tracer (PFT) gas decay method. After validation, in order to analyse the main parameters that affect ventilation, the model was used to predict the ventilation rates for a wide range of conditions. The purpose of this assessment was to find the best practices to improve natural ventilation. A simple analytical method to predict the ventilation flow rate of rooms is also presented. The method is based on the estimation of wind effect on the room through the evaluation of an average factor and on the assessment of relevant cross section of gaps and openings combined in series or in parallel. It is shown that it may be applied with acceptable accuracy for this type of buildings when ventilation is due essentially to wind action.

Numerical simulation of blast effects on fibre grout RC panels

The main purpose of the present dissertation is the simulation of the response of fibre grout strengthened RC panels when subjected to blast effects using the Applied Element Method, in order to validate and verify its applicability. Therefore, four experimental models, three of which were strengthened with a cement-based grout, each reinforced by one type of steel reinforcement, were tested against blast effects. After the calibration of the experimental set-up, it was possible to obtain and compare the response to the blast effects of the model without strengthening (reference model), and a fibre grout strengthened RC panel (strengthened model). Afterwards, a numerical model of the reference model was created in the commercial software Extreme Loading for Structures, which is based on the Applied Element Method, and calibrated to the obtained experimental results, namely to the residual displacement obtained by the experimental monitoring system. With the calibration verified, it is possible to assume that the numerical model correctly predicts the response of fibre grout RC panels when subjected to blast effects. In order to verify this assumption, the strengthened model was modelled and subjected to the blast effects of the corresponding experimental set-up. The comparison between the residual and maximum displacements and the bottom surface’s cracking obtained in the experimental and the numerical tests yields a difference of 4 % for the maximum displacements of the reference model, and a difference of 4 and 10 % for the residual and maximum displacements of the strengthened model, respectively. Additionally, the cracking on the bottom surface of the models was similar in both methods. Therefore, one can conclude that the Applied ElementMethod can correctly predict and simulate the response of fibre grout strengthened RC panels when subjected to blast effects.

Numerical simulation of galvanized rebars pullout

The usage of rebars in construction is the most common method for reinforcing plain concrete and thus bridging the tensile stresses along the concrete crack surfaces. Usually design codes for modelling the bond behaviour of rebars and concrete suggest a local bond stress – slip relationship that comprises distinct reinforcement mechanisms, such as adhesion, friction and mechanical anchorage. In this work, numerical simulations of pullout tests were performed using the finite element method framework. The interaction between rebar and concrete was modelled using cohesive elements. Distinct local bond laws were used and compared with ones proposed by the Model Code 2010. Finally an attempt was made to model the geometry of the rebar ribs in conjunction with a material damaged plasticity model for concrete.

Experimental and numerical approaches for structural assessment in new footbridge designs (SFRSCC–GFPR hybrid structure)

Within the civil engineering field, the use of the Finite Element Method has acquired a significant importance, since numerical simulations have been employed in a broad field, which encloses the design, analysis and prediction of the structural behaviour of constructions and infrastructures. Nevertheless, these mathematical simulations can only be useful if all the mechanical properties of the materials, boundary conditions and damages are properly modelled. Therefore, it is required not only experimental data (static and/or dynamic tests) to provide references parameters, but also robust calibration methods able to model damage or other special structural conditions. The present paper addresses the model calibration of a footbridge bridge tested with static loads and ambient vibrations. Damage assessment was also carried out based on a hybrid numerical procedure, which combines discrete damage functions with sets of piecewise linear damage functions. Results from the model calibration shows that the model reproduces with good accuracy the experimental behaviour of the bridge.

Numerical simulation of the punching shear behaviour of self-compacting fibre reinforced flat slabs

This paper presents the numerical simulations of the punching behaviour of centrally loaded steel fibre reinforced self-compacting concrete (SFRSCC) flat slabs. Eight half scaled slabs reinforced with different content of hooked-end steel fibres (0, 60, 75 and 90 kg/m3) and concrete strengths of 50 and 70 MPa were tested and numerically modelled. Moreover, a total of 54 three-point bending tests were carried out to assess the post-cracking flexural tensile strength. All the slabs had a relatively high conventional flexural reinforcement in order to promote the occurrence of punching failure mode. Neither of the slabs had any type of specific shear reinforcement rather than the contribution of the steel fibres. The numerical simulations were performed according to the Reissner-Mindlin theory under the finite element method framework. Regarding the classic formulation of the Reissner-Mindlin theory, in order to simulate the progressive damage induced by cracking, the shell element is discretized into layers, being assumed a plane stress state in each layer. The numerical results are, then, compared with the experimental ones and it is possible to notice that they accurately predict the experimental force-deflection relationship. The type of failure observed experimentally was also predicted in the numerical simulations.