Lack of exponential decay in viscoelastic materials with voids

We study the stability of solutions to a coupled evolution system associated with an isotropic porous and centrosymmetric viscoelastic solid with porous dissipation (a porous elastic system with history). With the help of the method of the semigroup theory and some novel observations, we prove successfully that the condition of equal wave-speed propagation is still necessary for exponential stability of the system in the case of Dirichlet–Dirichlet boundary conditions.     

Approximating thermo-viscoelastic heating of largely strained solid rubber components

Mechanically induced viscoelastic dissipation is difficult to compute when the constitutive model is defined by history integrals. The computation of the viscous energy dissipated is in the form of a double convolution integral. In this study, we present a method to approximate the dissipation for constitutive models in history integral form that represent Maxwell-like materials. The dissipation is obtained without directly computing the double convolution integral. The approximation requires that the total stress can be separated into elastic and viscous components, and that the relaxation form of the constitutive law is defined with a Prony series. A numerical approach often taken to approximate a history integral involves interpolating the history integral’s kernel across a time step. Integration then yields finite difference equations for the evolution of the viscous stresses in time. In the case when the material is modeled with a Prony series, the form of these finite difference equations is similar to the form of the finite difference equations for a Maxwell solid. Since the dissipation rate in a Maxwell solid can be easily computed from knowledge of its viscous stress and the Prony series constants (spring-dashpot constants), we computationally investigated employing a Maxwell solid’s dissipation function to couple thermal and large strain history integral based finite element models of solid rubber components. Numerical data is provided to support this analogy and to help understand its limitations. A rubber cylinder with an imbedded steel disk is dynamically loaded, and the non-uniform heating within the cylinder is computed.     

Forced convection in electro-osmotic/Poiseuille micro-channel flows of viscoelastic fluids: fully developed flow with imposed wall heat flux

The analytical solution for heat transfer in a dynamic and thermally fully developed channel flow of the simplified Phan-Thien–Tanner fluid induced by combined electro-osmosis and pressure gradient was obtained assuming that material properties are independent of temperature. The flow forcing was quantified by an appropriate dimensionless parameter and its effect and that of all other relevant dimensionless numbers is presented and discussed. Specifically, the forced convection occurs under conditions of constant wall heat flux and the solution includes the effects of Weissenberg number, electric double layer (EDL) thickness, forcing ratio parameter, viscous dissipation as well as of Joule heating due to the electric currents and was obtained under the simplifying Debye–Hückel approximation. Generally speaking, the Joule effect is stronger than the viscous dissipation except in very narrow channels, but these fall outside the validity of the Debye–Hückel conditions. For pure electro-osmosis, viscous dissipation is restricted to the near-wall region and virtually nonexistent elsewhere, so it is irrelevant for thin electric double layers and Joule heating is more relevant. As the EDL thickens and/or the pressure gradient contribution increases, the role of viscous dissipation grows and shear-thinning effects also appear more clearly on the Nusselt number. Generally speaking, an increase in internal heating results in lower Nusselt numbers and this effect is stronger than the effect of shear-thinning, which is responsible for a slight increase in the Nusselt number.     

A Continuum Mechanical Approach to Geodesics in Shape Space

In this paper concepts from continuum mechanics are used to define geodesic paths in the space of shapes, where shapes are implicitly described as boundary contours of objects. The proposed shape metric is derived from a continuum mechanical notion of viscous dissipation. A geodesic path is defined as the family of shapes such that the total amount of viscous dissipation caused by an optimal material transport along the path is minimized. The approach can easily be generalized to shapes given as segment contours of multi-labeled images and to geodesic paths between partially occluded objects. The proposed computational framework for finding such a minimizer is based on the time discretization of a geodesic path as a sequence of pairwise matching problems, which is strictly invariant with respect to rigid body motions and ensures a 1–1 correspondence along the induced flow in shape space. When decreasing the time step size, the proposed model leads to the minimization of the actual geodesic length, where the Hessian of the pairwise matching energy reflects the chosen Riemannian metric on the underlying shape space. If the constraint of pairwise shape correspondence is replaced by the volume of the shape mismatch as a penalty functional, one obtains for decreasing time step size an optical flow term controlling the transport of the shape by the underlying motion field. The method is implemented via a level set representation of shapes, and a finite element approximation is employed as spatial discretization both for the pairwise matching deformations and for the level set representations. The numerical relaxation of the energy is performed via an efficient multi-scale procedure in space and time. Various examples for 2D and 3D shapes underline the effectiveness and robustness of the proposed approach.     

