_{1}

^{*}

A newly discovered Ballistic Principle of the Property Balance in the Space (BPPBS) occupied by the gas is introduced to simplify and reduce computations in applications dealing with modeling of fluid dynamics problems. The integro-differential balance equations for mass, momentum, and energy, which were formulated by applying the BPPBS, are derived. The integro-differential balance equations for mass and momentum were further approximated for the collision-dominated flow regime. Then they were reduced to the corresponding vector differential equations by the method of vector differentiation with subsequent elimination of the terms belonging to the original equation. It was shown that in the collision-dominated flow regime, the derived vector differential equations of mass and momentum balance are identical to the corresponding Navier-Stokes equations. This finding validates the BPPBS and suggests that, in the collision-dominated flow regime, the formulated integro-differential forms of the balance are exact implicit solutions for corresponding Navier-Stokes equations. Six additional tests demonstrating the feasibility of the proposed method and validity of the BPPBS are presented here. The BPPBS and the methodology discussed here will be highly useful not only as the basis to solve the fluid dynamics problems, but also to model any dynamic system composed of presumably chaotically moving particles/elements, each carrying a specific amount of property/information.

Computational Fluid Dynamics (CFD) is widely used in many practical applications ranging from basic hydrodynamic and kinematic to fundamental cosmological applications. The fundamental basis of any CFD tool is a solver, which solves the Navier-Stokes equations that are a set of second-order partial differential equations. Navier-Stokes equations are based on the assumption that the fluid is a continuum, and they are formulated by considering the mass, momentum, and energy conservations for a control volume of any size. The flow is considered continuous and differentiable, allowing the mass, momentum, and energy balances to be expressed as partial differential equations. Scientists made further approximations and simplifications to the Navier-Stokes equation set until it can be solved [

In the most general form, the Navier-Stokes equations of mass and momentum conservation for compressible fluid are expressed as:

∂ ∂ t ρ + ∇ ⋅ ( ρ u → ) = 0 (1)

and

∂ ∂ t u → + ( u → ⋅ ∇ ) u → + 1 ρ ∇ p − g → = 1 ρ ∇ ⋅ Τ , (2)

respectively, for t > 0 and r → ∈ ℝ 3 . In the equations above, ρ is the density of the fluid, u → is mass flow velocity, ∇ p is the pressure gradient, ∇ ⋅ Τ represents stresses inside the fluid, and g → is the external force per unit mass. The right-hand term ∇ ⋅ Τ has generally unknown functionality and contains too many unknowns and up to date does not apply to practical problems.

The general major disadvantage of any existing “mathematical” approach in solving the fluid dynamics problem is based on using an infinitesimal fluid element viewed as a continuous medium, to which fundamental physical principles are applied. This approach contradicts the molecular or particle nature, thus providing a source of significant uncertainty in interpreting the results of modeling and calculations. In the flows at small-scale, when the characteristic hydrodynamic length scale approaches the fluid characteristic length scale, the Navier-Stokes description is expected to fail (see [

1) The ideal gas composed of a large number of identical particles (molecules or inert atoms).

2) The volume occupied by gas particles themselves is negligible compared to the volume occupied by the gas.

3) The particles obey Newton’s laws of motion, and they move in random motion.

4) The particles experience forces only during collisions; any collisions are wholly elastic and take a negligible amount of time.

Mathematicians and physicists consider the recent advancement of the method of the lattice Boltzmann (LB) equation as a significant alternative to standard computational fluid dynamics [

In this paper, we discuss a physical approach to solving fluid dynamics problems by a novel analytical molecular dynamics technique (NAMDT), which was initially presented in our patent [

In the following, Section 2 describes the physical principles of the fluid model and the general physical approach of forming balancing equations. Here we also introduce newly discovered Ballistic Principle of the Property Balance in the Space occupied by the gas, applying of which is expected to simplify and reduce computations in applications dealing with modeling of fluid dynamics problems. The section also provides an analytical representation of the general integro-differential forms of mass balance, momentum balance, and energy balance equations.

In Section 3, we provide integro-differential forms of mass balance and momentum balance equations, which were adapted to the collision-dominated flow regime.

In Section 4, we demonstrate seven validation tests aiming to prove the feasibility of the proposed method. The first two validation tests, 4.1 and 4.2, illustrate that interaction of the model gas and the ideal gas with the gas-solid interface are identical (compare the rate of collision per unit area and the pressure exerted on the surface); nevertheless, some the most essential properties characterizing the model gas and the ideal gas are different. In the third validation test, 4.3, we supported our approach by formulating the mass-balance and the momentum-balance in the three-dimensional unlimited incompressible gas space at the uniform temperature. Analogously, in the fourth validation test, 4.4.2, we supported our approach by formulating the mass-balance and the momentum-balance in one-dimensional incompressible gas space confined between two parallel plates at the uniform temperature.

Also, in the validation tests 4.4 and 4.5, we explain a method for obtaining an analytical solution describing the incompressible model gas flow at the uniform temperature in the channel, which is confined in the space between two infinite parallel plates. Validation test 4.4 demonstrates the model gas flow with diffuse particle scattering from the plates. Validation test 4.5 demonstrates the model gas flow in a case of mixed diffuse and specular particle scatterings from the plates being at rest.

Lastly, in the validation test 4.6, we supported the proposed approach by revealing that, in the collision-dominated flow regime, the differential equations, which were converted from the integro-differential mass and momentum balance equations, are identical to the corresponding Navier-Stokes equations.

Finally, in Section 5, we present the conclusions and highlight the prospectiveness of the proposed physical approach for developing a new generation of the CFD software based on the NAMDT.

We suggest an approach in which fluid flow in a fluid system is model gas flow in a model gas system, which is equivalent to the fluid system. The transport processes involve the exchange of properties such as mass, momentum, and energy between interacting particles. In a more abstract sense, all particle interactions and property randomizing are an exchange of property/information [

In the following, all references are made to absolute time, which is measured equally within the model gas system. Besides, when referencing an appropriate law of motion, Newton’s Second Law of motion is considered. Also, in the interests of simplicity, we analyzed the model gas flow at a uniform body force such as a gravitational field of force or the acceleration field if dealing with a particle of a unit mass.

We assign these unique properties to the model gas [

1) The model gas enables a distant transport of one or more properties, including one or more of mass, momentum, and energy by particles being in a constant state of mostly random motion and interaction by collisions.

2) Each of the particles of the model gas is assigned to travel by obeying a ballistic trajectory that is governed by a law of motion in free space. It overcomes a distance between any of the two points of the ballistic trajectory with certain survival probability.

3) Each of the particles is adapted to transport a combination of one or more properties, comprising mass, momentum, and energy between a point of initial collision and a point of ending collision.

4) Each point within the space occupied by the model gas is treated as a point of collisions for converging particles, each following a ballistic trajectory with the same ending point simultaneously.

5) Each point of collisions is treated as either a point source for diverging ballistic particles or a point sink for converging ballistic particles.

