Blow up of solutions of a generalized Boussinesq equation

Consider the Cauchy problem u_tt = (f(u))_xx + u_xxtt x R, t 0, u(x,0) = u₀(x), u_t(x,0) = u₁(x),7rcub; where f : R R C, f(0) = ). After treatment of the local existenceproblem, we show the blow up of the solution of the equation(1) under the folowing assumptions. Let > 0 be real, such that 2(l + 2)F(u) uf(u), (v₀, Pv₀)_l² + _- F(u₀)dx < 0 where P = 1 - d²/dx², F'(s), and v₀ is given by u₁(x,0) = (v₀(x))_x. Then we focus on various perturbations of the question. We alsostudy the vectorial case in the same way, and finally we giveexamples.     

On the behaviour of blow-up interfaces for an inhomogeneous filtration equation

We study the asymptotic behaviour of blow-up interfaces of thesolutions to the one-dimensional nonlinear filtration equationin inhomogeneous media where m>1 isa constant and (x) = |x|^– (for |x| 1, with > 2) isa bounded, positive, smooth, and symmetric function. The initialdata are assumed to be smooth, bounded, compactly supported,symmetric, and monotone. It is known that due to the fast decayof the density (x) as |x| the support of the solution increasesunboundedly in a finite time T. We prove that as tT^– theinterface behaves like O((T – t)^–b), where the exponentb > 0 (which depends on m and only) is given by a uniqueself-similar solution of the second kind satisfying the equation|x|^– u_t = (u^m)_xx. The corresponding rescaled profilesalso converge. We establish the stability of the self-similarsolution of the second kind for the exponential density (x)=e^–|x|for |x| 1. We give a formal asymptotic analysis of the blow-upbehaviour for the non-self-similar density (x) = e^–|x|2.Several exact self-similar solutions and their correspondingasymptotics are constructed.     

Discontinuous Galerkin finite element approximation of quasilinear elliptic boundary value problems I: the scalar case

^** Email: Paul.Houston{at}mcs.le.ac.uk^*** Email: Janice.Robson{at}comlab.ox.ac.uk^**** Email: Endre.Suli{at}comlab.ox.ac.uk We develop a one-parameter family of hp-version discontinuousGalerkin finite element methods, parameterised by [–1,1], for the numerical solution of quasilinear elliptic equationsin divergence form on a bounded open set ^d, d 2. In particular,we consider the analysis of the family for the equation –·{µ(x, |u|)u} = f(x) subject to mixed Dirichlet–Neumannboundary conditions on . It is assumed that µ is a real-valuedfunction, µ C( x [0, )), and thereexist positive constants m_µ and M_µ such that m_µ(t– s) µ(x, t)t – µ(x, s)s M_µ(t– s) for t s 0 and all x . Using a result from the theory of monotone operators for any valueof [–1, 1], the corresponding method is shown to havea unique solution u_DG in the finite element space. If u C¹() H^k(), k 2, then with discontinuous piecewise polynomials ofdegree p 1, the error between u and u_DG, measured in the brokenH¹()-norm, is (h^s–1/p^k–3/2), where 1 s min {p+ 1, k}.     

Numerical Evaluation of Fourier Integrals

A method is developed for evaluating Fourier integrals of theform A() = ¹_–1f(x) e^fax dx, 0. The method consists of expanding the function f in a seriesof Chebyshev polynomials and expressing the integral A() asa series of the Bessel functionsJ_r+(), r= 0, 1, 2,.... A partialsum A_N() of the series provides an approximant to A(). The principalfeature of the method is that one set of N+1 evaluations off(x) suffices for the calculation of A_N() for all , and alsothe truncation error A()–A_N() is essentially independentof . Numerical tests show that the method is accurate, economicaland reliable. An application to the inversion of Fourier andLaplace transforms is briefly described.     

On hearing the shape of a bounded domain with Robin boundary conditions

