Features of electron density and temperature in the 500–3500 km region of the plasmasphere—(1) Dawn and dusk sectors

Dawn-dusk features of the plasmasphere are examined for intervals in February and September 1969, using electrostatic probe data ofN _e andT _e from the ISIS-I satellite. Clear plasmatrough formation is seen in the vicinity of 70° geomagnetic latitude in both dawn and dusk sectors in the 1500–3500 km region, but the plasmatrough is absent in the altitude range 500–1500 km. The plasmatrough minimum near 70°φ exhibits no asymmetry between dawn and dusk sectors in its latitudinal position. TheT _e peak associated with the plasmatrough is more pronounced in the dawn sector. DawnN _e is less than duskN _e, but dawnT _e exceeds duskT _e. The influence of processes in the magnetosphere in causing these features is examined.     

Dirichlet forms for diffusion in ℝ<Superscript>2</Superscript> and jumps on fractals: The regularity problem

We define a Dirichlet form ɛ describing diffusion in ℝ ^d and jumps in a fractal Γ ⊂ ℝ ^d . The jump measure J is defined as an image of a jump measure j of a process in a non-Archimedean metric space. As the result the jump intensity depends on the hierarchical structure of Γ rather than the geometric distance in ℝ ^d . For a class of fractals in ℝ² we find a condition on the measure j so that the Dirichlet form ɛ is regular. The condition is given in terms of Hausdorff dimension of Γ.     

Closed-form solutions for a mode III radial matrix crack penetrating a circular inhomogeneity

In this research we address in detail a mode III radial matrix crack penetrating a circular inhomogeneity. One tip of the radial crack lies in the matrix, while the other tip of the radial crack lies in the circular inhomogeneity. In addition the two tips of the crack are mutually image points (or inverse points) with respect to the circular inhomogeneity-matrix interface. First we conformally map the crack onto a unit circle C_a in the new ζ-plane. Meanwhile the inhomogeneity-matrix interface is mapped onto C_b, a part of another circle in the ζ-plane. In addition C_a and C_b intersect at a vertex angle π/2. By using the method of image in the ζ-plane, closed-form solutions in terms of elementary functions are derived for three loading cases: (1) remote uniform antiplane shearing; (2) a screw dislocation located in the unbounded matrix; and (3) a radial Zener–Stroh crack.     

How fast are the particles of super-Brownian motion?

In this paper we investigate fast particles in the range and support ofsuper-Brownian motion in the historical setting. In this setting eachparticle of super-Brownian motion alive at time t is represented by apath w:[0,t]→ℝ ^d and the state of historical super-Brownian motionis a measure on the set of paths. Typical particles have Brownian paths,however in the uncountable collection of particles in the range of asuper-Brownian motion there are some which at exceptional times movefaster than Brownian motion. We determine the maximal speed of allparticles during a given time period E, which turns out to be afunction of the packing dimension of E. A path w in the support ofhistorical super-Brownian motion at time t is called a-fast if . Wecalculate the Hausdorff dimension of the set of a-fast paths in thesupport and the range of historical super-Brownian motion. A valuabletool in the proofs is a uniform dimension formula for the Browniansnake, which reduces dimension problems in the space of stopped paths to dimension problems on the line. Received: 27 January 2000 / Revised version: 28 August 2000 / Published online: 24 July 2001     

Covering a set with homothets of a convex body

We consider two problems mentioned in the book “Research Problems in Discrete Geometry” (Brass et al. in research problems in discrete geometry, vol xii+499. Springer, New York, pp ISBN: 978-0387-23815-8; 0-387-23815-8, 2005). First, let K and L be given convex bodies in \mathbbR^d{\mathbb{R}^{d}} . We prove that if the total volume of a family of positive homothets of K is sufficiently large then they permit a translative covering of L. This problem, in the case when K = L and the dimension is two, was originally posed by L. Fejes Tóth. The previously known bound (Januszewski in proc. of the International scientific conference on mathematics, pp 29–34. Žilina, 1998) on the total volume (in the case when K = L) was of order d ^d vol(K), we prove a bound that is exponential in the dimension. The second problem is the following: Find a condition, in terms of the coefficients of homothety, that is necessary for a family of positive homothets of K to cover K. The problem was phrased by V. Soltan, who conjectured that the sum of the coefficients is at least d. We confirm an asymptotic version of this conjecture.     

Singular Probability Distributions and Fractal Properties of Sets of Real Numbers Defined by the Asymptotic Frequencies of Their <Emphasis Type="Italic">s</Emphasis>-Adic Digits

