The Lie-group shooting method for solving nonlinear singularly perturbed boundary value problems

A new computational method for solving the second-order nonlinear singularly perturbed boundary value problems (SPBVPs) is provided in this paper. In order to overcome a highly singular behavior very near to the boundary as being not easy to treat by numerical method, we adopt a coordinate transformation from an x-domain to a t-domain via a rescaling technique, which can reduce the singularity within the boundary layer. Then, we construct a Lie-group shooting method (LGSM) to search a missing initial condition through the finding of a suitable value of a parameter r ∈ [0, 1]. Moreover, we can derive a closed-form formula to express the initial condition in terms of r, which can be determined properly by an accurate matching to the right-boundary condition. Numerical examples are examined, showing that the present approach is highly efficient and accurate.     

Free convection heat and mass transfer in a power law fluid past an inclined surface with thermophoresis

The problem of heat and mass transfer in a power law, two-dimensional, laminar, boundary layer flow of a viscous incompressible fluid over an inclined plate with heat generation and thermophoresis is investigated by the characteristic function method. The governing non-linear partial differential equations describing the flow and heat transfer problem are transformed into a set of coupled non-linear ordinary differential equation which was solved using Runge–Kutta shooting method. Exact solutions for the dimensionless temperature and concentration profiles, are presented graphically for different physical parameters and for the different power law exponents 0 < n < 0.5 and for n > 0.5.     

Solution of the Burgers equation using an implicit linearizing transformation

The initial-boundary-value problem on the semi-infinite interval and on a finite interval for the Burgers equation u_t = u_xx + 2u_xu is solved using a stream function ? and a linearizing transformation w = e^?. The transformation reduces the equation to a heat equation with appropriate initial and homogeneous time-dependent linear boundary conditions. One advantage of this method is that we never need to find an explicit expression for ? in our computations.     

Exact solutions for the rotational flow of a generalized Maxwell fluid between two circular cylinders

Unsteady flow of an incompressible generalized Maxwell fluid between two coaxial circular cylinders is studied by means of the Laplace and finite Hankel transforms. The motion of the fluid is produced by the rotation of cylinders around their common axis. The solutions that have been obtained, written in integral and series form in terms of the generalized G_a,b,c(·, t)-functions, are presented as a sum of the Newtonian solutions and the corresponding non-Newtonian contributions. They satisfy all imposed initial and boundary conditions and for λ → 0 reduce to the solutions corresponding to the Newtonian fluids performing the same solution. Furthermore, the corresponding solutions for ordinary Maxwell fluids are also obtained for β = 1. Finally, in order to reveal some relevant physical aspects of the obtained results, the diagrams of the velocity field ω(r, t) have been depicted against r and t for different values of the material and fractional parameters.     

Mathematical modeling and performance analysis of combed yarn production system: Based on few data

The paper describes the availability of combed sliver production system, a part of yarn production plant. The units under study are specialized single purpose machines. Performance analysis of the system is carried out to identify the key factors. The optimum value of ‘r’, where ‘r’ represent the number of repairman to repair the twelve carding machines (r ? 12), is calculated to maximizing the steady state availability of the system. The problem is formulated using probability consideration and supplementary variable technique. Probability considerations at various stages give differential-difference equations, which are solved using Lagrange method to obtain the state probabilities. The numerical analysis carried out helps in increasing the production rate by controlling the factors affecting the system i.e. availability optimization.     

Calculation of a static potential created by plane fractal cluster

In this paper we demonstrate new approach that can help in calculation of electrostatic potential of a fractal (self-similar) cluster that is created by a system of charged particles. For this purpose we used the simplified model of a plane dendrite cluster [1] that is generated by a system of the concentric charged rings located in some horizontal plane (see Fig. 2). The radiuses and charges of the system of concentric rings satisfy correspondingly to relationships: r_n = r₀ξⁿ and e_n = e₀bⁿ, where n determines the number of a current ring. The self-similar structure of the system considered allows to reduce the problem to consideration of the functional equation that similar to the conventional scaling equation. Its solution represents itself the sum of power-low terms of integer order and non-integer power-law term multiplied to a log-periodic function [5], [6]. The appearance of this term was confirmed numerically for internal region of the self-similar cluster (r₀ ≪ r ≪ r_N−1), where r₀, r_N−1 determine the smallest and the largest radiuses of the limiting rings correspondingly. The results were obtained for homogeneously (b > 0) and heterogeneously (b < 0) charged rings. We expect that this approach allows to consider more complex self-similar structures with different geometries of charge distributions.     

