A representation theorem for volume-preserving transformations

A general solution is presented for the partial differential equation ∂u/∂x=k(x), where u and x are n-vector fields, ∂u/∂x denotes the Jacobian of the transformation x→u and k(x) is a scalar-valued function. The solution for the case k(x)=1 is of special interest because it furnishes a representation theorem for volume-preserving transformations in an n-dimensional space. Such a representation for the case n=2 was obtained by Gauss. The solution for n=3, presented here, furnishes a representation for isochoric (volume-preserving) finite deformations, which are important in the mechanics of highly deformable incompressible solid materials.     

On steady plane flows in classical elasticity

A new class of boundary-value problems in mathematical elasticity is proposed, wherein the medium flows steadily relative to a non-embedded surface over which tractions or velocities are prescribed. Such flows are seen in metal forming operations where purely elastic streams enter and leave the working zone. The deformations are assumed here to be plane and isochoric. A general solution is formulated in terms of two complex potentials. Residual stress is accounted for in detail and a uniqueness theorem is proved. Some simple flows are examined, but it remains to develop a systematic procedure for matching the general solution to arbitrary boundary data.     

Similarity solutions for creeping flow and heat transfer in second grade fluids

The two-dimensional equations of motions for the slowly flowing and heat transfer in second grade fluid are written in cartesian coordinates neglecting the inertial terms. When the inertia terms are simply omitted from the equations of motions the resulting solutions are valid approximately for Re?1. This fact can also be deduced from the dimensionless form of the momentum and energy equations. By employing Lie group analysis, the symmetries of the equations are calculated. The Lie algebra consist of four finite parameter and one infinite parameter Lie group transformations, one being the scaling symmetry and the others being translations. Two different types of solutions are found using the symmetries. Using translations in x and y coordinates, an exponential type of exact solution is presented. For the scaling symmetry, the outcoming ordinary differential equations are more involved and only a series type of approximate solution is presented. Finally, some boundary value problems are discussed.     

A general Struble's technique for solving an nth order weakly non-linear differential system with damping

A general formula (based on the method of variation of parameters) has been presented for determining an approximate solution of an nth order n=2,3,… weakly non-linear differential system with several damping effects. The general solution covers the under-damped, undamped and over-damped cases. The formulation as well as determination of the solution is simple. The method is illustrated by several examples.     

On relations for general two-dimensional ideal plasticity problems under the full plasticity condition

Relations for two-dimensional ideal plasticity problems under the full plasticity condition are determined with material anisotropy, inhomogeneity, and compressibility properties taken into account. These properties are determined by the direction cosines of the principal stress, the coordinates of points in space, and the mean stress.For the yield strength we take a function of the form k = k(σ, n ₁, n ₂, n ₃, x, y, z). The desired relations are determined for the general plane ideal plasticity problem. The relations thus obtained are generalized to the cases of axisymmetric and spherical plasticity problems.     

Axisymmetric spreading of a thin liquid drop with suction or blowing at the horizontal base

The axisymmetric spreading under gravity of a thin liquid drop on a horizontal plane with suction or blowing of fluid at the base is considered. The thickness of the liquid drop satisfies a non-linear diffusion equation with a source term. For a group invariant solution to exist the normal component of the fluid velocity at the base, v_n, must satisfy a first-order quasi-linear partial differential equation. The general form of the group invariant solution for the thickness of the liquid drop and for v_n is derived. Two particular solutions are considered. Each solution depends essentially on only one parameter which can be varied to yield a range of models. In the first solution, v_n is proportional to the thickness of the liquid drop. The base radius always increases even for suction. In the second solution, v_n is proportional to the gradient of the thickness of the liquid drop. The thickness of the liquid drop always decreases even for blowing. A special case is the solution with no spreading or contraction at the base which may have application in ink-jet printing.     

Global Solution to the Generalized Riemann Problem in the Presence of a Boundary and Contact Discontinuities