Properties of the set T _s of “particularly nonnormal numbers” of the unit interval are studied in detail (T _s consists of real numbers x some of whose s-adic digits have the asymptotic frequencies in the nonterminating s-adic expansion of x, and some do not). It is proved that the set T _s is residual in the topological sense (i.e., it is of the first Baire category) and is generic in the sense of fractal geometry (T _s is a superfractal set, i.e., its Hausdorff-Besicovitch dimension is equal to 1). A topological and fractal classification of sets of real numbers via analysis of asymptotic frequencies of digits in their s-adic expansions is presented. Dedicated to V. S. Korolyuk on occasion of his 80th birthday __________ Published in Ukrains'kyi Matematychnyi Zhurnal, Vol. 57, No. 9, pp. 1163–1170, September, 2005.     

Stability of Riemann solutions with large oscillation for the relativistic Euler equations

We are concerned with entropy solutions of the 2×2 relativistic Euler equations for perfect fluids in special relativity. We establish the uniqueness of Riemann solutions in the class of entropy solutions in L^∞∩BV_loc with arbitrarily large oscillation. Our proof for solutions with large oscillation is based on a detailed analysis of global behavior of shock curves in the phase space and on special features of centered rarefaction waves in the physical plane for this system. The uniqueness result does not require specific reference to any particular method for constructing the entropy solutions. Then the uniqueness of Riemann solutions yields their inviscid large-time stability under arbitrarily largeL¹∩L^∞∩BV_loc perturbation of the Riemann initial data, as long as the corresponding solutions are in L^∞ and have local bounded total variation that allows the linear growth in time. We also extend our approach to deal with the uniqueness and stability of Riemann solutions containing vacuum in the class of entropy solutions in L^∞ with arbitrarily large oscillation.     

Feynman’s Operational Calculi: Disentangling Away from the Origin

In this paper we develop a method of forming functions of noncommuting operators (or disentangling) using functions that are not necessarily analytic at the origin in ℂ ⁿ . The method of disentangling follows Feynman’s heuristic rules from in (Feynman in Phys. Rev. 84:18–128, 1951) a mathematically rigorous fashion, generalizing the work of Jefferies and Johnson and the present author in (Jefferies and Johnson in Russ. J. Math. 8:153–181, 2001) and (Jefferies et al. in J. Korean Math. Soc. 38:193–226, 2001). In fact, the work in (Jefferies and Johnson in Russ. J. Math. 8:153–181, 2001) and (Jefferies et al. in J. Korean Math. Soc. 38:193–226, 2001) allow only functions analytic in a polydisk centered at the origin in ℂ ⁿ while the method introduced in this paper enable functions that are not analytic at the origin to be used. It is shown that the disentangling formalism introduced here reduces to that of (Jefferies and Johnson in Russ. J. Math. 8:153–181, 2001) and (Jefferies et al. in J. Korean Math. Soc. 38:193–226, 2001) under the appropriate assumptions. A basic commutativity theorem is also established.     

Time dependent neutron transport theory in multigroup-P L-approximation

Summary A method is developed for solving the time dependent neutron transport equation in multigroupP _L approximation for one-dimensional geometries. The partial differential equations in time and space are solved by means of a power series expansion in the spatial variable. The resulting ordinary differential equations are solved up to theNth order, and the last spatial coefficients are used to satisfy the boundary conditions. Integration of the purely time dependent differential equations is carried out by means of Lie series.Numerical oscillations, appearing for high ordersN, are avoided by subdividing each zone into smaller subzones. Even and odd spatial moments must be developed in opposite directions in each subzone, and the stationary solution representing the initial condition for the time dependent calculation must be developed in the same manner.Results of two calculations in spherical geometry are presented. One is the start-up of a small experimental reactor usingP ₃ theory, the other is a demonstration of neutron waves inP ₁ theory.

---

Zusammenfassung Es wird eine Methode zur Lösung der zeitabhängigen Neutronentransportgleichung in der Multigruppen-P _L-Näherung für eindimensionale Geometrien entwickelt. Für die Lösung der partiellen Differentialgleichungen in Ort und Zeit wird eine Potenzreihe im Ort angesetzt. Das sich ergebende gewöhnliche Differentialgleichungssystem wird bis zu einer gewählten OrdnungN erfüllt, und die letzten Koeffizienten der Ortsentwicklungen dienen zur Befriedigung der Randbedingungen. Die Integration in der Zeit erfolgt dann mit der Lie-Reihen-Methode.Zur Vermeidung der Oszillationen, die bei hohen OrdnungenN auftreten, werden die einzelnen Zonen in kleinere Abschnitte unterteilt. Die geraden und ungeraden Momente müssen in den Abschnitten in entgegengesetzter Richtung entwickelt werden. Die stationäre Lösung, die als Anfangsbedingung für die zeitabhängige Rechnung dient, muss mit demselben Entwicklungsschema ermittelt werden.Die Methode wird angewendet zur Lösung derP ₁- undP ₃-Gleichungen in sphärischen Reaktoren. Zwei Beispiele werden damit berechnet: Das Anfahren eines kleinen Versuchsreaktors inP ₃-Näherung und die Nachweisung von Neutronenwellen inP ₁-Theorie.

     

