Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “collision operators”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Kinetic theory of relativistic plasmas

The thermalization of particle kinetic motion by binary collisions is considered for a plasma with a Boltzmann constant-temperature product approximately equal to 10 to 100 times the product of the electron mass with the square of the speed of light. At this temperature, the principal mechanism for relaxation of electron motion is via radiationless electron-electron collisions (Moller scattering). Ions are nonrelativistic, but are energetic enough so that their Coulomb scattering can be treated in the Born approximation. Relaxation times are computed and Boltzmann-equation Fokker-Planck operators are derived for the various binary-collision processes. The expression for the rate of kinetic energy exchange between electron and ion gases is derived for the case where the gases are at different temperatures.

Gould, R. J.↗

Realistic Covariance Prediction For the Earth Science Constellations

Routine satellite operations for the Earth Science Constellations (ESC) include collision risk assessment between members of the constellations and other orbiting space objects. One component of the risk assessment process is computing the collision probability between two space objects. The collision probability is computed via Monte Carlo techniques as well as numerically integrating relative probability density functions. Each algorithm takes as inputs state vector and state vector uncertainty information for both objects. The state vector uncertainty information is expressed in terms of a covariance matrix. The collision probability computation is only as good as the inputs. Therefore, to obtain a collision calculation that is a useful decision-making metric, realistic covariance matrices must be used as inputs to the calculation. This paper describes the process used by NASA Goddard's Earth Science Mission Operations Project to generate realistic covariance predictions for three of the ESC satellites: Aqua, Aura, and Terra

Duncan, Matthew↗

USA Space Debris Environment, Operations, and Research Updates

Space Missions in 2017 Earth Satellite Population Collision Avoidance Maneuvers Post mission Disposal of U.S.A. Spacecraft Space Situational Awareness (SSA) and the Space Debris Sensor (SDS) A total of 86 space launches placed more than 400 spacecraft into Earth orbits during 2017, following the trend of increase over the past decade NASA has established conjunction assessment processes for its human spaceflight and uncrewed spacecraft to avoid accidental collisions with objects tracked by the U.S. Space Surveillance Network - NASA also assists other U.S. government spacecraft owners with conjunction assessments and subsequent maneuvers The ISS has conducted 25 debris collision avoidance maneuvers since 1999 - None in 2016-2017, but an ISS visiting vehicle had one collision avoidance maneuver in 2017 During 2017 NASA executed or assisted in the execution of 21 collision avoidance maneuvers by uncrewed spacecraft - Four maneuvers were conducted to avoid debris from Fengyun-1C - Two maneuvers were conducted to avoid debris from the collision of Cosmos 2251 and Iridium 33 - One maneuver was conducted to avoid the ISS NASA has established conjunction assessment processes for its human spaceflight and uncrewed spacecraft to avoid accidental collisions with objects tracked by the U.S. Space Surveillance Network - NASA also assists other U.S. government spacecraft owners with conjunction assessments and subsequent maneuvers The ISS has conducted 25 debris collision avoidance maneuvers since 1999 - None in 2016-2017, but an ISS visiting vehicle had one collision avoidance maneuver in 2017 During 2017 NASA executed or assisted in the execution of 21 collision avoidance maneuvers by uncrewed spacecraft - Four maneuvers were conducted to avoid debris from Fengyun-1C - Two maneuvers were conducted to avoid debris from the collision of Cosmos 2251 and Iridium 33 The 2014-15 NASA Engineering and Safety Center (NESC) study on the micrometeoroid and orbital debris (MMOD) assessment for the Joint Polar Satellite System (JPSS) provided the following findings - Millimeter-sized orbital debris pose the highest penetration risk to most operational spacecraft in LEO - The most effective means to collect direct measurement data on millimetersized debris above 600 km altitude is to conduct in situ measurements - There is currently no in situ data on such small debris above 600 km altitude Since the orbital debris population follows a power-law size distribution, there are many more millimeter-sized debris than the large tracked objects - Current conjunction assessments and collision avoidance maneuvers against the tracked objects (which are typically 10 cm and larger) only address a small fraction (<1%) of the mission-ending risk from orbital debris To address the millimeter-sized debris data gap above 600 km, NASA has recently developed an innovative in situ measurement instrument - the Space Debris Sensor (SDS) - One maneuver was conducted to avoid the ISS

