A three dimensional finite volume approach to the thermo-mechanical modelling of the shape casting of metals

This paper presents a three dimensional, thermos-mechanical modelling approach to the cooling and solidification phases associated with the shape casting of metals ei. Die, sand and investment casting. Novel vortex-based Finite Volume (FV) methods are described and employed with regard to the small strain, non-linear Computational Solid Mechanics (CSM) capabilities required to model shape casting. The CSM capabilities include the non-linear material phenomena of creep and thermo-elasto-visco-plasticity at high temperatures and thermo-elasto-visco-plasticity at low temperatures and also multi body deformable contact with which can occur between the metal casting of the mould. The vortex-based FV methods, which can be readily applied to unstructured meshes, are included within a comprehensive FV modelling framework, PHYSICA. The additional heat transfer, by conduction and convection, filling, porosity and solidification algorithms existing within PHYSICA for the complete modelling of all shape casting process employ cell-centred FV methods. The termo-mechanical coupling is performed in a staggered incremental fashion, which addresses the possible gap formation between the component and the mould, and is ultimately validated against a variety of shape casting benchmarks.

Evaluation of thermo-physical properties and drying kinetics of carrots in a convective hot air drying system

This paper presents an experimental study on the evolution of carrot properties along convective drying by hot air at different temperatures (50ºC, 60ºC and 70ºC). The thermo-physical properties calculated were: specific heat, thermal conductivity, diffusivity, enthalpy, heat and mass transfer coefficients. Furthermore, the data of drying kinetics were treated and adjusted according to the three empirical models: Page, Henderson & Pabis and Logarithmic. The sorption isotherms were also determined and fitted using the GAB model. The results showed that, generally, the thermo-physical properties presented a decline during the drying process, and the decrease was faster for the temperature of 70ºC. It was possible to verify that the Page model presented the best prediction ability for the representation of kinetics of the drying process. The GAB model used to fit the sorption isotherms showed a good prediction capacity and, at a given water activity, despite some variations, the amount of water sorbed increased with the decrease of drying temperature.

Le processus de thermo-érosion du pergélisol dans la zone de pergélisol continu

Les polygones à coin de glace sont très répandus dans la zone du pergélisol continu. Lorsque le ruissellement d’eau de fonte nivale s’infiltre de façon concentrée dans une cavité, il peut initier le processus de thermo-érosion du pergélisol (notamment des coins de glace) pouvant mener à la formation de ravins. Dans la vallée de Qalikturvik sur l’Ile Bylot (NU, Canada), le développement de ravins de thermo-érosion dans un milieu de polygones à coins de glace entraîne comme impact : i. la réorganisation des réseaux de drainage impliquant un assèchement des milieux humides en marge des chenaux d’érosion, ii. des variations dans le régime thermique et de l’humidité de proche-surface et iii. la prise en charge et le déplacement des sédiments vers l’extérieur du bassin-versant. L’objectif de cette thèse vise à approfondir les connaissances géomorphologiques propres au ravinement par thermo-érosion, d’examiner, caractériser et quantifier les impacts du ravinement (tel que sus-mentionné en i. ii. iii.) et le rôle de celui-ci dans une optique d’évolution du paysage périglaciaire à l’échelle temporelle de l’année à la décennie. Les ravins sont dynamiques : un ravin en particulier déclenché en 1999 et étudié depuis s’érodait à une vitesse de 38 à 50 m/a durant sa première décennie d’existence, pour atteindre une longueur totale de ~750 m et une surface érodée de ~25 000 m² en 2009. Des puits sont localisés près des zones de ravinement actives ; des levées alluviale, mares et polygones effondrés dans les zones stabilisées post-perturbation. Sur la terrasse de polygones recouvrant le plancher de la vallée au site à l’étude, 35 ravins furent identifiés et 1401 polygones furent perturbés avec 200 000 m³ de sols transportés. Une amélioration du drainage, une dégradation de la capacité de rétention de l’humidité, une transition d’un écoulement de ruissellement vers un écoulement canalisé caractérise les aires ravinées et leurs environs. Les polygones intacts sont homogènes d’un à l’autre et dans leurs centres ; les polygones perturbés ont une réponse hétérogène (flore, humidité et régime thermique). Les milieux érodés hétérogènes succèdent aux milieux homogènes et deviennent le nouvel état d’équilibre pour plusieurs décennies.

The SHIVir Survey: A Dynamical Catalogue of Virgo Cluster Galaxies and their Scaling Relations

