Irregular Behavior of Solutions for Fisher’s Equation

This paper is concerned with the irregular behavior of solutions for Fisher’s equation when initial data do not decay in a regular way at the spatial infinity. In the one-dimensional case, we show the existence of a solution whose profile and average speed are not convergent. In the higher-dimensional case, we show the existence of expanding fronts with arbitrarily prescribed profiles. We also show the existence of irregularly expanding fronts whose profile varies in time. Proofs are based on some estimate of the difference of two distinct solutions and a comparison technique. Dedicated to Professor Pavol Brunovsky on his 70th birthday.     

Nielsen方程的广义梯度表示及其稳定性分析

研究Nielsen方程的广义梯度表示以及方程零解稳定性．首先给出4类广义梯度系统及其性质．其次,给出完整系统和非完整系统的Nielsen方程转化成广义梯度系统的条件;将两类Nielsen方程分别化为广义梯度系统并研究方程的零解稳定性．最后,举例验证结果的应用并通过数值模拟验证结论的准确性．     

On Properties of the Generalized Wasserstein Distance

The Wasserstein distances W_p (p \({\geqq}\) 1), defined in terms of a solution to the Monge–Kantorovich problem, are known to be a useful tool to investigate transport equations. In particular, the Benamou–Brenier formula characterizes the square of the Wasserstein distance W₂ as the infimum of the kinetic energy, or action functional, of all vector fields transporting one measure to the other. Another important property of the Wasserstein distances is the Kantorovich–Rubinstein duality, stating the equality between the distance W₁(μ, ν) of two probability measures μ, ν and the supremum of the integrals in d(μ ?ν) of Lipschitz continuous functions with Lipschitz constant bounded by one. An intrinsic limitation of Wasserstein distances is the fact that they are defined only between measures having the same mass. To overcome such a limitation, we recently introduced the generalized Wasserstein distances \({W_p^{a,b}}\), defined in terms of both the classical Wasserstein distance W_p and the total variation (or L¹) distance, see (Piccoli and Rossi in Archive for Rational Mechanics and Analysis 211(1):335–358, 2014). Here p plays the same role as for the classic Wasserstein distance, while a and b are weights for the transport and the total variation term. In this paper we prove two important properties of the generalized Wasserstein distances: (1) a generalized Benamou–Brenier formula providing the equality between \({W_2^{a,b}}\) and the supremum of an action functional, which includes a transport term (kinetic energy) and a source term; (2) a duality à la Kantorovich–Rubinstein establishing the equality between \({W_1^{1,1}}\) and the flat metric.     

Crack Tip Equation of Motion in Dynamic Gradient Damage Models

We propose in this contribution to investigate the link between the dynamic gradient damage model and the classical Griffith’s theory of dynamic fracture during the crack propagation phase. To achieve this main objective, we first rigorously reformulate two-dimensional linear elastic dynamic fracture problems using variational methods and shape derivative techniques. The classical equation of motion governing a smoothly propagating crack tip follows by considering variations of a space-time action integral. We then give a variationally consistent framework of the dynamic gradient damage model. Owing to the analogies between the variational ingredients of these two models and under some basic assumptions concerning the damage band structuration, one obtains a generalized Griffith criterion which governs the crack tip evolution within the non-local damage model. Assuming further that the internal length is small compared to the dimension of the body, the previous criterion leads to the classical Griffith’s law through a separation of scales between the outer linear elastic domain and the inner damage process zone.     

Some Analytical Solutions for Groundwater Flow and Transport Equation

This paper presents the use of symmetry reduction method resulting in new exact solutions for the groundwater flow and transport equation. It is assumed that the radionuclides are transported by advection-diffusion in a single fracture and diffusion in the surrounding rock-matrix. The application of one-parameter group reduces the number of independent variables, and consequently the governing PDE of (1+2)-dimension reduces to set of ODEs which are solved analytically. This enables us to present some new exact time-dependent solutions of the advection-diffusion equation.     

Experimental Study of Nonlinear Flow in Micropores Under Low Pressure Gradient

As throat radius decrease to micro-nanoscale, seepage in unconventional reservoirs such as ultra-low permeability and tight reservoirs differs from conventional ones. Flow experiment in micropores is a promising approach to study characteristics of microflow. In this paper, a visual experimental device was established. Water flow through micropores with radius of 1.38–10.03 \(\upmu \hbox {m}\) was investigated, under 0.033–16 MPa/m. The results showed that in microscale, water flow did not agree with Poiseuille equation. Flow rate was lower than theoretical value and showed nonlinear characteristics. In the near wall area, due to the attraction of solid wall, a stagnant fluid layer was formed. It occupied flow space and thus lowered flow rate. Its thickness declined with pressure gradient increasing, which led to nonlinear flow characteristics. When the pressure gradient was very high, the thickness stopped declining and kept constant. Afterward, the flow transited to linear. In pores with smaller radius, the steady stagnant layer was thinner, but took a larger proportion of the flow space. For tubes of \(r = 1.38, 4.81, 10.03\,\upmu \hbox {m}\), the thickness of steady stagnant layer was 0.11, 0.23, 0.27 \(\upmu \hbox {m}\), respectively.     

