Integrability of the constrained rigid body

The integrability theory for the differential equations, which describe the motion of an unconstrained rigid body around a fixed point is well known. When there are constraints the theory of integrability is incomplete. The main objective of this paper is to analyze the integrability of the equations of motion of a constrained rigid body around a fixed point in a force field with potential U(γ)=U(γ ₁,γ ₂,γ ₃). This motion subject to the constraint 〈ν,ω〉=0 with ν is a constant vector is known as the Suslov problem, and when ν=γ is the known Veselova problem, here ω=(ω ₁,ω ₂,ω ₃) is the angular velocity and 〈?,?〉 is the inner product of $\mathbb{R}^{3}$ . We provide the following new integrable cases. (i) The Suslov’s problem is integrable under the assumption that ν is an eigenvector of the inertial tensor I and the potential is such that $$U=-\frac{1}{2I_1I_2}\bigl(I_1\mu^2_1+I_2 \mu^2_2\bigr), $$ where I ₁,I ₂, and I ₃ are the principal moments of inertia of the body, μ ₁ and μ ₂ are solutions of the first-order partial differential equation $$\gamma_3 \biggl(\frac{\partial\mu_1}{\partial\gamma_2}- \frac{\partial\mu_2}{\partial \gamma_1} \biggr)- \gamma_2\frac{\partial \mu_1}{\partial\gamma_3}+\gamma_1\frac{\partial\mu_2}{\partial \gamma_3}=0. $$ (ii) The Veselova problem is integrable for the potential $$U=-\frac{\varPsi^2_1+\varPsi^2_2}{2(I_1\gamma^2_2+I_2\gamma^2_1)}, $$ where Ψ ₁ and Ψ ₂ are the solutions of the first-order partial differential equation where $p=\sqrt{I_{1}I_{2}I_{3} (\frac{\gamma^{2}_{1}}{I_{1}}+\frac{\gamma^{2}_{2}}{I_{2}}+ \frac{\gamma^{2}_{3}}{I_{3}} )}$ . Also it is integrable when the potential U is a solution of the second-order partial differential equation where $\tau_{2}=I_{1}\gamma^{2}_{1}+I_{2}\gamma^{2}_{2}+I_{3}\gamma^{2}_{3}$ and $\tau_{3}=\frac{\gamma^{2}_{1}}{I_{1}}+\frac{\gamma^{2}_{2}}{I_{2}}+ \frac{\gamma^{2}_{3}}{I_{3}}$ . Moreover, we show that these integrable cases contain as a particular case the previous known results.     

Global Compensated Compactness Theorem for General Differential Operators of First Order

Let A ₁(x, D) and A ₂(x, D) be differential operators of the first order acting on l-vector functions ${u= (u_1, \ldots, u_l)}$ in a bounded domain ${\Omega \subset \mathbb{R}^{n}}$ with the smooth boundary ${\partial\Omega}$ . We assume that the H ¹-norm ${\|u\|_{H^{1}(\Omega)}}$ is equivalent to ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_1u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ and ${\sum_{i=1}^2\|A_iu\|_{L^2(\Omega)} + \|B_2u\|_{H^{\frac{1}{2}}(\partial\Omega)}}$ , where B _i  = B _i (x, ν) is the trace operator onto ${\partial\Omega}$ associated with A _i (x, D) for i = 1, 2 which is determined by the Stokes integral formula (ν: unit outer normal to ${\partial\Omega}$ ). Furthermore, we impose on A ₁ and A ₂ a cancellation property such as ${A_1A_2^{\prime}=0}$ and ${A_2A_1^{\prime}=0}$ , where ${A^{\prime}_i}$ is the formal adjoint differential operator of A _i (i = 1, 2). Suppose that ${\{u_m\}_{m=1}^{\infty}}$ and ${\{v_m\}_{m=1}^{\infty}}$ converge to u and v weakly in ${L^2(\Omega)}$ , respectively. Assume also that ${\{A_{1}u_m\}_{m=1}^{\infty}}$ and ${\{A_{2}v_{m}\}_{m=1}^{\infty}}$ are bounded in ${L^{2}(\Omega)}$ . If either ${\{B_{1}u_m\}_{m=1}^{\infty}}$ or ${\{B_{2}v_m\}_{m=1}^{\infty}}$ is bounded in ${H^{\frac{1}{2}}(\partial\Omega)}$ , then it holds that ${\int_{\Omega}u_m\cdot v_m \,{\rm d}x \to \int_{\Omega}u\cdot v \,{\rm d}x}$ . We also discuss a corresponding result on compact Riemannian manifolds with boundary.     

