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.     

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.     

On a Penalty-Perturbation Theory for Plate Problems

A penalty-perturbation method previously proposed by Westbrook(J. Inst. Maths Applics (1974) 14, 79–82) for the solutionsof static bending problems for elastic plates is analysed here.The method replaces the single fourth-order biharmonic equationby a system of three second-order equations which is "singularly"perturbed with respect to a small penalty parameter . The existenceof solutions of the perturbed problem for each > 0 is establishedand the behaviour of these solutions as 0 0 is studied. Inparticular, the results show that while these solutions arecontinuous in at = 0, analyticity in at = 0 is lost exceptin special cases.     

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.     

Optimal Parameter Choice in the {theta}-Method for Parabolic Equations

The general first-order method, known as the -method, is appliedto the semi-discrete form of a parabolic equation. It is shownthat to every required local accuracy there corresponds a valueof the parameter that is optimal in the sense of allowing thelargest step for which the error remains bounded below . Anasymptotic formula for in terms of is obtained, showing thatthe maximum step-size for the optimal -method is more than twiceas large as that for the Crank-Nicolson method. A numericalexample is given, showing good agreement between theory andpractice.     

The Critical Conditions in a Reactive Slab with Slowly Varying Surface Temperature

Liouville's non-linear partial differential equation is consideredfor an infinite rectangular strip domain with a slowly varyingboundary condition. The equation describes a layer of chemicallyreactive material under conditions where the resistance to surfaceheat transfer is negligible and the ambient temperature variesslowly along the surface. Symmetrical heating by a zero orderexothermic reaction is assumed. If is a small dimensionlesstemperature difference between regions where the surface temperatureis effectively constant, a perturbation series solution in may be determined provided the Frank-Kamenetskii parameter satisfies _c(). It is shown that a plausible value for thecritical parameter is _c() = _c(0) e^–e,where _c(0) = 0.878.The corresponding critical temperature distribution is shownto have a dependence on different from that for subcriticalcases.     

Existence of Solutions of a Two-Point Free-Boundary Problem Arising in the Theory of Porous Medium Combustion

The existence of solutions of a two-point free-boundary problemarising from the theory of travelling combustion waves in aporous medium is examined. The problem comprises a third-ordernonlinear ordinary differential equation posed on an unknowninterval of finite length; four boundary conditions are given,two at either end of the interval. The equations possess a trivialsolution for all values of the bifurcation parameter . A shootingtechnique is employed to prove the existence of a nontrivialsolution for 0 < < _c and nonexistence theorems are provedfor (0, _c).     

Numerical studies of the fourth Painlev equation

In this paper the authors investigate numerically solutionsof a special case of the fourth Painlev equation given by with a parameter, satisfying theboundary condition Equation (1a) arises asa symmetry reduction of the derivative nonlinear Schrdinger(DNLS) equation, which is a completely integrable soliton equationsolvable by inverse scattering techniques. Previous resultsconcerned with solutions of (1) are largely restricted to thecase when is an integer and very little has been proved when is a noninteger. Here a numerical approach to describing solutionsof (1) for noninteger is adopted, and information is obtainedcharacterizing connection formulae which describe how the asymptoticbehaviours of solutions as + relate to those as –.     

An analogue of Berry's conjecture for the phase in fractal obstacle scattering

In this paper, the authors consider the high-frequency asymptoticsof the phase s() of acoustic waves scattered by an obstacleRⁿ with fractal boundary. Under certain conditions, it is provedthat if is –Minkowski measurable with –Minkowskimeasure µ then there exists a positive constant C_n, dependingonlyon n and such that where     

A slender dielectric body embedded in an arbitrary external potential

This paper determines the electrostatic potential and field taking place both inside and outside a slender dielectric body embedded in a given potential ₀. This task is actually reducedto the determination of the occurring polarization surface-chargedensity qwhich depends on ₀, on the body shape but also onthe ratio = ₂/₁of the dielectic constants (₂outside and ₁inside the body). The adopted procedure consists in asymptoticallyexpanding and inverting (with respect to the small slendernessratio of the body) the well-known Fredholm boundary integralequation of the second kind governing the function q. The technicaldifficulties such an approach encounters are bypassed by employinga systematic formula in getting the asymptotic estimate ofcertain integrals depending upon a small parameter. Contraryto other works in the field, this method authorizes us to handlethe case of non-axisymmetric slender bodies. As an illustrationthe theory is applied to a body of elliptical cross-sectionand comparisons are presented for a slender dielectric ellipsoidembedded in a special potential ₀for which the exact densityqis obtained in a closed form.     

The large deformation of nonlinearly elastic rings in a two-dimensional compressible flow