Effect of variable density on hydromagnetic mixed convection flow of a non-Newtonian fluid past a moving vertical plate

The effect of temperature-dependent density on MHD mixed convection flow of power-law fluid past a moving semi-infinite vertical plate for high temperature differences between the plate and the ambient fluid is studied. The fluid density is assumed to decrease exponentially with temperature. The usual Boussinesq approximations are not considered due to the large temperature differences. The surface temperature of the moving plate was assumed to vary according to a power-law form, that is, T_w(x) = T_∞ + Ax^γ. The fluid is permeated by a uniform magnetic field imposed perpendicularly to the plate on the assumption of small magnetic Reynolds number. A numerical shooting algorithm for two unknown initial conditions with fourth-order Runge–Kutta integration scheme has been used to solve the coupled non-linear boundary value problem. The effects of various parameters on the velocity and temperature profiles as well as the local skin-friction coefficient and the local Nusselt number are presented graphically and in the tabular form. The results show that application of Boussinesq approximations in a non-Newtonian fluid subjected to high temperature differences gives a significant error in the values of the skin-friction coefficient and the application of an external magnetic field reduces this error markedly in the case of shear-thickening fluid.     

Special Functions on the Sphere with Applications to Minimal Surfaces

A function which is homogeneous in x, y, z of degree n and satisfies V_xx + V_yy + V_zz = 0 is called a spherical harmonic. In polar coordinates, the spherical harmonics take the form rⁿf_n, where f_n is a spherical surface harmonic of degree n. On a sphere, f_n satisfies ▵ f_n + n(n + 1)f_n = 0, where ▵ is the spherical Laplacian. Bounded spherical surface harmonics are well studied, but in certain instances, unbounded spherical surface harmonics may be of interest. For example, if X is a parameterization of a minimal surface and n is the corresponding unit normal, it is known that the support function, w = X · n, satisfies ▵w + 2w = 0 on a branched covering of a sphere with some points removed. While simple in form, the boundary value problem for the support function has a very rich solution set. We illustrate this by using spherical harmonics of degree one to construct a number of classical genus-zero minimal surfaces such as the catenoid, the helicoid, Enneper's surface, and Hennenberg's surface, and Riemann's family of singly periodic genus-one minimal surfaces.     

Eigenvalue problems for fractional ordinary differential equations

The eigenvalue problems are considered for the fractional ordinary differential equations with different classes of boundary conditions including the Dirichlet, Neumann, Robin boundary conditions and the periodic boundary condition. The eigenvalues and eigenfunctions are characterized in terms of the Mittag–Leffler functions. The eigenvalues of several specified boundary value problems are calculated by using MATLAB subroutine for the Mittag–Leffler functions. When the order is taken as the value 2, our results degenerate to the classical ones of the second-ordered differential equations. When the order α satisfies 1 < α < 2 the eigenvalues can be finitely many.     

Radical Rings with Soluble Adjoint Groups

An associative ring R, not necessarily with an identity, is called radical if it coincides with its Jacobson radical, which means that the set of all elements of R forms a group denoted by R^∘ under the circle operation r ∘ s = r + s + rs on R. It is proved that every radical ring R whose adjoint group R^∘ is soluble must be Lie-soluble. Moreover, if the commutator factor group of R^∘ has finite torsion-free rank, then R is locally nilpotent.     

Dynamics of landslide model with time delay and periodic parameter perturbations

In present paper, we analyze the dynamics of a single-block model on an inclined slope with Dieterich–Ruina friction law under the variation of two new introduced parameters: time delay T_d and initial shear stress μ. It is assumed that this phenomenological model qualitatively simulates the motion along the infinite creeping slope. The introduction of time delay is proposed to mimic the memory effect of the sliding surface and it is generally considered as a function of history of sliding. On the other hand, periodic perturbation of initial shear stress emulates external triggering effect of long-distant earthquakes or some non-natural vibration source. The effects of variation of a single observed parameter, T_d or μ, as well as their co-action, are estimated for three different sliding regimes: β < 1, β = 1 and β > 1, where β stands for the ratio of long-term to short-term stress changes. The results of standard local bifurcation analysis indicate the onset of complex dynamics for very low values of time delay. On the other side, numerical approach confirms an additional complexity that was not observed by local analysis, due to the possible effect of global bifurcations. The most complex dynamics is detected for β < 1, with a complete Ruelle–Takens–Newhouse route to chaos under the variation of T_d, or the co-action of both parameters T_d and μ. These results correspond well with the previous experimental observations on clay and siltstone with low clay fraction. In the same regime, the perturbation of only a single parameter, μ, renders the oscillatory motion of the block. Within the velocity-independent regime, β = 1, the inclusion and variation of T_d generates a transition to equilibrium state, whereas the small oscillations of μ induce oscillatory motion with decreasing amplitude. The co-action of both parameters, in the same regime, causes the decrease of block’s velocity. As for β > 1, highly-frequent, limit-amplitude oscillations of initial stress give rise to oscillatory motion. Also for β > 1, in case of perturbing only the initial shear stress, with smaller amplitude, velocity of the block changes exponentially fast. If the time delay is introduced, besides the stress perturbation, within the same regime, the co-action of T_d (T_d < 0.1) and small oscillations of μ induce the onset of deterministic chaos.     