Mixed Convection Boundary-Layer Flow Along a Vertical Cylinder Embedded in a Porous Medium Filled by a Nanofluid

The steady mixed convection boundary-layer flow on a vertical circular cylinder embedded in a porous medium filled by a nanofluid is studied for both cases of a heated and a cooled cylinder. The governing system of partial differential equations is reduced to ordinary differential equations by assuming that the surface temperature of the cylinder and the velocity of the external (inviscid) flow vary linearly with the axial distance x measured from the leading edge. Solutions of the resulting ordinary differential equations for the flow and heat transfer characteristics are evaluated numerically for various values of the governing parameters, namely the nanoparticle volume fraction ${\phi}$ , the mixed convection or buoyancy parameter ?? and the curvature parameter ??. Results are presented for the specific case of copper nanoparticles. A critical value ?? _c of ?? with ?? _c <?0 is found, with the values of | ?? _c| increasing as the curvature parameter ?? or nanoparticle volume fraction ${\phi}$ is increased. Dual solutions are seen for all values of ?? >??? _c for both aiding, ?? >?0 and opposing, ?? <?0, flows. Asymptotic solutions are also determined for both the free convection limit ${(\lambda \gg 1)}$ and for large curvature parameter ${(\gamma \gg 1)}$ .     

Vertical Motions of Heavy Inertial Particles Smaller than the Smallest Scale of the Turbulence in Strongly Stratified Turbulence

We study the statistics of the vertical motion of inertial particles in strongly stratified turbulence. We use Kinematic Simulation (KS) and Rapid Distortion Theory (RDT) to study the mean position and the root mean square (rms) of the position fluctuation in the vertical direction. We vary the strength of the stratification and the particle inertial characteristic time. The stratification is modelled using the Boussinesq equation and solved in the limit of RDT. The validity of the approximations used here requires that $ \sqrt{{L}/{g}} < {2\pi}/{\mathcal{N}} < \tau_{\eta} $ , where τ _η is the Kolmogorov time scale, g the gravitational acceleration, L the turbulence integral length scale and $\mathcal{N}$ the Brunt–Väisälä frequency. We introduce a drift Froude number $Fr_{d} = \tau_p g / \mathcal{N} L$ . When Fr _d ?<?1, the rms of the inertial particle displacement fluctuation is the same as for fluid elements, i.e. $\langle(\zeta_3 - \langle \zeta_3 \rangle)^2\rangle^{1/2} = 1.22\, u'/\mathcal{N} + \mbox{oscillations}$ . However, when Fr _d ?>?1, $\langle(\zeta_3 - \langle \zeta_3 \rangle)^2\rangle^{1/2} = 267 \, u' \tau_p$ . That is the level of the fluctuation is controlled by the particle inertia τ _p and not by the buoyancy frequency $\mathcal{N}$ . In other words it seems possible for inertial particles to retain the vertical capping while loosing the memory of the Brunt–Väisälä frequency.     

Statistical Analysis of Cross Scalar Dissipation Rate Transport in Turbulent Partially Premixed Flames: A Direct Numerical Simulation Study

