Lanchester's laws

Last updated

Lanchester's laws are mathematical formulae for calculating the relative strengths of military forces. The Lanchester equations are differential equations describing the time dependence of two armies' strengths A and B as a function of time, with the function depending only on A and B. [1] [2]

Contents

In 1915 and 1916 during World War I, M. Osipov [3] :vii–viii and Frederick Lanchester independently devised a series of differential equations to demonstrate the power relationships between opposing forces. [4] Among these are what is known as Lanchester's linear law (for ancient combat) and Lanchester's square law (for modern combat with long-range weapons such as firearms).

As of 2017 modified variations of the Lanchester equations continue to form the basis of analysis in many of the US Army’s combat simulations, [5] and in 2016 a RAND Corporation report examined by these laws the probable outcome in the event of a Russian invasion into the Baltic nations of Estonia, Latvia, and Lithuania. [6]

Lanchester's linear law

For ancient combat, between phalanxes of soldiers with spears for example, one soldier could only ever fight exactly one other soldier at a time. If each soldier kills, and is killed by, exactly one other, then the number of soldiers remaining at the end of the battle is simply the difference between the larger army and the smaller, assuming identical weapons.

The linear law also applies to unaimed fire into an enemy-occupied area. The rate of attrition depends on the density of the available targets in the target area as well as the number of weapons shooting. If two forces, occupying the same land area and using the same weapons, shoot randomly into the same target area, they will both suffer the same rate and number of casualties, until the smaller force is eventually eliminated: the greater probability of any one shot hitting the larger force is balanced by the greater number of shots directed at the smaller force.

Lanchester's square law

Lanchester's square law is also known as the N-square law.

Description

Idealized simulation of two forces damaging each other neglecting all other circumstances than the 1) Size of army 2) Rate of damaging. The picture illustrates the principle of Lanchester's square law. Damagerace.JPG
Idealized simulation of two forces damaging each other neglecting all other circumstances than the 1) Size of army 2) Rate of damaging. The picture illustrates the principle of Lanchester's square law.

With firearms engaging each other directly with aimed shooting from a distance, they can attack multiple targets and can receive fire from multiple directions. The rate of attrition now depends only on the number of weapons shooting. Lanchester determined that the power of such a force is proportional not to the number of units it has, but to the square of the number of units. This is known as Lanchester's square law.

More precisely, the law specifies the casualties a shooting force will inflict over a period of time, relative to those inflicted by the opposing force. In its basic form, the law is only useful to predict outcomes and casualties by attrition. It does not apply to whole armies, where tactical deployment means not all troops will be engaged all the time. It only works where each unit (soldier, ship, etc.) can kill only one equivalent unit at a time. For this reason, the law does not apply to machine guns, artillery with unguided munitions, or nuclear weapons. The law requires an assumption that casualties accumulate over time: it does not work in situations in which opposing troops kill each other instantly, either by shooting simultaneously or by one side getting off the first shot and inflicting multiple casualties.

Note that Lanchester's square law does not apply to technological force, only numerical force; so it requires an N-squared-fold increase in quality to compensate for an N-fold decrease in quantity.

Example equations

Suppose that two armies, Red and Blue, are engaging each other in combat. Red is shooting a continuous stream of bullets at Blue. Meanwhile, Blue is shooting a continuous stream of bullets at Red.

Let symbol A represent the number of soldiers in the Red force. Each one has offensive firepower α, which is the number of enemy soldiers it can incapacitate (e.g., kill or injure) per unit time. Likewise, Blue has B soldiers, each with offensive firepower β.

Lanchester's square law calculates the number of soldiers lost on each side using the following pair of equations. [7] Here, dA/dt represents the rate at which the number of Red soldiers is changing at a particular instant. A negative value indicates the loss of soldiers. Similarly, dB/dt represents the rate of change of the number of Blue soldiers.

The solution to these equations shows that:

The first three of these conclusions are obvious. The final one is the origin of the name "square law".

Relation to the salvo combat model

Lanchester's equations are related to the more recent salvo combat model equations, with two main differences.

