On the Temperature Slip Boundary Condition in a Mixed Convection Boundary-Layer Flow in a Porous Medium

A problem derived previously (Rohni et?al., Transp Porous Media 92:1?C14, 2012) for unsteady mixed convection flow in a porous medium involving a ??temperature slip?? boundary condition and fluid transfer through the boundary is considered. It is shown that the solution to this problem can be directly related to the solution of the corresponding problem for a prescribed surface temperature, involving a mixed convection parameter ??, an unsteadiness parameter A and transpiration parameter s. This latter problem is discussed in detail, particular attention being given to the steady analogue, A?=?0, allowing for fluid transfer through the surface, and to the unsteady problem, A?>?0, for an impermeable surface, s?=?0. Asymptotic results are obtained for large fluid transfer rates, ${s \gg 1}$ and ${s <0 , |s| \gg 1}$ and for large A. Particular attention is given to deriving asymptotic results for the critical points which determine the range of existence of solutions.     

Numerical Modelling of Steady-State Flow in 2D Cracked Anisotropic Porous Media by Singular Integral Equations Method

The equations governing plane steady-state flow in heterogeneous porous media containing curved-line intersecting cracks (Pouya and Ghabezloo in Transp Porous Media 84:511?C532, 2010) and the potential solution obtained for these equations are considered here. The theoretical results are first completed for the mass balance at crack intersections points. Then, a numerical procedure based on a singular integral equations method is described concretely to derive this solution for cracked materials. Closed-form expressions of elementary integrals for special choice of collocation points lead to a very quick and easy numerical method. It is shown that this method can be applied efficiently to the study of the steady-state flow in cracked materials with anisotropic matrix permeability and a dense distribution of curved-line intersecting cracks. Some applications of this method to the permeability of cracked materials are given.     

A Comment on the Flow and Heat Transfer Past a Permeable Stretching/Shrinking Surface in a Porous Medium: Brinkman Model

An analytical solution is presented for the boundary-layer flow and heat transfer over a permeable stretching/shrinking surface embedded in a porous medium using the Brinkman model. The problem is seen to be characterized by the Prandtl number $Pr$ , a mass flux parameter $s$ , with $s>0$ for suction, $s=0$ for an impermeable surface, and $s<0$ for blowing, a viscosity ratio parameter $M$ , the porous medium parameter $\Lambda $ and a wall velocity parameter $\lambda $ . The analytical solution identifies critical values which agree with those previously determined numerically (Bachok et al. Proceedings of the fifth International Conference on Applications of Porous Media, 2013) and shows that these critical values, and the consequent dual solutions, can arise only when there is suction through the wall, $s>0$ .     

Sensitivity Analysis for Dynamic Spontaneous Imbibition with Variable Inlet Saturation and Interfacial Coupling Effects

Dynamic spontaneous imbibition (DSI) plays an important role in oil reservoir characterization. Conventional equations used to characterize DSI consider neither interfacial coupling effects (ICE) nor variable inlet saturation (S*) for DSI. Yazzan et?al. (Transp Porous Media 87(1):309?C333, 2011a; 86(3):705?C717, 2011b) developed a set of equations, and a numerical solution scheme, to take into account ICE and variable S* for DSI. Based on these, a graphical user interface (GUI) simulator was built. A sensitivity analysis has been conducted to study the effect of the fluid and rock properties on DSI. The results reveal that including a variable S* has no significant impact; however, neglecting ICE results in an overestimation of the imbibition flow rate. Moreover, it is shown that the capillary and relative permeability curves determine the type of frontal advance, and that the imbibition recovery is proportional to the square root of time.     

On the Rigorous Derivation of the 3D Cubic Nonlinear Schrödinger Equation with a Quadratic Trap

We consider the dynamics of the three-dimensional N-body Schrödinger equation in the presence of a quadratic trap. We assume the pair interaction potential is N ^3β-1 V(N ^β x). We justify the mean-field approximation and offer a rigorous derivation of the three-dimensional cubic nonlinear Schrödinger equation (NLS) with a quadratic trap. We establish the space-time bound conjectured by Klainerman and Machedon (Commun Math Phys 279:169–185, 2008) for ${\beta \in (0, 2/7]}$ by adapting and simplifying an argument in Chen and Pavlovi? (Annales Henri Poincaré, 2013) which solves the problem for ${\beta \in (0, 1/4)}$ in the absence of a trap.     

