Influence of the Angular Momentum in Problems Continuum Mechanics
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova Influence of the Angular Momentum in Problems Continuum Mechanics EVELINA PROZOROVA Mathematical-Mechanical Department St. Petersburg State University Av. 28 , Peterhof, 198504 RUSSIA Abstract: - For continuum mechanics a model is proposed, that is built with consideration outside the integral term when deriving conservation laws using the Ostrogradsky-Gauss theorem. Performed analysis shows discrepancy between accepted classical conservation laws and classical theoretical mechanics and mathematics. As a result, the theory developed for potential flows was extended to flows with significant gradients of physical parameters. We have proposed a model that takes into account the joint implementation of the laws for balance of forces and angular momentums. It does not follow from the Boltzmann equation that the pressure in the Euler and Navier-Stokes equations is equal to one third of the sum the pressures on the corresponding coordinate axes. The vector definition of pressure is substantiated. It is shown that the symmetry condition for the stress tensor is one of the possible conditions for closing the problem. An example of solving the problem of the theory of elasticity is given. Keywords: angular moment, potential flow, circulation, Euler equation, Pascal’s law Received: September 26, 2020. Revised: February 13, 2021. Accepted: February 26, 2021. Published: March 4, 2021. 1. Introduction calculation of stresses in the rods [1,2] and In aeromechanics, shipbuilding, the physics in the calculations of viscoelastic problems of earth and atmosphere, there are no [3]. The probable reason for ignoring the understanding of some processes associated effect in the general theory is the proof of its with the interaction of flows with body absence, built on the basis of the theory, which are moving at high speeds; there are taking into account the interaction of no adequate models for the formation of neighboring elementary volumes only along vortex structures, etc. The objects of study the normal to the surface [4-6]. In this case, in continuum mechanics are solids, liquid, one can only prove the consistency of the gas, and plasma. In a mathematical whole theory. The influence of the rotational description, the basic laws are conservation part of the stress tensor can be traced by laws. This unites them. Therefore, the work analyzing the Hamel solution and its discusses general questions arising when generalizations [7]. Dividing the speed into writing conservation laws. The law of divergent and rotational components about equilibrium of forces was put as basis of the an axis passing through an arbitrary point classical theory, the law of conservation does not right, as according to theory, only angular momentum was considered as movement around the axis of inertia need to consequence of the fulfillment the law of consider. A model that is ignoring the balance forces. Each of laws was considered rotation of the elementary volume leads to independently. However, the gradient of the the symmetry of the stress tensor, which distributed moment leads to the formation of violates the "continuity" of the medium [8]. an additional force, that we need to take into When analyzing the derivation of account in the equations of motion. This is differential conservation laws in continuum evidenced by the works devoted to the mechanics, it turned out that the transition E-ISSN: 2224-3429 1 Volume 16, 2021
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova from the integral formulation to the version, use only the equilibrium part of the differential formulation (integration by Chapman-Enskog function, the rest of the parts) was performed without taking into terms included in the function do not account outside integral component. In contribute. Consequently, the full account of classical version, the Ostrogradsky-Gauss the viscous contribution is ensured by the theorem is formulated for a motionless body full velocity, in fact, by its corresponding without rotation [9]. Thus, the generalization first derivatives. However, the of the theorem to body that is rotating [10] macroparameters included in the equilibrium leads to a stronger influence of the no distribution function must be calculated not symmetry of the stress tensor in the entire according to zero approximation from the mechanics of a continuous medium. After Euler equations, but from the equations of integration by parts, the transition from the the first approximation (the Navier-Stokes integral formulation to the differential equations, if they are used). The equations formulation is difficult. However, one can of motion obtained from the Boltzmann take advantage of physical and geometric equation correspond to the zero order for the considerations and derive equations from Euler equations and the first order for the them. At present, the need to divide the Navier-Stokes equations. Using only stress tensor into two parts is not clear. The divergent part of velocity without rotation it speeds of various processes at the time of we obtain speed different from the initial writing the equations were relatively low one. The situation is saved by the fact that compared to modern one. Subsequently, the this circumstance affects the viscous theory that was developed for potential component of stresses, which is a first-order flows was used to flows with significant quantity. Hydrostatic pressure is a zero- gradients of physical parameters was order quantity. Using Pascal's law for expanded. The lack of symmetry leads to an equilibrium, the pressure is chosen equal to increase in the order of the derivatives the one third of the pressure on the coordinate Navier-Stokes equations. The role of areas. It does not follow from the Boltzmann boundary conditions increases accordingly. equation that the pressure in the Euler Their specific form is determined by the equation is equal to one third of the sum of problem under consideration. Since the the pressures on the corresponding original formulation obtained in experiments coordinate axes. There is no experimental is integral, then numerical methods using confirmation of this fact. However, the conservation laws in integral form are theory remains the same when determining preferable. These methods include the finite the different pressure on each of the sites, volume method. The implementation of the i.e. the use of one pressure is possible under method is achieved in various ways. equilibrium conditions (Pascal's law), but for Equations of an incompressible fluid are no equilibrium conditions the fact is not best approximated. For a compressible fluid, obvious. Attention is drawn to this in the the computational process becomes more textbook [11] and articles [12,13]. In the complicated. It is important to use the theory of elasticity, when considering the integral conservation laws obtained from the relationship between the components of the experiment. The transition from the deformation tensor and the stress tensor, the currently used equations to integral ones due experimental fact of the change in the to the rejection to investigate of the rotation components (stresses) normal to the effects does not allow the conservation laws elementary site is used in proportion to the to be fulfilled. The experimentally sum of the stresses from others components, confirmed Newton's law on the relationship and they all differ. If problems of the theory between the stress and strain tensor of elasticity was solving with a symmetric concerned the total velocity with the no tensor, it is necessary to satisfy the symmetric tensor. In the kinetic theory, the compatibility conditions [5,6]. To get coefficients of viscosity and thermal around these conditions when solving the conductivity, even in no equilibrium problem numerically, the function values in E-ISSN: 2224-3429 2 Volume 16, 2021
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova the cell are “averaged”. The result of the contribution to friction and pressure. calculation is then the values in the center. Measurements performed on the axes of This is not discussed in the articles, but in symmetry do not allow speaking about the private conversations the fact is admitted. degree of asymmetry of the stress tensor. Experimental fact of the no symmetry of the The existence of invariants higher than the stress tensor can be found in the works first and, accordingly, the principal axes of [14,15], theory in [16]. The existing theory the elementary volume in the proposed is connected with the fact that the derivation theory of elasticity becomes problematic. of conservation laws in the theory of The question of equilibrium in mechanics of elasticity excludes the contribution of the a continuous medium and the extension of distributed moment to the conditions of Pascal's law to the case of small equilibrium of forces. An interesting feature deformations, the need to replace the of the equations of hydrodynamics is the averaged pressure with its theoretical values absence of nesting of models equations into , , are discussed. The rejection each other. For example, the potential flow of pressure averaging makes it possible to cannot be obtained from the Euler equations. explain the formation of vortex bundles Static pressure, as follows from the kinetic when flowing around bodies in the vicinity theory, is a zero-order quantity, and the of the separation point. It is very difficult to terms associated with dissipative effects are judge the influence of the moment, limiting first-order quantities. As already noted, it ourselves to the results of calculating two- does not follow from the Boltzmann dimensional problems. If the tensor is equation that the pressure in the Euler asymmetric, then we have three equations equation is equal to one third of the sum of for four unknowns. The system is not the pressures on the corresponding closed. Therefore, additional assumptions coordinate axes. Inaccuracy in determining have to be made. One of these assumptions the velocities in the stress tensor does not is the assumption about the symmetry of the greatly affect the results. For an stress tensor. These issues are discussed in incompressible fluid, the speeds are the this work. same. The most important characteristic of an elementary volume is its position of the center of inertia and the correct projection of 2. Origin of the angular the boundaries of an arbitrary volume on the momentum coordinate axes. Due to the arbitrariness of Consider the interaction of three particles the directions of motion, the viscous interacting with each other. Since the components of the velocity must be taken angular momentum is equal to the force into account even in one-dimensional acting on the center of inertia, that is vector motion. With any restructuring of the flow, multiplied by the radius vector, we can the position of the center of inertia changes, a moment of force is created, the gradient of calculate the changes in the moment after a which is the force itself. Additional force is certain period of time added with stresses and change the direction of movement and the speed of the molecules. The result is a change in the equation of state even under equilibrium conditions and the formation of fluctuation phenomena. The equation of state (Newton's law) in classical mechanics ensures the symmetry of the stress tensor. Note again that the classical equations do not include all Fig. 1. Interaction of three particles the speeds that need to be considered. The vortex part is discarded. The experiment + + does not separate speeds, measuring the total = , + + E-ISSN: 2224-3429 3 Volume 16, 2021
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova ( , 1 , 2 , . . , 1 , 2 , … ) 1 … 1 … 1 ( 1 + 1̇ ∆ )+ 2 ( 2 + 2̇ ∆ )+ 3( 3 + 3̇ ∆ ) . To construct an equation directly for the +∆ = 1 + 2 + 3 velocity projection, we multiply the equation . by the velocity projection ξ on the At the next moment in time, the center of coordinate axis, for example, . Let us have inertia will shift. Therefore, an additional different pressure values 1 , 2 , 3 : force arises, since the gradient of the = 3 − , moment is a force. Thus, the moment creates = − 3 + 1 , a collective force for any movement of = 2 − 1 . particles with different speeds. The effect Really, works in the formation of fluctuations and must be taken into account when calculating the equation of state. In continuum [ ] mechanics, the integral conservation laws 1 2 3 obtained experimentally are written in differential form. Consequently, the impossibility of the same When passing to differential equations, the values of pressure p. Pressure is a vector like Ostrogradsky-Gauss theorem is used for the force. In this case, the Lamb equation [7] selected fixed volume, ignoring the possible has the form rotation of the elementary volume. The 2 discarded term represents the velocity + ( ) + × = − 2 1 circulation, its rotor part. Circulation is the ⃗ action of the moment. Consequently, the and there is a dependence of density as an rotational terms (moment) discarded in the average value over an elementary volume, construction of the Navier-Stokes equations then the Bernoulli equation makes sense. are formed not only in the case of large The finite volume method is often used to gradients, but also due to the additional numerically solve the Euler and Navier- circulation of the velocity around the Stokes equations. However, its formulation elementary volume [10]. is based on differential equations with a The influence of the moment is also symmetric stress tensor. In the mechanics of manifested through pressure. There are no liquid and gas, it is not customary to average experiments confirming Pascale's law in the calculation results over the calculated non-equilibrium cases. The kinetic theory cells. Consequently, since the mesh has (modified and classical) also does not follow finite dimensions, in old works, the no Pascal's law. symmetry of the tensor is partially taken into Consider the classical Boltzmann equation account in numerical calculations. In solid (the law of momentum does not hold) mechanics, averaging is offer performed 1 when solving problems by the finite element ( , ) = ∫ ( , , ) , method. With a large number of steps, this = ∫ ( , , ) , gives significant errors. The main argument ; = ∫ 2 ( , , ) , for choosing a symmetric tensor is the use of = − . the equilibrium equation for the force without the influence of the moment. ( + , + , + )drdξj = However, the momentum gradient is force. ( , , t ) + ( ) . Consequently, the requirement of distribution function, - radius vector; - simultaneous fulfillment of the laws of point coordinate; is the velocity of a point, balance of forces and moments of forces is the molecular weight, and, according to changes the equations. the definition of the distribution function , In general case the probability of finding the system at points ( , ) in the intervals is E-ISSN: 2224-3429 4 Volume 16, 2021
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova Here all designations are ( + + + ) = 1 standard, , , forces created by + + + + the moment, , , are external moments. ( + + + )= 3. Specific tasks + In the three-dimensional case, we have six 2 + + + unknowns and six equations for the no symmetric tensor, but used only three. In the flat version, there is one equation for the ( + + + )= moments and two for movement, unknown + four. Closing can be by the symmetry 3 + + + condition of the stress tensor. But you can do it in another way. Some problems with symmetric tensor cannot be solved. Let us ( + + + 3 ) − consider the problem when the stress distributions are given, with and without ( + + + 2 ) + taking into account the influence of the − + + = 0 angular momentum. The voltages can be the same, they can be different. Considered case of elongated plate. General view of ( + + + 2 ) − equations ( + + + 1 ) + + =0, + =0, − + + = 0 ( + ) ( + + + 1 ) − − ( + ) ( + + + 2 ) + + − = 0. − + + = 0 Fig. 2. Elongated plate We will assume that there is a point 1 = 1 , = 0 , = 01 is a where the stresses coincide. We will begin solution to the system of equations with the system of equations represents the and = . first two equations, instead of the third equation we use the symmetry condition of Thus, an additional assumption that closes the stress the systems of equations on the plane in classical mechanics is the condition of the ( + )− ( + )+ stress tensor symmetry. Another option is used in work [11], the assumption about the − = 0. Then distribution of stress. For a symmetric tensor, consider the solution: let = E-ISSN: 2224-3429 5 Volume 16, 2021
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova ( − ) + + ( ) or = + Later the theory was developed on the basis ( ). Then of Langevin and Fokker-Planck equations = − ∫0 ̇ = − ̇ + ( ) , 0 = [17-23]. The evolution of a Brownian 2 particle (fluctuation) is determined by its − ∫ (− ̈ + ̇ ) = ̈ 2 - ̇ + ( ). interaction with the environment, which is 1 Let us satisfy the boundary condition always collective. In the kinetic 0 ( = ) = 0. representation, the evolution of a system of 2 0 = ̈ ( − ) /2 - ̇ ( − ). Brownian particles is described by a = − ̇ ( − ) + ( ) nonlocal equation for the n-particle Equation for moment distribution function. The Langevin and (− ̇ + ̇ ) − ( ̈ ( − ) - ̇ − ̈ Fokker-Planck equations are obtained from ( − ) + ̇ ) = − = ( , ) = 0. the Liouville equation for specially selected The solution for the no symmetric tensor can models of integral kernels using be found as follows phenomenological conservation laws. In ( + ) recent years, the molecular modeling + =0, + =0, method has been widely used. The classical Langevin equation for one ( + )− ( particle ( + ) 1 + ) + = 0. = − ѷ V+ ( ) , where ( ) is a random force. At the first iteration, we consider equations A Markov Gaussian process is considered only for φ, leaving the coefficients from the with the condition that the average for an zero iteration. The first parenthesis vanishes. ensemble of particles < > = 0 . In our We represent the second parenthesis as version, this condition is fulfilled by virtue ( + ) + of the fulfillment of the theorem on the + = + conservation of the moment in a closed = 0 + volume. For equilibrium condition, this is . Then true. Langevin equation taking into account 1 1 − ∫ the influence of the angular momentum is + = 0. = 1 = ѷ 1 = − V + , where ѷ - 1 − | 1 = , = coefficient of friction of the selected = , 2 = 2 + ( ). particle, is the mass of the particle, is the moment of force acting on the particle, V Now you can find the corrected value of . – particle velocity. In this version, the tensor will be symmetric. The classical Langevin equation for one The tensor will not be symmetric if particle , = 0, ( ) = ( − ). In this 1 case, = 0, = 0. We get the results = − ѷ V+ ( ) , where ( ) is a of work random force. E. A. Bulanova [16]. If there are A more interesting situation arises in disturbances on both sides, then stitching is Coulomb plasma. The movement of the required at some point 1 , dependence on center of inertia entails a new distribution of can be arbitrary. The influence of the charges. The result is a change in the moment creates a collective effect. Let us strength of the electric field. In this case, the consider its effect in a rarefied gas. main force will be strength associated with Fluctuations are considered as random the movement of charges. deviations from the mean. This effect was first explained by A.Einstein. It consisted in 4. Conclusion the fact that the diffusion force should be Modern models of continuum mechanics do equal to the viscous Stokes drag force not give answers to many questions when np F s ∂np describing = DE cs ∂x E-ISSN: 2224-3429 6 Volume 16, 2021
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova the processes of turbulent flows, flow [7]. N.E. Kochin, I.A. Kibel, N.V. Rose, separation, in matters of destruction of Theoretical Hydromechanics, Moscow: Fiz- solids, etc. It is proposed to extend the Mat Literature, 1963, 583 suggested model to the case of large [8].A.A. Ilyushin, Asymmetry of strain and gradients and bring the model in line with stress tensors in the mechanics of a the main provisions of classical mechanics continuous medium. WEST. Moscow. and mathematics. The necessity of using no University. Ser. 1. Mathematics- symmetric stress tensor and a vector value of pressure in the conservation laws of mechanics.1996. pp. 6-14 continuum mechanics is substantiated. [9]. V.I. Smirnov, Higher mathematics Particular examples from the theory of the course. T.2. M .: Nauka, 1974 , 655. boundary layer, the theory of elasticity and [10]. E.V Prozorova . The Effect of Angular kinetic theory were given in previous works. Momentum and Ostrogradsky-Gauss We suggest an algorithm for taking into Theorem in the Equations of Mechanics. account the no symmetric stress tensor by Wseas Transaction on fluid DOI: the example of solving the problem of 10.37394/232013.2020.15.2 calculating a beam with a given distribution [11]. L.G. Loytsyansky, Mechanics of of one of the stresses. The ambiguity of liquid and gas. M .: Nauka. 1970, 904. using a symmetric stress tensor for closing [12]. Ouadie Koubaiti, Ahmed Elkhalfi, and the possibility of using the results with a Jaouad El-Mekkaoui, Nikos Mastorakis, symmetric tensor, supplementing the Solving the Problem of Constraints Due to solution with iteration, are shown. The Dirichlet Boundary Conditions in the proposed model can be applied to a wider Context of the Mini Element Method, range of problems in continuum mechanics International Journal of Mechanics, pp. 12- than the classical theory. 21, Volume 14, 2020. [13]. Tarik Chakkour, Fayssal Benkhaldoun, References Slurry Pipeline for Fluid Transients in Pressurized Conduits, International Journal [1]. Yu.N. Shevchenko, I. V. Prokhorenko, of Mechanics, pp. 1-11, Volume 14, 2020. The theory of elastic-plastic shells under [14]. Yu.G. Stepanov, Theoretical non-isothermal loading processes. Kiev. hydrodynamics, theory of elasticity, model Naukova Dumka. 1981, 296. of a continuous medium with moment [2]. R.A. Kayumov, Fundamentals of the stresses. Proceedings of the Krylov State theory of elasticity and elements of the Scientific Center. 2016. Issue 93 (377). theory of plates and shells: Textbook / R.A. pp.163-174 Kayumov - Kazan: Kazan Publishing House. state architect-builds un-that. 2016, 111 [15]. N.G. Kolbasnikov, Theory of metal [3]. P.F. Nedorezov, N.M. Sirotkina, pressure treatment. Deformation resistance Numerical methods for studying steady-state and ductility. St. Petersburg: St. Petersburg vibrations of viscoelastic rectangular plates State University Publishing House, 2000, and circular cylindrical shells. Saratov. 314. Saratov University Press. [16]. E.A. Bulanov, Moment stresses in the 1997, 72 mechanics of a solid, bulk and liquid body. [4]. A.I. Koshelev, M.A. Narbut, Lectures M.: University book. 2012, 140 on the mechanics of deformable solids. SPb: [17]. R. Balescu. Equilibrium and Publishing House of St. Petersburg nonequilibrium statistic mechanics. A University. 2003, 276. Wiley-Intersciences Publication John [5].G.Z. Sharafufutdinov, Some plane Willey and Sons. New-Yourk- problems of the theory of elasticity. London…1975, 256 Moscow: Scientific World. 2014, [18]. P. Resibois, M.De. Lener. Classical [6]. A.M. Katz, Elasticity theory. St. kinetic theory of fluids. John Wiley and Petersburg. 2002, 208 Sons. New-York, London,… 1977, 423 E-ISSN: 2224-3429 7 Volume 16, 2021
WSEAS TRANSACTIONS on APPLIED and THEORETICAL MECHANICS DOI: 10.37394/232011.2021.16.1 Evelina Prozorova [19]. N.G. Van Kampen, Stochastic [22]. V. Ya. Rudyak, Statistical processes in physics and chemistry. North- aerohydromechanics of homogeneous and Holland. 1984, 376 heterogeneous media. Vol.2 Novosibirsk: [20]. Physics of Simple Liquids. Edited by NGASU. 2005, 468 H.N.V. Temperley, I.S. Rowlinson, G.S. [23]. Yu.V. Klimontovich, Statistical theory Rushbrooke. Amsterdam, 1068, 308. of open systems. Vol.2, M.: Janus. K 2019, [21]. L. Boltzmann. Selected Works. 438 Moscow: Nauka, 1984, 146 17 D.N. Zubarev, Nonequilibrium Statistical Creative Commons Attribution License 4.0 Thermodynamics, Moscow: Nauka, 1971, (Attribution 4.0 International, CC BY 4.0) 416 This article is published under the terms of the Creative Commons Attribution License 4.0 https://creativecommons.org/licenses/by/4.0/deed.en_US E-ISSN: 2224-3429 8 Volume 16, 2021
You can also read