First, Lanchester's original equations form a continuous time model, whereas the basic salvo equations form a discrete time model. In a gun battle, bullets or shells are typically fired in large quantities. Each round has a relatively low chance of hitting its target, and does a relatively small amount of damage. Therefore, Lanchester's equations model gunfire as a stream of firepower that continuously weakens the enemy force over time.

By comparison, cruise missiles typically are fired in relatively small quantities. Each one has a high probability of hitting its target, and carries a relatively powerful warhead. Therefore, it makes more sense to model them as a discrete pulse (or salvo) of firepower in a discrete time model.

Second, Lanchester's equations include only offensive firepower, whereas the salvo equations also include defensive firepower. Given their small size and large number, it is not practical to intercept bullets and shells in a gun battle. By comparison, cruise missiles can be intercepted (shot down) by surface-to-air missiles and anti-aircraft guns. So it is important to include such active defenses in a missile combat model.

Lanchester's law in use

Lanchester's laws have been used to model historical battles for research purposes. Examples include Pickett's Charge of Confederate infantry against Union infantry during the 1863 Battle of Gettysburg, [8] the 1940 Battle of Britain between the British and German air forces, [9] and the Battle of Kursk. [10]

In modern warfare, to take into account that to some extent both linear and the square apply often, an exponent of 1.5 is used. [11] [12] [3] :7-5–7-8 Lanchester's laws have also been used to model guerrilla warfare. [13]

Attempts have been made to apply Lanchester's laws to conflicts between animal groups. [14] Examples include tests with chimpanzees [15] and ants. The chimpanzee application was relatively successful. A study of Australian meat ants and Argentine ants confirmed the square law, [16] a study of fire ants did not confirm the square law. [17]

Helmbold Parameters

The Helmbold Parameters provide quick, concise, exact numerical indices, soundly based on historical data, for comparing battles with respect to their bitterness and the degree to which side had the advantage. While their definition is modeled after a solution of the Lanchester Square Law's differential equations, their numerical values are based entirely on the initial and final strengths of the opponents and in no way depend upon the validity of Lanchester's Square Law as a model of attrition during the course of a battle.

The solution of Lanchester's Square Law used here can be written as:

where is the time since the battle began, and are the surviving fractions of the attacker's and defender's forces at time , is the Helmbold intensity parameter, is the Helmbold defender's advantage parameter, is the duration of the battle, and is the Helmbold bitterness parameter.

If the initial and final strengths of the two sides are known it is possible to solve for the parameters , , , and . If the battle duration is also known, then it is possible to solve for . [18] [19] [20]

If, as is normally the case, is small enough that the hyperbolic functions can, without any significant error, be replaced by their series expansion up to terms in the first power of , and if we adopt the following abbreviations for the casualty fractions

then the following approximate relations hold: [21]

That is a kind of "average" (specifically, the geometric mean) of the casualty fractions justifies using it as an index of the bitterness of the battle.

We note here that for statistical work it is better to use the natural logarithms of the Helmbold Parameters. We will call them, in an obvious notation, , , and .

Major findings

See Helmbold (2021):

  1. The Helmbold parameters and are statistically independent, i.e., they measure distinct features of a battle. [22]
  2. The probability that the defender wins, P(Dwins), is related to the defender's advantage parameter via the logistic function, P(Dwins) = 1 / (1 + exp(-z)), with z = -0.1794 + 5.8694 * logmu. [23] This logistic function is almost exactly skew-symmetric about logmu = 0, rising from P(Dwins) = 0.1 at logmu = -0.4, through P(DWins) = 0.5 at logmu = 0, to P(Dwins) = 0.9 at logmu = +0.4. Because the probability of victory depends on the Helmbold advantage parameter rather than the force ratio, it is clear that force ratio is an inferior and untrustworthy predictor of victory in battle.
  3. While the defender's advantage varies widely from one battle to the next, on average it has been practically constant since 1600CE. [24]
  4. Most of the other battle parameters (specifically the initial force strengths, initial force ratios, casualty numbers, casualty exchange ratios, battle durations, and distances advanced by the attacker) have changed so slowly since 1600CE that only the most acute observers would be likely to notice any change over their nominal 50-year military career. [25]
  5. Bitterness (), casualty fractions ( and in the above notation), and intensity () also changed slowly before 1939CE. But since then they have followed a startlingly steeper declining curve. [26]