A computational study of vortex-airfoil interaction using high-order finite difference methods

Simulations of the interaction between a vortex and a NACA0012 airfoil are performed with a stable, high-order accurate (in space and time), multi-block finite difference solver for the compressible Navier-Stokes equations.We begin by computing a benchmark test case to validate the code. Next, the flow with steady inflow conditions are computed on several different grids. The resolution of the boundary layer as well as the amount of the artificial dissipation is studied to establish the necessary resolution requirements. We propose an accuracy test based on the weak imposition of the boundary conditions that does not require a grid refinement.Finally, we compute the vortex-airfoil interaction and calculate the lift and drag coefficients. It is shown that the viscous terms add the effect of detailed small scale structures to the lift and drag coefficients.     

Asymptotic and finite element approximations for heat transfer in rotating compressible flow over an infinite porous disk

The laminar boundary layer equations for the compressible flow due to the finite difference in rotation and temperature rates are solved for the case of uniform suction through the disk. The effects of viscous dissipation on the incompressible flow are taken into account for any rotation rate, whereas for a compressible fluid they are considered only for a disk rotating in a stationary fluid. For the general case, the governing equations are solved numerically using a standard finite element scheme. Series solutions are developed for those cases where the suction effect is dominant. Based on the above analytical and numerical solutions, a new asymptotic finite element scheme is presented. By using this scheme one can significantly improve the pointwise accuracy of the standard finite element scheme.     

Classification of nasal inspiratory flow shapes by attributed finite automata.

In a significant proportion of individuals, the physiologic decrease of muscle tone during sleep results in increased collapsibility of the upper respiratory airway. At peak inspiratory flow, the pharyngeal soft tissues may collapse and cause airflow limitation or even complete occlusion of the upper airway (sleep apnea). While there are plenty of methods to detect sleep apnea, only a few can be used to monitor flow limitation in sleeping individuals. Nasal prongs connected to pressure sensor provide information of the nasal airflow over time. This paper documents a method to automatically classify each nasal inspiratory pressure profile into one without flow limitation or six flow-limited ones. The recognition of the sample signals consists of three phases: preprocessing, primitive extraction, and word parsing phases. In the last one, a sequence of signal primitives is treated as a word and we test its membership in the attribute grammars constructed to the signal categories. The method gave in practical tests surprisingly high performance. Classifying 94;pc of the inspiratory profiles in agreement with the visual judgment of an expert physician, the performance of the method was considered good enough to warrant further testing in well-defined patient populations to determine the pressure profile distributions of different subject classes.     

Effects of heat source/sink,radiation and work done by deformation on flow and heat transfer of a viscoelastic fluid over a stretching sheet

This paper presents a study of the flow and heat transfer of an incompressible homogeneous second-grade fluid over a non-isothermal stretching sheet. The governing partial differential equations are converted into ordinary differential equations by a similarity transformation. The effects of viscous dissipation, work due to deformation, internal heat generation/absorption and thermal radiation are considered in the energy equation, and the variations of dimensionless surface temperature and dimensionless surface temperature gradient as well as the heat transfer characteristics with various physical parameters are graphed and tabulated. Two cases are studied, namely, (i) a sheet with prescribed surface temperature (PST case) and (ii) a sheet with prescribed heat flux (PHF case).     

Shape design optimization of an expansion step in a channel with moving boundaries for a viscous flow

A shape design optimization problem for viscous flows has been investigated in the present study. An analytical shape design sensitivity expression has been derived for a general integral functional by using the adjoint variable method and the material derivative concept of optimization. A channel flow problem with a backward facing step and adversely moving boundary wall is taken as an example. The shape profile of the expansion step, represented by a fourth-degree polynomial, is optimized in order to minimize the total viscous dissipation in the flow field. Numerical discretizations of the primary (flow) and adjoint problems are achieved by using the Galerkin FEM method. A balancing upwinding technique is also used in the equations. Numerical results are provided in various graphical forms at relatively low Reynolds numbers. It is concluded that the proposed general method of solution for shape design optimization problems is applicable to physical systems described by nonlinear equations.     

