On the First Boundary Value Problem for Quasilinear Parabolic Equations with Two Independent Variables

This paper is concerned with the global solvability of the first initial boundary value problem for the quasilinear parabolic equations with two independent variables: a(t,x,u,u_xINF>)u_xxm u_t=f(t,x,u,u_xINF>). We investigate the case when the growth of [(|f(t,x,u,p)|)/(a(t,x,u,p))]{{|f(t,x,u,p)|}\over {a(t,x,u,p)}} with respect to p is faster than p² when |p|M X. Conditions which guarantee the global classical solvability of the problem are formulated.     

On the Riemann Problem for Non-Conservative Hyperbolic Systems

We consider the construction and the properties of the Riemann solver for the hyperbolic systemu_t + f(u)_x = 0, (0.1) u_t + f(u)_x = 0, (0.1) assuming only that Df is strictly hyperbolic. In the first part, we prove a general regularity theorem on the admissible curves T_i of the i-family, depending on the number of inflection points of f: namely, if there is only one inflection point, T_i is C^(1,1). If the i-th eigenvalue of Df is genuinely nonlinear, it is well known that T_i is C^(2,1). However, we give an example of an admissible curve T_i which is only Lipschitz continuous if f has two inflection points. In the second part, we show a general method for constructing the curves T_i, and we prove a stability result for the solution to the Riemann problem. In particular we prove the uniqueness of the admissible curves for (0.1). Finally we apply the construction to various approximations to (0.1): vanishing viscosity, relaxation schemes and the semidiscrete upwind scheme. In particular, when the system is in conservation form, we obtain the existence of smooth travelling profiles for all small admissible jumps of (0.1).     

On the damping function of shear relaxation modulus for entangled polymers

Published data of the damping function of the shear relaxation modulus, h(), are reviewed. This is the ratio of the relaxation modulus measured at a finite magnitude of shear, , to that at the limit of = 0. Majority of the data are in accord with the universal function of the Doi-Edwards tube model theory, in which the damping or the decrease of h() is attributed to the contraction along the tube of extended polymer chains. The weaker damping seems to be attributed to 1) comb-branching such as in LDPE; 2) lack of entanglement in too short chains; 3) bimodal molecular weight distribution. However, a star-branching does not cause a deviation from the tube model theory and a broadness of molecular weight distribution is not a major origin of a weaker damping. A star-branched polystyrene with 15 arms exhibits no strain dependence: h() = 1. For highly entangled systems with more than 50 entanglement points per molecule, the strain dependence is stronger than that of the Doi-Edwards theory. This could be due to a slip or an instability of deformation in the material.     

Effect of different initial conditions on a turbulent round free jet

Velocity measurements were made in two jet flows, the first exiting from a smooth contraction nozzle and the second from a long pipe with a fully developed pipe flow profile. The Reynolds number, based on nozzle diameter and exit bulk velocity, was the same (䏪,000) in each flow. The smooth contraction jet flow developed much more rapidly and approached self-preservation more rapidly than the pipe jet. These differences were associated with differences in the turbulence structure in both the near and far fields between the two jets. Throughout the shear layer for x<3d, the peak in the v spectrum occurred at a lower frequency in the pipe jet than in the contraction jet. For x́d, the peaks in the two jets appeared to be nearly at the same frequency. In the pipe jet, the near-field distributions of f(r) and g(r), the longitudinal and transverse velocity correlation functions, differed significantly from the contraction jet. The integral length scale L_u was greater in the pipe jet, whereas L_v was smaller. In the far field, the distributions of f(r) and g(r) were nearly similar in the two flows. The larger initial shear layer thickness of the pipe jet produced a dimensionally lower frequency instability, resulting in longer wavelength structures, which developed and paired at larger downstream distances. The regular vortex formation and pairing were disrupted in the shear layer of the pipe jet. The streamwise vortices, which enhance entrainment and turbulent mixing, were absent in the shear layer of the pipe jet. The formation of large-scale structures should occur much farther downstream in the pipe jet than in the contraction jet.     

Stabilizability of Two-Dimensional Navier—Stokes Equations with Help of a Boundary Feedback Control

