Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle tracking”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Analytical Dynamic Aperture Calculation Compared to Tracking Results

In studying dynamic aperture in the storage rings, it is becoming increasing beneficial to understand the growth of the fractal boundary of the dynamic aperture without tracking for millions of turns in a hadron ring. This note describes the J. Gao theory applied to the Electron Ion Collider (EIC) Hadron Storage Ring (HSR) and compares simulation results to theoretical predictions to show that the analytical method may prove useful in quick determinations of improved dynamic aperture based on sextupole and octupole configurations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Beam test performance of AstroPix sensor with 120 GeV protons

AstroPix is a High-Voltage CMOS Monolithic Active Pixel Sensor (HV-CMOS MAPS) developed for precision gamma-ray imaging and spectroscopy in the medium-energy regime, as well as for precise shower imaging and tracking in the Barrel Imaging Calorimeter (BIC) of the Electron Proton/Ion Collider (ePIC) detector at the future Electron–Ion Collider (EIC). We present beam test results of the AstroPix_v3 sensor using a 120 GeV proton beam at the Fermilab Test Beam Facility (FTBF), performed as part of the broader experimental campaign for the BIC prototype calorimeter. The sensor’s 500 µm pixel pitch enabled precise measurement of the beam profile, providing important information for the calorimeter performance studies. Using the measured 120 GeV proton data, we measure the energy deposit of minimum ionizing particles (MIP) and use them to extract the corresponding effective depletion depth at a single bias voltage of −150 V .

AstroPix↗
DOE OSTI2025

Intelligent experiments through real-time AI: Fast Data Processing and Autonomous Detector Control for sPHENIX and future EIC detectors

This R&D project, initiated by the DOE Nuclear Physics AI-Machine Learning initiative in 2022, leverages AI to address data processing challenges in high-energy nuclear experiments (RHIC, LHC, and future EIC). Our focus is on developing a demonstrator for real-time processing of high-rate data streams from sPHENIX experiment tracking detectors. The limitations of a 15 kHz maximum trigger rate imposed by the calorimeters can be negated by intelligent use of streaming technology in the tracking system. The approach efficiently identifies low momentum rare heavy flavor events in high-rate p+p collisions (3MHz), using Graph Neural Network (GNN) and High Level Synthesis for Machine Learning (hls4ml). Success at sPHENIX promises immediate benefits, minimizing resources and accelerating the heavy-flavor measurements. The approach is transferable to other fields. For the EIC, we develop a DIS-electron tagger using Artificial Intelligence - Machine Learning (AI-ML) algorithms for real-time identification, showcasing the transformative potential of AI and FPGA technologies in high-energy nuclear and particle experiments real-time data processing pipelines.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Optimisation of the Search for CP-symmetry Violation at the Deep Underground Neutrino Experiment

The Deep Underground Neutrino Experiment (DUNE) is a next-generation long baseline experiment, which will be situated in South Dakota. Its detectors will utilise liquid-argon time projection chamber technology, which is able to capture neutrino interactions with an incredible spatial and calorimetric resolution. With what will become the world’s most intense neutrino beam, a highly capable near detector, and four (10kt fducial mass) far detector modules, DUNE will be able to achieve an ambitious physics programme. Most notably, DUNE will determine whether charge-parity symmetry is broken in neutrino oscillations - a finding that would have significant implications for the understanding of the matter-antimatter asymmetry in our Universe.This thesis presents the optimisation of a CP-violation analysis at DUNE using thePandora pattern-recognition software. The analysis assumes a 3.5 year exposure (1.36 × 1023 protons on target) to a neutrino and an antineutrino beam (7 year total). Only the predicted data of the far detector modules is used; near detector samples are not included. The initial sensitivity to CP-violation is found to be 3.8σ+0.9σ−1.1σ in an estimate that includes oscillation parameter uncertainties, systematic uncertainties and statistical fluctuations, and assumes a normal-ordering of the neutrino mass hierarchy. The performance of the Pandora event reconstruction is linked to that of the analysis, which is found to be limited by the reconstruction of the initial track-like region of electrons and photons. The ShowerRefinement algorithm is developed in response to this and its implementation into the analysis workflow results in an improved sensitivity to CP-violation of 4.6σ+0.9σ−1.0σ. With a perfected neutrino interaction vertex placement, this is further increased to 5.1σ+1.0σ−1.1σ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

A thermal-driven graupel generation process to explain dry-season convective vigor over the Amazon