Matrix characterization of linear codes with arbitrary Hamming weight hierarchy

The support of an [n, k] linear code C over a finite field F_q is the set of all coordinate positions such that at least one codeword has a nonzero entry in each of these coordinate position. The rth generalized Hamming weight d_r(C), 1 ⩽ r ⩽ k, of C is defined as the minimum of the cardinalities of the supports of all [n, r] subcodes of C. The sequence (d₁(C), d₂(C), … , d_k(C)) is called the Hamming weight hierarchy (HWH) of C. The HWH, d_r(C) = n − k + r;  r = 1, 2 , …, k, characterizes maximum distance separable (MDS) codes. Therefore the matrix characterization of MDS codes is also the characterization of codes with the HWH d_r(C) = n − k + r; r = 1, 2, … , k. A linear code C with systematic check matrix [I∣P], where I is the (n − k) × (n − k) identity matrix and P is a (n − k) × k matrix, is MDS iff every square submatrix of P is nonsingular. In this paper we extend this characterization to linear codes with arbitrary HWH. Using this result, we characterize Near-MDS codes, Near-Near-MDS (N²-MDS) codes and A^μ-MDS codes. The MDS-rank of C is the smallest integer η such that d_η+1 = n − k + η + 1 and the defect vector of C with MDS-rank η is defined as the ordered set {μ₁(C), μ₂(C), μ₃(C), … , μ_η(C), μ_η+1(C)}, where μ_i(C) = n − k + i − d_i(C). We call C a dually defective code if the defect vector of the code and its dual are the same. We also discuss matrix characterization of dually defective codes. Further, the codes meeting the generalized Greismer bound are characterized in terms of their generator matrix. The HWH of dually defective codes meeting the generalized Greismer bound are also reported.     

Statistical properties of nucleotides in human chromosomes 21 and 22

In this paper the statistical properties of nucleotides in human chromosomes 21 and 22 are investigated. The n-tuple Zipf analysis with n = 3, 4, 5, 6, and 7 is used in our investigation. It is found that the most common n-tuples are those which consist only of adenine (A) and thymine (T), and the rarest n-tuples are those in which GC or CG pattern appears twice. With the n-tuples become more and more frequent, the double GC or CG pattern becomes a single GC or CG pattern. The percentage of four nucleotides in the rarest ten and the most common ten n-tuples are also considered in human chromosomes 21 and 22, and different behaviors are found in the percentage of four nucleotides. Frequency of appearance of n-tuple f(r) as a function of rank r is also examined. We find the n-tuple Zipf plot shows a power-law behavior for r < 4ⁿ⁻¹ and a rapid decrease for r > 4ⁿ⁻¹. In order to explore the interior statistical properties of human chromosomes 21 and 22 in detail, we divide the chromosome sequence into some moving windows and we discuss the percentage of ξη (ξ, η = A, C, G, T) pair in those moving windows. In some particular regions, there are some obvious changes in the percentage of ξη pair, and there maybe exist functional differences. The normalized number of repeats N₀(l) can be described by a power law: N₀(l) ∼ l^−μ. The distance distributions P₀(S) between two nucleotides in human chromosomes 21 and 22 are also discussed. A two-order polynomial fit exists in those distance distributions: log P₀(S) = a + bS + cS², and it is quite different from the random sequence.     

Multifractal scaling of crack images from pyroligneous acid dried sludge