Statistically planar turbulent partially premixed flames for different initial intensities of decaying turbulence have been simulated for global equivalence ratios <????> = 0.7 and <????> = 1.0 using three-dimensional simplified chemistry based Direct Numerical Simulations (DNS). The simulation parameters are chosen such that the combustion situation belongs to the thin reaction zones regime and a random bi-modal distribution of equivalence ratio ?? is introduced in the unburned gas ahead of the flame to account for mixture inhomogeneity. The DNS data has been used to analyse the statistical behaviour of the transport of the cross-scalar dissipation rate based on the fuel mass fraction Y _F and the mixture fraction ?? fluctuations $\,\tilde{\varepsilon}_{Y\xi}={\overline{\rho D\nabla Y_{F}^{\prime\prime}.\nabla \xi^{\prime\prime}} } \big/ {\bar {\rho }}$ (where $\bar{q}$ , $\tilde{q}={\overline{\rho q} } \big/ {\bar {\rho }}$ and $q^{\prime\prime} =q-\tilde {q}$ are Reynolds average, Favre mean and Favre fluctuation of a general quantity q) in the context of Reynolds Averaged Navier?CStokes simulations where ?? is the gas density and D is the gas diffusivity. The statistical behaviours of the unclosed terms in the $\tilde{\varepsilon }_{Y\xi } $ transport equation originating from turbulent transport T ₁, density variation T ₂, scalar?Cturbulence interaction T ₃, chemical reaction rate T ₄ and the molecular dissipation rate D ₂ have been analysed in detail. It has been observed that the contributions of T ₂, T ₃, T ₄ and D ₂ play important roles in the $\tilde{\varepsilon }_{Y\xi } $ transport for the globally stoichiometric cases, but in the globally fuel-lean cases the contributions of T ₂ and T ₄ become relatively weaker in comparison to the contributions of T ₃ and D ₂. The term T ₁ remains small in comparison to the leading order contributions of T ₃ and D ₂ for all cases, but the contribution of T ₁ plays a more important role in the low Damköhler combustion cases. The term T ₂ behaves as a sink term towards the unburned gas side but becomes a source term towards the burned gas side. The scalar?Cturbulence interaction term T ₃ has been found to be generally positive throughout the flame brush, but in globally stoichiometric cases the contribution of T ₃ becomes negative in regions of intense heat release. The combined contribution of (T ₄ ?C D ₂) remains mostly as a sink in all cases studied here. Models are proposed for the unclosed terms of the $\tilde{\varepsilon }_{Y\xi } $ transport equation in the context of Reynolds Averaged Navier?CStokes simulations, which are shown to satisfactorily predict the corresponding quantities extracted from the DNS data for all cases.     

A Concentration Phenomenon for Semilinear Elliptic Equations

For a domain ${\Omega \subset \mathbb{R}^{N}}$ we consider the equation $$-\Delta{u} + V(x)u = Q_n(x)|{u}|^{p-2}u$$ with zero Dirichlet boundary conditions and ${p\in(2, 2^*)}$ . Here ${V \geqq 0}$ and Q _n are bounded functions that are positive in a region contained in ${\Omega}$ and negative outside, and such that the sets {Q _n  > 0} shrink to a point ${x_0 \in \Omega}$ as ${n \to \infty}$ . We show that if u _n is a nontrivial solution corresponding to Q _n , then the sequence (u _n ) concentrates at x ₀ with respect to the H ¹ and certain L ^q -norms. We also show that if the sets {Q _n  > 0} shrink to two points and u _n are ground state solutions, then they concentrate at one of these points.     

Phragmén-Lindelöf theorems for nonlinear elliptic equations