For 2D Navier--Stokes equations defined in a bounded domain W \Omega we study stabilization of solution near a given steady-state flow [^(v)](x) \hat v(x) by means of feedback control defined on a part G \Gamma of boundary ?W \partial\Omega . New mathematical formalization of feedback notion is proposed. With its help for a prescribed number s > 0 \sigma > 0 and for an initial condition v₀(x) placed in a small neighbourhood of [^(v)](x) \hat v(x) a control u(t,x'), x￠ ? G x' \in \Gamma , is constructed such that solution v(t,x) of obtained boundary value problem for 2D Navier--Stokes equations satisfies the inequality: ||v(t,·)-[^(v)]||_H¹\leqslant ce^-st    for  t \geqslant 0 \|v(t,\cdot)-\hat v\|_{H^1}\leqslant ce^{-\sigma t}\quad {\rm for}\; t \geqslant 0 . To prove this result we firstly obtain analogous result on stabilization for 2D Oseen equations.     

Study of the flow around surface-mounted pyramids

The turbulent flow around square-based wall-mounted pyramids in thin and thick boundary layers was experimentally investigated as a function of the pyramid apex angle, ‘, and angle of attack, f, based on mean flow surface patterns, pressure and velocity measurements. For thin boundary layers, wake periodicity is observed. For slender pyramids (15°<‘<75°), the periodic formation and shedding of vortices is observed. The shedding frequency scales with the length scale L=w'(w'/h)^-1/4, where h is the pyramid height and w' is the frontal (projected) base width. For broad pyramids, wake periodicity exists but cannot be related to vortex shedding. Vortex shedding appears related to the existence of a double vortex structure along the side faces of the slender pyramids. The location of the separation point upstream and the attachment point downstream of the obstacle also scales with L. For thick boundary layers, no wake periodicity is observed and the mean flow structure in the wake differs from that seen for thin boundary layers. The location of the separation and attachment points scales only approximately with L.     

Effect of surfactant concentration on the formation and viscoelasticity of anionic wormlike micelle by the methods of rheology and freeze-fracture TEM

Wormlike micelles and network structures were obtained for the samples with the anionic surfactant, sodium dodecyl trioxyethylene sulfate (SDES), and multivalent counterion electrolyte AlCl₃ fixed in a proper concentration according to rheological measurements and the technique of freeze-fracture transmission electron microscopy (FF-TEM). The characteristics of a non-Newtonian fluid and the phenomenon of shear-induced structures transition were acquired from the steady-shear experiments. The TEM micrographs show the coexistence of wormlike micelles, entangled network structures, rod- and sphere- like micelles in the solutions. The Cole-Cole plot shows the rheological behavior of the non-linear viscoelasticity, the departure from the Maxwell model, and a relaxation time that is not simple. The plateau modulus G₀ and the relaxation time F were also investigated to express the rheological properties of the wormlike micellar solutions.     

Particle image velocimetry measurements of a backward-facing step flow

Particle image velocimetry (PIV) measurements were carried out on a backward-facing step flow at a Reynolds number of Re_h=U_Xh/9=4,660 (based on step height and freestream velocity). In-plane velocity, out-of-plane vorticity, Reynolds stress and turbulent kinetic energy production measurements in the x-y and x-z planes of the flow are presented. Proper orthogonal decomposition was performed on both the fluctuating velocity and vorticity fields of the x-y plane PIV data using the method of snapshots. Low-order representations of the instantaneous velocity fields were reconstructed using the velocity modes. These reconstructions provided insight into the contribution that the various length scales make to the spatial distribution of mean and turbulent flow quantities such as Reynolds stress and turbulent kinetic energy production. Large scales are found to contribute to the Reynolds stresses and turbulent kinetic energy production downstream of reattachment, while small scales contribute to the intense Reynolds stresses in the vicinity of reattachment.     

On multiple-path sonic anemometer measurement theory

In this paper, a model of the measuring process of sonic anemometers with more than one measuring path is presented. The main hypothesis of the work is that the time variation of the turbulent speed field during the sequence of pulses that produces a measure of the wind speed vector affects the measurement. Therefore, the previously considered frozen flow, or instantaneous averaging, condition is relaxed. This time variation, quantified by the mean Mach number of the flow and the time delay between consecutive pulses firings, in combination with both the full geometry of sensors (acoustic path location and orientation) and the incidence angles of the mean with speed vector, give rise to significant errors in the measurement of turbulence which are not considered by models based on the hypothesis of instantaneous line averaging. The additional corrections (relative to the ones proposed by instantaneous line-averaging models) are strongly dependent on the wave number component parallel to the mean wind speed, the time delay between consecutive pulses, the Mach number of the flow, the geometry of the sensor and the incidence angles of mean wind speed vector. Kaimal´s limit k_W1=1/l (where k_W1 is the wave number component parallel to mean wind speed and l is the path length) for the maximum wave numbers from which the sonic process affects the measurement of turbulence is here generalized as k_W1=C_l/l, where C_l is usually lesser than unity and depends on all the new parameters taken into account by the present model.     