This work is a continuation of our previous work. In the present paper we study the global structure stability of the Riemann solution $u=U(\frac{x}{t})$ containing only contact discontinuities for general n×n quasilinear hyperbolic systems of conservation laws in the presence of a boundary. We prove the existence and uniqueness of a global piecewise C ¹ solution containing only contact discontinuities to a class of the generalized Riemann problems for general n×n quasilinear hyperbolic systems of conservation laws in a half space. Our result indicates that this kind of Riemann solution $u=U(\frac{x}{t})$ mentioned above for general n×n quasilinear hyperbolic systems of conservation laws in the presence of a boundary possesses a global nonlinear structure stability. Some applications to quasilinear hyperbolic systems of conservation laws occurring in physics and other disciplines, particularly to the system describing the motion of the relativistic string in Minkowski space R ^1?+?n , are also given.     

Invariant linearization criteria for a three-dimensional dynamical system of second-order ordinary differential equations and applications

Second-order dynamical systems are of paramount importance as they arise in mechanics and many applications. It is essential to have workable explicit criteria in terms of the coefficients of the equations to effect reduction and solutions for such types of equations. One important aspect is linearization by invertible point transformations which enables one to reduce a non-linear system to a linear system. The solution of the linear system allows one to solve the non-linear system by use of the inverse of the point transformation. It was proved that the n-dimensional system of second-order ordinary differential equations obtained by projecting down the system of geodesics of a flat (n+1)-dimensional space can be converted to linear form by a point transformation. This is a generalization of the Lie linearization criteria for a scalar second-order equation. In this case it is of the maximally symmetric class for a system and the linearizing transformation as well as the solution can be directly written down. This was explicitly used for two-dimensional dynamical systems. The criteria were written down in terms of the coefficients and the linearizing transformation allowed for the general solution of the original system. Here the work is extended to a three-dimensional dynamical system and we find explicit criteria, including the linearization test given in terms of the coefficients of the cubic in the first derivatives of the system and the construction of the transformations, that result in linearization. Applications to equations of classical mechanics and relativity are given to illustrate our results.     

Determination of a dynamic model for urethane prosthetic compounds

A frequency-response testing technique for determining the dynamic behavior of urethane prosthetic compounds is discussed. Experimental preparation of strip specimens and test results are presented. Sinusoidal response data of the tested strips are compared with three computer synthesizers (viscoelastic, viscous and complex modulus) of the one-dimensional wave equation for deciding on a model which best represents the material and subsequently calculating the value of a dynamic loss factor. The closed form solutions for the three mathematical models subject to sinusoidal boundary conditions are expressed in terms of functions which can be easily programmed for machine computation in FORTRAN IV language involving complex arguments. Dynamicloss factors are required for experimental and finite-element studies of prosthetic left ventricle of animal and human hearts. For the test strips, the loss factor, δ, is found to be dependent on the frequency.f. The results presented show that a viscoelastic model with a frequency-dependent loss factor of the form δ=c/f ⁿ is an excellent representation for analyzing the dynamic behavior of urethane compounds subjected to a range of frequencies corresponding to heart-beat rates. Derivation of the model parametersc andn is explained in detail.     

Growth instabilities and folding in tubular organs: A variational method in non-linear elasticity

Morphoelastic theories have demonstrated that elastic instabilities can occur during the growth of soft materials, initiating the transition toward complex patterns. Within the framework of non-linear elasticity, the theory of incremental elastic deformations is classically employed for solving stability problems with finite strains. In this work, we define a variational method to study the bifurcation of growing cylinders with circular section. Accounting for a constant axial pre-stretch, we define a set of canonical transformations in mixed polar coordinates, providing a locally isochoric mapping. Introducing a generating function to derive an implicit gradient form of the mixed variables, the incompressibility constraint for the elastic deformation is solved exactly. The canonical representation allows to transform a generic boundary value problem, characterized by conservative body forces and surface traction loads, into a completely variational formulation. The proposed variational method gives a straightforward derivation of the linear stability analysis, which would otherwise require lengthy manipulations on the governing incremental equations. The definition of a generating function can also account for the presence of local singularities in the elastic solution. Bifurcation analysis is performed for few constrained growth problems of biomechanical interests, such as the mucosal folding of tubular tissues and surface instabilities in tumor growth. In a concluding section, the theoretical results are discussed for clarifying how anisotropy, residual strains and external constraints can affect the stability properties of soft tissues in growth and remodeling processes.     