We obtain theorems of Phragmén-Lindelöf type for the following classes of elliptic partial differential inequalities in an arbitrary unbounded domain \(\Omega \subseteq \mathbb{R}^n ,{\text{ }}n \geqq 2\) (A.1) $$\sum\limits_{i,j = 1}^n {\frac{\partial }{{\partial x_i }}\left( {a_{ij} 9(x)\frac{{\partial u}}{{\partial xj}}} \right)} + \sum\limits_{i = 1}^n {b_i (x,{\text{ }}u,{\text{ }}\nabla u)\frac{\partial }{{\partial x_i }}} \geqq f(x,{\text{ }}u)$$ where a _ij are elliptic in Ω and b _i ε L^∞(Ω) and where also a _ij are uniformly elliptic and Holder continuous at infinity and b _i = O(|x|⁺¹) as x → ∞; (A.2) $${\text{(A}}{\text{.2) }}\sum\limits_{i,j = 1}^n {a_{ij} (x,{\text{ }}u,{\text{ }}\nabla u)\frac{{\partial ^2 u}}{{\partial x_i \partial x_j }}} + \sum\limits_{i = 1}^n {b_i (x,{\text{ }}u,{\text{ }}\nabla u)\frac{\partial }{{\partial x_i }}} \geqq f(x,{\text{ }}u)$$ where a_ijare uniformly elliptic in Ω and b _iε L^∞(Ω); and finally (A.3) $${\text{div(}}\nabla u^p \nabla u {\text{)}} \geqq f{\text{(}}u{\text{), }}p > - 1,$$ where the operator on the left is the so-called P-Laplacian. The function f is always supposed positive and continuous. Moreover u is assumed throughout to be in the natural weak Sobolev space corresponding to the particular inequality under consideration, namely u ε. W _loc ^1,2 (Ω) ∩L _loc ^t8 (Ω) for (A.I), W _loc ^2,n(Ω) for (A.2), and W _loc ^1,p+2 (Ω) ∩ L _loc ^t8 (Ω) for (A.3). As a consequence of our results we obtain both non-existence and Liouville theorems, as well as existence theorems for (A.1).     

Modified dissipativity for a non-linear evolution equation arising in turbulence

We are concerned with the regularity properties for all times of the equation $$\frac{{\partial U}}{{\partial t}}\left( {t,x} \right) = - \frac{{\partial ^2 }}{{\partial x^2 }}\left[ {U\left( {t,{\text{0}}} \right) - U\left( {t,x} \right)} \right]^2 - v\left( { - \frac{{\partial ^2 }}{{\partial x^2 }}} \right)^\alpha U\left( {t,x} \right)$$ which arises, with α=1, in the theory of turbulence. Here U(t,·) is of positive type and the dissipativity α is a non-negative real number. It is shown that for arbitrary ν≧0 and ?>0, there exists a global solution in \(L^\infty [0,\infty ;H^{\tfrac{3}{2} - \varepsilon } (\mathbb{R})]\) . If ν>0 and \(\alpha > \alpha _{cr} = \tfrac{1}{2}\) , smoothness of initial data persists indefinitely. If 0≦α<α _cr, there exist positive constants ν₁(α) and ν₂(α), depending on the data, such that global regularity persists for ν>ν₁(α), whereas, for 0≦ν<ν₂(α), the second spatial derivative at the origin blows up after a finite time. It is conjectured that with a suitable choice of α _cr, similar results hold for the Navier-Stokes equation.     

Scaling of Conditional Lagrangian Time Correlation Functions of Velocity and Pressure Gradient Magnitudes in Isotropic Turbulence