The asymptotic expansions of the trace of the heat kernel (t)= [sum ]_j=1 exp(-t_j) for small positive t, where {_j} _j=1 arethe eigenvalues of the negative Laplacian -_n = -[sum ]ⁿ_k=1 (/x^k)²in Rⁿ (n = 2 or 3), are studied for a general multiply connectedbounded domain which is surrounded by simply connected boundeddomains _i with smooth boundaries _i (i = 1,...,m), where smoothfunctions Y_i (i = 1,...,m) are assuming the Robin boundary conditions(n_i + Y_i) = 0 on _i. Here /n_i denote differentiations along theinward-pointing normals to _i (i = 1,...,m). Some applicationsof an ideal gas enclosed in the multiply connected bounded containerwith Neumann or Robin boundary conditions are given.     

The Model Plasma Problem: Existence of a Solution Branch

We are interested in the model plasma problem –u = u⁺in ,u = –d on , _au⁺ dx=j where is a bounded domain in with boundary ; here, j isa given positive number, the function u and the positive number are the unknowns of the problem, and d is a real parameter.Using a variant of the implicit function theorem, we can provethe existence of a global solution branch parametrized by d.The method has the advantage that it can be used for analysingthe approximation of the above problem by a finite-element method.     

Elliptic Problems Involving Phase Boundaries Satisfying a Curvature Condition

Two theorems related to equilibrium free-boundary problems arepresented. One arises as a time-independent solution to thephase-field equations. The other is the relevant time-independentproblem for the Stefan model, modified for the surface tensioneffect. It also serves as a preliminary result for the phase-fieldformulation. Under appropriate conditions, we prove that, givenan appropriate positive constant and a smooth function u: R;,where is an annular domain in R², there exists a curve suchthat u(x)=—K(x) for all x , where K is the curvature.Using this result, we prove the existence of solutions to O=²+ ?(—³) + 2u that have a transition layer behaviour (from=—1 to =+1) for small and make the transition on thecurve . This proves there exist solutions to the phase fieldmodel that satisfy a Gibbs-Thompson relation.     

A Criticality Criterion for Thermal Explosions in Well-stirred Systems with Reactant Consumption

The autonomous differential equations for the temperature andreactant consumption in a first-order well-stirred exothermicreaction are considered. An examination of the phase-plane solutionsallows the qualitative behaviour of the Semenov number as afunction of maximum temperature rise ^* to be established. Inthe limit of infinite adiabatic temperature rise (B) and zeroactivation energy parameter ( = 0), the relationship between and stationary temperature _s is known to be e_¹ = _s. Criticalityarises at the maximum of (_s) and leads to the critical Semenovvalues (_s)_cr = 1, _cr = e^–1. For sufficiently large B,it is shown that the (^*) curve has a bifurcation at ^* = 1, withthe upper branch monotonically increasing and the lower branchmonotonically decreasing for ^* > 1. In the limit B thesebecome respectively the straight line = e^–1, _s 1 andthe unstable branch of = _se^–1, _s 1 and the unstablebranch of = _s e^–. Criticality for finite B is definedas occurring at the bifurcation, namely ^*_cr = 1, with _cr(B)the value of at this point. Values of these Semenoy numbersare obtainable from the numerical calculations of Boddingtonet al. [Proc. R. Soc. Lond. (1983), 390, 13–30]. The newcriterion is applied to an approximate phase-plane solution.The corresponding critical parameter is found to be _cr = e^–1[1+B^–(2–e^–1)+O(B^–1)].     

Experimental Observation of the Large-Amplitude Solutions of Duffing's and Related Equations

Experiments with a nonlinear electronic model show that certainsimple features of the solutions of where f(u) is an odd monotonic function of u for example u³,repeat in a regular pattern as either is decreased or U isincreased. For fixed U, the position of these features is periodicin 1/ and, when f(u) has the form u|u|^k–1 a quantitativerelation between the period in 1/ and U can be found. The occurrenceof large-amplitude chaotic solutions is found to depend notonly on the nonlinearity of f(u) for large U but also on itsbehaviour near u = 0. For the Duffing equation, which can bereduced to the range of parameters accessible to experiment is 0<1 and0<F5000.     

Injection of Ideal Fluid from a Slot into a Stream: Two Free Boundaries

A fluid is injected from a slot into a stream of another fluid.In a simple model this leads to a two-phase two-free-boundaryproblem with the jump relation |u^–|² – |u⁺|² = on the free boundary {u=0}, and |u| = ¹ on the free boundary{u > – Q}, where u is the stream function and Q isthe flux of the injected fluid. Using the variational theoryof Alt, Caffarelli & Friedman, we prove existence of (,₁, u) such that there is a smooth fit for both free boundaries.     