6) Each of the particles moving from the point source to the point sink is treated as a property carrier. The property carrier is created in the point source during the initial collision by obtaining one or more properties of specific values being intrinsic to the model gas surrounding the point source. It is ended in the point sink during ending collision by transferring one or more properties of specific values in the point sink.

7) The value of the property, which is delivered in the point sink, or the value of the property, which is taken away from the point source, is evaluated regarding whether the value of property carried by each of the particles is modified because of interaction with an external field.

8) The velocity of a point source equals the mass flow velocity of the model gas flow in a corresponding point of the initial collisions at the time of the initial collision.

One can note from the above that the model gas properties differ from the properties typically assigned to the ideal gas (see above in Introduction).

In this paper, we investigate the transport of properties that are conserved during the ballistic traveling time.

The schematic diagram above shows the model gas composed of identical randomly moving particles and positioned in the observer’s Cartesian coordinate system 100. Note that in the paper, the observer’s coordinate system is designated by index “100.” Here we consider an isotropic model, which requires that the coordinate system needs to be at rest. To detect the position of the event, the observer reads a space coordinate at the location of an event. Also, the clocks at any location within a system are synchronized. The observer allocates space coordinates and time by recording both the space coordinates and time at the clock nearest the event position [

In the point source 102 positioned in point r → ′ at time t ′ i and moving with velocity u → ( t ′ i , r → ′ ) , particle 101, as a properties carrier, obtains one or more of properties being intrinsic to the model gas surrounding the point of the initial collision at the time of the initial collision. In the point source, the particle acquires a thermal velocity component v T ( t ′ i , r → ′ ) of an arbitrary direction relatively to the point source. For certainty, the magnitude of the average thermal velocity component in three-dimensional configuration can be defined as:

v T = 3 k B T m , (3)

where k B is Boltzmann constant, T is the temperature, and m is the mass of a particle.

In the microscopic scale, the model gas flow is characterized by the group of particles of mass m, which move randomly and interact by collisions with effective collision cross-section σ c . In each of the points in space at a given time, the particle density n, the magnitude of thermal velocity v T , and the vector of mass flow velocity u → quantify the model gas. In the interests of simplicity, unless otherwise stated, the particles are considered to have a unit mass, which, in the presence of external force, are accelerated during ballistic traveling with acceleration g → . We have recognized that each point in space occupied by the model gas may serve as both a sink and a collector of property delivered by converging ballistic particles from the entire model gas system and a source or a disperser into the surrounding of the property taken away by diverging ballistic particles.

Here, we reasonably may expect maintenance of a general property balance in each of the points of collisions within the model gas system. We formulate the balance, illustrated as a word equation in

B i n Ψ _FS ( r → , t ) = B o u t Ψ _FS ( r → , t ) + ∂ ∂ t [ n ( t , r → ) Ψ ( t , r → ) ] . (4)

For identification, we call the quantitative relationship above as the Ballistic Principle of the Property Balance in the Space (BPPBS) occupied by the particles in presumably chaotic motion. The BPPBS applies in general to any gas system, including the model gas systems containing heterogeneous gas-solid interfaces. Also, for clarity, we call our model as the Ballistic Model (BM).

This conceptual relationship can be expanded to the infinite space, for example, in a hypothetical system with no gravitational force. Straight-line trajectories of the particles may start from the infinity.

In the interests of simplicity, we concentrate our further analysis on the homogeneous model gas flow in the three-dimensional space having uniform gas properties on its periphery. Therefore, Equation (4) is reduced to:

B i n Ψ _F ( r → , t ) = B o u t Ψ _FS ( r → , t ) + ∂ ∂ t [ n ( t , r → ) Ψ ( t , r → ) ] , (5)

where B i n Ψ _F is the net rate of property influx per unit volume, which is formed by the converging ballistic particles from the surrounding model gas in the given non-moving point r → at the given time t. Still, the space of the model gas system may be separated from the infinite space by defining, for example, non-uniform gas pressure over the surface confining the system. This situation is further discussed when analyzing the momentum balance.

Here we admit the virtual nature of the balance described by Equation (4) or Equation (5). We consider that the value of property/information carried by a particle ejected from a given point in space at a given time is not a result of preceding physical interactions by collisions of all virtual converging ballistic particles capable of targeting with a certain probability the given point in space at the given time, but the result of the expectation of that value because of the cumulative effect from the surrounding space, in which each point of space complies with the BPPBS. Each given point in space at a given time is a point of reality (present) or a pivot point of consuming the results of events from the gas space that occurred in the past and sending the result of consumption and balancing from the present into the future. The value of the property/information at the given point in space at the given time, which is to be sent in the future, can be determined by solving the balance equations shown above. The analytical tools needed to formulate the balance according to Equation (5) are described in more detail below.

We have recognized and explained afterward that there exists a combination of a specific direction of an initial instant vector of thermal velocity v → T ( t ′ i , r → ′ ) and a vector of mass flow velocity u → ( t ′ i , r → ′ ) , which allows each of the selected particles to arrive in the given non-moving point r → at the given time, t [

Parameters | Short description |
---|---|

∇ = i → ∂ ∂ x + j → ∂ ∂ y + k → ∂ ∂ z | the operator of vector differentiation |

t | given time |

t ′ i | the time of the initial collision of the converging particle |

r → | position of the ending point of the converging particle |

r → ′ | position of the starting point of the converging particle |

u → ( t ′ i , r → ′ ) | mass flow velocity in the point r → ′ at time t ′ i |

v T ( t ′ i , r → ′ ) | the average magnitude of the thermal velocity of converging particle in point r → ′ at time t ′ i |

Z V ( t ′ i , r → ′ ) | the rate of collisions per unit volume in the point of the collision r → ′ at the time t ′ i of the initial collision |

v → ( t ′ i , r → ′ , t , r → ) | velocity vector in the ending point r → at the given time t |

Q i ( t , t ′ i ) | the probability of free path traveling along the ballistic trajectory of the converging ballistic trajectory starting at time t ′ i and ending at time t |

Ψ i n ( t ′ i , r → ′ , t , r → ) | property content delivered by the converging ballistic particle in the ending point r → at the given time t |

φ i = t − t ′ i | traveling time between an initial and ending consecutive collisions or the ballistic traveling time |

n | particles density |

m | particle mass |

σ c | the cross-section of collisions |

P c = σ c n | the number of particles placed within a collision tube of a unit length |

σ c | the cross-section of collisions |

V | the volume of integration over space occupied by the model gas |

r → i c ( r → ′ , t ′ i , t ˜ ) | the position of a virtual ballistic particle at a time t ˜ , which has zero magnitude of thermal velocity in the starting point r → ′ at the time of the initial collision t ′ i |

v → i c ( r → ′ , t ′ i , t ˜ ) | the velocity vector of the virtual ballistic particle having a zero component of the thermal velocity at a time t ˜ |

n → i = r → − r → i c ( r → ′ , t ′ i , t ) | r → − r → i c ( r → ′ , t ′ i , t ) | | instant unit vector directing thermal velocity component, so a traveling particle targets point r → at time t |

r ˜ → = r → ( r → ′ , t ′ i , t ˜ ) | the position vector of a ballistic particle at time t ˜ |

v ˜ → = v → ( t ′ i , r → ′ , t ˜ ) | the velocity vector of the ballistic particle at time t ˜ |

g → | an external force that applies to a particle of a unit mass |

v r e l | the average magnitude of the velocity of the traveling particle with respect to a nearby passed particle |

We define the net rate of property influx from the model gas in the general non-moving point r → at the given time t by these six steps:

Step 1: Identifying the converging ballistic trajectory and trajectory characteristics

Step 1 includes:

1) Formulating position vector r ˜ → = r → ( r → ′ , t ′ i , t ˜ ) of particle 304 on trajectory 301 and velocity vector v ˜ → = v → ( t ˜ ) (not shown) at time t ˜ by applying Equations (6) and (7), respectively, given below:

r ˜ → = r → ( r → ′ , t ′ i , t ˜ ) = v T ( t ′ i , r → ′ ) ( t ˜ − t ′ i ) n → i + r → i c ( r → ′ , t ′ i , t ˜ ) (6)

and

v ˜ → = v → ( r → ′ , t ′ i , t ˜ ) = v T ( t ′ i , r → ′ ) n → i + v → i c ( r → ′ , t ′ i , t ˜ ) , (7)

where t ≥ t ˜ ≥ t ′ i , n → i is a unit vector defined in

Specifically, when Newton’s Laws of Motion govern the motion of the particles, then r → i c ( r → ′ , t ′ i , t ˜ ) and v → i c ( r → ′ , t ′ i , t ˜ ) are defined by Equations (8) and (9) below:

r → i c ( r → ′ , t ′ i , t ˜ ) = r → ′ + u → ( t ′ i , r → ′ ) ( t ˜ − t ′ i ) + 1 2 g → ( t ˜ − t ′ i ) 2 (8)

and

v → i c ( r → ′ , t ′ i , t ˜ ) = u → ( t ′ i , r → ′ ) + g → ( t ˜ − t ′ i ) (9)

2) Determining the time of the initial collision t ′ i in point r → ′ . It can be done by solving Equation (6) in which t ˜ = t . Where a model gas system is governed by Newton’s Law of Motion, it can be done by resolving, for each of the ballistic particles, the equation of projectile motion, given by Equation (10), with respect to the ballistic traveling time, φ i :