Applicability of the Linearized Governing Equation of Gas Flow in Porous Media

In order to make the non-linear gas flow Equation tractable, the linearization treatment has been commonly applied in many subsurface gas flow problems such as natural gas production, soil vapor extraction, barometric, and pneumatic pumping. In this study, the accuracies of two representative linearization methods denoted as the conventional and the Wu solutions (Wu et al. Transp. Porous Media 32(1):117–137, 1998), are investigated quantitatively based on a numerical solution. The conventional solution uses a linearized constant gas diffusivity, while the Wu solution employs a spatially averaged but time-dependent gas diffusivity. The numerical solution is obtained by implementing the stiff solver ODE15s in MATLAB to deal with the time derivative and using the finite-difference method to approximate the spatial derivative in the non-linear gas flow equation. Two scenarios, the one-dimensional gas flow with constant pressure difference between two boundaries and the one-dimensional radial gas flow with constant mass injection rate at the origin of the coordinate system, are considered. The percentage error, defined as the ratio of difference between the numerical solution and the linearization solution to the ambient pressure, is calculated. It is founded that the Wu solution generally provides more accurate pressure evaluation than the conventional solution. The conventional solution always underestimates the pressure, while the Wu solution generally underestimates the pressure near the higher pressure boundary and overestimates the pressure near the lower pressure boundary. The maximal percentage error of the conventional solution is insensitive to time. This observation can be explained through the property of the complementary error function involved in the convention solution. For the one-dimensional flow example, the maximal percentage error of the conventional solution is 1.7, 25.5, and 90% when the pressure at one boundary suddenly rises above the ambient pressure by 50, 200, and 400%, respectively. While for the same example, the maximal percentage error of the Wu solution is 1.1, 14, and 44%, respectively.     

Isentropic Gas Flow for the Compressible Euler Equation in a Nozzle

We study the motion of isentropic gas in a nozzle. Nozzles are used to increase the thrust of engines or to accelerate a flow from subsonic to supersonic. Nozzles are essential parts for jet engines, rocket engines and supersonicwind tunnels. In the present paper, we consider unsteady flow, which is governed by the compressible Euler equation, and prove the existence of global solutions for the Cauchy problem. For this problem, the existence theorem has already been obtained for initial data away from the sonic state, (Liu in Commun Math Phys 68:141–172, 1979). Here, we are interested in the transonic flow, which is essential for engineering and physics. Although the transonic flow has recently been studied (Tsuge in J Math Kyoto Univ 46:457–524, 2006; Lu in Nonlinear Anal Real World Appl 12:2802–2810, 2011), these papers assume monotonicity of the cross section area. Here, we consider the transonic flow in a nozzle with a general cross section area. When we prove global existence, the most difficult point is obtaining a bounded estimate for approximate solutions. To overcome this, we employ a new invariant region that depends on the space variable. Moreover, we introduce a modified Godunov scheme. The corresponding approximate solutions consist of piecewise steady-state solutions of an auxiliary equation, which yield a desired bounded estimate. In order to prove their convergence, we use the compensated compactness framework.     

The Analytical Solution of the Boussinesq Equation for Flow Induced by a Step Change of the Water Table Elevation Revisited

The classical problem of flow induced by a sudden change of the piezometic head in a semi-infinite aquifer is re-examined. A new analytical solution is derived, by combining an expression describing the water table elevation upstream, obtained by the Adomian’s decomposition approach, to an existing polynomial expression (Tolikas et al. in Water Resour Res 20:24–28, 1984), adequate for the downstream region; the parameters of both approximations are computed by matching the two solutions at the inflection point of the water table. Although several analytical solutions are available in the literature, we demonstrate that the expression we have developed in this issue is the most accurate for strong or moderate non-linear flows, where the degree of non-linearity is defined as the ratio of the piezometric head elevation at the origin to the initial water table elevation. For this type of flows the perturbation-series solution of Polubarinova-Kochina, characterized by previous studies as the best available analytical solution provides physically unacceptable results, while the analytical solution of Lockington (J Irrig Drain Eng 123:24–27, 1997), used to check the accuracy of numerical schemes, underestimates the penetration distance of the recharging front.     

Remark on Single Exponential Bound of the Vorticity Gradient for the Two-Dimensional Euler Flow Around a Corner