Some observers have noticed a similar post-WWII decline in casualties at the level of wars instead of battles. [27] [28] [29] [30]

See also

Related Research Articles

The propagation constant of a sinusoidal electromagnetic wave is a measure of the change undergone by the amplitude and phase of the wave as it propagates in a given direction. The quantity being measured can be the voltage, the current in a circuit, or a field vector such as electric field strength or flux density. The propagation constant itself measures the dimensionless change in magnitude or phase per unit length. In the context of two-port networks and their cascades, propagation constant measures the change undergone by the source quantity as it propagates from one port to the next.

<span class="mw-page-title-main">Astronomical coordinate systems</span> System for specifying positions of celestial objects

In astronomy, coordinate systems are used for specifying positions of celestial objects relative to a given reference frame, based on physical reference points available to a situated observer. Coordinate systems in astronomy can specify an object's position in three-dimensional space or plot merely its direction on a celestial sphere, if the object's distance is unknown or trivial.

<span class="mw-page-title-main">Noether's theorem</span> Statement relating differentiable symmetries to conserved quantities

Noether's theorem states that every continuous symmetry of the action of a physical system with conservative forces has a corresponding conservation law. This is the first of two theorems proven by mathematician Emmy Noether in 1915 and published in 1918. The action of a physical system is the integral over time of a Lagrangian function, from which the system's behavior can be determined by the principle of least action. This theorem only applies to continuous and smooth symmetries of physical space.

In statistics, the Gauss–Markov theorem states that the ordinary least squares (OLS) estimator has the lowest sampling variance within the class of linear unbiased estimators, if the errors in the linear regression model are uncorrelated, have equal variances and expectation value of zero. The errors do not need to be normal for the theorem to apply, nor do they need to be independent and identically distributed.

In statistics, ordinary least squares (OLS) is a type of linear least squares method for choosing the unknown parameters in a linear regression model by the principle of least squares: minimizing the sum of the squares of the differences between the observed dependent variable in the input dataset and the output of the (linear) function of the independent variable.

The Havriliak–Negami relaxation is an empirical modification of the Debye relaxation model in electromagnetism. Unlike the Debye model, the Havriliak–Negami relaxation accounts for the asymmetry and broadness of the dielectric dispersion curve. The model was first used to describe the dielectric relaxation of some polymers, by adding two exponential parameters to the Debye equation:

<span class="mw-page-title-main">Duffing equation</span> Non-linear second order differential equation and its attractor

The Duffing equation, named after Georg Duffing (1861–1944), is a non-linear second-order differential equation used to model certain damped and driven oscillators. The equation is given by

Difference in differences is a statistical technique used in econometrics and quantitative research in the social sciences that attempts to mimic an experimental research design using observational study data, by studying the differential effect of a treatment on a 'treatment group' versus a 'control group' in a natural experiment. It calculates the effect of a treatment on an outcome by comparing the average change over time in the outcome variable for the treatment group to the average change over time for the control group. Although it is intended to mitigate the effects of extraneous factors and selection bias, depending on how the treatment group is chosen, this method may still be subject to certain biases.

The Cole–Cole equation is a relaxation model that is often used to describe dielectric relaxation in polymers.

<span class="mw-page-title-main">Covariant formulation of classical electromagnetism</span> Ways of writing certain laws of physics

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

<span class="mw-page-title-main">Voigt effect</span>

