Ballistic Principle of the Property Balance in Space and Its Application to Modeling of Fluid Dynamics Problems

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.

tions 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 [1]. However, this intentional simplification of a fluid model may diminish the usefulness of the results of the computations. Also, the theoretical understanding of the solutions to these equations is still inadequate. Specifically, for three-dimensional Navier-Stokes equations and given initial conditions, mathematicians have not yet proved that smooth solutions always exist, and the solutions have limited energy per unit mass.
In the most general form, the Navier-Stokes equations of mass and momentum conservation for compressible fluid are expressed as: and ( ) respectively, for 0 t > and 3 ∈ r   . 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 [2] and [3]). In rarefied gases, the mean free path or the average distance traveled between intermolecular collisions is also considered a fail factor for the Navier-Stokes description [4]. Here we shall enumerate the most important properties of gases that often behave very close to ideal [5]: 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.
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.

Physical Principles of the Fluid Model
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 [14]. We have also realized that 1) transport phenomena include any situation that involves a net transfer of property/information between particles, which equals randomizing properties between interacting particles; 2) a randomized by property exchange between interacting particles physical/statistical property can be taken out into the surrounding gas [14]. However, we exclude interaction resulted in chemical reactions between colliding particles and fragmentation or fusion.
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.

Properties and Features of the Model Gas
We assign these unique properties to the model gas [12]: 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. 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. 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. Figure 1 illustrates a schematic of the ballistic trajectories of a particle between two subsequent collisions in a model gas system. 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 [12]. For clarity, the observer's Cartesian coordinate system 100 is orientated, so that y-axis is directed in the opposite direction of an applied field of external force 107, which provides, for each of particles, acceleration g  . Each of the particles (shown as black disks), particularly a particle 101 (shown as white disks), travels between its two consecutive collisions: an initial collision 102 and an ending collision 103 by following a trajectory 104 or 105 governed by applicable law of motion including Newton's laws of motion.
Note that in a lack of external force, all ballistic trajectories between consecutive collisions will be just straight lines, respectively, as indicated by trajectory 106.
Specifically, referring to Figure 1 where B k is Boltzmann constant, T is the temperature, and m is the mass of a particle.

Principles of Construction of the Property Balance in the Model Gas
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 T v , 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 Figure 2 as follows [12]. In a given non-moving point r  at a given time t, the net rate of property influx per unit 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. Journal of Applied Mathematics and Physics 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: where _F in B Ψ 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.

Defining a Net Rate of Total Property Influx Per Unit Volume in a General Non-Moving Point at a Given Time from the Surrounding Model Gas
We have recognized and explained afterward that there exists a combination of a specific direction of an initial instant vector of thermal velocity ( ) , which allows each of the selected particles to arrive in the given non-moving point r  at the given time, t [12]. These particles, which originate from the initial collisions within the whole model gas system, form the converging flux 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 influx per unit volume is presented in Table 1. Table 1. List of the model parameters associated with defining the net rate of total property influx per unit volume.

Parameters
Short description x y z the velocity vector of the ballistic particle at time t  g  an external force that applies to a particle of a unit mass rel v the average magnitude of the velocity of the traveling particle with respect to a nearby passed particle Figure 3 shows a perspective view of the model gas system for explaining a ballistic movement of the converging ballistic particle after the initial collision, which is affected by the external field of force. For clarity, the observer's Cartesian coordinate system 100 is oriented, so y-axis is along the negative direction of the applied acceleration field 306, which provides, for each particle, accelera- 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  (not shown) at time t  by applying Equations (6) and (7), respectively, given below: Table 1, and i t′ is defined in (2) below. In Equations (6) and (7) and which is obtained by substitution of Equation (8) in Equation (6) followed by the assignment of t t =  and substitution of i ϕ defined in Table 1 in Equation (6) and rearrangement of the terms. Upon resolving Equation ( 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 i t′ by presenting Equation (7) in the following form: where c i r  is defined as: In this, vector c i r  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 i t′ of the divergence, which is ob- and computing sp i R from Equation (14) as The velocity vector c i v  of the center of the expansion zone 302 at time t is computed as 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 which is obtained by assigning t t =  and substitution of (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: is the average magnitude of the velocity with respect to a nearby passed particle in the trajectory point ( )  Figure 4 shows a perspective view for explaining a method for determining the average magnitude of the instantaneous velocity of the ballistic particle with respect to a nearby passed particle in a three-dimensional configuration. For clarity, the observer's Cartesian coordinate system 100 is oriented, so y-axis is along the negative direction of the applied acceleration field, 401, which provides, for each particle, acceleration g  . Here, at time t  , particle P1 having a velocity 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: 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 ( ) (22) and executing the integration of the resulted equation: is the average magnitude of the velocity with respect to a nearby passed particle in the trajectory point, T v t r    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 ( ) in Equation (24): 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: where c is the speed of light. For 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 where v  is defined by Equation (17), i Q is a survival probability defined by Equation (18) 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: where , x y z 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  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 Ψ Step 5: Defining the rate of the property vector flux Ψ → FS r J   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: 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 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 by applying Equation (35) given below: 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. i t t ′ < .