Liou, J.-C.↗

Projection-Based Inter-Agent Collision Avoidance in Dual Agent Systems

Inter-agent collisions can occur in otherwise dynamically-stable (i.e., Lyapunov stable) dual-agent leaderfollower systems. These inter-agent collisions between the leader and the follower happen during the transient phase of the system’s evolution, although the steady-state behavior of the system is asymptotically/exponentially stable. Therefore, to avoid such inter-agent collisions, it is essential to control the relative error trajectory between the leader and the follower during the transient phase of the system’s evolution. In this paper, we introduce a novel projection operator based modelreference control architecture that can mitigate impending inter-agent collisions by modifying the transient dynamics of relative trajectories. This controller augments the follower’s baseline controller and consists of two essential components: a collision-free reference model based on the projection operator and a model reference tracking controller to guide the follower to follow the reference-model. This paper defines the concept of transient-instability in leader-follower systems, introduces collision mitigation controller architecture, and presents an illustrative example demonstrating its effectiveness.

Bayard, David S.↗

Helicopter Pilot Assessments of the Airborne Collision Avoidance System XR With Automated Maneuvering

The Airborne Collision Avoidance System X (ACAS X) is a next-generation collision avoidance tool developed to support different aircraft types and operations. The rotorcraft variant – referred to as ACAS XR – is designed to accommodate existing helicopter platforms as well as in-development, electric vertical takeoff and landing concepts, which are critical to the emerging concept of operations referred to as Advanced Air Mobility. The fundamental role of ACAS XR is to provide Detect and Avoid (DAA) and/or Collision Avoidance (CA) protection against airborne traffic. DAA alerting and guidance in the context of ACAS XR is caution-level and “suggestive,” and is to be used by the pilot if, and when, they decide to maneuver against an identified threat to DAA “well clear.” The CA alerting, by contrast, is warning-level and “directive,” with the pilot required to comply with the associated guidance to prevent a predicted Near Midair Collision (NMAC). The CA alerts generated by ACAS XR are referred to as Resolution Advisories (RAs) consistent with previous CA systems. Unlike earlier CA systems, ACAS XR issues RAs in the horizontal and vertical dimensions as well as multi-axis RAs (referred to as “Blended” RAs). According to the Minimal Operational Performance Standards of DAA systems for Unmanned Aircraft Systems, maneuvers to comply with RAs may be automated, whereas maneuvers based on DAA alerting assume a manual response. The current study was a human-in-the-loop simulation that presented rotorcraft pilots with ACAS XR alerts and guidance in a fixed-base eVTOL simulator with varying levels of automation. Objective results showed that pilots complied with all RAs within the expected 5-second time window and responded to DAA alerts quicker than in earlier studies with ACAS XU. Pilots often made larger horizontal deviations during Manual RAs, but often favored vertical and blended maneuvers. No NMACs occurred, and losses of well clear were mainly attributed to the obligation of the pilots and system to maneuver only after the CA phase of the encounter had begun. Other losses were due to pilots’ noncompliance or disregard for ACAS XR’s alerting and guidance. Noncompliance with RAs most frequently occurred when pilots determined they were too close to the terrain to continue to follow Descend RAs, performing vertical maneuvers instead of following Horizontal RAs, or rejecting Horizontal RA updates because they felt that enough maneuvering had been performed. Subjectively, pilots found the DAA and RA alerting and guidance intuitive and useful for VFR helicopter operations. Slightly more pilots preferred the Automated RA condition to the Manual RA condition. Lastly, they also felt that ACAS led to occasional unsafe Descend RAs. Caveats and future implications are discussed.

eVTOL↗

Helicopter Pilot Assessments of the Airborne Collision Avoidance System XR With Automated Maneuvering