By virtue of its proximity and richness, the Virgo galaxy cluster is a perfect testing ground to expand our understanding of structure formation in the Universe. Here, we present a comprehensive dynamical catalogue based on 190 Virgo cluster galaxies (VCGs) in the "Spectroscopy and H-band Imaging of the Virgo cluster" (SHIVir) survey, including kinematics and dynamical masses. Spectroscopy collected over a multi-year campaign on 4-8m telescopes was joined with optical and near-infrared imaging to create a cosmologically-representative overview of parameter distributions and scaling relations describing galaxy evolution in a rich cluster environment. The use of long-slit spectroscopy has allowed the extraction and systematic analysis of resolved kinematic profiles: Halpha rotation curves for late-type galaxies (LTGs), and velocity dispersion profiles for early-type galaxies (ETGs). The latter are shown to span a wide range of profile shapes which correlate with structural, morphological, and photometric parameters. A study of the distributions of surface brightnesses and circular velocities for ETGs and LTGs considered separately show them all to be strongly bimodal, hinting at the existence of dynamically unstable modes where the baryon and dark matter fractions may be comparable within the inner regions of galaxies. Both our Tully-Fisher relation for LTGs and Fundamental Plane analysis for ETGs exhibit the smallest scatter when a velocity metric probing the galaxy at larger radii (where the baryonic fraction becomes sub-dominant) is used: rotational velocity measured in the outer disc at the 23.5 i-mag arcsec^{-2} level, and velocity dispersion measured within an aperture of 2 effective radii, respectively. Dynamical estimates for gas-poor and gas-rich VCGs are merged into a joint analysis of the stellar-to-total mass relation (STMR), stellar TFR, and Mass-Size relation. These relations are all found to contain strong bimodalities or dichotomies between the ETG and LTG samples, alluding to a "mixed scenario'' evolutionary sequence between morphological/dynamical classes that involves both quenching and dry mergers. The unmistakable differentiation between these two galaxy classes appears robust against different classification schemes, and supports the notion that they are driven by different evolutionary histories. Future observations using integral field spectroscopy and including lower-mass galaxies should solidify this hypothesis.

José Rodrigues Santos de Sousa Ramos (1948-2007): The diversity of dynamical systems

José Rodrigues Santos de Sousa Ramos, mathematician of great merit, passed away on January 1st, 2007, in Lisbon. He was buried in Sobral da Adiça. His death was a huge loss for the development of mathematics in Portugal. The course of time will increase the dimension of this loss. Therefore, we dedicated this theme issue on Dynamical Systems to recall his memory and underline his work. We never forget you.

Dynamical modeling of biogas production supported by an experimental anaerobic digestion process

This paper presents the development of a combined experimental and numerical approach to study the anaerobic digestion of both the wastes produced in a biorefinery using yeast for biodiesel production and the wastes generated in the preceding microbial biomass production. The experimental results show that it is possible to valorise through anaerobic digestion all the tested residues. In the implementation of the numerical model for anaerobic digestion, a procedure for the identification of its parameters needs to be developed. A hybrid search Genetic Algorithm was used, followed by a direct search method. In order to test the procedure for estimation of parameters, first noise-free data was considered and a critical analysis of the results obtain so far was undertaken. As a demonstration of its application, the procedure was applied to experimental data.

A multiscale multisurface constitutive model for the thermo-plastic behavior of polyethylene

We present a multiscale model bridging length and time scales from molecular to continuum levels with the objective of predicting the yield behavior of amorphous glassy polyethylene (PE). Constitutive pa- rameters are obtained from molecular dynamics (MD) simulations, decreasing the requirement for ad- hoc experiments. Consequently, we achieve: (1) the identification of multisurface yield functions; (2) the high strain rate involved in MD simulations is upscaled to continuum via quasi-static simulations. Validation demonstrates that the entire multisurface yield functions can be scaled to quasi-static rates where the yield stresses are possibly predicted by a proposed scaling law; (3) a hierarchical multiscale model is constructed to predict temperature and strain rate dependent yield strength of the PE.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.

On the zeros of the Abelian integrals for a class of Liénard systems

A planar polynomial differential system has a finite number of limit cycles. However, finding the upper bound of the number of limit cycles is an open problem for the general nonlinear dynamical systems. In this paper, we investigated a class of Liénard systems of the form x'=y, y'=f(x)+y g(x) with deg f=5 and deg g=4. We proved that the related elliptic integrals of the Liénard systems have at most three zeros including multiple zeros, which implies that the number of limit cycles bifurcated from the periodic orbits of the unperturbed system is less than or equal to 3.

Implicit difference approximation for the time fractional diffusion equation

In this paper, we consider a time fractional diffusion equation on a finite domain. The equation is obtained from the standard diffusion equation by replacing the first-order time derivative by a fractional derivative (of order $0<\alpha<1$ ). We propose a computationally effective implicit difference approximation to solve the time fractional diffusion equation. Stability and convergence of the method are discussed. We prove that the implicit difference approximation (IDA) is unconditionally stable, and the IDA is convergent with $O(\tau+h^2)$, where $\tau$ and $h$ are time and space steps, respectively. Some numerical examples are presented to show the application of the present technique.