The Voigt effect is a magneto-optical phenomenon which rotates and elliptizes linearly polarised light sent into an optically active medium. The effect is named after the German scientist Woldemar Voigt who discovered it in vapors. Unlike many other magneto-optical effects such as the Kerr or Faraday effect which are linearly proportional to the magnetization, the Voigt effect is proportional to the square of the magnetization and can be seen experimentally at normal incidence. There are also other denominations for this effect, used interchangeably in the modern scientific literature: the Cotton–Mouton effect and magnetic-linear birefringence, with the latter reflecting the physical meaning of the effect.

The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members often asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.

In the Newman–Penrose (NP) formalism of general relativity, Weyl scalars refer to a set of five complex scalars which encode the ten independent components of the Weyl tensor of a four-dimensional spacetime.

<span class="mw-page-title-main">Viscoplasticity</span> Theory in continuum mechanics

Viscoplasticity is a theory in continuum mechanics that describes the rate-dependent inelastic behavior of solids. Rate-dependence in this context means that the deformation of the material depends on the rate at which loads are applied. The inelastic behavior that is the subject of viscoplasticity is plastic deformation which means that the material undergoes unrecoverable deformations when a load level is reached. Rate-dependent plasticity is important for transient plasticity calculations. The main difference between rate-independent plastic and viscoplastic material models is that the latter exhibit not only permanent deformations after the application of loads but continue to undergo a creep flow as a function of time under the influence of the applied load.

Free carrier absorption occurs when a material absorbs a photon, and a carrier is excited from an already-excited state to another, unoccupied state in the same band. This intraband absorption is different from interband absorption because the excited carrier is already in an excited band, such as an electron in the conduction band or a hole in the valence band, where it is free to move. In interband absorption, the carrier starts in a fixed, nonconducting band and is excited to a conducting one.

The salvo combat model provides a mathematical representation of anti-ship missile battles between modern warships. It was developed by Wayne Hughes at the U.S. Naval Postgraduate School in Monterey, California, and published in 1995. The salvo model describes the basic elements of modern missile combat in a very simple manner. This is similar to how Lanchester's square law provides a simple model of modern gun combat.

In the Newman–Penrose (NP) formalism of general relativity, independent components of the Ricci tensors of a four-dimensional spacetime are encoded into seven Ricci scalars which consist of three real scalars , three complex scalars and the NP curvature scalar . Physically, Ricci-NP scalars are related with the energy–momentum distribution of the spacetime due to Einstein's field equation.

<span class="mw-page-title-main">Relativistic Lagrangian mechanics</span> Mathematical formulation of special and general relativity

In theoretical physics, relativistic Lagrangian mechanics is Lagrangian mechanics applied in the context of special relativity and general relativity.

Response modeling methodology (RMM) is a general platform for statistical modeling of a linear/nonlinear relationship between a response variable (dependent variable) and a linear predictor (a linear combination of predictors/effects/factors/independent variables), often denoted the linear predictor function. It is generally assumed that the modeled relationship is monotone convex (delivering monotone convex function) or monotone concave (delivering monotone concave function). However, many non-monotone functions, like the quadratic equation, are special cases of the general model.

The porous medium equation, also called the nonlinear heat equation, is a nonlinear partial differential equation taking the form:

