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

Fractional pennes' bioheat equation: theoretical and numerical studies

In this work we provide a new mathematical model for the Pennes’ bioheat equation, assuming a fractional time derivative of single order. Alternative versions of the bioheat equation are studied and discussed, to take into account the temperature-dependent variability in the tissue perfusion, and both finite and infinite speed of heat propagation. The proposed bioheat model is solved numerically using an implicit finite difference scheme that we prove to be convergent and stable. The numerical method proposed can be applied to general reaction diffusion equations, with a variable diffusion coefficient. The results obtained with the single order fractional model, are compared with the original models that use classical derivatives.

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

Improvign mixing processes with the help of numerical tools

In several industrial applications, highly complex behaviour materials are used together with intricate mixing processes, which difficult the achievement of the desired properties for the produced materials. This is the case of the well-known dispersion of nano-sized fillers in a melt polymer matrix, used to improve the nanocomposite mechanical and/or electrical properties. This mixing is usually performed in twin-screw extruders, that promote complex flow patterns, and, since an in loco analysis of the material evolution and mixing is difficult to perform, numerical tools can be very useful to predict the evolution and behaviour of the material. This work presents a numerical based study to improve the understanding of mixing processes. Initial numerical studies were performed with generalized Newtonian fluids, but, due to the null relaxation time that characterize this type of fluids, the assumption of viscoelastic behavior was required. Therefore, the polymer melt was rheologically characterized, and, a six mode Phan-Thien-Tanner and Giesekus models were used to fit the rheological data. These viscoelastic rheological models were used to model the process. The conclusions obtained in this work provide additional and useful data to correlate the type and intensity of the deformation history promoted to the polymer nanocomposite and the quality of the mixing obtained.

Numerical modeling and assessment of the thermal environment in an industrial premise of automotive components

Dissertação de mestrado em Engenharia Industrial

Numerical simulation of the reflow soldering process

Dissertação de mestrado integrado em Mechanical Engineering

Numerical and experimental analysis of full scale arches reinforced with GFRP materials

In this contribution, original limit analysis numerical results are presented dealing with some reinforced masonry arches tested at the University of Minho-UMinho, PT. Twelve in-scale circular masonry arches were considered, reinforced in various ways at the intrados or at the extrados. GFRP reinforcements were applied either on undamaged or on previously damaged elements, in order to assess the role of external reinforcements even in repairing interventions. The experimental results were critically discussed at the light of limit analysis predictions, based on a 3D FE heterogeneous upper bound approach. Satisfactory agreement was found between experimental evidences and the numerical results, in terms of failure mechanisms and peak load.

Experimental and numerical studies of hybrid PCM embedded in plastering mortar for enhanced thermal behaviour of buildings

This paper proposes a methodology for improvement of energy efficiency in buildings through the innovative simultaneous incorporation of three distinct phase change materials (here termed as hybrid PCM) in plastering mortars for façade walls. The thermal performance of a hybrid PCM mortar was experimentally evaluated by comparing the behaviour of a prototype test cell (including hybrid PCM plastering mortar) subjected to realistic daily temperature profiles, with the behaviour of a similar prototype test cell, in which no PCM was added. A numerical simulation model was employed (using ANSYS-FLUENT) to validate the capacity of simulating temperature evolution within the prototype containing hybrid PCM, as well as to understand the contribution of hybrid PCM to energy efficiency. Incorporation of hybrid PCM into plastering mortars was found to have the potential to significantly reduce heating/cooling temperature demands for maintaining the interior temperature within comfort levels when compared to normal mortars (without PCM), or even mortars comprising a single type of PCM.

Pulse wave velocity distribution in a cohort study: from arterial stiffness to early vascular aging

BACKGROUND: By contrast with other southern European people, north Portuguese population registers an especially high prevalence of hypertension and stroke incidence. We designed a cohort study to identify individuals presenting accelerated and premature arterial aging in the Portuguese population. METHOD: Pulse wave velocity (PWV) was measured in randomly sampled population dwellers aged 18-96 years from northern Portugal, and used as a marker of early vascular aging (EVA). Of the 3038 individuals enrolled, 2542 completed the evaluation. RESULTS: Mean PWV value for the entire population was 8.4?m/s (men: 8.6?m/s; women: 8.2?m/s; P??10?m/s). Logistic regression models indicated gender differences concerning the risk of developing large artery damage, with women having the same odds of PWV above 10?m/s 10 years later than men. CONCLUSION: The population PWV values were higher than expected in a low cardiovascular risk area (Portugal). High prevalence rates of EVA and noteworthy large artery damage in young ages were found.

Numerical methods in multibody system dynamics

"Series: Solid mechanics and its applications, vol. 226"

Relation between the behaviors of P-wave and QT dispersions in elderly patients with heart failure

OBJECTIVE: To assess the relation between P-wave and QT dispersions in elderly patients with heart failure. METHODS: Forty-seven elderly patients (75.6±6 years) with stable heart failure in NYHA functional classes II or III and with ejection fractions of 37±6% underwent body surface mapping to analyze P-wave and QT dispersions. The degree of correlation between P-wave and QT dispersions was assessed, and P-wave dispersion values in patients with QT dispersion greater than and smaller than 100 ms were compared. RESULTS: The mean values of P-wave and QT dispersions were 54±14 ms and 68±27 ms, respectively. The correlation between the 2 variables was R=0.41 (p=0.04). In patients with QT dispersion values > 100 ms, P-wave dispersion was significantly greater than in those with QT dispersion values < 100 ms (58±16 vs 53±12 ms, p=0.04 ). CONCLUSION: Our results suggest that, in elderly patients with heart failure, a correlation between the values of P-wave and QT dispersions exists. These findings may have etiopathogenic, pathophysiologic, prognostic, and therapeutic implications, which should be investigated in other studies.