Entangled polymer orientation and stretch under large step shear deformations in primitive chain network simulations

Orientation and stretch of entangled polymers under large step shear deformations were investigated through primitive chain network simulations. In the simulations, entangled polymer dynamics is described by 3D motion of entanglements, 1D sliding of monomers along the chain, and creation/destruction of entanglements described by hooking/unhooking with surrounding chains at chain ends. In addition to the conventionally proposed relaxation mechanisms that are reptation, contour length fluctuations, and constraint release (both thermal and convective), the simulations also account for force balance among entanglement strands converging to an entanglement node, and nodes also fluctuate in space. Nonlinear step strain data for monodisperse polystyrene melts (Ferri and Greco, Macromolecules, 37:5931, 2006) were quantitatively reproduced by using the same two molecular-weight-independent parameters already adopted by Masubuchi et al. (J Non-Newt Fluid Mech, 149:87, 2008) to fit linear viscoelastic data of several monodisperse polystyrene melts. Analysis of the orientation tensor and of the chain stretch ratio indicates that the segment orientation and stretch realized in the simulation are quantitatively described by a simple three-chain model (Marrucci et al., Macromol Symp, 158:57, 2000a).     

Rheology of SiO2/(acrylic polymer/epoxy) suspensions. II. Nonlinear stress relaxation

Large deformation, nonlinear stress relaxation modulus G(t, γ) was examined for the SiO₂ suspensions in a blend of acrylic polymer (AP) and epoxy (EP) with various SiO₂ volume fractions (?) at various temperatures (T). The AP/EP contained 70 vol.% of EP. At ??≤?30 vol.%, the SiO₂/(AP/EP) suspensions behaved as a viscoelastic liquid, and the time-strain separability, G(t, γ)?=?G(t)h(γ), was applicable at long time. The h(γ) of the suspensions was more strongly dependent on γ than that of the matrix (AP/EP). At ??=?35 vol.% and T?=?100°C, and ??≥?40 vol.%, the time-strain separability was not applicable. The suspensions exhibited a critical gel behavior at ??=?35 vol.% and T?=?100°C characterized with a power law relationship between G(t) and t; G(t)?∝?t ^???n . The relaxation exponent n was estimated to be about 0.45, which was in good agreement with the result of linear dynamic viscoelasticity reported previously. G(t, γ) also could be approximately expressed by the relation $G(t,\gamma) \propto t^{-n^{\prime}}$ at ??=?40 vol.%. The exponent n ^′ increased with increasing γ. This nonlinear stress relaxation behavior is attributable to strain-induced disruption of the network structure formed by the SiO₂ particles therein.     

Effect of ionic interaction on linear and nonlinear viscoelastic properties of ethylene based ionomer melts

The effect of ionic interaction on linear and nonlinear viscoelastic properties was investigated using poly(ethylene-co-methacrylic acid) (E/MAA) and its ionomers which were partially neutralized by zinc or sodium. Dynamic shear viscosity and step-shear stress relaxation studies were performed. Stress relaxation moduli G(t, y) of the E/MAA and its sodium or zinc ionomers were factorized into linear relaxation moduli G°(t) and damping functions h(y). The relaxation modulus at the smallest strain in each ionomer agreed with the linear relaxation modulus calculated from storage modulus G and loss modulus G. In the linear region, the ionic interaction shifted the relaxation time longer with keeping the same relaxation time distribution as E/MAA. In the nonlinear region, the ionic interaction had no influence on h(y) when the ion content was low. At higher ion content, however, the ion bonding enhanced the strain softening of h(y).     

L1 Continuous Dependence Property¶for Systems of Conservation Laws

