Capillary waves past a flat plate in water of finite depth

Free-surface flow past a semi-infinite flat plate in a channelof finite depth is considered. The fluid is assumed to be inviscidand incompressible, and the flow to be two-dimensional and irrotational.Surface tension is included in the dynamic boundary conditionbut the effects of gravity are neglected. It is shown that thereis a three-parameter family of solutions with waves in the farfield and a discontinuity in slope at the separation point.This family includes as particular cases the solutions previouslycomputed by Osborn & Stump (2001, Phys. Fluids, 13, 616–623)and by Andersson & Vanden-Broeck (1996, Proc. R. Soc., 452,1985–1997).     

Stability of uniformly Morse-Smale gradient-like numerical methods for flows

In Garay (1996, Numer. Math., 72, 449–479) and Li (1997b,SIAM J. Math. Anal., 28, 381–388), it was shown that thequalitative properties of a Morse–Smale gradient-likeflow are preserved by its numerical approximations. In thispaper, we show that the qualitative properties of a family ofuniformly Morse–Smale gradient-like numerical methodsare preserved by the approximated flow. The techniques usedin the study of the structural stability theorem for diffeomorphismsare the main tools for this work.     

On the asymptotic behaviour for an electromagnetic system with a dissipative boundary condition

In this work we study some properties of solutions for the systemdescribing a three-dimensional non-homogeneous non-conductingdielectric with a general boundary condition with memory. Wefirst show the existence of the inverse of this boundary condition,which allows us to introduce a boundary free energy, similarto the one considered by Fabrizio & Morro (1996, Arch. Rat.Mech. Anal., 136, 359–381). Then, we prove existence anduniqueness theorems for weak and strong solutions of the evolutiveproblem in a finite time interval. Moreover, following Rivera& Olivera (1997, Boll. U.M.I., 11-A, 115–127), weexamine some dissipative properties of the boundary conditionand of its inverse and we give a useful energy estimate. Finally,when there is no memory in the boundary condition the exponentialdecay of the solution is proved.     

On the role of separation in compression wave generation by a train entering a tunnel hood with a window

An approximate analytical theory is proposed for calculatingthe compression wave generated when a train enters a tunnelfitted with an entrance hood with an open window. The pressurerise ahead of the entering train causes air to exhaust fromthe window in the form of a high-speed jet. The profile of thecompression wave transmitted into the tunnel is modified by theinteraction of the train nose with the window, by multiple reflectionsof wave energy between the window and the hood portal priorto transmission into the tunnel, and in addition by the productionof a pressure pulse by the jet. The wave generation problemcan be formulated in a quasi-one-dimensional manner, wherebythe pressure field generated in front of and to the sides ofthe train in the absence of the window is assumed to be scatteredby the window. A self-consistent solution is obtained by evaluatingthe jet flow from the window using a nonlinear empirical quationproposed and validated by Cummings (1984, Amer. Inst. Aeron.Astron. J., 22, 786–792; 1986, J. Acoust. Soc. Amer.,79, 942–951) for the velocity in the window-exit plane.Predictions are found to be in excellent agreement with measurementsof compression wave profiles obtained in model scale experimentsreported by Howe et al. (2003, J. Fluid Mech., 487, 211–243)at train speeds 350 km h^–1.     

Multiple Positive Solutions of Semilinear Differential Equations with Singularities

The existence of positive solutions of a second order differentialequation of the form z'+g(t)f(z)=0 (1.1) with the separated boundary conditions: z(0) – ßz'(0)= 0 and z(1)+z'(1) = 0 has proved to be important in physicsand applied mathematics. For example, the Thomas–Fermiequation, where f = z^3/2 and g = t^–1/2 (see [12, 13, 24]),so g has a singularity at 0, was developed in studies of atomicstructures (see for example, [24]) and atomic calculations [6].The separated boundary conditions are obtained from the usualThomas–Fermi boundary conditions by a change of variableand a normalization (see [22, 24]). The generalized Emden–Fowlerequation, where f = z^p, p > 0 and g is continuous (see [24,28]) arises in the fields of gas dynamics, nuclear physics,chemically reacting systems [28] and in the study of multipoletoroidal plasmas [4]. In most of these applications, the physicalinterest lies in the existence and uniqueness of positive solutions.     