Defining a Net Rate of Total Property Efflux Per Unit Volume from the General Non-Moving Point at the Given Time
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 Table 2. Figure 5 is a schematic shown to illustrate a method for analytical representation of the net rate of property efflux from a non-moving point in three-dimensional space occupied by the model gas. For clarity, the observer's Cartesian coordinate system 100 is oriented, so the y-axis is along the negative direction of the applied acceleration field 504. In Figure  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. Step and 2) Determining the time needed, for each particle diverging from a and substitution of a t 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 by applying Equation (45) given below: where a t t ϕ + ′ = + is the time of the particle positioning in point ′ r  , is the number of particles placed within a collision tube of a unit length in the is an average magnitude of the velocity with respect to a nearby passed particle in the trajectory point ( ) Step 3: Defining the vector/tensor field of property flux r J Ψ along with one of the ballistic trajectories of a diverging particle in point ′ r  .
Step 3 includes representing the property vector flux r J Ψ by applying Equation (46) as given below: a a n t Q t t t t t t Step 4: Defining the net rate of property efflux per unit volume _FS out B Ψ from the general non-moving point r  at the given time t.
Step 4 includes representing a a r t n t Q t t t t t t 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 → c a r r   .

Analytical Representation of a General Integro-Differential form of Property Balance Equation in the Three-Dimensional Model Gas System
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): 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 0 t ≥ as an implicit function of r  on 3  , i.e.
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 i t t ′ < (past). Whereas the diverging ballistic particles originated from collisions at time t (present) transport the balanced property/information into surrounding toward the future a t t ′ > .

Analytical Representation of a General Integro-Differential Form of Mass Balance Equation in the Three-Dimensional Model Gas System
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: Then, we obtain the following general integro-differential form of the mass balance equation

Analytical Representation of a General Integro-Differential Form of Momentum Balance Equation in the Three-Dimensional Model Gas System
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) , , 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.

Analytical Representation of a General Integro-Differential Form of Energy Balance Equation in the Three-Dimensional Model Gas 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 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 r  at time t, which may appear because of a non-uniform pressure applied to the system.

Integro-Differential Forms of Mass and Momentum Balance Equations in Collision-Dominated Flow Regime
We have recognized that, in the majority of real conditions on Earth and near-Earth space, first, the acceleration field g  applied to each of the model gas particles is approximately a constant, namely: which is expressed: 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 1 c P , which is expressed: We also note here that since point r  is excluded from integration in the domain of integration Ω in which ′ ≠ r r   , the operation of differentiation regarding a parameter r  is interchangeable with the operation of integration over some other variable ′ r  . Also, these functions or approximations for functions involving in calculations according to Equations (50), (53), and (56) are formulated: , , where 0 i n  is the unit vector of arbitrary direction from the point ′ r  , which is approximated from i n  as and 0 i ϕ is the ballistic traveling time, which is approximated from i ϕ as where 0 i t′ is the time of the particle positioning in point ′ r  , which is approximated from i t′ as 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.