The Mesa Problem: Diffusion Patterns for ut = {nabla} {middle dot} (um{nabla}u) as m -> +{infty}

In this paper we consider the limit m+ of solutions of the porous-mediumequation u_t = · (u^mu) (xR^N), with N > 1. We conjecturethat, for initial data with a unique maximum, the evolutionis characterized by the onset of a ‘mesa’ region,in which the solution is nearly spatially independent, surroundedby a region in which u is nearly equal to its initial value.The transition between these regions occurs near a surface whichis identified with the free boundary in a certain Stefan problemwhich can be studied using variational inequalities. Moreover,singular-perturbation theory can be used to describe the structureof the transition region.     

Global blow-up of separable solutions of the vorticity equation

In this paper we construct solutions to the equation on a finite interval in y which blow-up globallyin finite time. This equation arises in a number of physicalsituations and can be derived from the vorticity equation bylooking for stagnation-point type separable solutions for thetwo-dimensional streamfunction of the form xu(y, t). In theparticular application which has prompted the investigationreported in this paper, (*) is solved subject to boundary conditionsinvolving ²u/y². For this type of boundary condition the phenomenonof blow-up was first observed numerically by solving the initial-boundary-valueproblem for (*). These computations reveal that, depending onthe parameter combinations chosen, the solution to the initial-valueproblem may either blow-up globally in finite time or approacha steady state as t . Using the computations as a guide weconstruct the analytic behaviour of the solution close to theblow-up time using the methods of formal asymptotics.     

An Example on Sobolev Space Approximation

For each d2 we construct a connected open set R^d such that = int (clos()), and for each k 1 and each p [1, ), the subsetW^k, () fails to be dense in the Sobolev space W^{k, p}(), in thenorm of W^{k, p}(). 1991 Mathematics Subject Classification 46E35,46F05.     

Spreading speed and travelling wave solutions of a partially sedentary population

In this paper, we extend the population genetics model of Weinberger(1978, Asymptotic behavior of a model in population genetics.Nonlinear Partial Differential Equations and Applications (J.Chadam ed.). Lecture Notes in Mathematics, vol. 648. New York:Springer, pp. 47–98.) to the case where a fraction ofthe population does not migrate after the selection process.Mathematically, we study the asymptotic behaviour of solutionsto the recursion u_n+1 = Q_g[u_n], where In the above definition of Q_g, K is a probabilitydensity function and f behaves qualitatively like the Beverton–Holtfunction. Under some appropriate conditions on K and f, we showthat for each unit vector R^d, there exists a c^*_g() which hasan explicit formula and is the spreading speed of Q_g in thedirection . We also show that for each c c^*_g(), there existsa travelling wave solution in the direction which is continuousif gf '(0) 1.     

Convergence of a finite-element approximation of surfactant spreading on a thin film in the presence of van der Waals forces

We consider a fully practical finite-element approximationof the following system of nonlinear degenerate parabolic equations: (u)/(t) + . (u² [(v)]) - (1)/(3) .(u³ w)= 0, w = - c u - u^-+ a u^-3 , (v)/(t) + . (u v [(v)]) - v - .(u² v w) = 0. The above models a surfactant-driven thin-film flow in the presenceof both attractive, a>0, and repulsive, >0 with >3,van der Waals forces; where u is the height of the film, v isthe concentration of the insoluble surfactant monolayer and(v):=1-v is the typical surface tension. Here 0 and c>0 arethe inverses of the surface Peclet number and the modified capillarynumber. In addition to showing stability bounds for our approximation,we prove convergence, and hence existence of a solution to thisnonlinear degenerate parabolic system, (i) in one space dimensionwhen >0; and, moreover, (ii) in two space dimensions if inaddition 7. Furthermore, iterative schemes for solving the resultingnonlinear discrete system are discussed. Finally, some numericalexperiments are presented.     

On the Fredholm Alternative for the p-Laplacian in One Dimension