Desarrollo de Algoritmos Numéricos para el Estudio de la Propagación no Lineal de Ondas en Aplicaciones Aeroespaciales y Astrofísicas. Develop to Numerical Algorithms to Study non-Linear Waves Propagation in Aerospace and Astrophysics Applications.

En este proyecto se desarrollarán algoritmos numéricos para sistemas no lineales hiperbólicos-parabólicos de ecuaciones diferenciales en derivadas parciales. Dichos sistemas tienen aplicación en propagación de ondas en ámbitos aeroespaciales y astrofísicos.Objetivos generales: 1)Desarrollo y mejora de algoritmos numéricos con la finalidad de incrementar la calidad en la simulación de propagación e interacción de ondas gasdinámicas y magnetogasdinámicas no lineales. 2)Desarrollo de códigos computacionales con la finalidad de simular flujos gasdinámicos de elevada entalpía incluyendo cambios químicos, efectos dispersivos y difusivos.3)Desarrollo de códigos computacionales con la finalidad de simular flujos magnetogasdinámicos ideales y reales.4)Aplicación de los nuevos algoritmos y códigos computacionales a la solución del flujo aerotermodinámico alrededor de cuerpos que ingresan en la atmósfera terrestre. 5)Aplicación de los nuevos algoritmos y códigos computacionales a la simulación del comportamiento dinámico no lineal de arcos magnéticos en la corona solar. 6)Desarrollo de nuevos modelos para describir el comportamiento no lineal de arcos magnéticos en la corona solar.Este proyecto presenta como objetivo principal la introducción de mejoras en algoritmos numéricos para simular la propagación e interacción de ondas no lineales en dos medios gaseosos: aquellos que no poseen carga eléctrica libre (flujos gasdinámicos) y aquellos que tienen carga eléctrica libre (flujos magnetogasdinámicos). Al mismo tiempo se desarrollarán códigos computacionales que implementen las mejoras de las técnicas numéricas.Los algoritmos numéricos se aplicarán con la finalidad de incrementar el conocimiento en tópicos de interés en la ingeniería aeroespacial como es el cálculo del flujo de calor y fuerzas aerotermodinámicas que soportan objetos que ingresan a la atmósfera terrestre y en temas de astrofísica como la propagación e interacción de ondas, tanto para la transferencia de energía como para la generación de inestabilidades en arcos magnéticos de la corona solar. Estos dos temas poseen en común las técnicas y algoritmos numéricos con los que serán tratados. Las ecuaciones gasdinámicas y magnetogasdinámicas ideales conforman sistemas hiperbólicos de ecuaciones diferenciales y pueden ser solucionados utilizando "Riemann solvers" junto con el método de volúmenes finitos (Toro 1999; Udrea 1999; LeVeque 1992 y 2005). La inclusión de efectos difusivos genera que los sistemas de ecuaciones resulten hiperbólicos-parabólicos. La contribución parabólica puede ser considerada como términos fuentes y tratada adicionalmente tanto en forma explícita como implícita (Udrea 1999; LeVeque 2005).Para analizar el flujo alrededor de cuerpos que ingresan en la atmósfera se utilizarán las ecuaciones de Navier-Stokes químicamente activas, mientras la temperatura no supere los 6000K. Para mayores temperaturas es necesario considerar efectos de ionización (Anderson, 1989). Tanto los efectos difusivos como los cambios químicos serán considerados como términos fuentes en las ecuaciones de Euler. Para tratar la propagación de ondas, transferencia de energía e inestabilidades en arcos magnéticos de la corona solar se utilizarán las ecuaciones de la magnetogasdinámica ideal y real. En este caso será también conveniente implementar términos fuente para el tratamiento de fenómenos de transporte como el flujo de calor y el de radiación. Los códigos utilizarán la técnica de volúmenes finitos, junto con esquemas "Total Variation Disminishing - TVD" sobre mallas estructuradas y no estructuradas.

T-wave alternans: desempenho clínico, limitações, metodologias de análise

Reconhecer com precisão indivíduos com maior risco imediato de morte súbita cardíaca (MSC) ainda é uma questão em aberto. A natureza fortuita dos eventos cardiovasculares agudos não parece se adequar ao conhecido modelo de indução de taquicardia/fibrilação ventricular por um gatilho em sincronia a um substrato arritmogênico estático. Quanto ao mecanismo da MSC, uma instabilidade elétrica dinâmica explicaria melhor a raridade da associação simultânea de um gatilho certo a um substrato cardíaco apropriado. Diversos estudos tentaram medir essa instabilidade elétrica cardíaca (ou um equivalente válido) em uma sequência de batimentos cardíacos no ECG. Dentre os mecanismos possíveis podemos citar o prolongamento do QT, dispersão do QT, potenciais tardios, alternância de onda T ou T-wave alternans (TWA), e turbulência da frequência cardíaca. Este artigo se atém em particular ao papel da TWA no panorama atual da estratificação de risco cardíaco. Os achados sobre TWA ainda são heterogêneos, variando de um desempenho prognóstico muito bom até um quase nulo, dependendo da população clínica observada e protocolo clínico usado. Para preencher as atuais lacunas no conhecimento sobre TWA, profissionais médicos e pesquisadores devem explorar melhor as características técnicas das diversas tecnologias disponíveis para a avaliação de TWA e atentar ao fato de que os valores de TWA respondem a diversos outros fatores, além de medicamentos. Informações sobre mecanismos celulares e subcelulares da TWA estão fora do escopo deste artigo, mas são referenciados alguns dos principais trabalhos sobre este tópico, com o intuito de auxiliar no entendimento dos conceitos e fatos cobertos neste artigo.