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      , and in this decade, see  -  .
In the first two papers of the tensor-centric warfare (TCW) series   , 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, and , 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  ).
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  , with added various combinations of delta-strikes/missiles. We formulate this combat strike dynamics around the high-dimensional Dirac delta-function 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.
Figure 1. Narrow Gaussian approximation for the shifted Dirac delta-function: , corresponding to the values of (center) and (right).
Next, we introduce stochastic forces of arbitrary (non-delta) nature into our model by utilizing a simple Langevin representation of the discrete, continuous and bidirectional delta-strike spectrum between the Red and Blue aircraft. In the final instance we give the formal Hamiltonian formulation of this stochastic Langevin model. This will allow for a probabilistic interpretation of the battlefield and hence allow for a problem solving capability which is much closer to the reality of warfare.
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      , and  -  ).
2. 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― against Blue― combat dynamics in with a single strike (rocket) on each side:
where and are the Red and Blue “damping” coefficients, respectively, and and are the Red and Blue strike times; the terms with δ-functions represent single strikes/missiles.
Multi-strike. Now we generalize δ-functions in Equations (3) to include N multi-strikes for the Red forces (happening at times ) and M multi-strikes for the Blue forces (happening at times ), representing one-at-a-time (or, discrete time) multi-missile launchers:
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     ) and extend Equations (6) with random (non-delta) forces Red- and Blue- :
Hamilton-Langevin dynamics. Finally, we formulate the general Langevin dynamics in the Hamiltonian (symplectic geometric) framework to be able to extend the simple random forces Red- and Blue- in Equations (6) with the general Langevin forces derived from the hypothetical global Red-Blue combat Hamiltonian energy function, , defined in an Ising-Hamiltonian fashion as:
where is the Red-Blue interaction tensor, defined as the combined Red-Blue Combat tensor (see  ): , where is the random identity matrix (or, random Kronecker delta), and is the equilibrium probability distribution of the coupled -system. Besides, instead of the separate Poisson brackets for the Red forces and for the Blue forces (with the Boltzmann and temperature constants), we will replace them with the coupled Red-Blue Poisson bracket , which is the Lie bracket commutator (from Equations (1)): , where and 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     ) 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: and with constants ; and
• Feynman path integrals:
where represent the (Red, Blue) auxiliary-response variables and are their respective normalization constants.
3. 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 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  for determining the warfare symmetry) and additional random forces, since all implemented components already include random numbers.
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. 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. 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. 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. 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. 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.
Now, we discuss various simulated scenarios with delta-strikes as follows:
Single strike. The initial scenario consisted of nine Red and nine Blue aircraft which were both potential targets and missile strike forces. We have started with a single strike from each side, on the Red side, and 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  and the references therein).
Complete TCW equation. (with transparent terms) read:
where Hamilton-Langevin delta strikes, and , read:
4. 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  and  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.
The authors are grateful to Dr. Tim McKay and Dr. Brandon Pincombe, Joint and Operations Analysis Division, Defence Science & Technology Group, Australia―for their constructive comments which have improved the quality of this paper.
Appendix: Deterministic and Random Strikes
In this section we develop the necessary physical background of both deterministic (Newtonian) and random (Langevin) strikes, in the following seven steps (compare with  ):
Single strike. To start with, consider a one-dimensional (1D) damped motion of a ball of unit mass, moving with velocity v, within a viscous fluid with damping coefficient . The ball is pushed by an external striking force, defined in terms of the Dirac delta-function , with intensity/amplitude , occurring at time , as:
From the second Newton’s law we obtain the single-strike equation of motion:
where overdot represents the time derivative. The general solution of this ordinary differential equation (ODE) is given in terms of the Green function:
Multi-strike. Instead of a single strike given by the delta-force (9) we will now hit the ball by a series of N strikes with intensities , occurring at a sequence of times , so that the total striking force becomes:
If we replace this multi-strike force into the ODE (10), we obtain the multi-strike equation of motion:
The general solution of the ODE (12) is:
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 includes both the discrete spectrum and a continuous spectrum :
In this way we obtain the more general multi--strike equation of motion:
The general solution of Equation (14) is:
Bidirectional random strikes. Now, we introduce random strikes, by denoting the times at which strikes occur by 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: ) can be defined as:
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 at which the strikes happen are independent of the direction of the strikes, we can split (18) into the product:
Next, if the strikes happen in both directions with equal frequency, they simply cancell each other:
In order to characterize the strength of the force (16), consider a quadratic expression in f by calculating the correlation function for two times (see  ):
which simplifies into a single sum as:
since the ones for become 1 and for 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 denotes the mean time between any two strikes and represents the random fluctuation, the equation of motion driven solely by the random striking force becomes the simple Langevin equation:
Since both directions are equally possible, they cancel each other, and the average velocity vanishes. Therefore, from the integral solution (15) of Equation (14) we obtain
In the steady--state, Equation (22) reduces to:
which, for equal times t and , simplifies into:
Hamilton-Langevin equation. The simple 1D Langevin Equation (21) can be generalized in the framework of Hamiltonian dynamics (see  ) to the following nD form:
where and T denote the Boltzmann constant and the temperature, represents the Poisson bracket and H is the Hamiltonian of the system given as , where is the equilibrium probability distribution of . The Gaussian-distributed random forces have the correlation function:
which implies the Onsager reciprocity relations for the damping coefficients .
Probabilistic interpretations. The Hamilton-Langevin Equation (23) has two probabilistic interpretations (see   ):
• Fokker-Planck equation:
with a stationary solution given by the equilibrium distribution:
, with some constant c;
• Feynman path integral:
where represent the auxiliary-response variables and N is the normalization constant.
1Recall that the Dirac delta-function, , 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:
 Pincombe, A.H. and Pincombe, B.M. (2007) Markov Modelling on the Effectiveness of Sanctions: A Case Study of the Falklands War. ANZIAM Journal, 48, 527-541.
 Kosowski, L.R., Pincombe, A.H. and Pincombe, B.M. (2010) Irrelevance of the Fractal Dimension Term in the Fractal Attrition Equation, ANZIAM Journal, 52, 988-1008.
 Millikan, J., Wong, M. and Grieger, D. (2013) Suppression of Dismounted Soldiers: towards Improving Dynamic Threat Assessment in Closed Loop Combat Simulations. In: Piantadosi, J., Anderssen, R.S. and Boland, J., Eds., 20th International Congress on Modelling and Simulation (MODSIM 2013), Adelaide, 2013, 1054-1060.
 Bowden, F.D.J., Pincombe, B.M. and Williams, P.B. (2015) Feasible Scenario Spaces: A New Way of Measuring Capability Impacts. In: Weber, T., McPhee, M.J. and Anderssen, R.S., Eds, 21st International Congress on Modelling and Simulation (MODSIM2015), Gold Coast, 2015, 836-842.
 Ramirez, M., Papasimeon, M., Benke, L., Lipovetzky, N., Miller, T. and Pearce, A.R. (2017) Real-Time UAV Maneuvering via Automated Planning in Simulations. 26th International Joint Conference on Artificial Intelligence (IJCAI-17), Melbourne, 19-25 August 2017, 5243-5245.
 Ramirez, M., Papasimeon, Lipovetzky, N., Benke, L., Miller, T., Pearce, A.R., Scala, E. and Zamani, M. (2018) Integrated Hybrid Planning and Programmed Control for Real-Time UAV Maneuvering. Proceeding 17th International Conference on Autonomous Agents and Multiagent Systems (AAMAS 2018), Stockholm, 10-15 July 2018, 1318-1326.
 Kawasaki, K. (1973) Simple Derivations of Generalized Linear and Nonlinear Langevin Equations. Journal of Physics A: Mathematical, Nuclear and General, 6, 1289.
 Janssen, H.K. (1976) On a Lagrangean for Classical Field Dynamics and Renormalization Group Calculations of Dynamical Critical Properties. Zeitschrift für Physik B Condensed Matter, B, 23, 377-380.