In this paper, we consider the two–dimensional Euler flow under a simple symmetry condition, with hyperbolic structure in a unit square \({D = \{(x_1,x_2):0 < x_1+x_2 < \sqrt{2},0 < -x_1+x_2 < \sqrt{2}\}}\). It is shown that the Lipschitz estimate of the vorticity on the boundary is at most a single exponential growth near the stagnation point.     

Stability Analysis of a Vertical Curved Channel Flow with a Radial Temperature Gradient

The linear stability analysis of a Newtonian incompressible fluid in a vertical curved channel formed by two coaxial cylindrical surfaces with a radial temperature gradient and an azimuthal pressure gradient shows that critical modes are oscillatory and non-axisymmetric. We have derived a generalized Rayleigh discriminant which includes both the curvature and buoyancy effects. Centrifugal buoyancy induces weak asymmetry of the dependence of the control parameter critical values on the sign of the temperature gradient. The critical parameters depend on the temperature gradient, the radius ratio and the nature of the fluid. For a wide curvature channel flow, there are two critical modes: oscillatory Dean modes for small temperature gradients and oscillatory centrifugal-thermal modes for relatively large temperature gradients. Received 14 November 2001 and accepted 29 March 2002 Published online: 2 October 2002 Communicated by H.J.S. Fernando     

Steady Groundwater Flow to Drains on a Sloping Bed: Comparison of Solutions Based on Boussinesq Equation and Richards Equation

Analytical expressions for the scaling factor (A) in the Wooding and Chapman (J Geophys Res 71:2895–2902, 1966) solution for steady-state flow to drains on a sloping bed are presented. Otherwise A needs to be obtained by matching numerical and solutions. Corrections to various errors in other analytical solutions are given. The HYDRUS2D numerical model was used to generate results for steady-state flow to drains on a sloping bed which were compared to published Hele-Shaw cell results. The numerical results were used to compute both the pressure head on the bottom and the height of the phreatic surface. The numerical results for maximum water-table height are almost exactly the same as the published Hele-Shaw cell results and are greater than the numerical values for the maximum pressure heads on the sloping base. These HYDRUS2D model results were then compared with various analytical solutions, and it was found that Towner’s (Water Resour Res 11:144–147, 1975) solution gave the best results for both estimation of the maximum height of the phreatic surface and the position on the slope where this occurs.     

The Flow Structure Produced by Pulsed-jet Vortex Generators in a Turbulent Boundary Layer in an Adverse Pressure Gradient

The effect of pulsed jet vortex generators on the structure of an adverse pressure gradient turbulent boundary layer flow was investigated. Two geometrically optimised vortex generator configurations were used, co-rotating and counter-rotating. The duty cycle and pulse frequency were both varied and measurements of the skin friction (using hot films) and flow structure (using stereo PIV) were performed downstream of the actuators. The augmentation of the mean wall shear stress was found to be dependent on the net mass flow injected by the actuators. A quasi steady flow structure was found to develop far downstream of the injection location for the highest pulse frequency tested. The actuator near field flow structure was observed to respond very quickly to variations in the jet exit velocity.     

Solution of the Reynolds Equation for Viscous-Fluid Flow in the Clearance between a Rotating Ellipsoid and a Surface

Slow viscous-fluid flows in the narrow clearance (i) between a moving ellipsoid and a straight tube of elliptic cross section and (ii) between a rotating ellipsoid and a toroidal tube, including the case of an ellipsoid near a plane, are considered. A solution of the boundary-value problem for the Reynolds equation describing the flow in the clearance is found. The similarity of the pressure profiles in the “ellipsoid-plane” and “ cylinder-plane” systems is indicated.     

Finite Element Modelling of Flow Through a Porous Medium Between Two Parallel Plates Using The Brinkman Equation

Despite the widespread use of the Darcy equation to model porous flow, it is well known that this equation is inconsistent with commonly prescribed no slip conditions at flow domain walls or interfaces between different sections. Therefore, in cases where the wall effects on the flow regime are expected to be significant, the Darcy equation which is only consistent with perfect slip at solid boundaries, cannot predict velocity and pressure profiles properly and alternative models such as the Brinkman equation need to be considered. This paper is devoted to the study of the flow of a Newtonian fluid in a porous medium between two impermeable parallel walls at different Darcy parameters (Da). The flow regime is considered to be isothermal and steady. Three different flow regimes can be considered using the Brinkman equation: free flow (Da > 1), porous flow (high permeability, 1 > Da > 10⁻⁶) and porous flow (low permeability Da < 10⁻⁶). In the present work the described bench mark problem is used to study the effects of solid walls for a range of low to high Darcy parameters. Both no-slip and slip conditions are considered and the results of these two cases are compared. The range of the applicability of the Brinkman equation and simulated results for different cases are shown.     

