Annular flow of viscoelastic fluids: Analytical and numerical solutions

This work provides analytical and numerical solutions for the linear, quadratic and exponential Phan–Thien–Tanner (PTT) viscoelastic models, for axial and helical annular fully-developed flows under no slip and slip boundary conditions, the latter given by the linear and nonlinear Navier slip laws. The rheology of the three PTT model functions is discussed together with the influence of the slip velocity upon the flow velocity and stress fields. For the linear PTT model, full analytical solutions for the inverse problem (unknown velocity) are devised for the linear Navier slip law and two different slip exponents. For the linear PTT model with other values of the slip exponent and for the quadratic PTT model, the polynomial equation for the radial location (β) of the null shear stress must be solved numerically. For both models, the solution of the direct problem is given by an iterative procedure involving three nonlinear equations, one for β, other for the pressure gradient and another for the torque per unit length. For the exponential PTT model we devise a numerical procedure that can easily compute the numerical solution of the pure axial flow problem

Singularities of the dynamical structure factors of the spin-1/2 XXX chain at finite magnetic field

We study the longitudinal and transverse spin dynamical structure factors of the spin-1/2 XXX chain at finite magnetic field h, focusing in particular on the singularities at excitation energies in the vicinity of the lower thresholds. While the static properties of the model can be studied within a Fermi-liquid like description in terms of pseudoparticles, our derivation of the dynamical properties relies on the introduction of a form of the ‘pseudofermion dynamical theory’ (PDT) of the 1D Hubbard model suitably modified for the spin-only XXX chain and other models with two pseudoparticle Fermi points. Specifically, we derive the exact momentum and spin-density dependences of the exponents ζτ(k) controlling the singularities for both the longitudinal  and transverse (τ = t) dynamical structure factors for the whole momentum range  , in the thermodynamic limit. This requires the numerical solution of the integral equations that define the phase shifts in these exponents expressions. We discuss the relation to neutron scattering and suggest new experiments on spin-chain compounds using a carefully oriented crystal to test our predictions.

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

Exploring the potentialities of a new type of CFRP laminate for the simultaneous flexural and shear strengthening of RC beams: numerical research

This study aims to develop an innovative carbon fibre reinforced polymer (CFRP) laminate with a U configuration to address strengthening interventions, where the increment of both flexural and shear capacity of reinforced concrete (RC) elements is required. This strengthening solution combines the near surface mounted (NSM) and embedded through section (ETS) techniques in the same application, since these techniques have already evidenced high performance on flexural and shear strengthening of RC beams using FRP systems, respectively. In fact, the proposed hybrid technique aims to mobilize the advantages provided by these two strengthening techniques by using an innovative CFRP laminate. The strengthening efficacy of this new hybrid NSM/ETS technique was numerically assessed and compared to the corresponding efficiency of NSM and ETS techniques applied separately for the flexural and shear strengthening of RC beams, respectively. The numerical models are described and the main relevant results are presented and discussed.

Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed.

Local non-negative initial data scalar characterization of the Kerr solution

For any vacuum initial data set, we define a local, non-negative scalar quantity which vanishes at every point of the data hypersurface if and only if the data are Kerr initial data. Our scalar quantity only depends on the quantities used to construct the vacuum initial data set which are the Riemannian metric defined on the initial data hypersurface and a symmetric tensor which plays the role of the second fundamental form of the embedded initial data hypersurface. The dependency is algorithmic in the sense that given the initial data one can compute the scalar quantity by algebraic and differential manipulations, being thus suitable for an implementation in a numerical code. The scalar could also be useful in studies of the non-linear stability of the Kerr solution because it serves to measure the deviation of a vacuum initial data set from the Kerr initial data in a local and algorithmic way.

Numerical simulation of the reflow soldering process

Dissertação de mestrado integrado em Mechanical Engineering

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.

Phenanthroimidazoles as fluorimetric and colourimetric chemosensors in aqueous solution

2,4,5-Triaryl-imidazoles are versatile compounds with application in medicine, due to their biological activity, and materials sciences, for their interesting optical properties. These properties can be tuned by careful selection of substituents at positions 2, 4 and 5: replacement of the aryl group by an heterocyclic group results in larger π-conjugated systems with improved optical properties for application in nonlinear optics, OLEDs, DNA intercalators, and chemosensors. Moreover, it is expected that introducing more conjugation and rigidity into the resulting system will further improve its properties. The development of chromo/fluorescent probes that are capable of detecting ions with high sensitivity and selectivity in aqueous media is currently a topic of strong interest and the design of heteroditopic receptors that contain two or more different binding sites for the simultaneous complexation of cationic and anionic guests is a emerging field of supramolecular chemistry. In this communication, we report the synthesis of new phenanthroimidazoles substituted at position 2 with arylthienyl or arylfuryl moieties possessing substituents of different electronic character, in order to tune the chromo/fluoro response in the presence of relevant anions and metal cations. Their photophysical properties and chemosensory ability were studied in acetonitrile and mixtures of acetonitrile and water, and selective detection of cyanide was achieved in aqueous mixtures for some of the derivatives.

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.

High performance fiber reinforced concrete for the shear reinforcement: experimental and numerical research

High performance fiber reinforced concrete (HPFRC) is developing rapidly to a modern structural material with unique rheological and mechanical characteristics. Despite applying several methodologies to achieve self15 compacting requirements, some doubts still remain regarding the most convenient strategy for developing a HPFRC. In the present study, an innovative mix design method is proposed for the development of high17 performance concrete reinforced with a relatively high dosage of steel fibers. The material properties of the developed concrete are assessed, and the concrete structural behavior is characterized under compressive, flexural and shear loading. This study better clarifies the significant contribution of fibers for shear resistance of concrete elements. This paper further discusses a FEM-based simulation, aiming to address the possibility of calibrating the constitutive model parameters related to fracture modes I and II.

Flexural strengthening of RC beams using hybrid composite plate (HCP): experimental and analytical study

Hybrid Composite Plate (HCP) is a reliable recently proposed retrofitting solution for concrete structures, which is composed of a strain hardening cementitious composite (SHCC) plate reinforced with Carbon Fibre Reinforced Polymer (CFRP). This system benefits from the synergetic advantages of these two composites, namely the high ductility of SHCC and the high tensile strength of CFRPs. In the materialstructural of HCP, the ultra-ductile SHCC plate acts as a suitable medium for stress transfer between CFRP laminates (bonded into the pre-sawn grooves executed on the SHCC plate) and the concrete substrate by means of a connection system made by either chemical anchors, adhesive, or a combination thereof. In comparison with traditional applications of FRP systems, HCP is a retrofitting solution that (i) is less susceptible to the detrimental effect of the lack of strength and soundness of the concrete cover in the strengthening effectiveness; (ii) assures higher durability for the strengthened elements and higher protection to the FRP component in terms of high temperatures and vandalism; and (iii) delays, or even, prevents detachment of concrete substrate. This paper describes the experimental program carried out, and presents and discusses the relevant results obtained on the assessment of the performance of HCP strengthened reinforced concrete (RC) beams subjected to flexural loading. Moreover, an analytical approach to estimate the ultimate flexural capacity of these beams is presented, which was complemented with a numerical strategy for predicting their load-deflection behaviour. By attaching HCP to the beams’ soffit, a significant increase in the flexural capacity at service, at yield initiation of the tension steel bars and at failure of the beams can be achieved, while satisfactory deflection ductility is assured and a high tensile capacity of the CFRP laminates is mobilized. Both analytical and numerical approaches have predicted with satisfactory agreement, the load-deflection response of the reference beam and the strengthened ones tested experimentally.

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