References

  1. Lanchester F.W., Mathematics in Warfare in The World of Mathematics, Vol. 4 (1956) Ed. Newman, J.R., Simon and Schuster, 2138–2157; anthologised from Aircraft in Warfare (1916)
  2. Davis, Paul K. (1995). "Lanchester Equations and Scoring Systems". Aggregation, Disaggregation, and the 3:1 Rules in Ground Combat. Rand Corporation. doi:10.7249/MR638.
  3. 1 2 Osipov, M. (1991) [1915]. "The Influence of the Numerical Strength of Engaged Forces on Their Casualties" Влияние Численности Сражающихся Сторонъ На Ихъ Потери (PDF). Tsarist Russian Journal Military CollectionВоенный Сборник. Translated by Helmbold, Robert; Rehm, Allan. US Army Concepts Analysis Agency. Archived (PDF) from the original on 4 November 2021. Retrieved 23 January 2022.
  4. Wrigge, Staffan; Fransen, Ame; Wigg, Lars (September 1995). "The Lanchester Theory of Combat and Some Related Subjects" (PDF). FORSVARETS FORSKNINGSANSTALT.
  5. Christian, MAJ Joshua T. (23 May 2019). An Examination of Force Ratios (PDF). Fort Leavenworth, KS: US Army Command and General Staff College.PD-icon.svg This article incorporates public domain material from websites or documents of the United States Army .
  6. David A. Shlapak, and Michael W. Johnson, Reinforcing Deterrence on NATO’s Eastern Flank (Santa Monica, CA: RAND Corporation, 2016)
  7. Taylor JG. 1983. Lanchester Models of Warfare, volumes I & II. Operations Research Society of America.
  8. Armstrong MJ, Sodergren SE, 2015, Refighting Pickett's Charge: mathematical modeling of the Civil War battlefield, Social Science Quarterly.
  9. MacKay N, Price C, 2011, Safety in Numbers: Ideas of concentration in Royal Air Force fighter defence from Lanchester to the Battle of Britain, History 96, 304–325.
  10. Lucas, Thomas W.; Turkes, Turker (2004). "Fitting Lanchester equations to the battles of Kursk and Ardennes: Lanchester Equations for the Battles of Kursk and Ardennes". Naval Research Logistics (NRL). 51 (1): 95–116. doi:10.1002/nav.10101. S2CID   4809135.
  11. Race to the Swift: Thoughts on Twenty-First Century Warfare by Richard E. Simpkin
  12. FOWLER, CHARLES A. "BERT" (1 March 2006). "Asymmetric Warfare: A Primer".
  13. Deitchman, S. J. (1962). "A Lanchester Model of Guerrilla Warfare". Operations Research. 10 (6): 818–827. doi:10.1287/opre.10.6.818. ISSN   0030-364X. JSTOR   168104.
  14. Clifton, E. (2020). A Brief Review on the Application of Lanchester's Models of Combat in Nonhuman Animals. Ecological Psychology, 32, 181-191. doi:10.1080/10407413.2020.1846456
  15. Wilson, M. L., Britton, N. F., & Franks, N. R. (2002). Chimpanzees and the mathematics of battle. Proceedings of the Royal Society B: Biological Sciences, 269, 1107-1112. doi:10.1098/rspb.2001.1926
  16. Lymbery, Samuel J. (2023). "Complex battlefields favor strong soldiers over large armies in social animal warfare". PNAS. 120 (37): e2217973120. Bibcode:2023PNAS..12017973L. doi:10.1073/pnas.2217973120. PMC   10500280 . PMID   37639613 . Retrieved 18 September 2023.
  17. Plowes, N. J. R., & Adams, E. S. (2005). An empirical test of Lanchester's square law: mortality during battles of the fire ant Solenopsis invicta. Proceedings of the Royal Society B: Biological Sciences, 272, 1809-1814. doi:10.1098/rspb.2005.3162
  18. Helmbold 1961a.
  19. Helmbold 1961b.
  20. Helmbold 2021, pp. app A.
  21. Helmbold 2021, pp. 14–16, app A.
  22. Helmbold 2021, pp. 18–19.
  23. Helmbold 2021, pp. 17–18.
  24. Helmbold 2021, pp. 20, 68–69.
  25. Helmbold 2021, pp. 20, app C.
  26. Helmbold 2021, pp. 21, app C part 4.
  27. Lacina, Bethany & Nils Petter Gleditsch (2005) "Monitoring Trends in Flobal Combat: A New Dataset of Battle Deaths", Journal of Population (2005) 21:145-166
  28. Lacina, Bethany, Nils Petter Gleditsch, & Bruce Russett (2006) "The Declining Risk of Death in Battle", International Studies Quyarterly 50(3), 673-680
  29. Lacina, Bethany & Nils Petter Gleditsch, (2012) Journal of Conflict Resolution 57(6) 1109-1127
  30. Lacina, Bethany & Nils Petter Gleditsch, (2012) "The Waning of War Is Real: A Response to Gohdes and Price", Journal of Conflict Resolution

Bibliography