Tensor-Centric Warfare III : Combat Dynamics with Delta-Strikes

This paper is the third part of the complex combat dynamics series, called tensor-centric warfare (for the first two parts, see [1] [2]). In the present paper, we extend the tensor combat model from [1] and [2] to model the dynamics of delta-strikes/missiles, which are temporally confined strong kinetic effects. The scenarios analyzed here include both deterministic and random delta-strikes which mimic single, multiple and continuous-time missile attacks. We also look at the bidirectional random strike as well as the general Hamilton-Langevin dynamics framework and provide an interpretation of the results obtained through simulation.


Introduction
Various aspects of warfare modelling have a long and successful history in Defence Science & Technology Group, Australia.To get a glimpse of the combat modelling in the previous decade, see [3] [4] [5] [6] [7], and in this decade, see [8]- [15].
In the first two papers of the tensor-centric warfare (TCW) series [1] [2], we have developed combat tensor dynamics (generalized from the Lanchester-type equations) with included commutators (for analysis of warfare symmetry) and entropic Lie-dragging (for modeling warfare uncertainty), given by the following tensor equations: where the Red and Blue forces are defined as vector-fields, ( ) ( ) x , and the terms on the right-hand side respectively represent linear Lanchester-type terms, quadratic Lanchester-type terms, commutators and entropic Lie-dragging terms (for detailed explanation see description of Equation (7) in [2]).
In the present paper, we extend this smooth tensor combat dynamics, with various combinations of temporally confined strong kinetic effects, which we call delta-strikes/missiles, developed in a perturbative fashion (adding more-and-more complex terms at each step).To allow for analysis and comparison of the methodology we have adopted a similar scenario framework to the first TCW-paper [1], with added various combinations of delta-strikes/missiles.We formulate this combat strike dynamics around the high-dimensional Dirac delta-function ( ) t δ 1 (see Figure 1).Since the response of any (complex) system to the delta-function input is called the impulse response, in the simulations below we will look at various impulse responses of the Red and Blue forces to the delta-strikes from the opposite side.
The initial scenario consists of nine Red and nine Blue aircraft which can be potential targets and could also act as a missile strike force.This will allow us to analyze the combat dynamics of our delta-strikes.We will initially look at a single strike which will involve a missile attached by a Red entity targeting and striking a single Blue target and after a short delay Blue is responding in the same way.
Then we introduce a multi-strike scenario which will model multiple missile attacks by Red on a Blue target and Blue responding in the same way.The multi-strike scenario is further extended by continuous (rapid) fire by Red on specific Blue targets and Blue responding in the same way.This discrete + continuous delta-strike spectrum is then naturally extended with a bidirectional random strike scenario, modeled and simulated to depict a missile exchange between the two adversaries.Recall that the Dirac delta-function, ( ) t δ , in mathematics called delta-distribution and in control engineering called impulse function, represents the limit of the sequence of zero-mean Gaussian distributions: Its main properties are:  Intelligent Control and Automation Thus our aim is to develop sophisticated mathematical tools to aid problem formulation and solution rather than facilitate immediate predictive capability with unrealistically simplistic models (see [3] [4] [5] [6] [7], and [8]- [15]).