This paper is concerned with the uniqueness and L¹ continuous dependence of entropy solutions for nonlinear hyperbolic systems of conservation laws. We study first a class of linear hyperbolic systems with discontinuous coefficients: Each propagating shock wave may be a Lax shock, or a slow or fast undercompressive shock, or else a rarefaction shock. We establish a result of L¹ continuous dependence upon initial data in the case where the system does not contain rarefaction shocks. In the general case our estimate takes into account the total strength of rarefaction shocks. In the proof, a new time-decreasing, weighted L¹ functional is obtained via a step-by-step algorithm. To treat nonlinear systems, we introduce the concept of admissible averaging matrices which are shown to exist for solutions with small amplitude of genuinely nonlinear systems. Interestingly, for many systems of continuum mechanics, they also exist for solutions with arbitrary large amplitude. The key point is that an admissible averaging matrix does not exhibit rarefaction shocks. As a consequence, the L¹ continuous dependence estimate for linear systems can be extended to nonlinear hyperbolic systems using a wave-front tracking technique.     

An Existence Theorem¶for the Navier-Stokes Flow¶in the Exterior of a Rotating Obstacle

We consider the three-dimensional Navier-Stokes initial value problem in the exterior of a rotating obstacle. It is proved that a unique solution exists locally in time if the initial velocity possesses the regularity L^1/2. This regularity assumption is the same as that in the famous paper of Fujita &; Kato. An essential step for the proof is the deduction of a certain smoothing property together with estimates near t˸ of the semigroup, which is not an analytic one, generated by the operator \Cal Lu = -P[\De u+(\om×x)·\na u-\om×u]\Cal Lu= -P\left[\De u+(\om\times x)\cdot\na u-\om\times u\right] in the space L², where y stands for the angular velocity of the rotating obstacle and P denotes the projection associated with the Helmholtz decomposition.     

Entropy Solutions for Nonlinear Degenerate Problems

We consider a class of elliptic-hyperbolic degenerate equations g(u)-Db(u) +\divgf(u) = fg(u)-\Delta b(u) +\divg\phi (u) =f with Dirichlet homogeneous boundary conditions and a class of elliptic-parabolic-hyperbolic degenerate equations g(u)_t-Db(u) +\divgf(u) = fg(u)_t-\Delta b(u) +\divg\phi (u) =f with homogeneous Dirichlet conditions and initial conditions. Existence of entropy solutions for both problems is proved for nondecreasing continuous functions g and b vanishing at zero and for a continuous vectorial function J satisfying rather general conditions. Comparison and uniqueness of entropy solutions are proved for g and b continuous and nondecreasing and for J continuous.     

Experimental analysis of turbulent flow structure in a fully developed rib-roughened rectangular channel with PIV

An experimental study was carried out to investigate the turbulent water-flow structure over one-side micro-repeated ribs in a narrow two-dimensional rectangular channel by particle image velocimetry (PIV). Two rib pitch-to-height ratios p/k of 10 and 20 were investigated while the rib height was held constant at 4 mm. The rib height-to-channel equivalent diameter ratio k/D_e was 0.1. The streamwise mean velocity and turbulent kinetic energy fields in the fully developed flow region of the channel were calculated at three different positions of x=0, 3.8, and -6.2 mm, which corresponded to center, downstream, and upstream of the rib, respectively, and for two Reynolds numbers Re of 7,000 and 20,000. A large-scale turbulent eddy was generated by the rib promoter and then propagated into the mainstream flow, which led to the deformation of the velocity profile. Downstream of the rib, rotating and counter-rotating eddies were also generated by the rib promoter. The enhancement of the turbulent kinetic energy was not changed when the Reynolds number increased from 7,000 to 20,000 between p/k=20 and 10. The reattachment length L_R was measured from velocity vector fields in the developing, fully developed, and exit regions of the flow over the range Re=1,400-50,000. The results showed that the ratio p/k and the Reynolds number had no significant effect on the reattachment length beyond a critical value of Re=15,000, where L_R was found to be approximately 4 times the rib height.     

Propagation of shear-layer structures in the near-wall region of a turbulent boundary layer

Two non-intrusive techniques, namely laser Doppler anemometry (LDA) and the electrochemical method, have been used for simultaneous measurements of the instantaneous streamwise velocity (U) and longitudinal wall shear stress (S), evaluated in a zero pressure gradient turbulent boundary layer. The space-time correlation between the fluctuating velocity and shear stress suggests that the coherent flow structures are propagated (i) under a slight angle of 5° in the near-wall region and (ii) at an average angle of 15.5° for y⁺>30. It is shown that the time shift obtained from the correlation between the LDA and the electrochemical signals is due to the dynamic behaviour of the electrodiffusion probe, but also to the leaning character of the coherent structures.     

Experimental study on the capillary thinning of entangled polymer solutions