1 2 g → φ i 2 + [ v T ( t ′ i , r → ′ ) n → i + u → ( t ′ i , r → ′ ) ] φ i + r → ′ − r → = 0 , (10)

which is obtained by substitution of Equation (8) in Equation (6) followed by the assignment of t ˜ = t and substitution of φ i defined in

t ′ i = t − φ i (11)

for each of the converging ballistic particles.

3) Defining an instant unit vector directing thermal velocity component of each particle in point r → ′ at time t ′ i by presenting Equation (7) in the following form:

φ i v T n → i = r → − r → i c (12)

where r → i c is defined as:

r → i c ( r → ′ , t ′ i , t ) = r → ′ + u → ( t ′ i , r → ′ ) φ i + 1 2 g → ( φ i ) 2 (13)

In this, vector r → i c is interpreted as the location of the center of the expansion zone at time t or the location of a particle having zero magnitude of an arbitrary or thermal velocity in point r → ′ at the time t ′ i of the divergence, which is observed at time t.

4) Defining the size of the expansion zone R i s p , by executing scalar multiplication of Equation (12) on itself resulted and averaging as:

φ i 2 v T 2 = ( r → − r → i c ) 2 (14)

and computing R i s p from Equation (14) as

R i s p = φ i v T = | r → − r → i c | (15)

The velocity vector v → i c of the center of the expansion zone 302 at time t is computed as

v → i c ( r → ′ , t ′ i , t ) = u → ( t ′ i , r → ′ ) + g → φ i . (16)

The velocity v → of a particle reaching the general non-moving point r → at time t on any point of the control surface 302 is computed as

v → ( r → ′ , t ′ i , t ) = v T ( t ′ i , r → ′ ) n → i + u → ( t ′ i , r → ′ ) + g → φ i (17)

which is obtained by assigning t ˜ = t and substitution of φ i = t − t ′ i in Equation (7) given above and rearrangement of the terms.

Step 2: Defining the probability of free path traveling along the ballistic trajectory from the starting point to the ending point

Step 2 includes:

1) Defining the average magnitude of the velocity at a particular point of a trajectory of the ballistic particle with respect to nearby passed particles at a particular point of a trajectory.

2) Expressing the probability of traveling along the ballistic trajectory in the three-dimensional configuration by Equation (18) given below:

Q i ( t , t ′ i ) = exp ( − ∫ t ′ i t P c ( r ˜ → ( t ˜ ) ) v r e l ( r ˜ → ( t ˜ ) ) d t ˜ ) , (18)

where t ˜ is a parametric time t ′ i < t ˜ ≤ t , r ˜ → ( t ˜ ) is a point on the ballistic trajectory at the parametric time t ˜ , v r e l ( r ˜ → ( t ˜ ) ) is the average magnitude of the velocity with respect to a nearby passed particle in the trajectory point r ˜ → ( t ˜ ) , and P c ( r ˜ → ( t ˜ ) ) is the number of particles placed within a collision tube of a unit length in the trajectory point r ˜ → ( t ˜ ) .

Following sub-steps calculate the average magnitude of the velocity v r e l ( r ˜ → ( t ˜ ) ) of the ballistic particle with respect to nearby passed particles at a particular point of a trajectory B at a specified time t ˜ :

1) by defining an instant magnitude of the velocity of the converging ballistic particle in the trajectory point with respect to nearby particles in the trajectory point r ˜ → at time t ˜ as

v → r l ( r ˜ → ( t ˜ ) ) = v → 1 ( t ′ i , r → ′ , t ˜ ) − v → 2 ( r ˜ → ( t ˜ ) ) , (19)

which is formed by connecting the end of the instant velocity vector v → 2 ( t ˜ , r ˜ → ) with any point S on spherical surface 403 of radius v T ( r ˜ → ( t ˜ ) ) , where

v → 1 ( t ′ i , r → ′ , t ˜ ) = v T ( t ′ i , r → ′ ) n → 1 + u → ( t ′ i , r → ′ ) + g → ( t ˜ − t ′ i ) , (20)

where n → 1 is a unit vector having the point of origin r → ′ , v T ( t ′ i , r → ′ ) and u → ( t ′ i , r → ′ ) are thermal velocity and mass flow velocity components in the rest frame of the model gas in point r → ′ at time t ′ i , which are acquired by particle P1 because of a collision in this point, and where

v → 2 ( r ˜ → ( t ˜ ) ) = v T ( r ˜ → ( t ˜ ) ) n → i + u → ( r ˜ → ( t ˜ ) ) , (21)