The Airborne Collision Avoidance System X (ACAS X) is a next-generation collision avoidance tool developed to support different aircraft types and operations. The rotorcraft variant – referred to as ACAS XR – is designed to accommodate existing helicopter platforms as well as in-development, electric vertical takeoff and landing concepts, which are critical to the emerging concept of operations referred to as Advanced Air Mobility. The fundamental role of ACAS XR is to provide Detect and Avoid (DAA) and/or Collision Avoidance (CA) protection against airborne traffic. DAA alerting and guidance in the context of ACAS XR is caution-level and “suggestive,” and is to be used by the pilot if, and when, they decide to maneuver against an identified threat to DAA “well clear.” The CA alerting, by contrast, is warning-level and “directive,” with the pilot required to comply with the associated guidance to prevent a predicted Near Midair Collision (NMAC). The CA alerts generated by ACAS XR are referred to as Resolution Advisories (RAs) consistent with previous CA systems. Unlike earlier CA systems, ACAS XR issues RAs in the horizontal and vertical dimensions as well as multi-axis RAs (referred to as “Blended” RAs). According to the Minimal Operational Performance Standards of DAA systems for Unmanned Aircraft Systems, maneuvers to comply with RAs may be automated, whereas maneuvers based on DAA alerting assume a manual response. The current study was a human-in-the-loop simulation that presented rotorcraft pilots with ACAS XR alerts and guidance in a fixed-base eVTOL simulator with varying levels of automation. Objective results showed that pilots complied with all RAs within the expected 5-second time window and responded to DAA alerts quicker than in earlier studies with ACAS XU. Pilots often made larger horizontal deviations during Manual RAs, but often favored vertical and blended maneuvers. No NMACs occurred, and losses of well clear were mainly attributed to the obligation of the pilots and system to maneuver only after the CA phase of the encounter had begun. Other losses were due to pilots’ noncompliance or disregard for ACAS XR’s alerting and guidance. Noncompliance with RAs most frequently occurred when pilots determined they were too close to the terrain to continue to follow Descend RAs, performing vertical maneuvers instead of following Horizontal RAs, or rejecting Horizontal RA updates because they felt that enough maneuvering had been performed. Subjectively, pilots found the DAA and RA alerting and guidance intuitive and useful for VFR helicopter operations. Slightly more pilots preferred the Automated RA condition to the Manual RA condition. Lastly, they also felt that ACAS led to occasional unsafe Descend RAs. Caveats and future implications are discussed.

eVTOL↗

Realistic Covariance Prediction for the Earth Science Constellation

Routine satellite operations for the Earth Science Constellation (ESC) include collision risk assessment between members of the constellation and other orbiting space objects. One component of the risk assessment process is computing the collision probability between two space objects. The collision probability is computed using Monte Carlo techniques as well as by numerically integrating relative state probability density functions. Each algorithm takes as inputs state vector and state vector uncertainty information for both objects. The state vector uncertainty information is expressed in terms of a covariance matrix. The collision probability computation is only as good as the inputs. Therefore, to obtain a collision calculation that is a useful decision-making metric, realistic covariance matrices must be used as inputs to the calculation. This paper describes the process used by the NASA/Goddard Space Flight Center's Earth Science Mission Operations Project to generate realistic covariance predictions for three of the Earth Science Constellation satellites: Aqua, Aura and Terra.

Duncan, Matthew↗

Aerial Object Trajectory Classification by Training on Flight Controller Data and Testing on RADAR Generated Tracks

Onboard collision avoidance is needed to enable safe, autonomous flight operations for NASA projects such as Advanced Air Mobility (AAM), as well as many commercial applications. Real-time aerial object classification will improve onboard collision avoidance algorithm decision making and may reduce unnecessary activation of avoidance systems. This work trains an aircraft trajectory classifier using trajectories from flight controller logs and tests the classifier using RADAR collected trajectories during air to air experiments and ground to air experiments. In contrast to RADAR data, these flight controller logs are relatively abundant, which makes the possibility of substituting flight data for RADAR data an attractive, cost-effective option. The SVM model developed in this work achieved a 79.7% classification accuracy on the first second of radar trajectories of GA, multirotor sUAS, and fixed wing sUAS. Findings from this work indicate that it is feasible to classify sensor collected trajectories using a classifier trained on flight controller data.

Henry Holbrook↗

DataSet: Aerial Object Trajectory Classification by Training on Flight Controller Data and Testing on RADAR Generated