High-resolution finite difference schemes using curvilinear coordinate grids for DNS of compressible turbulent flow over wavy walls

Direct numerical simulations of compressible turbulent flow over wavy wall geometries have been carried out by solving N–S equations on general curvilinear coordinates. A 6th order WENO scheme with minimized dispersion and controllable dissipation is employed to compute the inviscid fluxes, a 4th order central difference scheme is applied to compute the viscous fluxes, and a 6th order conservative compact scheme is used for computing the geometrical metrics. An implicit LU-SGS method is used for time integration to improve computational efficiency over the explicit schemes such as the Runge–Kutta approach. The validity and applicability of the present algorithm is confirmed by comparing our results with laboratory experimental measurements and DNS results in the literature.     

Pseudospectral multi-domain method for incompressible viscous flow computation

We present a multi-domain pseudospectral method for the calculation of incompressible viscous flow. Governing equations are written in primitive variables formulation. Velocity components and pressure are discretized on the same grid of collocation points. A coupling algorithm for the Stokes problem is investigated and preliminary results are presented in the two-dimensional case.     

Patient-Specific Geometry Modeling and Mesh Generation for Simulating Obstructive Sleep Apnea Syndrome Cases by Maxillomandibular Advancement

The objective of this paper is the reconstruction of upper airway geometric models as hybrid meshes from clinically used Computed Tomography (CT) data sets in order to understand the dynamics and behaviors of the pre- and postoperative upper airway systems of Obstructive Sleep Apnea Syndrome (OSAS) patients by viscous Computational Fluid Dynamics (CFD) simulations. The selection criteria for OSAS cases studied are discussed because two reasonable pre- and postoperative upper airway models for CFD simulations may not be created for every case without a special protocol for CT scanning. The geometry extraction and manipulation methods are presented with technical barriers that must be overcome so that they can be used along with computational simulation software as a daily clinical evaluation tool. Eight cases are presented in this paper, and each case consists of pre- and postoperative configurations. The results of computational simulations of two cases are included in this paper as demonstration.     

Modified kinetic flux vector splitting (m-KFVS) method for compressible flows

To resolve many flow features accurately, like accurate capture of suction peak in subsonic flows and crisp shocks in flows with discontinuities, to minimise the loss in stagnation pressure in isentropic flows or even flow separation in viscous flows require an accurate and low dissipative numerical scheme. The first order kinetic flux vector splitting (KFVS) method has been found to be very robust but suffers from the problem of having much more numerical diffusion than required, resulting in inaccurate computation of the above flow features. However, numerical dissipation can be reduced by refining the grid or by using higher order kinetic schemes. In flows with strong shock waves, the higher order schemes require limiters, which reduce the local order of accuracy to first order, resulting in degradation of flow features in many cases. Further, these schemes require more points in the stencil and hence consume more computational time and memory. In this paper, we present a low dissipative modified KFVS (m-KFVS) method which leads to improved splitting of inviscid fluxes. The m-KFVS method captures the above flow features more accurately compared to first order KFVS and the results are comparable to second order accurate KFVS method, by still using the first order stencil.     

An optimal analysis of radiated nanomaterial flow with viscous dissipation and heat source

Microsystem Technologies - The present work explores the analytical study of nanofluid flow above a stretching medium with the heat source and viscous dissipation. Additional radiative effects are...     

Dynamical characteristics of fluid-conveying microbeams actuated by electrostatic force

In this paper, the dynamic characteristics and pull-in instability of electrostatically actuated microbeams which convey internal fluids are investigated. A theoretical model is developed by considering the elastic structure, laminar flow and electrostatic field to characterize the dynamic behavior. In addition, the energy dissipation induced by the fluid viscosity is studied through analyzing the fluid–structure interactions between the laminar fluid flow and oscillating microbeam by comprehensively considering the effects of velocity profile and fluid viscosity. The results indicate that the system is subjected to both the pull-in instability and the fluid-induced instability. It is demonstrated that as the flow velocity increases, both the static pull-in voltage and the dynamic one decrease for clamped–clamped microbeams while increase for clamped-free microbeams. It is also shown that the applied voltage and the steady flow can adjust the resonant frequency. The perturbation viscous flow caused by the vibration of microbeam is manifested to result in energy dissipation. The quality factor decreases with the increment of both the mode order and flow velocity. However, when the oscillating flow dominates, the flow velocity has no obvious effect.     