Bifurcation Corresponding to the Onset of Asymmetric Periodic Structures and a Self-Induced Pressure Gradient in a Modified Taylor Flow

With reference to the example of a modified Taylor flow, the bifurcation of the loss of flow symmetry with the onset of a self-induced pressure gradient is studied theoretically and numerically. A linear analysis shows that the bifurcation is supercritical. It is necessarily accompanied by the appearance of a longitudinal pressure gradient and takes place at values of the parameters for which the solution of the linear system for the perturbations satisfies the condition of zero mass flow. It is established that, as a result of the bifurcation, two asymmetric solutions with oppositely directed pressure gradients are simultaneously generated. In the supercritical region, the symmetric branch of the solutions is also retained but becomes unstable. Bifurcation of the loss of symmetry and a self-induced pressure gradient can occur only in nonlinear systems.     

Pore Network Modeling of the Effects of Viscosity Ratio and Pressure Gradient on Steady-State Incompressible Two-Phase Flow in Porous Media

We perform steady-state simulations with a dynamic pore network model, corresponding to a large span in viscosity ratios and capillary numbers. From these simulations, dimensionless steady-state time-averaged quantities such as relative permeabilities, residual saturations, mobility ratios and fractional flows are computed. These quantities are found to depend on three dimensionless variables, the wetting fluid saturation, the viscosity ratio and a dimensionless pressure gradient. Relative permeabilities and residual saturations show many of the same qualitative features observed in other experimental and modeling studies. The relative permeabilities do not approach straight lines at high capillary numbers for viscosity ratios different from 1. Our conclusion is that this is because the fluids are not in the highly miscible near-critical region. Instead they have a viscosity disparity and intermix rather than forming decoupled, similar flow channels. Ratios of average mobility to their high capillary number limit values are also considered. Roughly, these vary between 0 and 1, although values larger than 1 are also observed. For a given saturation, the mobilities are not always monotonically increasing with the pressure gradient. While increasing the pressure gradient mobilizes more fluid and activates more flow paths, when the mobilized fluid is more viscous, a reduction in average mobility may occur.

     

Nonlinear Multigrid Methods for Numerical Solution of the Unsaturated Flow Equation in Two Space Dimensions

Picard and Newton iterations are widely used to solve numerically the nonlinear Richards’ equation (RE) governing water flow in unsaturated porous media. When solving RE in two space dimensions, direct methods applied to the linearized problem in the Newton/Picard iterations are inefficient. The numerical solving of RE in 2D with a nonlinear multigrid (MG) method that avoids Picard/Newton iterations is the focus of this work. The numerical approach is based on an implicit, second-order accurate time discretization combined with a second-order accurate finite difference spatial discretization. The test problems simulate infiltration of water in 2D unsaturated soils with hydraulic properties described by Broadbridge–White and van Genuchten–Mualem models. The numerical results show that nonlinear MG deserves to be taken into consideration for numerical solving of RE.     

Hilbert-Huang变换方法在Duffing方程和圆柱尾迹流分析中的应用

本文采用N．E Huang最近提出的Hilbert-Huang变换方法分析了使用三阶Runge-Kutta方法数值求解的Duffing方程的性质和圆柱绕流尾流特性,分解得到了有限阶经验模态（Empirical Modes）,给出了相应的Hilbert谱。结果表明,Hilbert-Huang变换得到的主要含能模态物理意义清晰,相应的瞬时频率存在着明显的波内频率调制,很好地刻画了信号的时空局部特性。     

Development of an Intermittency Equation for the Modeling of the Supersonic/Hypersonic Boundary Layer Flow Transition

An intermittency transport equation is developed in this study to model the laminar-turbulence boundary layer transition at supersonic and hypersonic conditions. The model takes into account the effects of different instability modes associated with the variations in Mach numbers. The model equation is based on the intermittency factor γ concept and couples with the well-known SST k–ω eddy-viscosity model in the solution procedures. The particular features of the present model approach are that: (1) the fluctuating kinetic energy k includes the non-turbulent, as well as turbulent fluctuations; (2) the proposed transport equation for the intermittency factor γ triggers the transition onset through a source term; (3) through the introduction of a new length scale normal to wall, the present model employs the local variables only avoiding the use of the integral parameters, like the boundary layer thickness δ, which are often cost-ineffective with the modern CFD methods; (4) in the fully turbulent region, the model retreats to SST model. This model is validated with a number of available experiments on boundary layer transition including the incompressible, supersonic and hypersonic flows past flat plates, straight/flared cones at zero incidences, etc. It is demonstrated that the present model can be successfully applied to the engineering calculations of a variety of aerodynamic flow transition with a reasonably wide range of Mach numbers.