where v T ( r ˜ → ( t ˜ ) ) and u → ( r ˜ → ( t ˜ ) ) are thermal velocity and mass flow velocity components in the rest frame of the model gas in point r ˜ → at time t ˜ , which are acquired by particle P2 because of a collision in this point, and

2) by averaging the instant magnitude of the velocity overall directions of the thermal velocity component of one of the nearby particles in the trajectory point r ˜ → at time t ˜ . This is done by integrating the instant magnitude of the relative velocity of Equation (19) over the angle of ϑ from 0 to π and φ , which is the angle of rotation around axis OP from 0 to 2 π , and by normalizing by the solid angle of 4 π , which results in:

v r e l ( r ˜ → ( t ˜ ) ) = 1 2 ∫ 0 π | v → r ( r ˜ → ( t ˜ ) ) | 2 + [ v T ( r ˜ → ( t ˜ ) ) ] 2 − 2 v T ( r ˜ → ( t ˜ ) ) | v → r ( r ˜ → ( t ˜ ) ) | cos ( ϑ ) sin ( ϑ ) d ϑ , (22)

where

v → r ( r ˜ → ( t ˜ ) ) = v → 1 ( t ′ i , r → ′ , t ˜ ) − u → ( r ˜ → ( t ˜ ) ) . (23)

Note that, typically, the magnitude of the relative mass flow velocity or the mass flow velocity component of the passing particle P1 with respect to nearby passed particle P2 is insignificant compared to the magnitude of the thermal velocity of either passing particle P1 or nearby passed particle P2 or of both.

The average magnitude of the velocity of the traveling particle with respect to a nearby passed particle is calculated from Equation (24) given below, which is obtained by substitution of | v → r | ≅ v T ( t ′ i , r → ′ ) in Equation (22) and executing the integration of the resulted equation:

v r e l ( r ˜ → ( t ˜ ) ) = 1 6 v T ( t ′ i , r → ′ ) v T ( t ˜ , r ˜ → ) { ( v T ( t ′ i , r → ′ ) + v T ( r ˜ → ( t ˜ ) ) ) 3 − | v T ( t ′ i , r → ′ ) − v T ( r ˜ → ( t ˜ ) ) | 3 } , (24)

where v r e l ( r ˜ → ( t ˜ ) ) is the average magnitude of the velocity with respect to a nearby passed particle in the trajectory point, v T ( t ′ i , r → ′ ) is the magnitude of the thermal velocity in the starting point of the ballistic trajectory, and v T ( r ˜ → ( t ˜ ) ) is the thermal velocity of a nearby passed particle in the trajectory point.

Also, typically, the magnitudes of the thermal velocity of nearby particles are approximately identical. For non-relativistic particles, the average magnitude of the velocity with respect to each particle moving in an arbitrary direction is calculated from Equation (25) given below, which is obtained by substitution of v T ( r ˜ → ( t ˜ ) ) = v T ( t ′ i , r → ′ ) = v T in Equation (24):

v r e l = 4 3 v T . (25)

Analogously, for relativistic particles, the average magnitude of the velocity with respect to each particle moving in an arbitrary direction is calculated from Equation (25) given below:

v r e l = v T 2 ∫ 0 π 2 − 2 cos ( ϑ ) − [ 1 − cos 2 ( ϑ ) ] v T 2 / c 2 1 − v T 2 cos ( ϑ ) / c 2 sin ( ϑ ) d ϑ (26)

where c is the speed of light. For v T c ≪ 1 , integrating Equation (26) will yield Equation (25). For v T c ≅ 1 , integrating Equation (26) will yield:

v r e l ≅ c . (27)

Step 3: Defining the net rate of particle efflux per unit volume from a point source positioned in a point of the initial collisions and moving with the mass flow velocity of the model gas in that point.

Step 3 includes the following sub-steps:

1) Defining the particle flux J → r ′ N along the ballistic trajectory in a point of the space r → surrounding the point of the initial collision, the step that includes representing J → r ′ N by applying Equation (28) as given below:

J → r ′ N = 1 2 n ( t ′ i , r → ′ ) Q i ( t , t ′ i ) v → ( t ′ i , r → ′ , t , r → ) , (28)

where v → is defined by Equation (17), Q i is a survival probability defined by Equation (18), and n ( t ′ i , r → ′ ) is particle density at a specific point r → ′ at time t ′ i .

2) Representing, in a coordinate system associated with the point of origin r → i c that moves with velocity v → i c , the vector field of the particle flux J → r ′ N through the in the control surface 302 of

J → C S N = 1 2 n ( t ′ i , r → ′ ) Q i ( t , t ′ i ) [ v → ( t ′ i , r → ′ , t , r → ) − v → i c ] = 1 2 n ( t ′ i , r → ′ ) Q i ( t , t ′ i ) v T ( t ′ i , r → ′ ) n → i . (29)

3) Applying and executing the divergence operator ∇ ′ ⋅ to the vector field of Equation (29) followed by shrinking the volume of the auxiliary control volume to infinitely small volume, i.e., r → → r → ′ , which, in formula form, is expressed as:

Z V ( t ′ i , r → ′ ) = { ∇ ′ ⋅ [ J → C S N ] } r → → r → ′ = 1 2 { ∇ ′ ⋅ [ n ( t ′ i , r → ′ ) Q i ( t , t ′ i ) v T ( t ′ i , r → ′ ) n → i ] } r → → r → ′ , (30)

where

∇ ′ = i → ∂ ∂ x ′ + j → ∂ ∂ y ′ + k → ∂ ∂ z ′ , (31)

and includes representing the particle flux production rate, or the net rate of particle efflux per unit volume, or the rate of collisions per unit volume Z V ( t ′ i , r → ′ ) , in a point of the initial collision moving with the mass flow velocity u ( t ′ i , r → ′ ) at time t ′ i by following Equation (32) given below:

Z V ( t ′ i , r → ′ ) = 1 2 n ( t ′ i , r → ′ ) P c ( t ′ i , r → ′ ) v r e l ( t ′ i , r → ′ ) , (32)