Large-eddy simulations (LESs) are conducted for each day of the intensive observation periods (IOPs) of the Green Ocean Amazon (GoAmazon) field campaign to characterize the updrafts and microphysics within deep convective cores while contrasting those properties between Amazon wet and dry seasons. Mean Doppler velocity (V dop ) values simulated using LESs are compared with 2-year measurements from a radar wind profiler (RWP) as viewed by statistical composites separated according to wet- and dry-season conditions. In the observed RWP and simulated LES V dop composites, we find more intense low-level updraft velocity, vigorous graupel generation, and intense surface rain during the dry periods compared with the wet periods. To investigate coupled updraft–microphysical processes further, single-day golden cases are selected from the wet and dry periods to conduct detailed cumulus thermal tracking analysis. Tracking analysis reveals that simulated dry-season environments generate more droplet-loaded low-level thermals than wet-season environments. This tendency correlates with seasonal contrasts in buoyancy and vertical moisture advection profiles in large-scale forcing. Employing a normalized time series of mean thermal microphysics, the simulated cumulus thermals appear to be the primary generator of cloud droplets. When subsequent thermals penetrate the ice crystal layer, droplets within the thermals interact with entrained ice crystals, which enhances riming in the thermals. This appears to be a production pathway of graupel/hail particles within simulated deep convective cores. In addition, less-diluted dry-case thermals tend to be elevated higher, and graupel grows further during sedimentation after spilling out from thermals. Therefore, greater concentrations of low-level moist thermals likely result in more graupel/hail production and associated dry-season convective vigor.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Impact of Time-Dependent Reactor and Sensor Physics on Core Power Synthesis (Rev.1)

Online synthesis of power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins and to provide essential knowledge associated with the burnup of the fuel. In light-water reactors, power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report examines the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a pressurized water reactor (PWR), considering the typical low-enriched uranium (LEU, 3%-5% enrichment) fuel cycle as well as the higher enrichment LEU+ (5%-8% enrichment) fuel cycle. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative PWR LEU core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used as well to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori power distribution assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters and the extended nature of the LEU+ fuel cycle in comparison with LEU. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly underapproximated or overapproximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumptions leads to severe levels of error, approaching 25% at maximum (for LEU). The authors also discovered that V SPNDs are extraordinarily robust in both the LEU and LEU+ fuel cycles considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Search for dark matter produced in association with a Standard Model Higgs boson decaying into b -quarks using the full Run 2 dataset from the ATLAS detector

The production of dark matter in association with Higgs bosons is predicted in several extensions of the Standard Model. An exploration of such scenarios is presented, considering final states with missing transverse momentum and b-tagged jets consistent with a Higgs boson. The analysis uses proton-proton collision data at a centre-of-mass energy of 13 TeV recorded by the ATLAS experiment at the LHC during Run 2, amounting to an integrated luminosity of 139 fb -1 . The analysis, when compared with previous searches, benefits from a larger dataset, but also has further improvements providing sensitivity to a wider spectrum of signal scenarios. These improvements include both an optimised event selection and advances in the object identification, such as the use of the likelihood-based significance of the missing transverse momentum and variable-radius track-jets. No significant deviation from Standard Model expectations is observed. Limits are set, at 95% confidence level, in two benchmark models with two Higgs doublets extended by either a heavy vector boson Z' or a pseudoscalar singlet a and which both provide a dark matter candidate χ. In the case of the two-Higgs-doublet model with an additional vector boson Z', the observed limits extend up to a Z' mass of 3 TeV for a mass of 100 GeV for the dark matter candidate. The two-Higgs-doublet model with a dark matter particle mass of 10 GeV and an additional pseudoscalar a is excluded for masses of the a up to 520 GeV and 240 GeV for tan β = 1 and tan β = 10 respectively. Limits on the visible cross-sections are set and range from to 0.05 fb to 3.26 fb, depending on the missing transverse momentum and b-quark jet multiplicity requirements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Search for Light Pseudoscalar Bosons, Pair-Produced in Higgs Boson Decays in the Four-Electron Final State in Proton-Proton Collisions at $\sqrt{s}=13$ TeV

A search for pairs of light neutral pseudoscalar bosons (𝐴) resulting from the decay of a Higgs boson is performed. The search is conducted using LHC proton-proton collision data at $\sqrt{s}=13$ TeV, collected with the CMS detector in 2016–2018 and corresponding to an integrated luminosity of 138 fb −1 . The 𝐴 boson decays into a highly collimated electron-positron pair. A novel multivariate algorithm using tracks and calorimeter information is developed to identify these distinctive signatures, and events are selected with two such merged electron-positron pairs. No significant excess above the standard model background predictions is observed. Upper limits on the branching fraction for 𝐻 → 𝐴⁢𝐴 → 4⁢𝑒 are set at 95% confidence level, for masses between 10 and 100 MeV and proper decay lengths below 100 μ⁢m, reaching branching fraction sensitivities as low as 10 −5 . This is the first search for Higgs boson decays to four electrons via light pseudoscalars at the LHC. It significantly improves the experimental sensitivity to axionlike particles with masses below 100 MeV.