Efficient Algorithms for Finding the Maximum Number of Disjoint Paths in Grids

In a rectangular grid, given two sets of nodes, (sources) and (sinks), of size each, the disjoint paths (DP) problem is to connect as many nodes in to the nodes in using a set of “disjoint” paths. (Both edge-disjoint and vertex-disjoint cases are considered in this paper.) Note that in this DP problem, a node in can be connected to any node in . Although in general the sizes of and do not have to be the same, algorithms presented in this paper can also find the maximum number of disjoint paths pairing nodes in and . We use the network flow approach to solve this DP problem. By exploiting all the properties of the network, such as planarity and regularity of a grid, integral flow, and unit capacity source/sink/flow, we can optimally compress the size of the working grid (to be defined) from O(N²) to O(N^1.5) and solve the problem in O(N^2.5) time for both the edge-disjoint and vertex-disjoint cases, an improvement over the straightforward approach which takes O(N³) time.     

A note on the C2/G/1 queue and the C2/G/1 loss system

In this note, we consider the steady-state probability of delay (PW) in the C₂/G/1 queue and the steady-state probability of loss (p_loss) in the C₂/G/1 loss system, in both of which the interarrival time has a two-phase Coxian distribution. We show that, for c_X ²<1, where c_X is the coefficient of variation of the interarrival time, both p_loss and PW are increasing in β(s), the Laplace–Stieltjes transform of the general service-time distribution. This generalises earlier results for the GE₂/G/1 queue and the GE₂/G/1 loss system. The practical significance of this is that, for c_X ²<1, p_loss in the C₂/G/1 loss system and PW in the C₂/G/1 queue are both increasing in the variability of the service time. This revised version was published online in June 2006 with corrections to the Cover Date.     

A Problem in Linear Matrix Approximation

Let A be a normal operator in ??(H), H a complex Hilbert space, and let ? _A = ? {AX - XA:X ∈ ??(H)} be the commutator subspace of ??(H) associated with A. If B in ??(H) commutes with A, then B is orthogonal to ?_A with respect to the spectral norm; i.e., the null operator is an element of best approximation of B in ? _A. This was proved by J. Anderson in 1973 and extended by P. J. Maher with respect to the Schatten p-norm recently. We take a look at their result from a more approximation theoretical point of view in the finite dimensional setting; in particular, we characterize all elements of best approximation of B in R_A and prove that the metric projection of H onto ?_A is continuous.     

Pullback attractors for a non-autonomous incompressible non-Newtonian fluid

This paper studies the pullback asymptotic behavior of solutions for a non-autonomous incompressible non-Newtonian fluid in two-dimensional (2D) bounded domains. We first prove the existence of pullback attractors AV in space V (has H²-regularity, see notation in Section 2) and AH in space H (has L²-regularity) for the cocycle corresponding to the solutions of the fluid. Then we verify the regularity of the pullback attractors by showing AV=AH, which implies the pullback asymptotic smoothing effect of the fluid in the sense that the solutions become eventually more regular than the initial data.     

Homogenization of the parabolic cauchy problem in the Sobolev class H<Superscript>1</Superscript>(R<Superscript>d</Superscript>)

Homogenization in the small period limit for the solution ue of the Cauchy problem for a parabolic equation in R^d is studied. The coefficients are assumed to be periodic in R^d with respect to the lattice ɛG. As ɛ → 0, the solution u ɛ converges in L2(R^d) to the solution u0 of the effective problem with constant coefficients. The solution u ɛis approximated in the norm of the Sobolev space H ¹(R^d) with error O( ɛ); this approximation is uniform with respect to the L2-norm of the initial data and contains a corrector term of order ɛ. The dependence of the constant in the error estimate on time t is given. Also, an approximation in H ¹(R^d) for the solution of the Cauchy problem for a nonhomogeneous parabolic equation is obtained.     

Random Flights in Higher Spaces

We consider in this paper random flights in ℝ ^d performed by a particle changing direction of motion at Poisson times. Directions are uniformly distributed on hyperspheres S ₁ ^d . We obtain the conditional characteristic function of the position of the particle after n changes of direction. From this characteristic function we extract the conditional distributions in terms of (n+1)−fold integrals of products of Bessel functions. These integrals can be worked out in simple terms for spaces of dimension d=2 and d=4. In these two cases also the unconditional distribution is determined in explicit form. Some distributions connected with random flights in ℝ³ are discussed and in some special cases are analyzed in full detail. We point out that a strict connection between these types of motions with infinite directions and the equation of damped waves holds only for d=2. Related motions with random velocity in spaces of lower dimension are analyzed and their distributions derived.     