where n ( t ′ i , r → ′ ) is particle density, P c ( t ′ i , r → ′ ) is the number of particles placed within a collision tube of a unit length in the corresponding point of the initial collisions at the time of the initial collision, and v r e l ( t ′ i , r → ′ ) is the average magnitude of the velocity with respect to a nearby passed particle in the corresponding point of initial collisions at the time of the initial collision. The above is obtained upon acknowledgment that the control volume 305, which is confined by inflated control surface 302, is isolated (see

Step 4: Defining property flux in a given non-moving point at a given time from one of the point sources of the model gas.

Step 4 includes representing the property vector flux J → r → ′ → r → Ψ originated from the point source of the initial collisions in point r → ′ at time t ′ i , which moves in the space of the model gas with a mass-flow velocity u → ( t ′ i , r → ′ ) , and being detected by a point sink positioned in a being at the rest point r → at the given time t:

J → r → ′ → r → Ψ ( t , r → ) = 1 4 π 1 | r → − r → i c | 2 Q i ( t , t ′ i ) Z V ( t ′ i , r → ′ ) v → ( t ′ i , r → ′ , t , r → ) v T ( t ′ i , r → ′ ) Ψ i n ( t ′ i , r → ′ , t , r → ) d V ′ . (33)

Step 5: Defining the rate of the property vector flux J → F S → r → Ψ in point r → at the given time t, which is originated from initial collisions within entire space occupied by the model gas.

Step 5 includes applying Equation (34) given below, which is obtained by integrating Equation (33) over the volume of the model gas system:

J → F S → r → Ψ = ∭ V 1 4 π 1 | r → − r → i c | 2 Q i ( t , t ′ i ) Z V ( t ′ i , r → ′ ) v → ( t ′ i , r → ′ , t , r → ) v T ( t ′ i , r → ′ ) Ψ i n ( t ′ i , r → ′ , t , r → ) d V ′ . (34)

here point r → is excluded from integration in the equation above because we are interested in calculating the total rate of the property flux in the point sink at r → , which is originated from the surrounding point sources of the initial collisions at r → ′ .

Step 6: Defining the net rate of property influx per unit volume B i n Ψ _FS formed by the flow of ballistic particles and converging from the gas space in the general non-moving point at the given time.

Step 6 further includes a step of representing B i n Ψ _FS by applying Equation (35) given below:

B i n Ψ _FS ( r → , t ) = − ∇ ⋅ ∭ V 1 4 π 1 | r → − r → i c | 2 Q i ( t , t ′ i ) Z V ( t ′ i , r → ′ ) v → ( t ′ i , r → ′ , t , r → ) v T ( t ′ i , r → ′ ) Ψ i n ( t ′ i , r → ′ , t , r → ) d V ′ . (35)

From the equation above, one can conclude that the net rate of property influx in the general non-moving point r → at the given time t is resulted from the impact of the flow of ballistic particles converging from the gas space. The equation above calculates the impact in point r → at the given time t from all initial collisions of converging ballistic particles having trajectories allowing them to target point r → at the given time t, the initial collisions taking place during all preceding dynamic history of the system preceding the given time, i.e. t ′ i < t .

To define the net rate of property efflux per unit volume into surroundings from the general non-moving point r → at the given time, t, the linear dimensions of the main control volume surrounding point r → are selected to be sufficiently small for preventing two and more consecutive collisions of the same particle within the main control volume. We anticipate that the net rate of property efflux per unit volume is formed by diverging particles. Each of the diverging particles is selected from all available particles by the ballistic trajectory having the starting point in the given non-moving point at the given time. The table of the model parameters associated with defining the net rate of total property efflux per unit volume is presented in

We define the net rate of property efflux from the general non-moving point r → at the given time t by these four steps:

Step 1: Identifying a trajectory and trajectory characteristics, for each particle diverging from the general non-moving point at the given time.

Parameters | Short description |
---|---|

t | given time |

r → | position of the starting point of a ballistic trajectory of the diverging particle |

r → ′ | position of the ending point of a ballistic trajectory of the diverging particle |

t ′ a | time of positioning the ending point of the diverging particle |

n ( t , r → ) | particle density in the starting point r → at the given time t |

u → ( r → , t ) | mass flow velocity in point r → at time t |

v T ( r → , t ) | the average magnitude of the thermal velocity of a diverging particle in point r → at time t |

v → + ( t , r → , t ′ a ) | the velocity vector of the diverging particle at the time of positioning in point r → ′ |

r → a c ( t , r → , t ˜ ) | the position of a virtual ballistic particle at time t ˜ , which has zero magnitude of thermal velocity in the starting point r → at the time of initial collision t |

v → a c ( t , r → , t ˜ ) | the velocity vector of the virtual ballistic particle at time t ˜ |

Q + ( t ′ a , t ) | the probability of free path traveling along the ballistic trajectory starting at time t and ending at time t ′ a |

Ψ ( t , r → , t ′ a , r → ′ ) | property content carried by the diverging particle at the time t ′ a of crossing in point r → ′ enclosing control surface 503 |

Step 1 includes:

1) Formulating position vector r ˜ → ′ (506) of the particle on trajectory 505 and velocity vector v ˜ → ′ (not shown) at time t ˜ by applying Equations (36) and (37), respectively, given below:

r ˜ → ′ = r → ′ ( t , r → , t ˜ ) = v T ( t , r → ) n → + c ( t ˜ − t ) + r → a c ( t , r → , t ˜ ) , (36)

and

v ˜ → ′ + ( t ˜ ) = v → ′ + ( t , r → , t ˜ ) = v T ( t , r → ) n → + c + v → a c ( t , r → , t ˜ ) , (37)

where n → + c is a unit vector with the initial point of origin r → at time t and t ′ a ≥ t ˜ ≥ t . In Equations (36) and (37) above, r → a c ( t , r → , t ˜ ) and v → a c ( t , r → , t ˜ ) are defined by an appropriate law of motion.

Specifically, when Newton’s Law of Motion governs the model gas system, r → a c ( t , r → , t ˜ ) and v → a c ( t , r → , t ˜ ) are defined by Equations (38) and (39), respectively, given below:

r → a c ( t , r → , t ˜ ) = r → + u → ( t , r → ) ( t ˜ − t ) + 1 2 g → ( t ˜ − t ) 2 (38)

and

v → a c ( t , r → , t ˜ ) = u → ( t , r → ) + g → ( t ˜ − t ) . (39)

2) Determining the time needed, for each particle diverging from a non-moving point r → at the given time t to cross a control surface enclosing point r → in point r → ′ of the control surface. It can be done by solving Equation (36) in which t ˜ = t ′ a . Specifically, when Newton’s Laws of Motion govern the motion of the particles, then it can be done by resolving the equation of projectile motion, Equation (40) given below, with respect to the ballistic traveling time φ +

r → ′ = r → + [ v T ( t , r → ) n → + c + u → ( t , r → ) ] φ + + 1 2 g → φ + 2 (40)

which is obtained by assigning t ˜ = t ′ a and substitution of φ + = t ′ a − t in Equation (36) and rearrangement of the terms. Traveling time φ + is needed to indicate a time t ′ a at which the departing particle has reached a point r → ′ on control surface 503 enclosing point r → .

3) Defining an instant unit vector directing thermal velocity component along the diverging ballistic trajectory includes:

presenting Equation (40) in the following form:

φ + v T n → + c = r → ′ − r → a c , (41)

where n → + c is a unit vector and r → a c (not shown) is the location of a virtual ballistic particle leaving point r → at time t, which would have zero magnitude of the thermal velocity and is observed at time t ′ a . In a case, where a model gas system is governed by Newton’s Laws of Motion, r → a c is expressed as:

r → a c = r → + u → ( t , r → ) φ + + 1 2 g → φ + 2 (42)

and deriving n → + c from Equation (41), which is given as

