Numerical Method to Predict Void Formation during The Liquid Composite Molding Process

Void formation during the injection phase of the liquid composite molding process can be explained as a consequence of the non-uniformity of the flow front progression. This is due to the dual porosity within the fiber perform (spacing between the fiber tows is much larger than between the fibers within in a tow) and therefore the best explanation can be provided by a mesolevel analysis, where the characteristic dimension is given by the fiber tow diameter of the order of millimeters. In mesolevel analysis, liquid impregnation along two different scales; inside fiber tows and within the open spaces between the fiber tows must be considered and the coupling between the flow regimes must be addressed. In such cases, it is extremely important to account correctly for the surface tension effects, which can be modeled as capillary pressure applied at the flow front. Numerical implementation of such boundary conditions leads to illposing of the problem, in terms of the weak classical as well as stabilized formulation. As a consequence, there is an error in mass conservation accumulated especially along the free flow front. A numerical procedure was formulated and is implemented in an existing Free Boundary Program to reduce this error significantly.

Numerical methods and tangible interfaces for pollutant dispersion simulation

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia do Ambiente, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

MASS CONSERVATION ENHANCEMENT OF FREE BOUNDARY FLOWS IN COMPOSITES MANUFACTURING

Undesirable void formation during the injection phase of the liquid composite molding process can be understood as a consequence of the non-uniformity of the flow front progression, caused by the dual porosity of the fiber perform. Therefore the best examination of the void formation physics can be provided by a mesolevel analysis, where the characteristic dimension is given by the fiber tow diameter. In mesolevel analysis, liquid impregnation along two different scales; inside fiber tows and within the spaces between them; must be considered and the coupling between these flow regimes must be addressed. In such case, it is extremely important to account correctly for the surface tension effects, which can be modeled as capillary pressure applied at the flow front. When continues Galerkin method is used, exploiting elements with velocity components and pressure as nodal variables, strong numerical implementation of such boundary conditions leads to ill-posing of the problem, in terms of the weak classical as well as stabilized formulation. As a consequence, there is an error in mass conservation accumulated especially along the free flow front. This article presents a numerical procedure, which was formulated and implemented in the existing Free Boundary Program in order to significantly reduce this error.

Mass Conservation Enhancement of Free Boundary Mesolevel Flows during LCM Processes of Composites Manufacturing

Undesirable void formation during the injection phase of the liquid composite moulding process can be understood as a consequence of the non-uniformity of the flow front progression, caused by the dual porosity of the fibre perform. Therefore the best examination of the void formation physics can be provided by a mesolevel analysis, where the characteristic dimension is given by the fibre tow diameter. In mesolevel analysis, liquid impregnation along two different scales; inside fibre tows and within the open spaces between them; must be considered and the coupling between these flow regimes must be addressed. In such case, it is extremely important to account correctly for the surface tension effects, which can be modelled as capillary pressure applied at the flow front. Numerical implementation of such boundary conditions leads to ill-posing of the problem, in terms of the weak classical as well as stabilized formulation. As a consequence, there is an error in mass conservation accumulated especially along the free flow front. This contribution presents a numerical procedure, which was formulated and implemented in the existing Free Boundary Program in order to significantly reduce this error.

Geomorphic dam-break flows. Part II: numerical simulation

Proceedings of the Institution of Civil Engineers - Water Management 163 Issue WM6

Framework development for providing accessibility to qualitative spatial calculi

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies.

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

Validation of NASA rotor 67 with openFOAM's transonic density-based solver

Dissertação para obtenção do Grau de Mestre em Engenharia Mecânica

Predictive analytical and numerical modeling for orthogonal cutting

