Self-Similarity and the Evaluation of Long-Time Nonhydrostatic Effects in Compositionally Driven Gravity Flows in Deep Surroundings

In the study of compositionally driven gravity currents involving one or more homogeneous fluid layers, it has been customary to adopt the hydrostatic assumption for the pressure field in each layer which, in turn, leads to a depth‐independent horizontal velocity field in each of these layers and significant simplifications to the governing equations. Under this hydrostatic paradigm, each layer will then have its motion governed by the well‐known reduced dimension shallow‐water equations. For the so‐called ‐layer or reduced gravity shallow‐water equations, similarity solutions for fixed volume gravity currents released in rectangular geometry have been found. Very few attempts have been made to evaluate contributions arising from the possible loss of hydrostatic balance in the context of the problems treated using the classic shallow‐water approach. Where such attempts have been pursued, they have usually been carried out in a time‐independent context or using layer‐averaged equations and very small amplitude disturbances. The vast majority of these studies into nonhydrostatic effects do not include any relevant numerical work to assess these effects. In this paper, we develop an approach for evaluating nonhydrostatic contributions to the flow field for bottom gravity currents in deep surroundings and rectangular geometry. Our approach makes no assumptions on the amplitudes of the disturbances and does not depend on layer‐averaging in the governing equations. We seek asymptotic expansions of the solutions to the Euler equations for a shallow fluid by using the small parameter δ², where δ is the aspect ratio of the flow regime. At leading order the equations enforce hydrostatic balance while those obtained at first order retain certain nonhydrostatic effects which we evaluate. Our method for evaluation of these first‐order contributions employs the self‐similar nature of the solution to the leading‐order equations in the new first‐order equations without any vertical averaging procedures being employed.     

Analytical solution for the space fractional diffusion equation by two-step Adomian Decomposition Method

Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are increasingly used in modeling practical superdiffusive problems in fluid flow, finance and other areas of application. This paper presents the analytical solutions of the space fractional diffusion equations by two-step Adomian Decomposition Method (TSADM). By using initial conditions, the explicit solutions of the equations have been presented in the closed form and then their solutions have been represented graphically. Two examples, the first one is one-dimensional and the second one is two-dimensional fractional diffusion equation, are presented to show the application of the present technique. The solutions obtained by the standard decomposition method have been numerically evaluated and presented in the form of tables and then compared with those obtained by TSADM. The present TSADM performs extremely well in terms of efficiency and simplicity.     

Implementación eficiente de una formulación semirrecursiva para la dinámica de sistemas multicuerpo de gran tamaño

This article shows an efficient implementation of a dynamic semi-recursive formulation for large and complex multibody system simulations, with interesting applications in the automotive field and especially with industrial vehicles. These systems tend to have a huge amount of kinematic constraints, becoming usual the presence of redundant but compatible systems of equations. The maths involved in the solution of these problems have a high computational cost, making very challenging to achieve real-time simulations.In this article, two implementations to increase the efficiency of these computations will be shown. The difference between them is the way they consider the Jacobian matrix of the constraint equations. The first one treats this matrix as a dense one, using the BLAS functions to solve the system of equations. The second one takes into account the sparse pattern of the Jacobian matrix, introducing the sparse function MA48 from Harwell.Both methodologies have been applied on two multibody system models with different sizes. The first model is a vehicle IVECO DAILY 35C15 with 17 degrees of freedom. The second one is a semi-trailer truck with 40 degrees of freedom. Taking as a reference the standard C/C + + implementation, the efficiency improvements that have been achieved using dense matrices (BLAS) have been of 15% and 50% respectively. The results in the first model have not improved significantly by using sparse matrices, but in the second one, the times with sparse matrices have been reduced 8% with respect to the BLAS ones.     

Homoclinic Bifurcations for Partial Differential Equations in Unbounded Domains

The connection between low-dimensional chaos in ordinary differential equations, and turbulence in fluids and other systems governed by partial differential equations, is one that is in many circumstances not clear. We discuss some examples of turbulent fluid flow, and consider ways in which they may be related to much simpler sets of ordinary differential equations, whose behavior can be reasonably well understood. (We are not advocating drastic Fourier truncation.) The generation of aperiodic solutions through the occurrence of homoclinic orbits is briefly analysed for ordinary differential equations, and the same kind of heuristic analysis is sketched for partial differential equations (in one space dimension). It is suggested that such an analysis can explain certain features of chaos, which have been observed in real fluids.     