Optimal error estimates for the hp-version interior penalty discontinuous Galerkin finite element method

We consider the hp-version interior penalty discontinuous Galerkinfinite-element method (hp-DGFEM) for second-order linear reaction–diffusionequations. To the best of our knowledge, the sharpest knownerror bounds for the hp-DGFEM are due to Rivière et al.(1999,Comput. Geosci., 3, 337–360) and Houston et al.(2002,SIAM J. Numer. Anal., 99, 2133–2163). These are optimalwith respect to the meshsize h but suboptimal with respect tothe polynomial degree p by half an order of p. We present improvederror bounds in the energy norm, by introducing a new functionspace framework. More specifically, assuming that the solutionsbelong element-wise to an augmented Sobolev space, we deducefully hp-optimal error bounds.     

Shape-preserving interpolation by G1 and G2 PH quintic splines

The interpolation of a planar sequence of points p₀, ..., p_Nby shape-preserving G¹ or G² PH quintic splines with specifiedend conditions is considered. The shape-preservation propertyis secured by adjusting ‘tension’ parameters thatarise upon relaxing parametric continuity to geometric continuity.In the G² case, the PH spline construction is based on applyingNewton–Raphson iterations to a global system of equations,commencing with a suitable initialization strategy—thisgeneralizes the construction described previously in NumericalAlgorithms 27, 35–60 (2001). As a simpler and cheaperalternative, a shape-preserving G¹ PH quintic spline schemeis also introduced. Although the order of continuity is lower,this has the advantage of allowing construction through purelylocal equations.     

On the global stability of a temporal discretization scheme for the Navier-Stokes equations

The time discretization by a linear backward Euler scheme forthe non-stationary viscous incompressible Navier–Stokesequations with a non-zero external force in a bounded 2D domainwith no-slip boundary condition or periodic boundary conditionis studied. Improved global stability results are obtained. The boundedness of the solution sequence in V and D(A) normsuniform with respect to &t for t [0, ) is proved. A similarresult in the V norm was previously obtained by (Geveci, 1989Math. Comp., 53, 43–53) for the non-forced system. A differentapproach is used here. As a corollary, the global attractorfor the approximation scheme is proved to exist, which is boundedin both V and D(A) spaces, thus compact in both H and V spaces.Applying the same techniques developed here, we are able toimprove the main result of (Hill and Süli 2000 IMA J. Numer.Anal., 20, 633–667) by showing that besides the existenceof a global attractor, the whole solution sequence is uniformlybounded in V as well, which is of significance from the pointof view of computing. As a corollary of local convergence results,upper semi-continuity of the attractor with respect to the numericalperturbation induced by the linear scheme is also establishedin both H and V spaces. Finally, some preliminary estimates,which are to our knowledge the first of their kind, on the dimensionsof the attractors in H and V spaces are also obtained.     

Modular Form Congruences and Selmer Groups

Motivated by Cremona and Mazur's notion of visibility of elementsin Shafarevich–Tate groups [6, 27], there have been anumber of recent works which test its compatibility with theBirch and Swinnerton–Dyer conjecture and the Bloch–Katoconjecture. These conjectures provide formulas for the ordersof Shafarevich–Tate groups in terms of values of L-functions.For example, one may see recent work of Agashe, Dummigan, Steinand Watkins [1, 2, 10, 11]. In their examples, they find thatthe presence of visible elements agrees with the expected divisibilityproperties of the relevant L-values.     

A Functional Relation for Accessory Parameters for Genus 2 Algebraic Curves with an Order 4 Automorphism