Existence and Stability of Solutions for Steady Flows of Fibre Suspension Flows

We establish existence, uniqueness, convergence and stability of solutions to the equations of steady flows of fibre suspension flows. The existence of a unique steady solution is proven by using an iterative scheme. One of the restrictions imposed on the data confirms a well known fact proven in Galdi and Reddy (J Non-Newtonian Fluid Mech 83:205–230, 1999), Munganga and Reddy (Math Models Methods Appl Sci 12:1177–1203, 2002) and Munganga et al. (J Non-Newtonian fluid Mech 92:135–150, 2000) that the particle number N _p must be less than 35/2. Exact solutions are calculated for Couette and Poiseuille flows. Solutions of Poiseuille flows are shown to be more accurate than those of Couette flow when a time perturbation is considered.     

Measurement and Correlation of Pressure Drop Characteristics for Air Flow Through Sintered Metal Porous Media

Sintered metal porous media are currently used to replace conventional orifices as restrictors in air-bearing systems. The flow properties in porous media are generally approximated by Darcy and Forchheimer regimes in different flow regions. In this study, an ISO expanded expression is proven defective when it is used to represent flow properties through porous media under slight pressure drops ( ${<}10$  kPa). A modified Forchheimer equation is therefore developed to correlate the pressure drop with flow rate, including compressibility and inertial effects. Experimental and theoretical investigations on pressure drop characteristics are conducted with a series of metal-sintered porous media. Permeability is first determined in a strict Darcy region with $Re<0.1$ , followed by the inertial coefficient with $Re>0.1$ , rather than determining these two simultaneously. The theoretical mass flow rate in terms of the modified Forchheimer equation provides close approximations to the experimental data.     

On the Formulation of Continuum Thermodynamic Models for Solids as General Equations for Non-equilibrium Reversible-Irreversible Coupling

The purpose of this work is the comparison of some aspects of the formulation of material models in the context of continuum thermodynamics (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997) with their formulation in the form of a General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC: e.g., Grmela and Öttinger in Phys. Rev. E 56: 6620–6632, 1997; Öttinger and Grmela in Phys. Rev. E 56: 6633–6655, 1997; Öttinger in Beyond equilibrium thermodynamics, Wiley, New York, 2005; Grmela in J. Non-Newton. Fluid Mech. 165: 980–998, 2010). A GENERIC represents a generalization of the Ginzburg-Landau model for the approach of non-equilibrium systems to thermodynamic equilibrium. Originally developed to formulate non-equilibrium thermodynamic models for complex fluids, it has recently been applied to anisotropic inelastic solids in a Eulerian setting (Hütter and Tervoort in J. Non-Newton. Fluid Mech. 152: 45–52, 2008; 53–65, 2008; Adv. Appl. Mech. 42: 254–317, 2009) as well as to damage mechanics (Hütter and Tervoort in Acta Mech. 201: 297–312, 2008). In the current work, attention is focused for simplicity on the case of thermoelastic solids with heat conduction and viscosity in a Lagrangian setting (e.g., ?ilhavý in The mechanics and thermodynamics of continuous media, Springer, Berlin, 1997, Chaps. 9–12). In the process, the relation of the two formulations to each other is investigated in detail. A particular point in this regard is the concept of dissipation and its model representation in both contexts.     

Analytical Validation of a Continuum Model for Epitaxial Growth with Elasticity on Vicinal Surfaces

Within the context of heteroepitaxial growth of a film onto a substrate, terraces and steps self-organize according to misfit elasticity forces. Discrete models of this behavior were developed by Duport et al. (J Phys I 5:1317–1350, 1995) and Tersoff et al. (Phys Rev Lett 75:2730–2733, 1995). A continuum limit of these was in turn derived by Xiang (SIAM J Appl Math 63:241–258, 2002) (see also the work of Xiang and Weinan Phys Rev B 69:035409-1–035409-16, 2004; Xu and Xiang SIAM J Appl Math 69:1393–1414, 2009). In this paper we formulate a notion of weak solution to Xiang’s continuum model in terms of a variational inequality that is satisfied by strong solutions. Then we prove the existence of a weak solution.     

Thin Shear Layer Structures in High Reynolds Number Turbulence