Onboard collision avoidance is needed to enable safe, autonomous flight operations for NASA projects such as Advanced Air Mobility (AAM), as well as many commercial applications. Real-time aerial object classification will improve onboard collision avoidance algorithm decision making and may reduce unnecessary activation of avoidance systems. This work trains an aircraft trajectory classifier using trajectories from flight controller logs and tests the classifier using RADAR collected trajectories during air to air experiments and ground to air experiments. In contrast to RADAR data, these flight controller logs are relatively abundant, which makes the possibility of substituting flight data for RADAR data an attractive, cost-effective option. The SVM model developed in this work achieved a 79.7% classification accuracy on the first second of radar trajectories of GA, multirotor sUAS, and fixed wing sUAS. Findings from this e that it is feasible to classify sensor collected trajectories using a classifier trained flight controller data.

Chester V Dolph↗

Deep-Space Conjunction Assessment: Recent Developments and Future Evolution

The Multi-mission Automated Deep-space Conjunction Assessment Process (MADCAP) is a NASA Jet Propulsion Laboratory (JPL) capability used to perform conjunction assessment in shared deep-space environments. MADCAP began performing conjunction assessment at Mars and the Moon in 2011, with the Sun/Earth libration points added to its functionality in 2020. There has been an increasing number of missions operating in these environments in recent years, leading to an elevated frequency of close conjunction events, especially in the Lunar orbital environment. MADCAP provides this service not only to NASA missions, but to any operator who is willing to share ephemerides. Since there is no space surveillance network for deep space environments, ephemeris sharing is the only way in which spacecraft operators can ensure the safety of their spacecraft from collision in these orbit regimes. NASA published a set of conjunction assessment best practices in 2020 that cover the MADCAP process. This paper details recent MADCAP operational experience in the deep space environments, including statistics and process improvements. Updates to the MADCAP software and automation framework implemented to handle the recent growth in the number of deep space missions are also discussed. Future enhancements planned in anticipation of increasingly crowded deep-space environments, such as non-standard runs based on exploratory scenarios, are also discussed.

conjunction assessment↗

Development of a Nonlinear Probability of Collision Tool for the Earth Observing System

The Earth Observing System (EOS) spacecraft Terra, Aqua, and Aura fly in constellation with several other spacecraft in 705-kilometer mean altitude sun-synchronous orbits. All three spacecraft are operated by the Earth Science Mission Operations (ESMO) Project at Goddard Space Flight Center (GSFC). In 2004, the ESMO project began assessing the probability of collision of the EOS spacecraft with other space objects. In addition to conjunctions with high relative velocities, the collision assessment method for the EOS spacecraft must address conjunctions with low relative velocities during potential collisions between constellation members. Probability of Collision algorithms that are based on assumptions of high relative velocities and linear relative trajectories are not suitable for these situations; therefore an algorithm for handling the nonlinear relative trajectories was developed. This paper describes this algorithm and presents results from its validation for operational use. The probability of collision is typically calculated by integrating a Gaussian probability distribution over the volume swept out by a sphere representing the size of the space objects involved in the conjunction. This sphere is defined as the Hard Body Radius. With the assumption of linear relative trajectories, this volume is a cylinder, which translates into simple limits of integration for the probability calculation. For the case of nonlinear relative trajectories, the volume becomes a complex geometry. However, with an appropriate choice of coordinate systems, the new algorithm breaks down the complex geometry into a series of simple cylinders that have simple limits of integration. This nonlinear algorithm will be discussed in detail in the paper. The nonlinear Probability of Collision algorithm was first verified by showing that, when used in high relative velocity cases, it yields similar answers to existing high relative velocity linear relative trajectory algorithms. The comparison with the existing high velocity/linear theory will also be used to determine at what relative velocity the analysis should use the new nonlinear theory in place of the existing linear theory. The nonlinear algorithm was also compared to a known exact solution for the probability of collision between two objects when the relative motion is strictly circular and the error covariance is spherically symmetric. Figure I shows preliminary results from this comparison by plotting the probabilities calculated from the new algorithm and those from the exact solution versus the Hard Body Radius to Covariance ratio. These results show about 5% error when the Hard Body Radius is equal to one half the spherical covariance magnitude. The algorithm was then combined with a high fidelity orbit state and error covariance propagator into a useful tool for analyzing low relative velocity nonlinear relative trajectories. The high fidelity propagator is capable of using atmospheric drag, central body gravitational, solar radiation, and third body forces to provide accurate prediction of the relative trajectories and covariance evolution. The covariance propagator also includes a process noise model to ensure realistic evolutions of the error covariance. This paper will describe the integration of the nonlinear probability algorithm and the propagators into a useful collision assessment tool. Finally, a hypothetical case study involving a low relative velocity conjunction between members of the Earth Observation System constellation will be presented.