Generalized Electrodynamics With Ternary Internal Structure

In Refs. [2]–[7] we suggested generalized dynamic equations of motion of relativistic charged particles inside electromagnetic fields. The dynamic equations had been formulated in terms of external as well as internal momenta. Evolution equations for external momenta, the Lorentz-force equations, had been derived from evolution equations for internal momenta. In this paper, along with relativistic dynamics we generalize electromagnetic fields within the scope of ternary algebras. The full theory is constructed in 4D euclidean space. This space possesses an advantage to build ternary mappings from three vectors onto one. The dynamics is given by non-linear evolution equations with cubic characteristic polynomial. In polar representation the internal momenta obey the Jacobi equations whereas external momenta obey the Weierstrass equations for elliptic functions. The generalized electromagnetic fields are defined by the triple fields where the first one has properties of the electric field and the other two have properties of the magnetic field. The field equations for the triple fields analogous to the Maxwell equations are suggested.     

Group theoretic method for unsteady free convection flow of a micropolar fluid along a vertical plate in a thermally stratified medium

The group theoretic approach is applied for solving the problem of unsteady natural convection flow of micropolar fluid along a vertical flat plate in a thermally stratified medium. The application of two-parameter transformation group reduces the number of independent variables in the governing system consisting of partial differential equations and a set of auxiliary conditions from three to only one independent variable, and consequently the system of governing partial differential equations with boundary conditions reduces to a system of ordinary differential equations with appropriate boundary conditions. Numerical solution of the velocity, microrotation and heat transfer have been obtained. The possible forms of the ambient temperature variation with position and time are derived.     

Lie group transformations for self‐similar shocks in a gas with dust particles

In this paper, Lie group of transformation method is used to investigate the self‐similar solutions for the system of partial differential equations describing one‐dimensional unsteady plane flow of an inviscid gas with large number of small dust particles. The forms of the drag force D and the heat transfer rate Q experienced by the particle not in equilibrium with the gas have been derived for which the system of equations admits self‐similar solutions. A particular solution to the problem in one case have been found out and is used to study the effect of the dust particles on the similarity exponent. Copyright © 2009 John Wiley & Sons, Ltd.     

Stability in distribution of neutral stochastic functional differential equations with Markovian switching

Stability in distribution of stochastic differential equations with Markovian switching and stochastic differential delay equations with Markovian switching have been studied by several authors and this kind of stability is an important property for stochastic systems. There are several papers which study this stability for stochastic differential equations with Markovian switching and stochastic differential delay equations with Markovian switching technically. In our paper, we are concerned with the general neutral stochastic functional differential equations with Markovian switching and we derive the sufficient conditions for stability in distribution. At the end of our paper, one example is established to illustrate the theory of our work.     

Multiple level finite element analysis and its applications to tidal current flow in Tokyo Bay

A finite element method for the analysis of a one level and a multiple level current flow is presented. The basic equations can be derived from the three-dimensional Navier-Stokes equations under the shallow water assumptions. The standard finite element method has been introduced using the linear interpolation function based on a triangular finite element. For each level, the finite element subdivisions are not required to be coincident. To integrate the discretized equations numerically in time, an improved two step explicit scheme is employed. The multiple level finite element method is applied to a tidal flow analysis of Tokyo Bay.The multiple level tidal flow analysis is performed at the entrance channel of Tokyo Bay. The density of water is assumed to be constant for each level. The vertical profiles of the numerical velocity are compared with those of the observed velocity. The flow directions and the order of velocity are both well in agreement with the observed data. The tidal flow pattern in Tokyo Bay has been shown to be expressed by the multiple level flow assuming that the density of seawater is levelwise constant.The numerical tidal flow computation of Tokyo Bay carried out using a one level model is compared with observed data. The one level numerical values will be used to specify the boundary conditions for the multiple level analysis. Both numerical and observed results correspond extremely well in this computation. The two dominant circulated residual flows have been computed, and they coincide with the observed facts.     

Relationship between the Udwadia–Kalaba equations and the generalized Lagrange and Maggi equations