n → + c = r → ′ − r → a c | r → ′ − r → a c | , (43)

where n → + c is a unit vector of origin r → a c at time t ′ a , directing the vector of the thermal velocity component and where t ′ a is the time of positioning in point r → ′ , which is on the control surface 503.

4) Defining the velocity vector v → + , (507) of each particle at the moment of crossing the control surface in point r → ′ includes, where a model gas system is governed by Newton’s Laws of Motion, representing v → + by applying Equation (44) as given below:

v → + ( t , r → , t ′ a , r → ′ ) = v T ( t , r → ) n → + c + u → ( t , r → ) + g → φ + , (44)

which is obtained by assigning t ˜ = t ′ a and substitution of t ′ a − t = φ + in Equation (37).

Step 2: Defining the probability of traveling along the ballistic trajectory from the general non-moving point r → at time t to one of the points in space surrounding the general non-moving point.

Step 2 includes representing the probability Q + ( t ′ a , t ) by applying Equation (45) given below:

Q + ( t ′ a , t ) = Q + ( 0 , φ i + ) = exp ( − ∫ t t ′ a P c ( r ˜ → ′ ( t ˜ ) ) v r e l ( r ˜ → ′ ( t ˜ ) ) d t ˜ ) , (45)

where t ′ a = t + φ + is the time of the particle positioning in point r → ′ , P c ( r ˜ → ′ ( t ˜ ) ) is the number of particles placed within a collision tube of a unit length in the trajectory point r ˜ → ′ ( t ˜ ) , v r e l ( r ˜ → ′ ( t ˜ ) ) is an average magnitude of the velocity with respect to a nearby passed particle in the trajectory point r ˜ → ′ ( t ˜ ) .

Step 3: Defining the vector/tensor field of property flux J r Ψ along with one of the ballistic trajectories of a diverging particle in point r → ′ .

Step 3 includes representing the property vector flux J r Ψ by applying Equation (46) as given below:

J r Ψ = 1 2 n ( t , r → ) Q + ( t ′ a , t ) v → + ( t , r → , t ′ a , r → ′ ) Ψ ( t , r → , t ′ a , r → ′ ) . (46)

Step 4: Defining the net rate of property efflux per unit volume B o u t Ψ _FS from the general non-moving point r → at the given time t.

Step 4 includes representing B o u t Ψ _FS by applying Equation (47) as given below:

B o u t Ψ _FS ( r → , t ) = 1 2 { ∇ ⋅ [ n ( t , r → ) Q + ( t ′ a , t ) v → + ( t , r → , t ′ a , r → ′ ) Ψ ( t , r → , t ′ a , r → ′ ) ] } r → ′ → r → , (47)

which is obtained by executing the divergence operator ∇ ⋅ to the vector field of Equation (46) and by shrinking the volume of control volume 508 confined by the control surface 503 to an infinitely small volume, i.e., r → ′ → r → , which also leads to the limit r → a c → r → .

The integro-differential form of property balance equation is formulated by Equation (48) given below, which is obtained by substitution of Equations (35) and (47) in Equation (5):

∂ ∂ t [ n ( t , r → ) Ψ ( t , r → ) ] + 1 2 { ∇ ⋅ [ n ( t , r → ) Q + ( t ′ a , t ) v → + ( t , r → , t ′ a , r → ′ ) Ψ ( t , r → , t ′ a , r → ′ ) ] } r → ′ → r → = − ∇ ⋅ ∭ V 1 4 π 1 | r → − r → i c | 2 Q i ( t , t ′ i ) Z V ( t ′ i , r → ′ ) v → ( t ′ i , r → ′ , t , r → ) v T ( t ′ i , r → ′ ) Ψ i n ( t ′ i , r → ′ , t , r → ) d V ′ , (48)

here point r → is not included in integration for converging ballistic particles. It implies that any singularity in the right-hand of the equation above is excluded, so Equation (48) defines Ψ at t ≥ 0 as an implicit function of r → on ℝ 3 , i.e. r → ∈ ℝ 3 .

The analytical representation of the general integro-differential form of property balance equation shown by Equation (48) is generally valid for any homogeneous fluid system with any configuration of the external field of force.

Remark that the integro-differential property balance equation needs to be formed for each unknown property/variable so that the number of equations in a system of balance equations is sufficient to determine each of the unknown properties characterizing the model gas flow. In the following, we provide general governing integro-differential forms of mass balance, momentum balance, and energy balance equations.

In continuation of the discussion at the end of Section 2.2, we would like to highlight that the equation above is a general symbolic representation of the rule that shall be obeyed at any given time t in any given point r → of space occupied by presumably chaotically moving particles experiencing random collisions. Equation (48) also shows that the balance at time t in any given point r → is formed by the exhaustive combination of converging ballistic particles from the surrounding, which can target with a certain probability the given point in space r → at the given time t. The converging in point r → at time t ballistic particles are originated from preceding collisions at times t ′ i < t (past). Whereas the diverging ballistic particles originated from collisions at time t (present) transport the balanced property/information into surrounding toward the future t ′ a > t .

To formulate a general integro-differential form of mass balance equation in a given non-moving point of space occupied by the model at a given time, we will modify Equation (48) by assigning:

Ψ = Ψ i n = 1. (49)

Then, we obtain the following general integro-differential form of the mass balance equation

∂ ∂ t [ n ( t , r → ) ] + 1 2 { ∇ ⋅ [ n ( t , r → ) Q + ( t ′ a , t ) v → + ( t , r → , t ′ a , r → ′ ) ] } r → ′ → r → = − ∇ ⋅ ∭ V 1 4 π 1 | r → − r → i c | 2 Q i ( t , t ′ i ) Z V ( t ′ i , r → ′ ) v → ( t ′ i , r → ′ , t , r → ) v T ( t ′ i , r → ′ ) d V ′ . (50)

To formulate a general integro-differential form of momentum balance equation in a given non-moving point of space occupied by the model at a given time, we will modify Equation (48) by assigning:

Ψ ( t , r → , t ′ a , r → ′ ) = v → + ( t , r → , t ′ a , r → ′ ) (51)

in the left-hand of the equation and

Ψ i n ( t ′ i , r → ′ , t , r → ) = v → ( t ′ i , r → ′ , t , r → ) (52)

in the right-hand of the equation. Then we obtain:

{ m ∂ ∂ t [ n ( t , r → ) v → + ( t , r → , t ′ a , r → ′ ) ] } r → ′ → r → + 1 2 { ∇ ⋅ [ m n ( t , r → ) Q + ( t ′ a , t ) v → + ( t , r → , t ′ a , r → ′ ) v → + ( t , r → , t ′ a , r → ′ ) ] } r → ' → r → + ∇ p ( t , r → ) = − m ∇ ⋅ ∭ V 1 4 π 1 | r → − r → i c | 2 Q i ( t , t ′ i ) Z V ( t ′ i , r → ′ ) v → ( t ′ i , r → ′ , t , r → ) v → + ( t , r → , t ′ a , r → ′ ) v T ( t ′ i , r → ′ ) d V ′ . (53)