Analytical Representation of an Integro-Differential Form of the Mass Balance Equation in the Collision-Dominated Flow Regime
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 r  at given time t and originates from a preceding collision within the model gas system at a time earlier than time t. Such a combination of converging and diverging ballistic particles capable of targeting or escaping point r  at given time t is treated as an exhaustive combination. Therefore, from a physical viewpoint, the solution ( ) , t u r   is unique. 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 ′ → r r   in the second left-hand term of the resulting integro-differential mass balance equation; 3) remembering that point r  is not included in integration for converging ballistic particle; and 4) neglecting terms containing 1 c P rel v .
Upon executing the above, Equation (50) is reduced to: ,  ,  2  1  1  1  ,  ,  ,  ,  d  ,  4   1  ,  ,  ,  ,  , , d , the rate of collisions per unit volume in the given non-moving point r  at time t, which is expressed as

Analytical Representation of an Integro-Differential Form of the Momentum Balance Equation in the Collision-Dominated Flow Regime
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 r  at given time t and originates from a preceding collision within the model gas system at a time earlier than time t. Such a combination of converging and diverging ballistic particles capable of targeting or escaping point r  at given time t is treated as an exhaustive combination.
Therefore, from a physical viewpoint, the solution ( ) , t u r   is unique. 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 2) substituting approximations for functions involving in Equation (48) It implies that Equation (80) defines u  at 0 t ≥ as an implicit function of

Validation Tests
The following validation tests are aimed to demonstrate the feasibility of the proposed method. Figure 6 is a perspective view for explaining a method for determining the total rate of collisions per unit area on a surface being in contact with the gas. In a system with no gravitational force, the ballistic particles have straight-line trajectories. Here we limit further our consideration to an incompressible model gas at the uniform temperature in a steady-state condition. In the semi-sphere filled with the model gas over the being at rest surface s A having a directional vector n  , each particle having an instantaneous randomly directed vector of the thermal velocity of magnitude, T v may have the instant vector-velocity component directing a particle toward the surface s A . In Figure 6  , the total rate of collisions per unit area on the surface s A is given as:

Determining the Total Rate of Collisions Per Unit Area on a Surface Being in Contact with the Gas
and V Z is the rate of collisions per unit volume, which, upon substitution of Equation (25) in Equation (77) is expressed as where n is particle density, T v is the magnitude of the thermal velocity, and c P is the number of particles placed within a collision tube of a unit length.
Finally, using the geometry illustrated in Figure 6, considering that ( ) and Figure 6. The perspective view of the geometry and coordinate system for determining the total rate of collisions per unit area and the pressure exerted on the surface from the gas volume. , we obtain: 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.

Determining the Pressure Exerted on the Surface from the Entire Gas Volume Being in Contact with the Surface
Referring to the previous section and Figure 6 we recognize that each particle of the model gas carries the momentum and the momentum component delivered by the particle to the surface s A is 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 Figure 7 and assigning 2 = ∆ in y p Ψ , the pressure P exerted on the surface s A is given as: 3 cos 3 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) [15]. Journal of Applied Mathematics and Physics

Direct Validation of the BPPBS in the Three-Dimensional Unlimited Incompressible Gas Space at the Uniform Temperature
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, T v , are constant. In the model gas system being at rest and characterized by the above conditions, when the applied to the model gas system external pressure is uniform, the mass-flow velocity, u  , is expected to be zero. Now we may formulate the mass-balance equation by adopting Equations (76) as 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 r  and transforming the Cartesian coordinate system into the spherical coordinate system, then, after defining Analogously, after substitution of Equations (74), (83), and (94) in Equation (93) and taking into account that point 0 ρ′ =   is not included into integration, we calculate the right-hand of Equations (93) as follows: here c P is the number of particles placed within a collision tube of a unit length, V Z is the rate of collisions per unit volume, which is defined by Equation (84), and ρ ρ ′ ′   is a unit vector with the point of origin at 0 ρ′ =   .
That the results of integration provided by Equations (95)

Incompressible Model Gas Flow between Two Infinite Parallel Plates at the Uniform Temperature in a Case of Diffuse Particle Scatterings
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. Figure 7 shows schematically a one-dimensional 2 , where u is the gas slip velocity near the wall, w u is the tangential velocity of the wall, and f λ is the mean free path.
The method described below uses none of the Maxwell-type assumptions.

Properties and Features of the Gas-Solid Interface
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 [13]: 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 σ in a corresponding point of diffuse particle.
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 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 0 T mv + n , is zero.
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.