We investigate the existence of a weak solution u to the quasilineartwo-point boundary value problem We assume that 1 < p < p ¬ = 2, 0 < a < , andthat f L¹(0,a) is a given function. The number _k stands forthe k-th eigenvalue of the one-dimensional p-Laplacian. Let_{p p} x/a) denote the eigenfunction associated with ₁; then _p(k_p x/a) is the eigenfunction associated with _k. We show the existenceof solutions to (P) in the following cases. (i) When k=1 and f satisfies the orthogonality condition the set of solutions is bounded. (ii) If k=1 and f_t L¹(0,a) is a continuous family parametrizedby t [0,1], with then there exists some t^* [0,1] such that (P) has a solutionfor f = f_t^*. Moreover, an appropriate choice of t_* yields asolution u with an arbitrarily large L¹(0,a)-norm which meansthat such f cannot be orthogonal to _pp x/a. (iii) When k 2 and f satisfies a set of orthogonality conditionsto _p(k p x/a) on the subintervals , again, the set of solutions is bounded. is a continuous family satisfying either or another related condition, then there exists some t^* [0,1]such that (P) has a solution for f = f_t*. Prüfer's transformation plays the key role in our proofs.2000 Mathematical Subject Classification: primary 34B16, 47J10;secondary 34L40, 47H30.     

On a Comparison Principle and the Critical Fujita Exponents for Solutions of Semilinear Parabolic Inequalities

The paper examines blow-up phenomena for the inequality u_t–Lu–|u|^q–1u_t–L–||^q–1 (*) in the half-space x Rⁿ, n 1, where L is a linear second-order partial differential operatorin divergence form. The paper studies weak solutions of (*) that belong only locallyto the corresponding Sobolev spaces in the half-space x Rⁿ. It also requires no conditionsfor the behavior of solutions of (*) on the hyperplane t = 0. The existence of critical blow-up exponents is obtained forsolutions of (*) as a special case of a comparison principlefor the corresponding solutions of (*). For example, the well-knownFujita result is a consequence of the comparison principle. The approach developed in the paper is directly applicable tothe study of analogous problems involving nonlinear differentialoperators. Its elliptic analogue has been recently developedby the authors.     

Concavity and B-Concavity of Solutions of Quasilinear Filtration Equations

Spatial concavity properties of non-negative weak solutionsof the filtration equations with absorption u_t = ((u))_xx–(u)in Q = Rx(0, ), '0, 0 are studied. Under certain assumptionson the coefficients , it is proved that concavity of the pressurefunction is a consequence of a ‘weak’ convexityof travelling-wave solutions of the form V(x, t) = (x–t+a).It is established that the global structure of a so-called properset B = {V} of such particular solutions determines a propertyof B-concavity for more general solutions which is preservedin time. For the filtration equation u_t = ((u))_xx a semiconcavityestimate for the pressure, v_xx(t+)^–1'(), due to the B-concavityof the solution to the subset B of the explicit self-similarsolutions (x/t+)) is proved. The analysis is based on the intersection comparison based onthe Sturmian argument of the general solution u(x, t) with subsetsB of particular solutions. Also studied are other aspects ofthe B-concavity/convexity with respect to different subsetsof explicit solutions.     

Periodic solutions and oscillation of discrete non-linear delay population dynamics model with external force

^** Email: emelabbasy{at}mans.edu.eg^*** Email: shsaker{at}mans.edu.eg In this paper, we consider the discrete non-linear delay populationdynamics model [graphic: see PDF] where m is a positive integer, p(n), Q(n) and (n) are positiveperiodic sequences of period . By the method that involves theapplication of the Gaines and Mawhins coincidence degree theory,we prove that there exists a positive -periodic solution (n). We prove that every positive solutionof (*) which does not oscillate about (n)satisfies lim_t[y(n)–(n)]=0.We establish some necessary and sufficient conditions for theoscillation of every positive solution about (n), and finally, we establish the lower and upperbounds of the oscillatory solutions.     

Logistic Type Equations on RN by a Squeezing Method Involving Boundary Blow-Up Solutions

We study, on the entire space R^N(N 1), the diffusive logisticequation u_t–u=u–u^p, u0 (1.1) and its generalizations. Here p > 1 is a constant. Problem(1.1) plays an important role in understanding various populationmodels and some other problems in applied mathematics. When = 1 and p = 2, it is also known as the Fisher equation andKPP equation, due to the pioneering works of Fisher [8] andKolmogoroff, Petrovsky and Piscounoff [18].