Combat Dynamics with Delta-Strikes
In this section we will perturbatively extend tensor Lie-Lanchester Equations (1) so to develop general combat tensor dynamics with both deterministic and random delta-strikes.We formulate this combat δ-strike dynamics in the following seven steps (using Einstein's summation convention over repeated indices), corresponding to the physical background steps given in the Appendix.
Single strike.We start with the pure deterministic framework, with the basic nD Red-a R against Blue-a B combat dynamics in n  with a single strike (rocket) on each side: V. ( ) ( ) Red : , Blue : .
Discrete and continuous striking spectra.Now we extend Equations ( 4) by adding continuous striking spectrum, representing many-at-a-time (or, continuous time) rocket launchers on both sides: where δ-terms with sums correspond to discrete-time launchers (one rocket at a time) and δ-terms with integrals correspond to continuous-time launchers (many rockets at a time).
Bidirectional random strikes.Now we extend Equations ( 5) by including bidirectional (±) random strikes (i.e., random bidirectional striking at each other), with Red-and Blue-related random numbers: ( Simple Langevin equations.Now we move to the simple stochastic Langevin dynamics framework (see [16] V. where ab J is the Red-Blue interaction tensor, defined as the combined Red-Blue Combat tensor (see [1]): , where R L and B L denote the Lie derivatives in the direction of the Red and Blue vector-fields, respectively.In this way, we come to our final Red-Blue combat model in the form of tensor Hamilton-Langevin equations: ( This is our final tensor Lie-Lanchester combat model with δ-strikes/missiles. For completeness purpose, we finish this section with two probabilistic interpretations (compare with [16] [17] [18] [19]) of the random parts of the Hamilton-Langevin Equations ( 7):  Fokker-Planck equations: with the Red and Blue stationary solutions given by the equilibrium distributions: ( ) ( ) ( ) where ( )

Simulations and Discussion
Simulations of the combat strike-dynamics are performed in Mathematica @ as follows.For the calibration/comparing purpose, we firstly show the simulation of the tensor Lie-Lanchester Equations (1) given in Figure 2. Also, for better simulation results, we have approximated the Dirac δ-function by the narrow Gaussian (see Figure 1) implemented in Mathematica as: . We perform the complex impulse-response simulations as follows: 1) We start with the single strike combat simulation (see Figure 3), given by Equations ( 3), with the following parameters: , where 1 t tFin = is the final simulation time (= 10 units), ( ) , , , ρ ρ ρ ρ are real pseudo-random numbers generated from the [0,1]-interval.
2) Next, we simulate the multi-strike combat (see Figure 4), given by Equations ( 4), with the following new parameters: 3) Then, we simulate both discrete and continuous strike-spectra (see Figure 5), given by Equations (5).
4) Then, we simulate bidirectional random strikes (see Figure 6), given by Equations (6).5) As a final simulation (see Figure 7), we have implemented the Hamilton-Langevin Equations ( 8), where we have simplified γ and χ tensors and both skipped the computation-heavy Lie bracket terms (used in [2] for determining the warfare symmetry) and additional random forces, since all implemented components already include random numbers.the Blue side, which involved a missile attached by a Red entity targeting and striking a single Blue target and after a short delay Blue was responding in the same way.We can see from Figure 3 that since the Red (attacker) strike comes first and the Blue (defender) responds later, the Red-Blue phase-plot already shows quite different and more complex dynamics than the smooth dynamics on the initial Figure 2. From control-theoretic perspective, in the time plots of the Red and Blue forces we can see the expected effect of the impulse control input followed by the transients.The effect of the single impulse on each side is the convergence of the trajectories in the phase plane, forming the first "combat signature" for this specific type of attack.Multi-strike.Next, we have introduced a discrete-time multi-missile attacks by Red, ( ) , on a Blue target and Blue, ( ) , responding in the same way.We can see from Figure 4 that the amplitudes are now several times higher than in the single strike, for the Red the amplitude is all positive and for the Blue it starts negative and then grows positive.In contrast to the previous Figure, we can see now that the sequence of impulses leads to general irregularity, determined by the rate of the sequence of strikes, which cumulatively results in the increased complexity of the combat signature in the Red-Blue phase-plot.
Discrete and continuous striking spectra.The above multi-strike scenario has been further extended by continuous-time fire by Red, ( ) ( ) , on specific Blue targets and Blue, ( ) ( ) , responding in the same way.We can see from Figure 5 that the continuous spectrum defined by the integral of the high-dimensional impulse function has the "low-pass filtering" effect, which drastically reduces the overall complexity.This is related to the global resonance of the Red-Blue combat system.In other words, tuning the time intervals between the strikes produces the resonance effect (like in the radio tuning to a station), giving the significant increase in the amplitudes for both forces-in a half time (five units) the amplitudes are about six times bigger.However, this reduction of the overall complexity has the highly-potent psychological effect, as was well-known as "Katyusha-effect'' in the WW2: the continuous sound and vibration reduces morale and effects the combat readiness of the soldiers and commanders.
Bidirectional random strikes.The above discrete + continuous delta-strike spectrum has been naturally extended with a bidirectional random strike scenario, ( )( ) by the Red forces and by the Blue forces, simulated to depict a missile exchange between the two adversaries.We can see from Figure 6 that the overall combat complexity again increases but with obvious filtering effect from the continuous fire (many at a time).The combat signature in the Red-Blue phase-plot reappears, but is much smoother (low-pass filtering effect is still present).
Hamilton-Langevin dynamics.In the final instance we have given the formal Hamilton-Langevin stochastic generalization of the previous models.This stochastic formulation allows for a probabilistic interpretation of the battlefield and hence allows for a problem solving capability which is much closer to the reality of warfare.We can see from Figure 7 that adding the Hamilton-Langevin stochastic (non-delta) forces to the combat scenario from the previous Figure again adds the complexity to the combat signature in the Red-Blue phase-plane.
Basically, the more noise is added to dynamics, the dynamics becomes more complex.Introduction of random fluctuations into a complex dynamical/control system like our Red-Blue combat system, where the parts of the overall system are interacting (Red and Blue are engaging each other), results in the significant increase of the overall complexity/nonlinearity.However, we can see from the combat signature that the overall Red-Blue combat system is still performing within the edge of chaos (but not over the edge), where it is supposed to perform at its best (in the AI commutinty, it is generally assumed that neural networks have their highest performance at the edge of chaos, see [20] and the references therein).