Tannery effluent (sludge, wastewater) is treated by natural way. The waste sludge has been taken for two treatment process. The alkali chemicals are neutralized by pyroligneous acid which is obtained by pyrolysis process of wood. This sludge is sent out for drying. The dried sludge contains some crack pattern formation. Photographs were used to record two sludge cracking surfaces. Experiment has been performed to study the fracture pattern formation in thin film sludge. We studied changes of crack surface of a sludge by image analysis and also assessed applicability of fractal scaling to crack surfaces. The calculated crack surface dimension shows that the fracture surface exhibit fractal structure. Image size was 256 × 256 pixels. MFA (multifractal analysis) was carried out to the method of moments, i.e., the probability distribution was estimated for moments ranging from ?150 < q < 150 and the generalized dimension were calculated from the log/log slope of the probability distribution for the respective moments over box sizes. Generalized dimension D(q) were attained for this box size range, which are capable of characterizing heterogeneous spatial crack structure. Multifractal spectra analyzed two fracture surface of the image and results were indicated that the width of spectra increases due to pyroligneous acid. Multifractal method is sensitive enough to measure the fracture distribution and can be seen as a different approach for changing research of crack images of manure sludge drying.     

Homotopy analysis of unsteady boundary layer flow adjacent to permeable stretching surface in a porous medium

This communication deals with the unsteady boundary layer flow of a viscous fluid in porous medium started due to the impulsively stretching of the plane wall. The wall is assumed to be porous so that suction or injection is possible. Complete analytic solution which is uniformly valid for all the dimensionless times 0 ⩽ τ < 0 in the whole spatial region 0 ⩽ η < ∞ is obtained by a purely analytic technique, namely the homotopy analysis method. Results are discussed through graphs.     

On the stagnation-point flow towards a stretching sheet with homogeneous–heterogeneous reactions effects

The effects of homogeneous–heterogeneous reactions on the steady boundary layer flow near the stagnation point on a stretching surface is studied. The possible steady-states of this system are analyzed in the case when the diffusion coefficients of both reactant and auto catalyst are equal. The strength of this effect is represented by the dimensionless parameter K and K_s. It is shown that for a fluid of small kinematic viscosity, a boundary layer is formed when the stretching velocity is less than the free stream velocity and an inverted boundary layer is formed when the stretching velocity exceeds the free stream velocity. The uniqueness of this problem lies on the fact that the solutions are possible for all values of λ > 0 (stretching surface), while for λ < 0 (shrinking surface), solutions are possible only for its limited range.     

Inverse problems for parabolic equations 2

Let u_t − u_xx = h(t) in 0 ⩽ x ⩽ π, t ⩾ 0. Assume that u(0, t) = v(t), u(π, t) = 0, and u(x, 0) = g(t). The problem is: what extra data determine the three unknown functions {h, v, g} uniquely? This question is answered and an analytical method for recovery of the above three functions is proposed.     

Electrostatics in fractal geometry: Fractional calculus approach

An electric field in a composite dielectric with a fractal charge distribution is obtained in the spherical symmetry case. The method is based on the splitting of a composite volume into a fractal volume V_d ～ r^d with the fractal dimension d and a complementary host volume V_h = V₃ ? V_d. Integrations over these fractal volumes correspond to the convolution integrals that eventually lead to the employment of the fractional integro-differentiation.     

Conceptual blending: Student reasoning when proving “conditional implies conditional” statements

Conceptual blending describes how humans condense information, combining it in novel ways. The blending process may create global insight or new detailed connections, but it may also result in a loss of information, causing confusion. In this paper, we describe the proof writing process of a group of four students in a university geometry course proving a statement of the form conditional implies conditional, i.e., (p → q) ⇒ (r → s). We use blending theory to provide insight into three diverse questions relevant for proof writing: (1) Where do key ideas for proofs come from?, (2) How do students structure their proofs and combine those structures with their more intuitive ideas?, and (3) How are students reasoning when they fail to keep track of the implication structure of the statements that they are using? We also use blending theory to describe the evolution of the students’ proof writing process through four episodes each described by a primary blend.     

An extended rational interpolation method

To interpolate function, f(x), a ? x ? b, when we have some information about the values of f(x) and their derivatives in separate points on {x₀, x₁, … , x_n} ? [a, b], the Hermit interpolation method is usually used. Here, to solve this kind of problems, extended rational interpolation method is presented and it is shown that the suggested method is more efficient and suitable than the Hermit interpolation method, especially when the function f(x) has singular points in interval [a, b]. Also for implementing the extended rational interpolation method, the direct method and the inverse differences method are presented, and with some examples these arguments are examined numerically.