Three-dimensional tomographic time dependent PIV measurements of high Reynolds number (Re) laboratory turbulence are presented which show the existence of long-lived, highly sheared thin layer eddy structures with thickness of the order of the Taylor microscale and internal fluctuations. Highly sheared layer structures are also observed in direct numerical simulations of homogeneous turbulence at higher values of Re (Ishihara et al., Annu Rev Fluid Mech 41:165–180, 2009). But in the latter simulation, where the fluctuations are more intense, the layer thickness is greater. A rapid distortion model describes the structure and spectra for the velocity fluctuations outside and within ‘significant’ layers; their spectra are similar to the Kolmogorov (C R Acad Sci URSS 30:299–303, 1941) and Obukhov (Dokl Akad Nauk SSSR 32:22–24, 1941) statistical model (KO) for the whole flow. As larger-scale eddy motions are blocked by the shear layers, they distort smaller-scale eddies leading to local zones of down-scale and up-scale transfer of energy. Thence the energy spectrum for high wave number k is $E_X (k)\sim Bk^{-2p}$ . The exponent p depends on the forms of the large eddies. The non-linear interactions between the distorted inhomogeneous eddies produce a steady local structure, which implies that 2p?=?5/3 and a flux of energy into the thin-layers balancing the intense dissipation, which is much greater than the mean $\left<\epsilon\right>$ . Thence $B\sim\left<\epsilon\right>^{2/3}$ as in KO. Within the thin layers the inward flux energises extended vortices whose thickness and spacing are comparable with the viscous microscale. Although peak values of vorticity and velocity of these vortices greatly exceed those based on the KO scaling, the form of the viscous range spectrum is consistent with their model.     

On the Energy Dissipation Rate of Solutions to the Compressible Isentropic Euler System

In this paper we extend and complement the results in Chiodaroli et al. (Global ill-posedness of the isentropic system of gas dynamics, 2014) on the well-posedness issue for weak solutions of the compressible isentropic Euler system in 2 space dimensions with pressure law p(ρ) = ρ ^γ , γ ≥ 1. First we show that every Riemann problem whose one-dimensional self-similar solution consists of two shocks admits also infinitely many two-dimensional admissible bounded weak solutions (not containing vacuum) generated by the method of De Lellis and Székelyhidi (Ann Math 170:1417–1436, 2009), (Arch Ration Mech Anal 195:225–260, 2010). Moreover we prove that for some of these Riemann problems and for 1 ≤ γ < 3 such solutions have a greater energy dissipation rate than the self-similar solution emanating from the same Riemann data. We therefore show that the maximal dissipation criterion proposed by Dafermos in (J Diff Equ 14:202–212, 1973) does not favour the classical self-similar solutions.     

Analytic derivation of the Cox�CMerz rule using the MLD ��toy�� model for polydisperse linear polymers

A general constitutive formalism, the ??na?ve?? polydisperse MLD model, has been developed by Mead et al. (Macromolecules 31:7895?C7914, 1998) and Mead (Rheol Acta 46:369?C395, 2007) at both the tube coordinate level and the mathematically simplified ??toy?? level independent of the tube coordinate. The model includes constraint release generated by convection-driven chain retraction (which is equivalent to ??convective constraint release?? (CCR)), reptation, and tube contour length fluctuations. The properties of the mathematically simplified na?ve polydisperse ??toy?? MLD model are explored in linear and nonlinear steady shear flows where we analytically derive the Cox?CMerz rule relating the steady shear viscosity to the modulus of the linear viscoelastic dynamic viscosity. The Cox?CMerz rule relating the linear viscoelastic material properties and the nonlinear material properties is shown to be a direct consequence of convective constraint release. The specific feature of CCR that leads to this result is that the relaxation rate due to convective constraint release is proportional to the shear rate, $\dot{{\gamma }}$ , independent of molecular weight. The viability of this well-known empirical relationship is a direct consequence of a coincidence in the mathematical structure of the linear viscoelastic material properties and convective constraint release. There is no physical analogy or relationship between the molecular relaxation mechanisms operative in linear (diffusive relaxation) and nonlinear (convective relaxation) flow regimes. The polydisperse MLD model predictions of the individual molecular weight component contributions to the flow curve, and interpretations thereof, are effectively identical to those first postulated by Bersted (J Appl Polym Sci 19:2167?C2177, 1975, J Appl Polym Sci 20:2705?C2714, 1976). Following the theoretical developments, a limited experimental study is executed with a commercial polydisperse polystyrene melt. Nearly quantitative agreement between the polydisperse MLD theory and experimental measurements of steady-shear viscosity and dynamic moduli is achieved over a wide range of shear rates.     