Determining the limiting solutions of nonstationary axisymmetric Hele-Shaw problems

Numerical solution of the Hele-Shaw problem reduces to solution of three boundary-value problems of determining analytic functions of a complex variable in each time step: conformal mapping of the range of the parametric variable to the physical plane, the Dirichlet problems for determining the electric-field strength, and the Riemann-Hilbert problem for calculating partial time derivatives of the coordinates of points of the interelectrode space (the images of the points on the boundary of the parametric plane are fixed). Unlike in the two-dimensional problem, the electric-field strength is determined using integral transformations of an analytic function. Approximation by spline function is performed, and more accurate and steady (than the well-known ones) general solution algorithms for the nonstationary axisymmetric problems are described. Results of a numerical study of the formation of stationary and self-similar configurations are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 87–99, July–August, 2009.     

The free energy of mixing for n-variant martensitic phase transformations using quasi-convex analysis

The construction of effective models for materials that undergo martensitic phase transformations requires usable and accurate functional representations for the free energy density. The general representation of this energy is known to be highly non-convex; it even lacks the property of quasi-convexity. A quasi-convex relaxation, however, does permit one to make certain estimates and powerful conclusions regarding phase transformation. The general expression for the relaxed free energy is however not known in the n-variant case. Analytic solutions are known only for up to 3 variants, whereas cases of practical interests involve 7-13 variants. In this study we examine the n-variant case utilizing relaxation theory and produce a seemingly obvious but very powerful observation regarding a lower bound to the quasi-convex relaxation that makes practical evolutionary computations possible. We also examine in detail the 4-variant case where we explicitly show the relation between three different forms of the free energy of mixing: upper bound by lamination, the Reuß lower bound, and a lower estimate of the -measure bound. A discussion of the bounds and their utility is provided; sample computations are presented for illustrative purposes.     

Direct numerical simulation of supersonic pipe flow at moderate Reynolds number

We study compressible turbulent flow in a circular pipe at computationally high Reynolds number. Classical related issues are addressed and discussed in light of the DNS data, including validity of compressibility transformations, velocity/temperature relations, passive scalar statistics, and size of turbulent eddies. Regarding velocity statistics, we find that Huang’s transformation yields excellent universality of the scaled Reynolds stresses distributions, whereas the transformation proposed by Trettel and Larsson (2016) yields better representation of the effects of strong variation of density and viscosity occurring in the buffer layer on the mean velocity distribution. A clear logarithmic layer is recovered in terms of transformed velocity and wall distance coordinates at the higher Reynolds number under scrutiny (Re_τ ≈ 1000), whereas the core part of the flow is found to be characterized by a universal parabolic velocity profile. Based on formal similarity between the streamwise velocity and the passive scalar transport equations, we further propose an extension of the above compressibility transformations to also achieve universality of passive scalar statistics. Analysis of the velocity/temperature relationship provides evidence for quadratic dependence which is very well approximated by the thermal analogy proposed by Zhang et al. (2014). The azimuthal velocity and scalar spectra show an organization very similar to canonical incompressible flow, with a bump-shaped distribution across the flow scales, whose peak increases with the wall distance. We find that the size growth effect is well accounted for through an effective length scale accounting for the local friction velocity and for the local mean shear.     

Smooth positon solutions of the focusing modified Korteweg–de Vries equation

The n-fold Darboux transformation \(T_{n}\) of the focusing real modified Korteweg–de Vries (mKdV) equation is expressed in terms of the determinant representation. Using this representation, the n-soliton solutions of the mKdV equation are also expressed by determinants whose elements consist of the eigenvalues \(\lambda _{j}\) and the corresponding eigenfunctions of the associated Lax equation. The nonsingular n-positon solutions of the focusing mKdV equation are obtained in the special limit \(\lambda _{j}\rightarrow \lambda _{1}\), from the corresponding n-soliton solutions and by using the associated higher-order Taylor expansion. Furthermore, the decomposition method of the n-positon solution into n single-soliton solutions, the trajectories, and the corresponding “phase shifts” of the multi-positons are also investigated.     