Conclusions and Future Work
The purpose of this paper, the third part of the tensor-centric warfare (TCW) series, was to introduce a new framework for modeling complex warfighting (MCW), by combining continuous tensor quantities with temporally confined strong kinetic effects, which we call delta-strikes/missiles.The verification and/or validation of the proposed framework against the real military examples will be the subject of the future study.
We have presented various combat Red-Blue dynamics of delta-strikes, both deterministic and random-together with the corresponding impulse responses on the opposite side.We have shown that the TCW methodology presented in our previous papers [1] and [2] can be applied to analyzing the dynamic behavior of both deterministic and random delta-strike which mimic single, multiple and continuous-time missile attacks.
Also, Mathematica simulation of the Red-Blue pair of Fokker-Planck equations and the corresponding Red-Blue pair of Feynman path integrals is planned for the future study.Also, we foresee further development of the TCW framework to include tensors for learning and adaptations using deep learning.This will also involve developing the "Entropy Battle" software environment as an experimental testbed for future TCW development.The Entropy Battle will form the backbone of a future Battle Simulation Environment which is hoped to be utilized within the context of military planning.0 for , where e for j N j Discrete and continuous striking spectra.Now, we can imagine that, within a certain time interval, the strikes are continuously exerted on the ball, in a rapid-fire fashion, so that the total striking force ( ) tot f t includes both the discrete spectrum ( ) dis f t and a continuous spectrum ( ) In this way we obtain the more general multi--strike equation of motion: The general solution of Equation ( 14) is: where ( ) Bidirectional random strikes.Now, we introduce random strikes, by denoting the times at which strikes occur by j τ and indicate their 1D direction by ( ) , where the choice of the plus or minus sign is random.Assuming (for simplicity) the same intensity α for all random strikes, the random striking force (with random numbers: . Inclusion of the random striking force (16) into Equation ( 14) gives the randomized equation of motion: The general solution of Equation ( 17) is: Average random strikes.If we now observe many similar scenarios, we can perform an average ... over all of them, to obtain the average random striking force: Assuming that the times j τ at which the strikes happen are independent of In order to characterize the strength of the force ( 16), consider a quadratic expression in f by calculating the correlation function for two times , t t′ (see [22]): which simplifies into a single sum as: since the ones for j k = become 1 and for j k ≠ cancel each other.This correlation function is usually evaluated by assuming the Poisson process for the times of the strikes.
Simple Langevin equation.Now, if we assume: 1) that the average random striking force (18) vanishes, and 2) that the correlation function ( 19) can be further simplified as: where 0 t denotes the mean time between any two strikes and   , , where B k and T denote the Boltzmann constant and the temperature, Probabilistic interpretations.The Hamilton-Langevin Equation ( 23) has two probabilistic interpretations (see [24] [25]):  Fokker-Planck equation: with a stationary solution given by the equilibrium distribution: ( ) ( ) where i v  represent the auxiliary-response variables and N is the normalization constant.

(
Blue forces (with the Boltzmann and temperature constants), we will replace them with the coupled Red-Blue Poisson bracket , is the Lie bracket commutator (from Equations (1)): and  Feynman path integrals: the (Red, Blue) auxiliary-response variables and ( ) , N M are their respective normalization constants.

Figure 2 .
Figure 2. Sample simulation of the combat tensor Lie-Lanchester Equations (1)-for the comparison/calibration purpose: Red amplitude (left) is losing, Blue amplitude (center) is winning and Red-Blue phase plot (right) shows nine points of conflict; the amplitudes range from −50 to 50.

Figure 3 .
Figure 3. Sample Red-Blue combat simulation for 10 time units of the single strike Equations (3) with the parameters given in the text.We can see that the Red (attacker) strike comes first (left) and the Blue (defender) responds later (center).The amplitudes range up to 150.The Red-Blue phase-plot (right) already shows complex combat dynamics.

Figure 4 .
Figure 4. Sample Red-Blue combat simulation for 10 time units of the multi-strike Equations (4) with the parameters given in the text.Here the amplitudes are several times higher than in the single strike, for the Red (left) the amplitude is all positive and for the Blue (center) it starts negative and then grows positive.The Red-Blue phase-plot (right) is even more complex.

Figure 5 .
Figure 5. Sample Red-Blue combat simulation for 5 time units with both discrete and continuous strike-spectra given by Equations.(5) with the parameters given in the text.The amplitudes are similar as before, with the roles reversed: the Blue (center) the amplitude is all positive and for the Red (left) it starts negative and then grows positive.The Red-Blue phase-plot (right) is somewhat less complex, showing the filtering capacity of the integral in the continuous spectrum.

Figure 6 .
Figure 6.Sample Red-Blue combat simulation for 3.3 time units with added bidirectional random strikes given by Eqs.(6) with the parameters given in the text.Since complexity increases, we had to shorten simulation time.A characteristic signature is shown in the Red-Blue phase plane.

Figure 7 .
Figure 7. Sample Red-Blue combat simulation for 3.3 time units with added Hamilton-Langevin strikes given by Equations (8) with the parameters given in the text.Here we have even higher complexity with a new signature in the Red-Blue phase plane.
Since both directions are equally possible, they cancel each other, and the average velocity vanishes.Therefore, from the integral solution(15) equal times t and τ , simplifies into: equation.The simple 1D Langevin Equation (21) can be generalized in the framework of Hamiltonian dynamics (see[23]) to the V.Ivancevic et al.
bracket and H is the Hamiltonian of the system given as Ivancevic et al.