Galerkin discontinuous approximation of the transport equation and viscoelastic fluid flow on quadrilaterals

Numerical simulation of industrial processes involving viscoelastic liquids is often based on finite element methods on quadrilateral meshes. However, numerical analysis of these methods has so far been limited to triangular meshes. In this work, we consider quadrilateral meshes. We first study the approximation of the transport equation by a Galerkin discontinuous method and prove an 𝒪(h^k+1/2) error estimates for the Q_k finite element. Then we study a differential model for viscoelastic flow with unknowns u the velocity, p the pressure, and σ the viscoelastic part of the extra-stress tensor. The approximations are ((Q₁)² transforms of) Q_k+1 continuous for u, Q_k discontinuous for σ, and P_k discontinuous for p, with k ≥ 1. Upwinding for σ is obtained by the Galerkin discontinuous method. We show that an error estimate of order 𝒪(h^k+1/2) is valid in the energy norm for the three unknowns. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 97–114, 1998     

Global strong solutions for a class of compressible non‐Newtonian fluids with vacuum

We study a class of compressible non‐Newtonian fluids in one space dimension. We prove, by using iterative method, the global time existence and uniqueness of strong solutions provided that the initial data satisfy a compatibility condition and the initial density is small in its H¹‐norm. The main difficulty is due to the strong nonlinearity of the system and the initial vacuum. Copyright © 2010 John Wiley & Sons, Ltd.     

Waves on the surface of a falling power-law fluid film

Waves that occur at the surface of a falling film of thin power-law fluid on a vertical plane are investigated. Using the method of integral relations an evolution equation is derived for two types of waves equation which are possible under long wave approximation. This equation reveals the presence of both kinematic and dynamic wave processes which may either act together or singularly dominate the wave field depending on the order of different parameters. It is shown that, at a small flow rate, kinematic waves dominate the flow field and the energy is acquired from the mean flow during the interaction of the waves, while, for high flow rate, inertial waves dominate and the energy comes from the kinematic waves. It is also found that this exchange of energy between kinematic and inertial waves strongly depends on the power-law index n. Linear stability analysis predicts the contribution of different terms in the wave mechanism. Further, it is found that the surface tension plays a double role: for a kinematic wave process, it exerts dissipative effects so that a finite amplitude case may be established, but for a dynamic wave process it yields dispersion. Further, it is shown that the non-Newtonian character n plays a vital role in controlling the role of the term that contains surface tension in the above processes.     

广义二阶流体管内轴向流动

在流体的本构关系中引入分数阶导数运算，对于介于粘性与弹性之间的流体的描述更具有合理性。本文将这种关系引入二阶流体，研究其管内轴向流动。我们先求出了１／２阶导数的解析解，用以验证Ｌａｐｌａｃｅ数值反演的ＣＲＵＭＰ方法的有效性。然后用ＣＲＵＭＰ法分析二阶流体管内轴向流动的特征。分析表明粘弹性特征越明显的流体，其速度与应力对分数导数的阶数越具有敏感性。     

A continuum model of gas flows with localized density variations

We discuss the kinetic representation of gases and the derivation of macroscopic equations governing the thermomechanical behavior of a dilute gas viewed at the macroscopic level as a continuous medium. We introduce an approach to kinetic theory where spatial distributions of the molecules are incorporated through a mean-free-volume argument. The new kinetic equation derived contains an extra term involving the evolution of this volume, which we attribute to changes in the thermodynamic properties of the medium. Our kinetic equation leads to a macroscopic set of continuum equations in which the gradients of thermodynamic properties, in particular density gradients, impact on diffusive fluxes. New transport terms bearing both convective and diffusive natures arise and are interpreted as purely macroscopic expansion or compression. Our new model is useful for describing gas flows that display non-local-thermodynamic-equilibrium (rarefied gas flows), flows with relatively large variations of macroscopic properties, and/or highly compressible fluid flows.     

Evaporation of an emulsion trapped in a yield stress fluid

The present work deals with emulsions of volatile alkanes in an aqueous clay suspension, Laponite, which forms a yield stress fluid. For a large enough yield stress (i.e. Laponite concentration), the oil droplets are prevented from creaming and the emulsions are thus mechanically stabilized. We have studied the evaporation kinetics of the oil phase of those emulsions in contact with the atmosphere. We show that the evaporation process is characterized by the formation of a sharp front separating the emulsion from a droplet-free Laponite phase, and that the displacement of the front vs. time follows a diffusion law. Experimental data are confronted to a diffusion-controlled model, in the case where the limiting step is the diffusion of the dissolved oil through the aqueous phase. The nature of the alkane, as well as its volume fraction in the emulsion, has been varied. Quantitative agreement with the model is achieved without any adjustable parameter and we describe the mechanism leading to the formation of a front.     

A fluctuating lattice Boltzmann scheme for the one-dimensional KPZ equation

Based on the well-known mapping between the Burgers equation with noise and the Kardar–Parisi–Zhang (KPZ) equation for fluctuating interfaces, we develop a fluctuating lattice Boltzmann (LB) scheme for growth phenomena, as described by the KPZ formalism. A very simple LB-KPZ scheme is demonstrated in 1+1 spacetime dimensions, and is shown to reproduce the scaling exponents characterizing the growth of one-dimensional fluctuating interfaces.     

A Mixture Model of Truncated Zeta Distributions with Applications to Scientific Collaboration Networks

The degree distribution has attracted considerable attention from network scientists in the last few decades to have knowledge of the topological structure of networks. It is widely acknowledged that many real networks have power-law degree distributions. However, the deviation from such a behavior often appears when the range of degrees is small. Even worse, the conventional employment of the continuous power-law distribution usually causes an inaccurate inference as the degree should be discrete-valued. To remedy these obstacles, we propose a finite mixture model of truncated zeta distributions for a broad range of degrees that disobeys a power-law behavior in the range of small degrees while maintaining the scale-free behavior. The maximum likelihood algorithm alongside the model selection method is presented to estimate model parameters and the number of mixture components. The validity of the suggested algorithm is evidenced by Monte Carlo simulations. We apply our method to five disciplines of scientific collaboration networks with remarkable interpretations. The proposed model outperforms the other alternatives in terms of the goodness-of-fit.     

Nonlinear Exact Solutions of the 2-Dimensional Rotational Euler Equations for the Incompressible Fluid

In this paper, the Clarkson–Kruskal direct approach is employed to investigate the exact solutions of the2-dimensional rotational Euler equations for the incompressible fluid. The application of the method leads to a system of completely solvable ordinary differential equations. Several special cases are discussed and novel nonlinear exact solutions with respect to variables x and y are obtained. It is of interest to notice that the pressure p is obtained by the second kind of curvilinear integral and the coefficients of the nonlinear solutions are solitary wave type functions like tanh(kt/2)and sech(kt/2) due to the rotational parameter k = 0. Such phenomenon never appear in the classical Euler equations wherein the Coriolis force arising from the gravity and Earth's rotation is ignored. Finally, illustrative numerical figures are attached to show the behaviors that the exact solutions may exhibit.     

Cavitation in superfluid helium-4 at low temperature

We have studied the nucleation of bubbles in pure superfluid helium-4 at temperatures down to 65 mK. We have found that the nucleation is a stochastic process, and that at temperatures below 600 mK the nucleation rate is independent of temperature. These results are consistent with the assumption that the nucleation takes place via quantum tunneling. Received: 15 November 1997 / Received in final form: 19 December 1997 / Accepted: 22 January 1998