Non-Newtonian pseudoplastic fluids: Analytical results and exact solutions

A one layer model of laminar non-Newtonian fluids (Ostwald-de Waele model) past a semi-infinite flat plate is revisited. The stretching and the suction/injection velocities are assumed to be proportional to x^1/(1−2n) and x⁻¹, respectively, where n is the power-law index which is taken in the interval . It is shown that the boundary-layer equations display both similarity and pseudosimilarity reductions according to a parameter γ, which can be identified as suction/injection velocity. Interestingly, it is found that there is a unique similarity solution, which is given in a closed form, if and only if γ=0 (impermeable surface). For γ≠0 (permeable surface) we obtain a unique pseudosimilarity solution for any 0≠γ≥−((n+1)/3n(1−2n))^n/(n+1). Moreover, we explicitly show that any pseudosimilarity solution exhibits similarity behavior and it is, in fact, similarity solution to a modified boundary-layer problem for an impermeable surface. In addition, the exact similarity solution of the original boundary-layer problem is used, via suitable transverse translations, to construct new explicit solutions describing boundary-layer flows induced by permeable surfaces.     

Multiple-zone sliding contact with friction on an anisotropic thermoelastic half-space

The paper studies a class of multiple-zone sliding contact problems. This class is general enough to include frictional and thermal effects, and anisotropic response of the indented material. In particular, a rigid die (indenter) slides with Coulomb friction and at constant speed over the surface of a deformable and conducting body in the form of a 2D half-space. The body is assumed to behave as a thermoelastic transversely isotropic material. Thermoelasticity of the Green–Lindsay type is assumed to govern. The solution method is based on integral transforms and singular integral equations. First, an exact transform solution for the auxiliary problem of multiple-zone (integer n > 1) surface tractions is obtained. Then, an asymptotic form for this auxiliary problem is extracted. This form can be inverted analytically, and the result applied to sliding contacts with multiple zones. For illustration, detailed calculations are provided for the case of two (n = 2) contact zones. The solution yields the contact zone width and location in terms of sliding speed, friction, die profile, and also the force exerted. Calculations for the hexagonal material zinc illustrate effects of speed, friction and line of action of the die force on relative contact zone size, location of maximal values for the temperature and the compressive stress, and the maximum temperature for a given maximum stress. Finally, from our general results, a single contact zone solution follows as a simple limit.     

Extrema of Young’s modulus for cubic and transversely isotropic solids

For a homogeneous anisotropic and linearly elastic solid, the general expression of Young’s modulus E(n), embracing all classes that characterize the anisotropy, is given. A constrained extremum problem is then formulated for the evaluation of those directions n at which E(n) attains stationary values. Cubic and transversely isotropic symmetry classes are dealt with, and explicit solutions for such directions n are provided. For each case, relevant properties of these directions and corresponding values of the modulus are discussed as well. Results are shown in terms of suitable combinations of elements of the elastic tensor that embody the discrepancy from isotropy. On the basis of such material parameters, for cubic symmetry two classes of behavior can be distinguished and, in the case of transversely isotropic solids, the classes are found to be four. For both symmetries and for each class of behavior, some examples for real materials are shown and graphical representations of the dependence of Young’s modulus on direction n are given as well.     

On the Canonical Elastic Moduli of Linear Plane Anisotropic Elasticity

The linear, planar, anisotropic elastic equilibrium equations are transformed to canonical form, through linear transformations of both coordinates and unknown displacement functions, together with a linear combination of equations. Correspondingly, the six original material moduli are replaced by two canonical elastic moduli. Similar results have been reached by Olver in 1988. However, the method demonstrated in this paper is more concise and direct. As an example, the general solution to the canonical equations is obtained in the case of a pair of double roots.