McKinley, David P.↗

Module Testing Procedures and Results for the CMS Phase-2 Outer Tracker Upgrade at Fermilab

The High-Luminosity LHC (HL-LHC) will operate at significantly increased luminosities and is expected to deliver about 3000 $fb^{-1}$ of proton-proton collision data at $\sqrt{s}=14$ TeV over a decade of operation. To maintain efficient tracking and triggering performance under high pileup and radiation conditions, the CMS experiment is upgrading its tracking detector for Phase-2 operations. The upgraded Outer Tracker (OT) will consist of Pixel-Strip (PS) and Strip-Strip (2S) silicon modules capable of providing tracking information to the Level-1 trigger at 40 MHz. Production and qualification of OT modules have been ongoing for about one year across several assembly and testing centers in the US, Europe, India, and Pakistan, requiring extensive testing to ensure stable operation in the HL-LHC environment. Fermilab is responsible for the production and testing of a significant fraction of the OT modules. The testing activities include IV characterization of silicon sensors, noise and pedestal measurements, verification of communication between module components, and burn-in studies using cold-box systems operated under controlled thermal conditions. Results from module testing and qualification studies performed during production will be presented.

Baradia, Sweta [UC, Davis (main)]↗

Capabilities of an Acoustic Camera to Inform Fish Collision Risk with Current Energy Converter Turbines

A diversified energy portfolio may include marine energy in the form of current energy converters (CECs) such as tidal or in-river turbines. New technology development in the research stage typically requires monitoring for environmental effects. A significant environmental effect of concern for CECs is the risk of moving parts (e.g., turbine blades) colliding with animals such as fishes. CECs are installed in energetic locations in which it is difficult to operate sensors to fulfill monitoring requirements for informing collision risk. Collecting data (i.e., about blade strikes or near-misses) that inform interactions of fishes with CECs is usually attempted using active acoustic sensors or video cameras (VCs). Limitations of low-light conditions or water turbidity that preclude effective use of VCs are overcome by using high-resolution multibeam echosounders (or acoustic cameras (ACs)). We used an AC at two sites to test its ability to detect artificial and real fish targets and determine if strike, near-miss, and near-field behavior could be observed. Interactions with fish and artificial targets with turbines have been documented but strike confirmation with an AC is novel. The first site was in a tidal estuary with a 25 kW turbine and water clarity sufficient to allow VC data to be collected concurrently with AC data showing turbine blade strike on tethered artificial fish targets. The second site was a turbid, debris-laden river with a 5 kW turbine where only AC data were collected due to high water turbidity. Data collection at the second site coincided with downstream Pacific salmon (Oncorhynchus spp.) smolt migration. Physical fish capture downstream of the turbine was performed with an incline plane trap (IPT) to provide context for the AC observations, by comparing fish catches. Discrimination between debris and fishes in the AC data was not possible, because active movement of fishes was not discernable. Nineteen fishes were released upstream of the turbine to provide known times of possible fish/turbine interactions, but detection was difficult to confirm in the AC data. ACs have been used extensively in past studies to count large migratory fish such as Pacific salmon, but their application for small fish targets has been limited. The results from these two field campaigns demonstrate the ability of ACs to detect targets in turbid water and observe blade strikes, as well as their limitations such as the difficulty of distinguishing small fishes from debris in a high-energy turbid river. Recommendations are presented for future applications associated with CEC device testing.

16 TIDAL AND WAVE POWER↗

Measurements of Born cross sections for e + e − → Λ c + Λ ¯ c ( 2595 ) − + c . c . and e + e − → Λ c + Λ ¯ c ( 2625 ) − + c . c . at s = 4918.0 and 4950.9 MeV