The nonlinear nonlocal system of the equilibrium equations ofan elastic ring under the action of an external two-dimensionaluniformly subsonic potential barotropic steady-state gas flowis considered. The configurations of the elastic ring are identifiedby a pair of functions (, ). The simple curve represents theshape of the ring and the real-valued function identifies theorientation of the material sections of the ring. The pressurefield on the ring depends nonlocally on , and on two parametersU and P which represent the pressure and the velocity at infinity.The system is shown to be equivalent to a fixed-point problem,which is then treated with continuation methods. It is shownthat the solution branch ensuing from certain equilibrium states((₀, ₀), 0, P₀) in the solution-parameter space of ((₀, ₀),0, P₀) either approaches the boundary of the admissible ((,), U,p)'s in a well-defined sense, or is unbounded, or is homotopicallynontrivial in the sense that there exists a continuous map from the branch to a two-dimensional sphere which is not homotopicin the sphere to a constant, while restricted to the branchminus ((₀, ₀), 0, P₀) is homotopic to a constant in the sphere.Furthermore, by fixing the pressure parameter at P₀ and by consideringthe one-parameter problem in ((, ), U), the following holds.Every hyperplane in the solution-parameter space of the ((,), U)'s which contains the equilibrium state ((₀, ₀), 0) anddoes not include a welldetermined one-dimensional subspace intersectsthe solution branch above at a point different from ((₀, ₀),0).     

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.     

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.     

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)].     

Exact Remainders for Asymptotic Expansions of Fractional Integrals

This is a continuation of work begun in an earlier paper inwhich we used the theory of distributions to derive explicitexpressions for the remainder terms associated with the asymptoticexpansions of the Stieltjes transform. In this paper similarresults are obtained for the fractional integral of order definedby 1f(x)=1/f()^x_o(x-)^x-1 f(t)dt, >. Heref(t) is a locally integrable function on [0, ) and satisfies f(t) a_st-5–0(ó >0), s=0 as      

Asymptotic Convergences in a One-Phase Continuous-Casting Stefan Problem with High Extraction Velocity

Asymptotic convergence for variational solutions and for thefree boundaries in a one-phase continuous-casting problem ofthe Stefan type are derived when the Péclet number ß+. This corresponds, in particular, to a model with high extractionvelocity, where an ultraparabolic problem is obtained from aparabolic one and, in the stabilized case (i.e. when t +), aparabolic problem results from an elliptic one. The asymptoticstabilization in time as t + is also discussed in the limitß = +.     

Computer extended series solution for unsteady flow in a contracting or expanding pipe

Unsteady flow in a semi-infinite contracting or expanding pipeis reinvestigated using long series analysis. The proposed seriesmethod is useful in analysing the problem for a moderately largeconstant ( = aa/, where a = a(t), the radius of the pipe isa function of time, a(t) is the velocity of the wall, and iskinematic viscosity). For positive values of (expansion ofthe pipe) accuracy of the series representing shear stress andpressure gradient is increased from = 2.89 to = 6.0 by extractingthe singularity followed by completion of the series. For negativevalues of (contraction of the pipe), we revert the series whichresults into the increase of the region of validity of the transposedseries from = -25.0 to = -2.89. Later we use Padé approximantsfor summing them. Also, the asymptotic solution for large valuesof is obtained and it agrees closely with pure numerical valuesof shear stress at the wall and pressure gradient.     

Short-time Asymptotic Solutions of the Heat Conduction Equation with Spatially Varying Coefficients

An asymptotic solution of the heat conduction equation withspatially varying coefficients is developed for small times.The method followed consists of an application of the Laplacetransformation and use of the Liouville-Green approximationin the subdominant solution of the resulting second-order differentialequation. The approximate solution is inverted by contour integration.The resulting asymptotic expression has a time-dependence identicalto that applying to the case of constant properties providedthat an appropriately averaged value of the thermal diffusivity is used, namely, The error term is of the order of t where is a measure of spatialvariability of the coefficients in the heat equation.     

Fitted and Unfitted Finite-Element Methods for Elliptic Equations with Smooth Interfaces

This paper considers the finite-element approximation of theelliptic interface problem: -?(u) + cu = f in Rⁿ (n = 2 or3), with u = 0 on , where is discontinuous across a smoothsurface in the interior of . First we show that, if the meshis isoparametrically fitted to using simplicial elements ofdegree k - 1, with k 2, then the standard Galerkin method achievesthe optimal rate of convergence in the H¹ and L² norms overthe approximations _l⁴ of _l where _{l 2}. Second, since itmay be computationally inconvenient to fit the mesh to , weanalyse a fully practical piecewise linear approximation ofa related penalized problem, as introduced by Babuska (1970),based on a mesh that is independent of . We show that, by choosingthe penalty parameter appropriately, this approximation convergesto u at the optimal rate in the H¹ norm over _l⁴ and in the L²norm over any interior domain _l^* satisfying _l^* _l^** _l⁴ for somedomain _l^**. Present address: School of Mathematical and Physical Sciences,University of Sussex, Brighton BN1 9QH     

A Mathematical Theory of Full-wave Rectification

The results of a systematic asymptotic analysis of two basicfull-wave rectifier circuits are presented in this paper. Thephysical circuits are modelled by two systems of nonlinear differentialequations which are both singularly perturbed in the limit RC.Here is the angular frequency of the source, R is the resistanceof the load, and C is the capacitance of the circuit. The methodof multiple scales is used to construct the asymptotic behaviourof the solutions as RC. The results contain transient, ripple,and steady-state information. The latter two components simplifyand yield the standard results when the diodes' forward impedanceis small. The analysis of more complicated full-wave rectifiersfollows from the procedure described.