Hayrapetyan, A. [Yerevan Physics Institute]↗
DOE OSTI2021

Tonopah Test Range Air Monitoring: CY2020 Meteorological, Radiological, and Wind Transported Particulate Observations

In 1963, the U.S. Department of Energy (DOE) (formerly the Atomic Energy Commission [AEC]), implemented Operation Roller Coaster on the Tonopah Test Range (TTR) and an adjacent area of the Nevada Test and Training Range (NTTR) (formerly the Nellis Air Force Range). This operation resulted in radionuclide-contaminated soils at the Double Tracks site and Clean Slate I, II, and III sites. This report documents observations made during ongoing monitoring of radiological, meteorological, and dust conditions at stations installed adjacent to the Clean Slate sites and at the TTR Sandia National Laboratories (SNL) Range Operations Center (ROC). The monitoring effort supports the DOE Office of Legacy Management in its stewardship of the Clean Slate sites as part of their mission to protect human health and the environment at sites of former nuclear device development and testing. Specific monitoring objectives include determining if wind blowing across the Clean Slate sites is transporting particles of radionuclide-contaminated soil beyond the physical and administrative boundaries of the sites, and providing information for planning long-term stewardship of the sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Unravelling crystal growth of nanoparticles

Crystal growth from nanoscale constituents is a ubiquitous phenomenon in biology, geology and materials science. Numerous studies have focused on understanding the onset of nucleation and on producing high-quality crystals by empirically sampling constituents with different attributes and varying the growth conditions. However, the kinetics of post-nucleation growth processes, an important determinant of crystal morphology and properties, have remained underexplored due to experimental challenges associated with real-space imaging at the nanoscale. Here we report the imaging of the crystal growth of nanoparticles of different shapes using liquid-phase transmission electron microscopy, resolving both lateral and perpendicular growth of crystal layers by tracking individual nanoparticles. We observe that these nanoscale systems exhibit layer-by-layer growth, typical of atomic crystallization, as well as rough growth prevalent in colloidal systems. Surprisingly, the lateral and perpendicular growth modes can be independently controlled, resulting in two mixed crystallization modes that, until now, have received only scant attention. Combining analytical considerations with molecular dynamics and kinetic Monte Carlo simulations, we develop a comprehensive framework for our observations, which are fundamentally determined by the size and shape of the building blocks. Furthermore, these insights unify the understanding of crystal growth across four orders of magnitude in particle size and suggest novel pathways to crystal engineering.

Imaging techniques↗
DOE OSTI2025

Wire Plane Transparency: Commissioning the Time Projection Chambers for the Short-Baseline Near Detector

The Short Baseline Near Detector (SBND) serves as the near detector for Fermilab's Short Baseline Neutrino (SBN) programme. It is a 112-ton Liquid Argon Time Projection Chamber (LArTPC) designed to study neutrino-argon interactions and search for new physics phenomena such as sterile neutrinos. Situated just 110 m from the Booster Neutrino Beam (BNB), SBND just began its first physics run and is expected to record over two million neutrino interactions annually. Commissioning the time projection chambers of SBND is crucial for optimising detector performance. Each of the two anodes is composed of three charge-sensing wire planes which detect drifting ionisation electrons from charged particles traversing the detector. The wire planes, arranged with different 2D orientations, facilitate the 3D reconstruction of neutrino interactions. Electron transport induces signals on the first two planes, with the signal being collected by the rear plane. The first two induction planes require transparency to ensure all drift electrons reach the collection plane. This transparency depends on the bias voltages applied to the wires and the spacing between them. Fine-tuning these bias voltages ensures the optimisation of wire plane transparency. In this poster, I will demonstrate how transparency can be assessed by analysing the waveforms from cosmic-ray muon tracks traversing the detector at specific orientations.

McCusker, Bethany [Unlisted, US]↗
DOE OSTI2022

A Forward Feedback Control Scheme for a Solar Thermochemical Moving Bed Counter-Current Flow Reactor