In this thesis, a predictive analytical and numerical modeling approach for the orthogonal cutting process is proposed to calculate temperature distributions and subsequently, forces and stress distributions. The models proposed include a constitutive model for the material being cut based on the work of Weber, a model for the shear plane based on Merchants model, a model describing the contribution of friction based on Zorev’s approach, a model for the effect of wear on the tool based on the work of Waldorf, and a thermal model based on the works of Komanduri and Hou, with a fraction heat partition for a non-uniform distribution of the heat in the interfaces, but extended to encompass a set of contributions to the global temperature rise of chip, tool and work piece. The models proposed in this work, try to avoid from experimental based values or expressions, and simplifying assumptions or suppositions, as much as possible. On a thermo-physical point of view, the results were affected not only by the mechanical or cutting parameters chosen, but also by their coupling effects, instead of the simplifying way of modeling which is to contemplate only the direct effect of the variation of a parameter. The implementation of these models was performed using the MATLAB environment. Since it was possible to find in the literature all the parameters for AISI 1045 and AISI O2, these materials were used to run the simulations in order to avoid arbitrary assumption.

Mathematical and numerical methods for inverse problems using fundamental solutions techniques

Fundação para a Ciência e a Tecnologia - SFRH/BD/27914/2006

Numerical prediction of diffusion and electric field-induced iron nanoparticle transport

Zero valent iron nanoparticles (nZVI) are considered very promising for the remediation of contaminated soils and groundwaters. However, an important issue related to their limited mobility remains unsolved. Direct current can be used to enhance the nanoparticles transport, based on the same principles of electrokinetic remediation. In this work, a generalized physicochemical model was developed and solved numerically to describe the nZVI transport through porous media under electric field, and with different electrolytes (with different ionic strengths). The model consists of the Nernst–Planck coupled system of equations, which accounts for the mass balance of ionic species in a fluid medium, when both the diffusion and electromigration of the ions are considered. The diffusion and electrophoretic transport of the negatively charged nZVI particles were also considered in the system. The contribution of electroosmotic flow to the overall mass transport was included in the model for all cases. The nZVI effective mobility values in the porous medium are very low (10−7–10−4 cm2 V−1 s−1), due to the counterbalance between the positive electroosmotic flow and the electrophoretic transport of the negatively charged nanoparticles. The higher the nZVI concentration is in the matrix, the higher the aggregation; therefore, low concentration of nZVI suspensions must be used for successful field application.

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

Tracking error of exchange-traded funds: Evidence from the UK

This paper is mainly concerned with the tracking accuracy of Exchange Traded Funds (ETFs) listed on the London Stock Exchange (LSE) but also evaluates their performance and pricing efficiency. The findings show that ETFs offer virtually the same return but exhibit higher volatility than their benchmark. It seems that the pricing efficiency, which should come from the creation and redemption process, does not fully hold as equity ETFs show consistent price premiums. The tracking error of the funds is generally small and is decreasing over time. The risk of the ETF, daily price volatility and the total expense ratio explain a large part of the tracking error. Trading volume, fund size, bid-ask spread and average price premium or discount did not have an impact on the tracking error. Finally, it is concluded that market volatility and the tracking error are positively correlated.

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.

Damage propagation in composite Materials meso-mechanical models

Composite materials have a complex behavior, which is difficult to predict under different types of loads. In the course of this dissertation a methodology was developed to predict failure and damage propagation of composite material specimens. This methodology uses finite element numerical models created with Ansys and Matlab softwares. The methodology is able to perform an incremental-iterative analysis, which increases, gradually, the load applied to the specimen. Several structural failure phenomena are considered, such as fiber and/or matrix failure, delamination or shear plasticity. Failure criteria based on element stresses were implemented and a procedure to reduce the stiffness of the failed elements was prepared. The material used in this dissertation consist of a spread tow carbon fabric with a 0°/90° arrangement and the main numerical model analyzed is a 26-plies specimen under compression loads. Numerical results were compared with the results of specimens tested experimentally, whose mechanical properties are unknown, knowing only the geometry of the specimen. The material properties of the numerical model were adjusted in the course of this dissertation, in order to find the lowest difference between the numerical and experimental results with an error lower than 5% (it was performed the numerical model identification based on the experimental results).