Let C be a genus 2 algebraic curve defined by an equation ofthe form y² = x(x² – 1)(x – a)(x – 1/a). Asis well known, the five accessory parameters for such an equationcan all be expressed in terms of a and the accessory parameter b corresponding to a. The main result of the paper is thatif a' = 1 – a², which in general yields a non-isomorphiccurve C', then b'a'(a'² – 1) = – – ba(a²– 1). This is proven by it being shown how the uniformizing functionfrom the unit disk to C' can be explicitly described in termsof the uniformizing function for C.     

On the Number of Solutions of an Equation Over a Finite Field

The Stöhr–Voloch approach is used to obtain a newbound for the number of solutions in (F_q)² of an equation f(X,Y) = 0, where f(X, Y) is an absolutely irreducible polynomialwith coefficients in a finite field F_q.     

Semi-discrete central-upwind schemes with reduced dissipation for Hamilton-Jacobi equations

We introduce a new family of Godunov-type semi-discrete centralschemes for multidimensional Hamilton–Jacobi equations.These schemes are a less dissipative generalization of the central-upwindschemes that have been recently proposed in Kurganov, Noelleand Petrova (2001, SIAM J. Sci. Comput., 23, pp. 707–740).We provide the details of the new family of methods in one,two, and three space dimensions, and then verify their expectedlow-dissipative property in a variety of examples.     

Rayleigh waves having generalised lateral dependence

It is shown that the surface-guided elastic waves found by Kiselevfor isotropic materials and having displacements depending linearlyupon the Cartesian coordinate orthogonal to the sagittal planemay be generalised in many ways. For surface waves on any anisotropichalf-space, a simple procedure applied to the displacementswithin the standard surface wave having dependence eⁱ, where k · x – t and k is the (surface) wave vector,yields displacements depending linearly upon the surface cartesiancoordinate orthogonal to the group velocity vector. Moreover,repeated application of this (differentiation) procedure yieldsa hierarchy of waves having algebraic dependence of successivelyincreasing degree. For isotropic materials, substantial simplificationand generalization are possible. Solutions of all algebraicdegrees have identical depth dependence. This allows the solutionsto be constructed iteratively and motivates a search for generalsolutions having depth dependence of the normal displacementthe same as in the standard surface wave. The procedure givesa new derivation of the solutions found by Achenbach havingamplitude of the normal displacement of the surface given byany solution to the two-dimensional Helmholtz equation. Furthermore,exploiting the scale invariance (a property of surface waveson any homogeneous half-space) shows that in every surface-guideddisturbance of an elastic half-space, the elevation of the freesurface is a solution of the wave equation in two dimensions(the membrane equation). Using the paraxial approximation tothe membrane equation, high-frequency Rayleigh waves propagatingas narrow beams are described in terms of a scalar Gaussianbeam.     

The effects of a complexing agent on travelling waves in autocatalytic systems with applied electric fields

The effects of electric fields on the reaction fronts that arisein a system governed by an autocatalytic reaction and a complexationreaction between the autocatalyst and a complexing agent areconsidered. The complexation reaction is assumed to be fastrelative to the autocatalytic reaction and the equations forthis limit are derived. The corresponding travelling waves arediscussed, the case of quadratic autocatalysis being treatedin detail. The existence of minimum speed waves is examined,being dependent on the ratio of diffusion coefficients D, theconcentration S₀ and equilibrium constant K of the complexationreaction as well as the electric field strength E. It is seenthat, for some parameter values, minimum speed waves have negativeautocatalayst concentrations, and waves which have the lowestspeed consistent with non-negative concentrations are also obtained.Numerical integrations of the initial-value problem are performedfor representative parameter values. These show the developmentof the appropriate travelling wave (when it exists) as the largetime behaviour of the system, and, in cases where no travellingwave exists, the numerical integrations show the electrophoreticseparation of substrate and autocatalyst.     

Local and global behaviour of steady-state solutions of the Sel'kov model

In this paper we discuss steady-state solutions of the systemof reaction-diffusion equations known as the Sel'kov model.This model has been the subject of much discussion; in particular,analytical and numerical results have been discussed by Lopez-Gomezet al. (1992, IMA J. Num. Anal. 12, 405–28). We show thata simple analysis of the bifurcation function associated withthe system can explain many of the numerical observations, suchas the formation and development of loops of nontrivial solutions,in a simpler and more complete manner than the analysis of Lopez-Gomezet al. This allows for a clearer understanding of the qualitativebehaviour of the set of nontrivial solutions and hence of thebifurcation diagram.     

Locally Nilpotent p-Groups whose Proper Subgroups are Hypercentral or Nilpotent-by-Chernikov

Let G be a group and P be a property of groups. If every propersubgroup of G satisfies P but G itself does not satisfy it,then G is called a minimal non-P group. In this work we studylocally nilpotent minimal non-P groups, where P stands for ‘hypercentral’or ‘nilpotent-by-Chernikov’. In the first case weshow that if G is a minimal non-hypercentral Fitting group inwhich every proper subgroup is solvable, then G is solvable(see Theorem 1.1 below). This result generalizes [3, Theorem1]. In the second case we show that if every proper subgroupof G is nilpotent-by-Chernikov, then G is nilpotent-by-Chernikov(see Theorem 1.3 below). This settles a question which was consideredin [1–3, 10]. Recently in [9], the non-periodic case ofthe above question has been settled but the same work containsan assertion without proof about the periodic case. The main results of this paper are given below (see also [13]).     

Error expansion for a first-order upwind difference scheme applied to a model convection-diffusion problem

A singularly perturbed convection–diffusion problem isconsidered. The problem is discretized using a simple first-orderupwind difference scheme on general meshes. We derive an expansionof the error of the scheme that enables uniform error boundswith respect to the perturbation parameter in the discrete maximumnorm for both a defect correction method and the Richardsonextrapolation technique. This generalizes and simplifies resultsobtained in earlier publications by Fröhner et al.(2001,Numer. Algorithms, 26, 281–299) and by Natividad &Stynes (2003, Appl. Numer. Math., 45, 315–329). Numericalexperiments complement our theoretical results.     

Optimal control of parabolic equation with an infinite number of variables for non-standard functional and time delay

A distributed control problem for the parabolic operator withan infinite number of variables and time delay is considered.The performance index has an integral form. Constraints on controlsare imposed. To obtain optimality conditions for the Neumannproblem, the generalization of the Dubovitskii–Milyutintheorem given by Walczak in WALCZAK, S. Folia Mathematics, 1,187–196 and WALCZAK, S. J. Optim. Theory Appl., 42, 561–582was applied.     

Semidiscrete finite element Galerkin approximations to the equations of motion arising in the Oldroyd model

In this paper, a semidiscrete finite element Galerkin methodfor the equations of motion arising in the 2D Oldroyd modelof viscoelastic fluids with zero forcing function is analysed.Some new a priori bounds for the exact solutions are derivedunder realistically assumed conditions on the data. Moreover,the long-time behaviour of the solution is established. By introducinga Stokes–Volterra projection, optimal error bounds forthe velocity in the L(L²) as well as in the L(H¹)-norms andfor the pressure in the L(L²)-norm are derived which are validuniformly in time t > 0.     

Funk-Hecke Formula for Orthogonal Polynomials on Spheres and on Balls

Analogues of the Funk–Hecke formula for spherical harmonicsare proved for Dunkl's h-harmonics associated to the reflectiongroups, and for orthogonal polynomials related to h-harmonicson the unit ball. In particular, an analogue and its applicationare discussed for the weight function (1–|x|²)^µ–1/2on the unit ball in R^d. 2000 Mathematics Subject Classification33C50, 33C55, 42C10.