Abstract Pelletized thermochemical energy storage media has a potential for long-duration energy storage. Production of solid-state energy storage media can be done within a cavity chemical reactor that captures concentrated solar radiation from a solar thermal field. The temperature stability of a solar reactor is directly influenced by the solar flux intercepted. This paper presents a low-order physical model to simulate the dynamic response of temperature inside a tubular plug-flow reactor prototype. Solid granular particles are fed to the reactor from the top whereas a counter-current flowing gas enters the reactor from the bottom. An in-house code was developed to model transient heat transfer of the reactor wall, gas, and moving particles. The model was preliminarily validated with packed beds for different temperature ranges and two gas flowrates. Dynamic response of the reactor temperature is simulated for different input power and gas/particle flowrates. The results show that the system response can be controlled efficiently by utilizing input power (solar flux) as a control parameter. A conventional proportional integral (PI) controller is designed to control the temperature inside the reactor and to maintain it during the solar flux intermittency. The controller parameters are tuned using the Ziegler–Nichols method to ensure optimal system response. The results show that the feedback control model is successful in tracking different reference reactor temperatures within a reasonable settling time of 30 min and eliminated overshoot. This study can be extended to include a hybrid reactor with a multi-input, multi-output variable system.

Energy & Fuels↗
DOE OSTI2024

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗
DOE OSTI2021

Beyond 4D tracking: using cluster shapes for track seeding

Tracking is one of the most time consuming aspects of event reconstruction at the Large Hadron Collider (LHC) and its high-luminosity upgrade (HL-LHC). Innovative detector technologies extend tracking to four-dimensions by including timing in the pattern recognition and parameter estimation. However, present and future hardware already have additional information that is largely unused by existing track seeding algorithms. The shape of pixel-clusters provides an additional dimension for track seeding that can significantly reduce the combinatorial challenge of track finding. We use neural networks to show that cluster shapes can reduce significantly the rate of fake combinatorical backgrounds while preserving a high efficiency. We demonstrate this using the information in cluster singlets, doublets and triplets. Numerical results are presented with simulations from the TrackML challenge.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Ship-based measurements of ice nuclei concentrations over the Arctic, Atlantic, Pacific and Southern oceans

Ambient concentrations of ice-forming particles measured during ship expeditions are collected and summarised with the aim of determining the spatial distribution and variability in ice nuclei in oceanic regions. The presented data from literature and previously unpublished data from over 23 months of ship-based measurements stretch from the Arctic to the Southern Ocean and include a circumnavigation of Antarctica. In comparison to continental observations, ship-based measurements of ambient ice nuclei show 1 to 2 orders of magnitude lower mean concentrations. To quantify the geographical variability in oceanic areas, the concentration range of potential ice nuclei in different climate zones is analysed by meridionally dividing the expedition tracks into tropical, temperate and polar climate zones. We find that concentrations of ice nuclei in these meridional zones follow temperature spectra with similar slopes but vary in absolute concentration. Typically, the frequency with which specific concentrations of ice nuclei are observed at a certain temperature follows a log-normal distribution. A consequence of the log-normal distribution is that the mean concentration is higher than the most frequently measured concentration. Finally, the potential contribution of ship exhaust to the measured ice nuclei concentration on board research vessels is analysed as function of temperature. We find a sharp onset of the influence at approximately -36 °C but none at warmer temperatures that could bias ship-based measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

H⁰ Stark Stripping and Component Irradiation in Fermilab Booster

In foil stripping of H⁻ some fraction of the emerging neutral H⁰ will be in excited states, which can then strip through the Stark effect in the magnetic field of the downstream orbit bump magnet. The resultant H⁺ will experience a depleted net kick compared to protons emerging from the foil and will track on trajectories different from the nominal circulating beam. This will lead to irradiation of downstream machine components. An analysis of these processes is of particular importance looking forward to the much higher beam power of the Fermilab PIP-II era. This study investigates where these errant protons will be lost, how much power is deposited, and whether this will be a shielding concern.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

SSAPy - Space Situational Awareness for Python

SSAPy is a fast and flexible orbit modeling and analysis tool for orbits spanning from low-Earth into the cislunar regime. Orbits can be flexibly specified from common input formats such as Keplerian elements or two-line element (TLE) data files. SSAPy allows users to model satellites and specify parameters such as satellite area, mass, and drag coefficients. SSAPy includes a customizable force-propagation with a range of Earth, Lunar, radiation, atmospheric, and maneuvering models. SSAPy makes use of various community integration methods and can calculate time-evolved orbital quantities, including satellite magnitudes and state vectors. Users can specify various space- and ground-based observation models with support for multiple coordinate and reference frames. SSAPy also supports orbit analysis and propagation methods such as multiple hypothesis tracking and has built-in uncertainty quantification. The majority of SSAPy’s methods are vectorized and parallelizable, allowing for effective use of high-performance computer (HPC) systems. Finally, SSAPy has plotting functionality, allowing users to visualize orbits and trajectories. Examples are shown in Figure 1 and Figure 2.

97 MATHEMATICS AND COMPUTING↗