A finite volume method to solve the 3D Navier–Stokes equations on unstructured collocated meshes

A new method to solve the Navier–Stokes equations for incompressible viscous flows and the transport of a scalar quantity is proposed. This method is based upon a fractional time step scheme and the finite volume method on unstructured meshes. The governing equations are discretized using a collocated, cell-centered arrangement of velocity and pressure. The solution variables are stored at the cell-circumcenters. Theoretical results and numerical properties of the scheme are provided. Predictions of lid-driven cavity flow, flows past a cylinder and heat transport in a cylinder are performed to validate the method.     

Analytical solution for suction and injection flow of a viscoplastic Casson fluid past a stretching surface in the presence of viscous dissipation

In this study, steady two-dimensional flow of a viscoplastic Casson fluid past a stretching surface is considered under the effects of thermal radiation and viscous dissipation. Both suction and injection flows situations are considered. The partial differential governing equations are transformed into ordinary differential equations and solved analytical. Analytical solutions for velocity and temperature are obtained in terms of hypergeometric function and discussed graphically. Moreover, numerical results are also obtained by Runge–Kutta–Fehlberg fourth–fifth-order (RKF45) method and compared with the analytical results. The results showed that the injection and suction parameter can be used to control the direction and strength of flow. The effects of Casson parameter on the temperature and velocity are quite opposite. The effects of thermal radiation on the temperature are much more stronger in case of injection. The heat transfer coefficient shows higher value for Casson fluid while for Newtonian fluid is the lowest.

     

Large stencil viscous flux linearization for the simulation of 3D compressible turbulent flows with backward-Euler schemes

The purpose of this article is to study different approximate linearizations of the RANS equations viscous fluxes, for numerical simulations of compressible turbulent flows with backward-Euler schemes. The explicit convective flux under consideration is centred with artificial dissipation. The discrete viscous flux, calculated from cell-centred evaluation of the gradients, needs less computations and memory storage than other discretizations. But, in other respects, the balance of this numerical flux has a large stencil, which is not coherent with the 3-point per mesh direction stencil of classical implicit stages. Therefore 3-point and 5-point per mesh direction approximate linearizations are built from the thin layer flux formula. The stability condition of the corresponding backward-Euler schemes is given for a scalar linear equation (for the basic non-factored version of scheme and with LU-relaxation). Multigrid and monogrid computations of turbulent flow around two external configurations are performed with Wilcox’s k-ω turbulence model. The 5-point per mesh direction linearizations, coherent with the differential of the fluxes balance of thin layer approximation of explicit viscous fluxes, leads to the most efficient implicit stages.     

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov–Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics.     

Inertial thin film flow on planar surfaces featuring topography

A range of problems is investigated, involving the gravity-driven inertial flow of a thin viscous liquid film over an inclined planar surface containing topographical features, modelled via a depth-averaged form of the governing unsteady Navier-Stokes equations. The discrete analogue of the resulting coupled equation set, employing a staggered mesh arrangement for the dependent variables, is solved accurately using an efficient full approximation storage (FAS) algorithm and a full multigrid (FMG) technique; together with error-controlled automatic adaptive time-stepping and proper treatment of the associated nonlinear convective terms. An extensive set of results is presented for flow over both one- and two-dimensional topographical features, and errors quantified via detailed comparisons drawn with complementary experimental data and predictions from finite element analyses where they exist. In the case of one-dimensional (spanwise) topography, moderate Reynolds numbers and shallow/short topographical features, the results obtained are in close agreement with corresponding finite element solutions of the full free-surface problem. For the case of flow over two-dimensional (localised) topography, it is shown that the free-surface disturbance is influenced significantly by the presence of inertia leading, as in the case of spanwise topography, to an increase in the magnitude and severity of the resulting capillary ridge and trough formations: the effect of inclination angle and topography aspect ratio are similarly explored.