Weak Solutions of the Navier–Stokes Equations for Compressible Flows in a Half-Space with No-Slip Boundary Conditions

We consider the Navier–Stokes equations for the motion of compressible, viscous flows in a half-space ${\mathbb{R}^n_+,}$ n =  2,  3, with the no-slip boundary conditions. We prove the existence of a global weak solution when the initial data are close to a static equilibrium. The density of the weak solution is uniformly bounded and does not contain a vacuum, the velocity is Hölder continuous in (x, t) and the material acceleration is weakly differentiable. The weak solutions of this type were introduced by D. Hoff in Arch Ration Mech Anal 114(1):15–46, (1991), Commun Pure and Appl Math 55(11):1365–1407, (2002) for the initial-boundary value problem in ${\Omega = \mathbb{R}^n}$ and for the problem in ${\Omega = \mathbb{R}^n_+}$ with the Navier boundary conditions.     

A Quantitative Modulus of Continuity for the Two-Phase Stefan Problem

We derive the quantitative modulus of continuity $$\omega(r)=\left[ p+\ln \left( \frac{r_0}{r}\right)\right]^{-\alpha (n, p)},$$ which we conjecture to be optimal for solutions of the p-degenerate two-phase Stefan problem. Even in the classical case p = 2, this represents a twofold improvement with respect to the early 1980’s state-of-the-art results by Caffarelli– Evans (Arch Rational Mech Anal 81(3):199–220, 1983) and DiBenedetto (Ann Mat Pura Appl 103(4):131–176, 1982), in the sense that we discard one logarithm iteration and obtain an explicit value for the exponent α(n, p).     

Hölder Continuity and Injectivity of Optimal Maps

Consider transportation of one distribution of mass onto another, chosen to optimize the total expected cost, where cost per unit mass transported from x to y is given by a smooth function c(x, y). If the source density f ⁺(x) is bounded away from zero and infinity in an open region ${U' \subset \mathbf{R}^n}$ , and the target density f ^?(y) is bounded away from zero and infinity on its support ${\overline{V} \subset \mathbf{R}^n}$ , which is strongly c-convex with respect to U′, and the transportation cost c satisfies the ${\mathbf{A3}_{\rm w}}$ condition of Trudinger and Wang (Ann Sc Norm Super Pisa Cl Sci 5, 8(1):143–174, 2009), we deduce the local Hölder continuity and injectivity of the optimal map inside U′ (so that the associated potential u belongs to ${C^{1,\alpha}_{loc}(U')}$ ). Here the exponent α > 0 depends only on the dimension and the bounds on the densities, but not on c. Our result provides a crucial step in the low/interior regularity setting: in a sequel (Figalli et al., J Eur Math Soc (JEMS), 1131–1166, 2013), we use it to establish regularity of optimal maps with respect to the Riemannian distance squared on arbitrary products of spheres. Three key tools are introduced in the present paper. Namely, we first find a transformation that under ${\mathbf{A3}_{\rm w}}$ makes c-convex functions level-set convex (as was also obtained independently from us by Liu (Calc Var Partial Diff Eq 34:435–451, 2009)). We then derive new Alexandrov type estimates for the level-set convex c-convex functions, and a topological lemma showing that optimal maps do not mix the interior with the boundary. This topological lemma, which does not require ${\mathbf{A3}_{\rm w}}$ , is needed by Figalli and Loeper (Calc Var Partial Diff Eq 35:537–550, 2009) to conclude the continuity of optimal maps in two dimensions. In higher dimensions, if the densities f ^± are Hölder continuous, our result permits continuous differentiability of the map inside U′ (in fact, ${C^{2,\alpha}_{loc}}$ regularity of the associated potential) to be deduced from the work of Liu et al. (Comm Partial Diff Eq 35(1):165–184, 2010).     

Some Indentities for the Gradient of the Principal Invariants, Traces and Determinants via Grassmann Calculus

Let A be a second order tensor in a finite dimensional space. In this work we determine the gradient of the principal invariants of A and obtain some trace and determinant identities using only some standard rigorous statements concerning Grassmann calculus. We recover some of the results of Dui et al. (J. Elast. 75:193–196, 2004) and of Truesdell and Noll (The Non-linear Field Theories of Mechanics, Springer, Berlin, 2002) and solve an old problem proposed in SIAM Review concerning a determinant identity from a new perspective in a concise and simple manner.