In their paper “A New Perspective on Constrained Motion,” F. E. Udwadia and R. E. Kalaba propose a new form of matrix equations of motion for nonholonomic systems subject to linear nonholonomic second-order constraints. These equations contain all of the generalized coordinates of the mechanical system in question and, at the same time, they do not involve the forces of constraint. The equations under study have been shown to follow naturally from the generalized Lagrange and Maggi equations; they can be also obtained using the contravariant form of the motion equations of a mechanical system subjected to nonholonomic linear constraints of second order. It has been noted that a similar method of eliminating the forces of constraint from differential equations is usually useful for practical purposes in the study of motion of mechanical systems subjected to holonomic or classical nonholonomic constraints of first order. As a result, one obtains motion equations that involve only generalized coordinates of a mechanical system, which corresponds to the equations in the Udwadia–Kalaba form.     

Dissipative mechanism of a semilinear higher order parabolic equation in RN

It is known that the concept of dissipativeness is fundamental for understanding the asymptotic behavior of solutions to evolutionary problems. In this paper we investigate the dissipative mechanism for some semilinear fourth-order parabolic equations in the spaces of Bessel potentials and discuss some weak conditions that lead to the existence of a compact global attractor. While for second-order reaction–diffusion equations the dissipativeness mechanism has already been satisfactorily understood (see Arrieta et al. (2004), doi:10.1142/S0218202504003234 [7]), for higher order problems in unbounded domains it has not yet been fully developed. As shown throughout the paper, one of the main differences from the case of reaction–diffusion equations stems from the lack of a maximum principle. Thus we have to rely here on suitable energy estimates for the solutions. As in the case of second-order reaction–diffusion equations, we show here that both linear and nonlinear terms have to collaborate in order to produce dissipativeness. Thus, the dissipative mechanisms in second-order and fourth-order equations are similar, although the lack of a maximum principle makes the proofs more difficult and the results not as complete.Finally, we make essential use of the sharp results of Cholewa and Rodriguez-Bernal (2012), doi:10.1016/j.na.2011.08.022 [12], on linear fourth-order equations with a very large class of linear potentials.     

Spatial Lanchester models

Lanchester equations have been widely used to model combat for many years, nevertheless, one of their most important limitations has been their failure to model the spatial dimension of the problems. Despite the fact that some efforts have been made in order to overcome this drawback, mainly through the use of Reaction–Diffusion equations, there is not yet a consistently clear theoretical framework linking Lanchester equations with these physical systems, apart from similarity. In this paper, a spatial modeling of Lanchester equations is conceptualized on the basis of explicit movement dynamics and balance of forces, ensuring stability and theoretical consistency with the original model. This formulation allows a better understanding and interpretation of the problem, thus improving the current treatment, modeling and comprehension of warfare applications. Finally, as a numerical illustration, a new spatial model of responsive movement is developed, confirming that location influences the results of modeling attrition conflict between two opposite forces.     

Separation of the elasticity theory equations with radial inhomogeneity : PMM vol. 38, n 6, 1974, pp. 1139–1144

The separation of a system of three elasticity theory equations in the static case to a system of two equations and one independent equation for a space with a radial inhomogeneity is presented in a spherical coordinate system. These equations are solved by separation of variables for specific kinds of radial inhomogeneity. In particular, solutions are found for the Lamé coefficients μ = const, λ (ifr) is an arbitrary function, μ = μ_or^β, λ = λ_or^β.While methods of solving problems associated with the equilibrium of an elastic homogeneous sphere have been studied sufficiently [1], problems with spherical symmetry of the boundary conditions have mainly been solved for an inhomogeneous sphere [2, 3],For a particular kind of inhomogeneity dependent on one Cartesian coordinate, the equations have been separated completely in [4], A system of three equations with a radial inhomogeneity in a spherical coordinate system is separated below by a method analogous to [4].     

The new exact solutions of variant types of time fractional coupled Schr\"{o}dinger equations in plasma physics

In the present article, the new exact solutions of fractional coupled Schr\"{o}dinger type equations have been studied by using a new reliable analytical method. We applied a relatively new method for finding some new exact solutions of time fractional coupled equations viz. time fractional coupled Schr\"{o}dinger--KdV and coupled Schr\"{o}dinger--Boussinesq equations. The fractional complex transform have been used here along with the property of local fractional calculus for reduction of fractional partial differential equations (FPDE) to ordinary differential equations (ODE). The obtained results have been plotted here for demonstrating the nature of the solutions.     

The diffraction of plane elastic unsteady waves by a delaminated inclusion in the case of smooth contact in the delamination region