Using e + e − collision data collected with the BESIII detector operating at the BEPCII collider, the Born cross sections of e + e − → Λ c + Λ ¯ c ( 2595 ) − + c . c . and e + e − → Λ c + Λ ¯ c ( 2625 ) − + c . c . are measured for the first time at center-of-mass energies of s = 4918.0 and 4950.9 MeV. Nonzero cross sections are observed very close to the production threshold. The measured Born cross sections of e + e − → Λ c + Λ ¯ c ( 2625 ) − + c . c . are about 2–3 times greater than those of e + e − → Λ c + Λ ¯ c ( 2595 ) − + c . c . , providing the similar behavior as semileptonic decays of Λ b 0 , but different behavior from that in the hadronic decays of Λ b 0 . The Born cross sections are 15.6 ± 3.1 ± 0.9 pb and 29.4 ± 3.7 ± 2.7 pb for e + e − → Λ c + Λ ¯ c ( 2595 ) − + c . c . , and are 43.4 ± 4.0 ± 4.1 pb and 76.8 ± 6.5 ± 4.2 pb for e + e − → Λ c + Λ ¯ c ( 2625 ) − + c . c . at s = 4918.0 and 4950.9 MeV, respectively. Based on the polar angle distributions of the Λ ¯ c ( 2625 ) − and Λ c ( 2625 ) + , the form-factor ratios | G E | 2 + 3 | G M | 2 / | G C | are determined for e + e − → Λ c + Λ ¯ c ( 2625 ) − + c . c . for the first time, which are 5.95 ± 4.07 ± 0.15 and 0.94 ± 0.32 ± 0.02 at s = 4918.0 and 4950.9 MeV, respectively. All of these first uncertainties are statistical and second systematic. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Search for C -even states decaying to D s ± D s * ∓ with masses between 4.08 and 4.32 GeV / c 2

Six C -even states, denoted as X , with quantum numbers J P C = 0 − + , 1 ± + , or 2 ± + , are searched for via the e + e − → γ D s ± D s * ∓ process using ( 1667.39 ± 8.84 ) pb − 1 of e + e − collision data collected with the BESIII detector operating at the BEPCII storage ring at center-of-mass energy of s = ( 4681.92 ± 0.30 ) MeV . No statistically significant signal is observed in the mass range from 4.08 GeV / c 2 to 4.32 GeV / c 2 . The upper limits of σ [ e + e − → γ X ] · B [ X → D s ± D s * ∓ ] at a 90% confidence level are determined. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Precise Measurement of the e + e − → D s + D s − Cross Section at Center-of-Mass Energies from Threshold to 4.95 GeV

Using the e + e − collision data collected with the BESIII detector operating at the BEPCII collider, at center-of-mass energies from the threshold to 4.95 GeV, we present precise measurements of the cross section for the process e + e − → D s + D s − using a single-tag method. The resulting cross section line shape exhibits several new structures, thereby offering an input for a future coupled-channel analysis and model tests, which are critical to understand vector charmonium-like states with masses between 4 and 5 GeV. Published by the American Physical Society 2024

Ablikim, M.↗

SNOWMASS2021 (AF7): Sub-scale models and infrastructure for accelerator magnet development

The LHC is undoubtfully the most important instrument the HEP community has on its disposal and is widely considered an amazing technological achievement. At the core of this instrument lay accelerator magnets based on NbTi technology – a technology well developed and widely used in the commercial sector too. Yet, a decade after the accelerator started it still did not reach its design collision energy of 14 TeV and it operated for years well below it. The reason for this is simple – magnet training (even though only a subset of the magnets is mostly affected).

43 PARTICLE ACCELERATORS↗

The generalized diffusion-convection equation

Starting from the Boltzmann equation, a transport equation is derived for energetic particles in a moving magnetized plasma in which the scattering centers that keep the particles quasi-isotropic are moving with a velocity that is not necessarily the same as that of the plasma. The scattering is characterized by three very loose constraints: (1) there is a rest frame for each scatterer in which the particles scatter elastically; (2) in this frame the scattering will not disturb an isotropic distribution; and (3) the momentum transfer in an average collision may be described by a tensor operating on the particles original momentum. Since the strength of the scattering is not specified, the derivation should be as valid for plasma microturbulence as for hard-sphere scattering. The results show clearly which phenomena are responsible for tying the particles to the plasma in the transport equation.

Jones, Frank C.↗