Identification of memory kernels in thermo-viscoelasticity

This paper discusses the possibility of identifying the shear and structure relaxation kernels in glassy materials by means of a single, simple and non-intrusive experiment. The material should be thermorheologically simple and the kernels expressed in the form of Prony series. The experiment considered consists in measuring the thickness variations over time of a flat sample cooled symmetrically from both side from a temperature above the glass transition temperature down to room temperature. The comparison of experimental observations with theoretically calculated responses allows the identification of the coefficients of the Prony series for the shear and structure relaxation kernels using a least-square type method. This paper illustrates the success of the method with `artificially' created experimental observations and with up to two exponential terms in the Prony series for the shear and structure relaxation kernels.Comment: 4 Pages, Proceedings of ESAFORM 200

Fixed Point Iterative Schemes for Initial Shape Identification

The question of interest in the present study is ”given a work-piece subject to prescribed loads, how to define its initial shape such that the work-piece matches a prescribed geometry after deformation?”. This question is particularly relevant in forming processes where the tolerated mismatch between the deformed and desired geometries may be lower than a Micron. To tackle this optimal shape design problem, a range of fixed point iterative schemes is proposed, i.e. the next initial shape is deduced from the previous one by subtracting the error in the previous final shape possibly corrected by an additional term. The required form of this additional, corrective term is revealed through a convergence analysis of the schemes. The schemes are applied to a test problem and their performance compared. The problem consists in designing the hole in a membrane such that its contour matches a prescribed shape when the membrane is stretched by a given load

Inertial and dimensional effects on the instability of a thin film

We consider here the effects of inertia on the instability of a flat liquid film under the effects of capillary and intermolecular forces (van der Waals interaction). Firstly, we perform the linear stability analysis within the long wave approximation, which shows that the inclusion of inertia does not produce new regions of instability other than the one previously known from the usual lubrication case. The wavelength, $\lambda_m$, corresponding to he maximum growth, $\omega_m$, and the critical (marginal) wavelength do not change at all. The most affected feature of the instability under an increase of the Laplace number is the noticeable decrease of the growth rates of the unstable modes. In order to put in evidence the effects of the bidimensional aspects of the flow (neglected in the long wave approximation), we also calculate the dispersion relation of the instability from the linearized version of the complete Navier-Stokes (N-S) equation. Unlike the long wave approximation, the bidimensional model shows that $\lambda_m$ can vary significantly with inertia when the aspect ratio of the film is not sufficiently small. We also perform numerical simulations of the nonlinear N-S equations and analyze to which extent the linear predictions can be applied depending on both the amount of inertia involved and the aspect ratio of the film

Efficient and accurate time adaptive multigrid simulations of droplet spreading

An efficient full approximation storage (FAS) Multigrid algorithm is used to solve a range of droplet spreading flows modelled as a coupled set of non-linear lubrication equations. The algorithm is fully implicit and has embedded within it an adaptive time-stepping scheme that enables the same to be optimized in a controlled manner subject to a specific error tolerance. The method is first validated against a range of analytical and existing numerical predictions commensurate with droplet spreading and then used to simulate a series of new, three-dimensional flows consisting of droplet motion on substrates containing topographic and wetting heterogeneities. The latter are of particular interest and reveal how droplets can be made to spread preferentially on substrates owing to an interplay between different topographic and surface wetting characteristics

Gravity-driven flow of continuous thin liquid films on non-porous substrates with topography

A range of two- and three-dimensional problems is explored featuring the gravity-driven flow of a continuous thin liquid film over a non-porous inclined flat surface containing well-defined topography. These are analysed principally within the framework of the lubrication approximation, where accurate numerical solution of the governing nonlinear equations is achieved using an efficient multigrid solver. Results for flow over one-dimensional steep-sided topographies are shown to be in very good agreement with previously reported data. The accuracy of the lubrication approximation in the context of such topographies is assessed and quantified by comparison with finite element solutions of the full Navier–Stokes equations, and results support the consensus that lubrication theory provides an accurate description of these flows even when its inherent assumptions are not strictly satisfied. The Navier–Stokes solutions also illustrate the effect of inertia on the capillary ridge/trough and the two-dimensional flow structures caused by steep topography. Solutions obtained for flow over localized topography are shown to be in excellent agreement with the recent experimental results of Decré & Baret (2003) for the motion of thin water films over finite trenches. The spread of the ‘bow wave’, as measured by the positions of spanwise local extrema in free-surface height, is shown to be well-represented both upstream and downstream of the topography by an inverse hyperbolic cosine function. An explanation, in terms of local flow rate, is given for the presence of the ‘downstream surge’ following square trenches, and its evolution as trench aspect ratio is increased is discussed. Unlike the upstream capillary ridge, this feature cannot be completely suppressed by increasing the normal component of gravity. The linearity of free-surface response to topographies is explored by superposition of the free surfaces corresponding to two ‘equal-but-opposite’ topographies. Results confirm the findings of Decré & Baret (2003) that, under the conditions considered, the responses behave in a near-linear fashion

Flow of evaporating, gravity-driven thin liquid films over topography

The effect of topography on the free surface and solvent concentration profiles of an evaporating thin film of liquid flowing down an inclined plane is considered. The liquid is assumed to be composed of a resin dissolved in a volatile solvent with the associated solvent concentration equation derived on the basis of the well-mixed approximation. The dynamics of the film is formulated as a lubrication approximation and the effect of a composition-dependent viscosity is included in the model. The resulting time-dependent, nonlinear, coupled set of governing equations is solved using a full approximation storage multigrid method. The approach is first validated against a closed-form analytical solution for the case of a gravity-driven, evaporating thin film flowing down a flat substrate. Analysis of the results for a range of topography shapes reveal that although a full-width, spanwise topography such as a step-up or a step-down does not affect the composition of the film, the same is no longer true for the case of localized topography, such as a peak or a trough, for which clear nonuniformities of the solvent concentration profile can be observed in the wake of the topography

Slug Self-Propulsion in a Capillary Tube Mathematical Modeling and Numerical Simulation

A composite droplet made of two miscible fluids in a narrow tube generally moves under the action of capillarity until complete mixture is attained. This physical situation is analysed here on a combined theoretical and numerical analysis. The mathematical framework consists of the two-phase flow phase-field equation set, an advection-diffusion chemical concentration equation, and closure relationships relating the surface tensions to the chemical concentration. The numerical framework is composed of the COMSOL Laminar two-phase flow phase-field method coupled with an advection-diffusion chemical concentration equation. Through transient studies, we show that the penetrating length of the bidroplet system into the capillary tube is linear at early-time regime and exponential at late-time regime. Through parametric studies, we show that the rate of penetration of the bidroplet system into the capillary tube is proportional to a time-dependent exponential function. We also show that this speed obeys the Poiseuille law at the early-time regime. A series of position, speed-versus-property graphs are included to support the analysis. Finally, the overall results are contrasted with available experimental data, grouped together to settle a general mathematical description of the phenomenon, and explained and concluded on this basis