On the tensor rank of the multiplication in the finite fields

First, we prove the existence of certain types of non-special divisors of degree g−1 in the algebraic function fields of genus g defined over Fq. Then, it enables us to obtain upper bounds of the tensor rank of the multiplication in any extension of quadratic finite fields Fq by using Shimura and modular curves defined over Fq. From the preceding results, we obtain upper bounds of the tensor rank of the multiplication in any extension of certain non-quadratic finite fields Fq, notably in the case of F₂. These upper bounds attain the best asymptotic upper bounds of Shparlinski-Tsfasman-Vladut [I.E. Shparlinski, M.A. Tsfasman, S.G. Vladut, Curves with many points and multiplication in finite fields, in: Lecture Notes in Math., vol. 1518, Springer-Verlag, Berlin, 1992, pp. 145-169].     

On the Fiedler vectors of graphs that arise from trees by Schur complementation of the Laplacian

The utility of Fiedler vectors in interrogating the structure of graphs has generated intense interest and motivated the pursuit of further theoretical results. This paper focuses on how the Fiedler vectors of one graph reveal structure in a second graph that is related to the first. Specifically, we consider a point of articulation r in the graph G whose Laplacian matrix is L and derive a related graph G_{r} whose Laplacian is the matrix obtained by taking the Schur complement with respect to r in L. We show how Fiedler vectors of G_{r} relate to the structure of G and we provide bounds for the algebraic connectivity of G_{r} in terms of the connected components at r in G. In the case where G is a tree with points of articulation rR, we further consider the graph G_R derived from G by taking the Schur complement with respect to R in L. We show that Fiedler vectors of G_R valuate the pendent vertices of G in a manner consistent with the structure of the tree.     

On the modelling of distributed outflow in one-dimensional models of arterial blood flow

Summary One-dimensional mathematical models of the flow of blood in arteries commonly make use of the concept of distributed outflow to simulate the loss of blood via side branches. It is shown that the usual approach leads to physically unrealistic results in certain special cases involving flow in rigid tubes, unless we introduce in the momentum equation a suitable additional term depending on the outflow. The influence of such a term on the pressure- and velocity waves calculated from an existing model of the aorta is investigated.

---

Zusammenfassung Eindimensionale mathematische Modelle des Blutstromes in Arterien verwenden meist das Konzept von distribuiertem Abfluß um den Blutverlust über Verzweigungen zu simulieren. Es wird gezeigt, daß die übliche Behandlungsweise in speziellen Fällen von Strömungen in starren Röhren zu physikalisch unrealistischen Resultaten führt, es sei denn, daß in der Impulsgleichung ein passendes zusätzliches Glied eingeführt wird, welches vom Abfluß abhängt. Der Einfluß dieses Gliedes auf die Druck- und Geschwindigkeitswellen wird mit Hilfe eines bestehenden Modelles der Aorta untersucht.

     

The lifetimes of positrons in oxides of arsenic and antimony

The lifetimes of positrons in crystalline powders of dimeric oxides of arsenic and antimony, and in trioxymethylene and polyoxymethylene are reported. In As₄O₆ and polyoxymethylene a longer componentτ ₄ is observed in addition to the usually observedτ ₂ component. The variation ofτ ₄ with temperature and pressure is studied in As₄O₆. On the basis of these results and the X-ray diffraction pictures, it is suggested that ortho-positronium atoms quenched in intercrystallite regions in As₄O₆ gives rise toτ ₄.     

Analysis of a typical railway turnout sleeper system using grillage beam analogy

A simplified grillage beam analogy was performed to investigate the behaviour of railway turnout sleeper system with a low value of elastic modulus on different support moduli. This study aimed at determining an optimum modulus of elasticity for an emerging technology for railway turnout application - fibre composites sleeper. The finite element simulation suggests that the changes in modulus of elasticity of sleeper, E_sleeper and the sleeper support modulus, U_s have a significant influence on the behaviour of turnout sleepers. The increase in U_s from 10 to 40 MPa resulted in a 15% reduction in the bending moment while the increase in E_sleeper from 1 GPa to 10 GPa has resulted in almost 75% increase in the bending moment. The shear forces in turnout sleepers is not sensitive to both the changes of the E_sleeper and U_s while the sleeper with low E_sleeper tend to undergo greater settlement into the ballast. An E_sleeper of 4 GPa was found optimal for an alternative fibre composite turnout sleeper provided that the U_s is at least 20 MPa from the consideration of sleeper ballast pressure and maximum vertical deflection. It was established that the turnout sleeper has a maximum bending moment of 19 kN-m and a shear force of 158 kN under service conditions.