In the equation above, we also introduce a term of the pressure force exerted on its surroundings in point r → at time t, which may appear because of a non-uniform pressure applied to a bounded system.

To formulate a general integro-differential form of energy balance equation in a given non-moving point of space occupied by the model at a given time, we will modify Equation (48) by assigning:

in the left-hand of the equation and

in the right-hand of the equation. Then we obtain:

In the second left-term of the equation above, we also introduced a term of the pressure force work on its surroundings in point

We have recognized that, in the majority of real conditions on Earth and near-Earth space, first, the acceleration field

second, for each of the model gas particles, the average magnitude of the thermal velocity component,

third, in the collision-dominated flow regime, the magnitude of velocity gained or lost because of interaction with the acceleration field during the ballistic traveling time

fourth, in the collision-dominated flow regime, the relative change of any property value or any parameter characterizing the model gas

and fifth, in the collision-dominated flow regime, the relative change of any property value or any parameter characterizing the model gas is insignificant on the length scale of the average distance between the model gas particles

We also note here that since point

1) Approximating

where

and

2) Approximating

where

3) For

4) Approximating

where

and

5) Approximating

where

6) Approximating

7) Approximating

8) Representing

which has the following property:

Considering that governing equations require satisfaction of the balance of any model gas property at any time, these analytical representations of approximations for mass balance, momentum balance, and energy balance are provided.

Here the mass balance equation is formed by considering a unique combination of ballistic particles converging from the entire model gas system in a given point at a given time and the diverging ballistic particles from the given point at the given time. Besides, each of the converging ballistic particles can target point

To formulate an integro-differential form of mass balance equation in the collision-dominated flow regime, which applies to a given non-moving point of space occupied by the model at a given time, we will modify Equation (50) by following steps of:

1) substituting approximations from Equation (62) to Equation (73) in Equation (50) and executing vector differentiation;

2) executing limit

3) remembering that point

4) neglecting terms containing

Upon executing the above, Equation (50) is reduced to:

where

Equation (76) defines function

Here we again should note that the momentum balance equation is formed by considering a unique combination of ballistic particles converging from the entire model gas system in a given point at a given time and the diverging ballistic particles from the given point at the given time. Besides, each of the converging ballistic particles can target point

To formulate an integro-differential form of momentum balance equation in a given non-moving point of space occupied by the model at a given time, we will modify Equation (53) as follows:

1) assigning

and

2) substituting approximations for functions involving in Equation (48) by equations from Equation (62) to Equation (73) and executing vector differentiation;

3) executing limit

4) remembering that point

5) neglecting terms containing

Upon executing the above, Equation (56) is reduced to:

It implies that Equation (80) defines

The following validation tests are aimed to demonstrate the feasibility of the proposed method.

Adopting Equation (34) to the conditions above and assigning

where

and

where n is particle density,

Finally, using the geometry illustrated in

and

and substituting Equations (82), (83), (84), (85), and (86) in Equation (81), in which

That the result of derivation above according to our method applied to the model gas is identical to the result of the derivation of the rate of collisions per unit area of an ideal gas, which one can find in any course of the kinetic theory of gases, supports the Ballistic Model.

Referring to the previous section and

and the momentum component delivered by the particle to the surface

where m is mass of the particle/molecule.

If the particle undergoes an elastic collision with the surface, in such a collision, the momentum passed on the surface is steady-state. Adopting Equation (34) to the conditions listed initially referring to

Now, using the geometry illustrated in

Again, that the result of derivation above according to our method, which considers an impact on the surface of unlimited number particles of the model gas (integration along y direction from zero to infinity), is identical to the result of the derivation of the pressure of an ideal gas according to the kinetic theory of gases, corroborates with the Ballistic Model.

Besides, analyzing validation tests 4.1 and 4.2 above, one may note that the proposed method provides the possibility of quantifying impacts of the limited number of particles on the gas-solid interface. This, from the practical viewpoint, is important in many applications dealing with MEMS technology (capacitive sensing, electrostatic actuation mechanisms, a block of sensing mass in micro gyroscopes, accelerometers, switches, mirrors, pressure sensors and so on) [

Recognizing that in a steady-state model incompressible gas flow at the uniform temperature with no gravitational force, all variables describing flow depend on the position in space and are not dependent on time, so that the particle density n, the mass m, and the magnitude of the thermal velocity of the particle/molecule,

and the momentum-balance equation by adopting Equations (80) as

We need to verify these two equalities above by integration of the right-hands of Equations (92) and (93).

Placing the basis of the coordinate system in point

and substituting Equations (83) and (94) in Equation (92), in which

Analogously, after substitution of Equations (74), (83), and (94) in Equation (93) and taking into account that point

here

That the results of integration provided by Equations (95) and (96) are identical to the left-hands of Equations (92) and (93), respectively, evidently supports the validity of the BPPBS in three-dimensional infinite space.

Note: A short communication about these two validation tests below is present on the website and can be accessed on [

The velocity profile generated in the model gas due to the pressure gradient along the channel is analyzed by an analytical method based on the proposed model gas flow described above.

where u is the gas slip velocity near the wall,

The method described below uses none of the Maxwell-type assumptions.

Here we assign the following additional unique properties to the model gas being in contact with a gas-solid interface revealing mixed diffuse and specular scattering of particles [

1) Each collision on a gas-solid interface of the model gas system, which has resulted in the scattering of the diffuse particles from the gas-solid interface, is treated as an act of interaction involving a property transport from the gas-solid interface to the scattered particle.

2) Each point of the diffuse particle scattering on the gas-solid interface is treated as a heterogeneous point source for each of the scattered particles.

3) The velocity of each of the heterogeneous point sources on the gas-solid interface equals the velocity of the gas-solid interface of corresponding points of diffuse particle scattering.

4) The point source strength of the heterogeneous point sources on the gas-solid interface is directly proportional to a property accommodation coefficient

Note: Diffuse scattering is an act of interaction involving property transfer from a gas-solid interface to a scattered particle. Specular scattering does not involve property exchange between the gas-solid interface and a scattered particle.