A solution of the problem of the diffraction of unsteady elastic waves by a thin strip-like delaminated rigid inclusion in an unbounded elastic medium under conditions of planer strain is proposed. We have in mind an inclusion, one side of which is completely bonded with the medium while, the other side is delaminated and conditions of smooth contact are satisfied on it. The method of solution is based on the use of discontinuous solutions of the Lamé equations of motion under conditions of planer strain, which have been constructed earlier in the space of Laplace transforms. As a result, the problem reduces to solving a system of three singular integral equations for the transforms of the unknown discontinuities. The inverse transforms are found by a numerical method, based on the replacement of a Mellin integral by a Fourier series.     

Geometrically non-linear equations in the theory of momentless shells with applications to problems on the non-classical forms of loss of stability of a cylinder

Geometrically non-linear and linearized equations in the theory of momentless shells are set up based on the kinematic relations in [Paimushin VN, Shalashilin VI. Relations of the theory of deformations in the quadratic approximation and problems of constructing refined versions of the geometrically non-linear theory of laminar structural components. Prikl Mat Mekh 2005; 69(5): 861–81]. The use of these equations, unlike in the case of the well-known equations, enables one to avoid the occurrence of spurious bifurcation points in solving real problems. Non-classical problems of the stability of cylindrical shells under an external pressure, axial compression and torsion are considered, which can be formulated on the basis of the derived equations of the theory of momentless shells. Their exact analytical solutions are found and enable one to estimate the quality of the previously obtained relations [Paimushin VN, Shalashilin VI. Relations of the theory of deformations in the quadratic approximation and problems of constructing refined versions of the geometrically non-linear theory of laminar structural components. Prikl Mat Mekh 2005; 69(5): 861–81] and the richness of content of the equations which have been constructed compared with well-known equations in the mechanics of thin shells. It is established that the majority of the new forms of loss of stability of cylindrical shells which are revealed relate to a number of shear forms, the onset of which is possible before the flexural forms which have been well studied up to now, in the case of small values of the shear modulus of a shell material with a very highly pronounced anisotropy in its properties.     

Implementation of a finite‐volume method for the determination of effective parameters in fissured porous media

Many studies have proposed one‐equation models to represent transport processes in heterogeneous porous media. This approach is based on the assumption that dependent variables such as pressure, temperature, or concentration can be expressed in terms of a single large‐scale averaged quantity in regions having very different chemical and/or mechanical properties. However, one can also develop large‐scale averaged equations that apply to the distinct regions that make up a heterogeneous porous medium. This approach leads to region‐averaged equations that contain traditional convective and dispersive terms, in addition to exchange terms that account for the transfer between the different media. In our approach, the fissures represent one region, and the porous media blocks represent the second region. The analysis leads to upscaled equations having a domain of validity that is clearly identified in terms of time and length‐scale constraints. Closure problems are developed that lead to the prediction of the effective coefficients that appear in the region averaged equations, and the main purpose of this article is to provide solutions to those closure problems. The method of solution makes use of an unstructured grid and a joint element method in order to take care of the special characteristics of the fissured network. This new numerical method uses the theory developed by Quintard and Whitaker and is applied on considerably more complex geometries than previously published results. It has been tested for several special cases such as stratified systems and “sugarbox” media, and we have compared our calculations with other computational methods. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 237–263, 2000     

Variational iteration method for solving sixth-order boundary value problems

In this paper, we have shown that sixth-order boundary value problems can be transformed into a system of integral equations, which can be solved by using variational iteration method. The analytical results of the equations have been obtained in terms of convergent series with easily computable components. Several examples are given to illustrate the efficiency and implementation of the method. It is observed that the proposed technique is more useful and is easier to implement because one does not need to calculate the Adomian’s polynomials which is itself a difficult task.     

Effect of Space Dimensions on Equilibrium Solutions of Cahn-Hilliard and Conservative Allen-Cahn Equations

In this study, we investigate the effect of space dimensions on the equilibrium solutions of the Cahn-Hilliard (CH) and conservative Allen-Cahn (CAC) equations in one, two, and three dimensions. The CH and CAC equations are fourth-order parabolic partial and second-order integro-partial differential equations, respectively. The former is used to model phase separation in binary mixtures, and the latter is used to model mean curvature flow with conserved mass. Both equations have been used for modeling various interface problems. To study the space-dimension effect on both the equations, we consider the equilibrium solution profiles for symmetric, radially symmetric, and spherically symmetric drop shapes. We highlight the different dynamics obtained from the CH and CAC equations. In particular, we find that there is a large difference between the solutions obtained from these equations in three-dimensional space.