Analytical Derivation of the Velocity Profile Induced in the Model Gas Due to the Pressure Gradient along the Channel between Two Parallel Plates
Step 1: Formulating an integral form of the x u -momentum balance equation Considering the above, we obtain the following integral form of the x u momentum balance equation in a steady-state model gas flow [12]: Step 2: Obtaining a differential form matching to the corresponding integral form of the x u momentum balance equation.
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: 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 0 y = and y H = for the balance establishment.
Step 5: Solving a system of the obtained in Step 4 algebraic equations.
We obtained the following functional relationship of fluid velocity x u from other properties and geometry parameters characterizing the model gas system: The equation above reveals there is a minimum in the normalized mass flow rate (at about 0.3 Kn ≅ ), which is called the Knudsen paradox in the literature [18]. It is explained by the fact that [w]ith increasing Knudsen number, the additional diffusive flux becomes significant and deviates from the no-slip solution. Beyond a critical limit ( 0.3 Kn > ), the diffusive flux dominates the convective flux [18]. Figure 8 compares the non-dimensional flow rates calculated by Equation (110) for the proposed model (Ballistic Model) and a selection of slip models proposed by various authors. The non-dimensional flow rate was proposed in the form [19] ( ) ( ) where corresponding slip coefficients 1 C and 2 C are taken from Table 6 in [20] The results of the linearized Boltzmann solution of [21] are plotted from Figure 2 of [17].
In Figure 8, the flow rate according to Equations (110) is rescaled by 2 π factor for the comparison purposes.
We compared the non-dimensional flow rates calculated by Equation (110) for the Ballistic Model and a selection of slip models proposed by various au-  Table 6 in [20]).

Incompressible Model Gas Flow at the Uniform Temperature in the Space between Two Being at Rest Infinite Parallel Plates in a Case of Mixed Diffuse and Specular Particle Scatterings
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 σ is the momentum accommodation coefficient, which is the probability, for an incident particle, to accommodate momentum from the gas-solid interface and to scatter back in the model gas as a diffuse particle. Comparative analysis of the integral equations describing model gas flow confined between the parallel plates with purely diffuse scattering (see Equation (103)) and mixed diffuse and specular scattering (see Equation (112)) results in finding they have similar forms if where xH u and 0 x u represent coefficients of Equation (103), and the right-hand of the equation above represents the similar coefficient in Equation (112). The analysis also leads to the conclusion that the velocity profile can be described by Equation (105), where coefficient A is expressed by Equation (106). (107), we obtain:
Finally, substituting in Equation (113) and substituting the resulting equation in Equation (108), then executing corresponding integrations and algebra operations, we obtain: Remarkably, the derived tangential slip velocity coefficient C does contain the term being proportional to 2 σ σ − . Its appearance is the result of the application of the BM but not the result of usage of the semi-empirical Maxwell-type assumptions as of Equation (97). (105)

Substitution in Equation
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

Reducing Integro-Differential Forms of Mass and Momentum Balance Equations into Corresponding Differential Forms
According to Tenenbaum and Pollard [26], to test whether an implicit function defined by the relation ( ) , 0 f x y = is a solution of a given differential equation; there is a need to show that the function does satisfy a given differential equation for every x on I.
The standard procedure in calculus to prove that an implicit function defined by relationship ( ) , 0 f x y = (a) is a solution of a given differential equation , f x y defines as an implicit function of x on the same interval I: a x b < < . 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 [27]. Specifically, we will use this technique for reducing vector integro-differential balance equations into corresponding vector differential balance equations. Here we need also note that the point r  should be included when we apply the operator of vector differentiation at this point. In the domain of integration Ω in which ′ ≠ r r   , operation of differentiation regarding a parameter r  is interchangeable with the operation of the integration over some other variable ′ r  . Even if domain Γ includes singularity point r  , one can easy determine that the integrals containing the term are zeroed. The "fate" of other integrals containing singularity point r  because of the Green function,

Reducing the Integro-Differential Form of Mass Balance Equation into a Corresponding Vector Differential Form
Deriving a vector differential form of mass balance equation is shown: Step  Step 3. Comparing the first three right-hand terms in the equation above and Equation (76), we may suggest a modification of Equation (76)      Step 2. Applying divergence operator ∇ ⋅ to the left-and right-hands of the equation above and neglecting insignificant terms by the procedures, which is similar to the described in Section III:  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.

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