The transient elongation behavior of entangled polymer and wormlike micelles (WLM) solutions has been investigated using capillary breakup extensional rheometry (CaBER). The transient force ratio X = 0.713 reveals the existence of an intermediate Newtonian thinning region for polystyrene and WLM solutions prior to the viscoelastic thinning. The exponential decay of X(t) in the first period of thinning defines an elongational relaxation time λ _x which is equal to elongational relaxation time λ _e obtained from exponential diameter decay D(t) indicating that the initial stress decay is controlled by the same molecular relaxation process as the strain hardening observed in the terminal regime of filament thinning. Deviations in true and apparent elongational viscosity are discussed in terms of X(t). A minimum Trouton ratio is observed which decreases exponentially with increasing polymer concentration leveling off at Tr_min = 3 for the solutions exhibiting intermediate Newtonian thinning and Tr_min ≈ 10 otherwise. The relaxation time ratio λ _e/ λ _s, where λ _s is the terminal shear relaxation time, decreases exponentially with increasing polymer concentration and the data for all investigated solutions collapse onto a master curve irrespective of polymer molecular weight or solvent viscosity when plotted versus the reduced concentration c[ η], with [ η] being the intrinsic viscosity. This confirms the strong effect of the nonlinear deformation in CaBER experiments on entangled polymer solutions as suggested earlier. On the other hand, λ _e ≈λ _s is found for all WLM solutions clearly indicating that these nonlinear deformations do not affect the capillary thinning process of these living polymer systems.     

Rheological properties of branched polystyrenes: nonlinear shear and extensional behavior

Nonlinear shear and uniaxial extensional measurements on a series of graft-polystyrenes with varying average numbers and molar masses of grafted side chains are presented. Step-strain measurements are performed to evaluate the damping functions of the melts in shear. The damping functions show a decreasing dependence on strain with an increase in mass fraction of grafted side chains. Extensional viscosities of the melts of graft-polystyrenes exhibit a growing strain hardening with increasing average number of grafted side chains as long as the side branches have a sufficient molar mass to be entangled. Graft-polystyrenes with side arms below the critical molar mass M _c for entanglements of linear polystyrene but above the entanglement molar mass M _e of linear polystyrene (M _e < M _w,br < M _c) still show a distinct strain hardening. With decreasing molar mass of the grafted side chains (M _w,br < M _e) the nonlinear-viscoelastic properties of the graft-polystyrene melts approach the behavior for a linear polystyrene with comparable polydispersity.Electronic Supplementary Material Supplementary material is available for this article at     

Compliance of actin filament networks measured by particle-tracking microrheology and diffusing wave spectroscopy

We monitor the time-dependent shear compliance of a solution of semi-flexible polymers, using diffusing wave spectroscopy (DWS) and video-enhanced single-particle-tracking (SPT) microrheology. These two techniques use the small thermally excited motion of probing microspheres to interrogate the local properties of polymer solutions. The solutions consist of networks of actin filaments which are long semi-flexible polymers. We establish a relationship between the mean square displacement (MSD) of microspheres imbedded in the solution and the time-dependent creep compliance of the solution, <Δr ²(t)>=(k _B T/πa)J(t). Here, J(t) is the creep compliance, <Δr ²(t)> is the mean-square displacement, and a is the radius of the microsphere chosen to be larger than the mesh size of the polymer network. DWS allows us to measure mean square displacements with microsecond temporal resolution and Ångström spatial resolution. At short times, the mean square displacement of a 0.96μm diameter sphere in a concentrated actin solution displays sub-diffusion. <Δr ²(t)>∝t ^{, with a characteristic exponent =0.78±0.05, which reflects the finite rigidity of actin. At long times, the MSD reaches a plateau, with a magnitude that decreases with concentration. The creep compliance is shown to be a weak function of polymer concentration and scales as J _p ∝c –1.2±0.3. This exponent is correctly described by a recent model describing the viscoelasticity of semi-flexible polymer solutions. The DWS and video-enhanced SPT measurements of the compliance plateau agree quantitatively with compliance measured independently using classical mechanical rheometry for a viscous oil and for a solution of flexible polymers. This paper extends the use of DWS and single-particle-tracking to probe the local mechanical properties of polymer networks, shows for the first time the proportionality between mean square displacement and local creep compliance, and therefore presents a new, direct way to extract the viscoelastic properties of polymer systems and complex fluids.}