We study Lagrangian statistics of the magnitudes of velocity and pressure gradients in isotropic turbulence by quantifying their correlation functions and their characteristic time scales. In a recent work (Yu and Meneveau, Phys Rev Lett 104:084502, 2010), it has been found that the Lagrangian time-correlations of the velocity and pressure gradient tensor and vector elements scale with the locally-defined Kolmogorov time scale, evaluated from the locally-averaged dissipation-rate (? _r ) and viscosity (ν) according to $\tau_{K,r}=\sqrt{\nu/\epsilon_r}$ . In this work, we study the Lagrangian time-correlations of the absolute values of velocity and pressure gradients. It has long been known that such correlations display longer memories into the inertial-range as well as possible intermittency effects. We explore the appropriate temporal scales with the aim to achieve collapse of the correlation functions. The data used in this study are sampled from the web-services accessible public turbulence database (http://turbulence.pha.jhu.edu). The database archives a 1024⁴ (space+time) pseudo-spectral direct numerical simulation of forced isotropic turbulence with Taylor-scale Reynolds number Re _λ ?=?433, and supports spatial differentiation and spatial/temporal interpolation inside the database. The analysis shows that the temporal auto-correlations of the absolute values extend deep into the inertial range where they are determined not by the local Kolmogorov time-scale but by the local eddy-turnover time scale defined as $\tau_{e,r}= r^{2/3}\epsilon_r^{-1/3}$ . However, considerable scatter remains and appears to be reduced only after a further (intermittency) correction factor of the form of (r/L) ^χ is introduced, where L is the turbulence integral scale. The exponent χ varies for different variables. The collapse of the correlation functions for absolute values is, however, less satisfactory than the collapse observed for the more rapidly decaying strain-rate tensor element correlation functions in the viscous range.     

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.     

Finite-dimensional asymptotic behavior of some semilinear damped hyperbolic problems

We prove that the solution semigroup $$S_t \left[ {u_0 ,v_0 } \right] = \left[ {u(t),u_t (t)} \right]$$ generated by the evolutionary problem $$\left\{ P \right\}\left\{ \begin{gathered} u_{tt} + g(u_t ) + Lu + f(u) = 0, t \geqslant 0 \hfill \\ u(0) = u_0 , u_t (0) = \upsilon _0 \hfill \\ \end{gathered} \right.$$ possesses a global attractorA in the energy spaceE _o=V×L ²(Ω). Moreover,A is contained in a finite-dimensional inertial setA attracting bounded subsets ofE ₁=D(L)×V exponentially with growing time.     

Flow behaviour of highly concentrated emulsions of supersaturated aqueous solution in oil

A set of highly concentrated water-in-oil emulsions with supersaturated dispersed phase were investigated in this work to verify and/or develop the models that have been presented both in the literature and in this work. The material used to form emulsions consisted of supersaturated oxidiser solution, hydrocarbon oil and PIBSA-based surfactants. The interfacial characteristics for different surfactant types were first examined. Then, the rheology of samples was studied, and different scaling methods and fitting of experimental data were studied. On the basis of flow curve measurements and observed $\tau _\emph{v} \sim \dot {\gamma }^{1/2}$ scaling, a modified version of Windhab model was suggested which showed excellent fitting of experimental results. The linear dependences of ?? _y0/?? versus 1/d ₃₂ for studied emulsions showed non-zero intercept which implies a non-linear dependence (resulting from interdroplet interaction) to fulfil the zero-intercept requirement. It was established that the zero intercept condition was fulfilled in the $\tau _{y0} \sim \sigma /d_{32}^2 $ scaling, although the experimental results for different surfactants were not superimposed.     

Biaxial gage factors for piezoresistive strain gages

Theoretical considerations of piezoresistive strain gages show that the change in electrical resistivity depends on the biaxial state of strain at the surface of the specimen to which the gage is bonded. In particular, whenV is the initial voltage across the gage and ( \( \in _{11} , \in _{22} , \in _{12} \) ) is the surface-strain state at the point of attachment, the gage-voltage change ΔV is given by \(\frac{{\Delta V}}{V} = G_{11} \in _{11} + G_{22} \in _{22} + G_{12} \in _{12} \) whereG ₁₁,G ₂₂ andG ₁₂ are the biaxial gage factors. Experiments were conducted on a nominally one-dimensional gage. Kulite type DLP-120-500, bonded to a standard ASTM flat tensile specimen of CR 1018 steel. For this gage, typical values were found to beG ₁₁?26,G ₂₂??1.4 andG ₁₂??1.1. SinceG ₂₂ andG ₁₂ are less than 6 percent ofG ₁₁, it is concluded that contributions from these two factors (called transverse and shear sensitivities) will be significant only when the gage is oriented such that \( \in _{11}<< \left( { \in _{22} , \in _{22} } \right)\) . However, in the interest of completeness and accuracy, all biaxial gage factors should be reported.     

C 1 Solutions for Semi-Implicit Systems of Differential Equations

The present note is a continuation of the author??s effort to study the existence of continuously differentiable solutions to the semi-implicit system of differential equations (1) $$f(x^{\prime}(t)) = g(t, x(t))$$ (2) $$\quad x(0) = x_0,$$ where

${\quad\Omega_g \subseteq \mathbb{R} \times\mathbb{R}^n}$ is an open set containing (0, x ₀) and ${g:\Omega_g \rightarrow\mathbb{R}^n}$ is a continuous function,

${\quad\Omega_f \subseteq \mathbb{R}^n}$ is an open set and ${f:\Omega_f\rightarrow\mathbb{R}^n}$ is a continuous function.

The transformation of (1)?C(2) into a solvable explicit system of differential equations is trivial if f is locally injective around an element ${\gamma\in \Omega_f\cap f^{-1}(g(0,x_0))}$ . In this paper, we study (1)?C(2) when such a translation is not possible because of the inherent multivalued nature of f ^?1.     

The Problems of the Turbulent Burning Velocity

The paper reviews the practical problems in measuring a turbulent burning velocity that gives the mass rate of burning. These largely centre on identifying an appropriate flame surface to associate with the turbulent burning velocity, u _t , and the density of the unburned mixture. Such a flame surface has been identified, in terms of the mean reaction progress variable, $\bar {c}$ , for explosive flame propagation in a fan-stirred bomb. Measurement of $\bar {c}$ makes possible an estimation of the flame surface density, ??, from the relationship ${\it \Sigma} =k\bar {c}\left( {1-\bar {c}} \right)$ . It is shown that in such explosions, mass rates of burning derived from the measured total flame surface area agreed well with those found from the measured turbulent burning velocity. Flamelet considerations identify appropriate dimensionless correlating parameters for u _t . As a result, correlations of turbulent burning velocity divided by the effective rms turbulent velocity, are plotted against the turbulent Karlovitz stretch factor, K, for different values of the Markstein number for flame strain rate, Ma_sr. These plots cover a wide range of variables, including pressure and fuels, and are indicative of different regimes of turbulent combustion. At the lower values of K, there is some evidence of increases in u _t and k due to high-frequency flame surface wrinkling arising from flame instabilities. These increase as Ma_sr becomes more negative. It is found from the developed value of the mean flame surface density throughout the flame brush that, to a first approximation, an increase in u _t for a given mixture is accompanied by a proportional increase in the volume of the brush. The analysis shows that the volume fraction of the turbulent flame brush that is reacting is quite small.     

On the stability or instability of the singular solution of the semilinear heat equation with exponential reaction term

We study questions of existence, uniqueness and asymptotic behaviour for the solutions of u(x, t) of the problem $$\begin{gathered} {\text{ }}u_t - \Delta u = \lambda e^u ,{\text{ }}\lambda {\text{ > 0, }}t > 0,{\text{ }}x{\text{ }}\varepsilon B, \hfill \\ (P){\text{ }}u(x,0) = u_0 (x),{\text{ }}x{\text{ }}\varepsilon B, \hfill \\ {\text{ }}u(x,t) = 0{\text{ }}on{\text{ }}\partial B \times (0,\infty ), \hfill \\ \end{gathered} $$ where B is the unit ball $\{ x\varepsilon R^N :|x|{\text{ }} \leqq {\text{ }}1\} {\text{ and }}N \geqq 3$ . Our interest is focused on the parameter λ ₀=2(N?2) for which (P) admits a singular stationary solution of the form $$S(x) = - 2log|x|$$ . We study the dynamical stability or instability of S, which depends on the dimension. In particular, there exists a minimal bounded stationary solution u which is stable if $3 \leqq N \leqq 9$ , while S is unstable. For $N \geqq 10$ there is no bounded minimal solution and S is an attractor from below but not from above. In fact, solutions larger than S cannot exist in any time interval (there is instantaneous blow-up), and this happens for all dimensions.     

Maximal and Minimal Spreading Speeds for Reaction Diffusion Equations in Nonperiodic Slowly Varying Media

This paper investigates the asymptotic behavior of the solutions of the Fisher-KPP equation in a heterogeneous medium, $$\partial_t u = \partial_{xx} u + f(x,u),$$ associated with a compactly supported initial datum. A typical nonlinearity we consider is ${f(x,u) = \mu_0 (\phi (x)) u(1-u)}$ , where??? ₀ is a 1-periodic function and ${\phi}$ is a ${\mathcal{C}^1}$ increasing function that satisfies ${\lim_{x \to+\infty}\phi (x) = +\infty}$ and ${\lim_{x \to +\infty}\phi' (x) =0}$ . Although quite specific, the choice of such a reaction term is motivated by its highly heterogeneous nature. We exhibit two different behaviors for u for large times, depending on the speed of the convergence of ${\phi}$ at infinity. If ${\phi}$ grows sufficiently slowly, then we prove that the spreading speed of u oscillates between two distinct values. If ${\phi}$ grows rapidly, then we compute explicitly a unique and well determined speed of propagation w _??, arising from the limiting problem of an infinite period. We give a heuristic interpretation for these two behaviors.     

A well posed problem in singular Fickian diffusion

We prove that the problem of solving $$u_t = (u^{m - 1} u_x )_x {\text{ for }} - 1< m \leqq 0$$ with initial conditionu(x, 0)=φ(x) and flux conditions at infinity \(\mathop {\lim }\limits_{x \to \infty } u^{m - 1} u_x = - f(t),\mathop {\lim }\limits_{x \to - \infty } u^{m - 1} u_x = g(t)\) , admits a unique solution \(u \in C^\infty \{ - \infty< x< \infty ,0< t< T\} \) for every φεL¹(R), φ≧0, φ≡0 and every pair of nonnegative flux functionsf, g ε L _loc ^∞ [0, ∞) The maximal existence time is given by $$T = \sup \left\{ {t:\smallint \phi (x)dx > \int\limits_0^t {[f} (s) + g(s)]ds} \right\}$$ This mixed problem is ill posed for anym outside the above specified range.     

Improved theoretical model for wavy filmwise condensation

This paper presents a numerical solution for wavy laminar film-wise condensation on vertical walls. Integral method is achieved based on the recently developed simple wave equations. Solutions are obtained for ranges of dimensionless groups as follows: $$1.5 \leqslant \left( {Pr = \frac{{^{\mu C} p}}{k}} \right) \leqslant 6.0$$ $$10 \leqslant \left( {G = \frac{{^h fg}}{{^{C_p \Delta T} }}} \right) \leqslant 400$$ $$100 \leqslant \left( {S = \left( {\frac{{\sigma ^2 \rho }}{{g_\rho \mu ^4 }}} \right)^{{1 \mathord{\left/ {\vphantom {1 5}} \right. \kern-\nulldelimiterspace} 5}} } \right) \leqslant 400$$ $$1000 \leqslant \left( {L = \frac{{{\rm H}_t }}{{^\delta cr}}} \right) \leqslant 10000$$ . Such ranges cover the expected situations in industrial applications. It is found that the Reynolds number (Re=h_LΔTH_t/h_fg) is a linear function of L on the log-log plane. It is also relatively insensitive to small variations of Pr at high values of this number. At situations where G less than 200 the Re appears to be dependent on S. Agreement with experimental observation is improved over that obtained from previous analytical theories.