In the model gas system confined by gas-solid interfaces with the purely diffuse scattering of particles, each particle initiated from the preceding diffuse scatterings from the interfaces delivers in the CV some property obtained from the location of the initial diffuse scattering (ballistic trajectory 703 from plate 2 and ballistic trajectory 704 from plate 1 of

Recognizing that in a steady-state model incompressible gas flow at the uniform temperature, all variables describing flow depend on the position in space and are not dependent on time so that the particle density n, the mass m and the magnitude of the thermal velocity of the particle/molecule,

for mass-balance, and

for y-momentum-balance. Here

and

The left-hand of Equation (99) has zero value because of the average momentum of the diverging particle, which is measured by the instant momentum

Solving Equation (98) will cause finding that

Substituting Equation (102) in Equation (99), we may verify that the mass and momentum balance in any point of space between plate 1 and plate 2 along axis are concerved, which suggests the valididty of the BPPBS in the one-dimensional configuration.

Step 1: Formulating an integral form of the

Considering the above, we obtain the following integral form of the

Step 2: Obtaining a differential form matching to the corresponding integral form of the

Here we use the method of differentiation (twice). Each step of the differentiation is followed by the step of subsequent elimination of the integral terms by using the original equation, which is given as:

Step 3: Obtaining a general solution having arbitrary coefficients.

Integrating twice the equation above, we obtain:

where

and B and C are the arbitrary coefficients.

Step 4: Determining the values of each of the arbitrary coefficients.

Introducing Equation (105) in Equation (103) (for a specific number of points (two) within the model gas system), we determine the values of each of the arbitrary coefficients. For certainty, we selected points

Step 5: Solving a system of the obtained in Step 4 algebraic equations.

We obtained the following functional relationship of fluid velocity

and

where Plate 1 and Plate 2 are at rest, i.e.

where

The equation above reveals there is a minimum in the normalized mass flow rate (at about

where corresponding slip coefficients

In

We compared the non-dimensional flow rates calculated by Equation (110) for the Ballistic Model and a selection of slip models proposed by various authors [

Here we assume the symmetricity of the expected solution (because of the zero velocity of the confining plates). We also consider that the influx of

where

where

Finally, substituting

Remarkably, the derived tangential slip velocity coefficient C does contain the term being proportional to

Substitution in Equation (105) of Equations (106), (114), and (115) will yield:

That the tangential slip velocity coefficient C is analytically defined from the continuum through the slip and transition to free-molecule flow regimes gave us confidence that our approach is valid. Again, our method uses none of the Maxwell-type assumptions. Still, coefficient C contains the terms proportional to

According to Tenenbaum and Pollard [

Definition 3.6 [

on an interval I:

1) it defines y as an implicit function of x on I, i.e., if there exists a function

2)

for every x on I.

The standard procedure in calculus to prove that an implicit function defined by relationship

To reduce an equation to an ordinary differential equation, we apply the method of differentiation for integral equations (ones, twice, and so on) with subsequent elimination of the terms belonging to the original equation [

Deriving a vector differential form of mass balance equation is shown:

Step 1. Normalizing Equation (76) by

Step 2. Applying divergence operator

Step 3. Comparing the first three right-hand terms in the equation above and Equation (76), we may suggest a modification of Equation (76) by introducing a coefficient of proportionality to the mass,

Substituting the first three right-hand terms of Equation (120) by the left-hand terms of Equation (121) and rearranging terms, we finally obtain:

Rearranging terms in the equation above, we derived the following reduced vector-differential form of the mass balance equation:

In the collision-dominated flow regime, the relative change of any property value or any parameter characterizing the model gas is insignificant on the length scale of the average distance between the gas particles or on the time scale of traveling time between consecutive collisions. Therefore, for a high frequency of collisions quantified by

which is identical to Equation (1). Since the vector differentiation of Equation (76) yields Equation (123), which eventually is reduced to Equation (1) and Equation (76) defines

Deriving a vector differential form of mass balance equation is shown:

Step 1. Normalizing Equation (80) by

Step 2. Applying divergence operator

Step 3. Comparing the first four right-hand terms in the equation above and Equation (76), we may suggest a modification of Equation (80) by introducing a coefficient of proportionality to the mass,

Substituting the first five right-hand terms of Equation (126) by the left-hand terms of Equation (127) and rearranging terms, we finally obtain:

Considering the conservation of mass by subtracting Equation (122) multiplied by

Comparing the first, the second, the third, and the fourth left-hand terms of Equations (2) with the first, the third, the fourth, and the fifth left-hand terms of Equation (129), respectively, suggests that each pair of the corresponding compared terms is identical. The second left-hand term of the equation above is identified as a normalized self-diffusion force. Additional comparison of the right-hand terms of Equations (2) and (129) may suggest that stress inside the fluid has the following functional dependence:

Since the vector differentiation of Equation (80) yields Equation (129), which resembles Equation (2) and Equation (80) defines

1) Modeling of fluid dynamics problems by the NAMDT is based on the recognition that each particle composing the model gas travels with a probability between any of two points in space occupied by the model gas while following a ballistic trajectory governed by a law of motion in free space. Each ballistic particle is treated as a property carrier transporting one or more of mass, momentum, and energy between the points of consecutive collisions and each point in space occupied by the model gas is both a sink accumulating property delivered by converging ballistic particles from the entire model gas system and a source dispersing property by diverging ballistic particles.

2) Based on the proposed model gas properties, we formulated the Ballistic Principle of the Property Balance in the Space occupied by the gas, application of which may simplify and reduce computations in applications dealing with modeling of fluid dynamics problems.

3) Following the above principles, a general integro-differential form of the property balance equation is proposed. The general integro-differential form of the property balance is further modified to derive the integro-differential forms of mass balance, momentum balance, and energy balance.

4) The following two direct tests validate the BPPBS:

a) demonstration of the mass-balance and momentum-balance conservation in a given non-moving point in three-dimensional unlimited incompressible gas space with a lack of gravitational force at the uniform temperature by analytical verification of the balance between the rate of mass and momentum influx and the rate of mass and momentum efflux, respectively, in the given point of the gas space;

b) demonstration of the mass-balance and momentum-balance conservation in a given non-moving point in one-dimensional incompressible gas space confined between two parallel plates at the uniform temperature by analytical verification of the balance between the rate of mass and momentum influx and the rate of mass and momentum efflux, respectively, in the given point of the gas space.

5) The BPPBS is also validated by demonstrating that, in the collision dominated flow regime, the differential equations, which we converted from the derived integro-differential mass and momentum balance equations, are identical to the corresponding Navier-Stokes equations. This finding supports the assumption that, in the collision-dominated flow regime, the formulated integro-differential forms of the balance are exact implicit solutions for corresponding Navier-Stokes equations.

6) The analytical solution for determining the velocity profile induced in the model gas flow due to the pressure gradient along the channel is demonstrated. That the analytical solution is valid to explain velocity profiles in the wide range of gas pressure from the continuum through the slip and transition to free-molecule flow regimes gives additional confidence that the BM is valid.

7) From the practical viewpoint, the proposed method can be a fundamental base of a new generation of the CFD software in which the solver does not use the Navier-Stokes or Boltzmann equations. Using the solver operating on the exact implicit solutions of the balance equations will eliminate uncertainty, improve predictability, and shorten the computational time.

8) Although the BM is formulated to solve the CFD problems, it can be eventually used to model any dynamic system composed of presumably chaotically moving particles/elements, each carrying a specific amount of property/information.

The author declares no conflicts of interest regarding the publication of this paper.

Kislov, N. (2020) Ballistic Principle of the Property Balance in Space and Its Application to Modeling of Fluid Dynamics Problems. Journal of Applied Mathematics and Physics, 8, 1081-1122. https://doi.org/10.4236/jamp.2020.86084