Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “consistent scheme”

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,261 records · Page 70

Numerical investigations of low-density nozzle flow by solving the Boltzmann equation

A two-dimensional finite-difference code to solve the BGK-Boltzmann equation has been developed. The solution procedure consists of three steps: (1) transforming the BGK-Boltzmann equation into two simultaneous partial differential equations by taking moments of the distribution function with respect to the molecular velocity u(sub z), with weighting factors 1 and u(sub z)(sup 2); (2) solving the transformed equations in the physical space based on the time-marching technique and the four-stage Runge-Kutta time integration, for a given discrete-ordinate. The Roe's second-order upwind difference scheme is used to discretize the convective terms and the collision terms are treated as source terms; and (3) using the newly calculated distribution functions at each point in the physical space to calculate the macroscopic flow parameters by the modified Gaussian quadrature formula. Repeating steps 2 and 3, the time-marching procedure stops when the convergent criteria is reached. A low-density nozzle flow field has been calculated by this newly developed code. The BGK Boltzmann solution and experimental data show excellent agreement. It demonstrated that numerical solutions of the BGK-Boltzmann equation are ready to be experimentally validated.

Deng, Zheng-Tao↗

Impact of Additive Manufacturing on Internal Cooling Channels With Varying Diameters and Build Directions

Abstract The use of additive manufacturing (AM) processes, such as direct metal laser sintering, provides the design freedom required to incorporate complex cooling schemes in gas turbine components. Additively manufactured turbine components have a range of cooling feature sizes and, because of the inherent three-dimensionality, a wide range of build angles. Previous studies have shown that AM built directions influence internal channel surface roughness that, in turn, augment heat transfer and pressure loss. This study investigates the impact of AM on channel feature size and builds direction relative to tolerance, surface roughness, pressure losses, and convective cooling. Multiple AM coupons were built from Inconel 718 consisting of channels with different diameters and a variety of build directions. An experimental rig was used to measure pressure drop to calculate friction factor and was used to impose a constant surface temperature boundary condition to collect Nusselt number over a range of Reynolds numbers. Significant variations in surface roughness and geometric deviations from the design intent were observed for distinct build directions and channel sizes. These differences led to notable impacts in friction factor and Nusselt number augmentations, which were a strong function of build angle.

Engineering↗

Pre-Ionization Controlled Laser Plasma Formation for Ignition Applications (Final Scientific/Technical Report)

This report consists of three grand funded papers. Abstracts for each follow. 1) In this work, we present experimental and modeling studies of air pre-ionization using ultraviolet (UV) laser pulses and its effect on laser breakdown of an overlapped near-infrared (NIR) pulse. Experimental studies are conducted with a 266 nm beam (fourth harmonic of Nd:YAG) for UV pre-ionization and an overlapped 1064 nm NIR beam (fundamental of Nd:YAG), both having pulse duration of 10 ns. Results show that the UV beam produces a pre-ionized volume which assists in breakdown of the NIR beam, leading to reduction in NIR breakdown threshold by factor of >2. Numerical modeling is performed to examine the ionization and breakdown of both beams. The modeled breakdown threshold of the NIR, including assist by pre-ionization, is in reasonable agreement with the experimental results. 2) The present contribution compares the energy absorption, optical emission, temperature, and fluid dynamics of ultraviolet (UV) λ = 266 nm and near infrared (NIR) λ = 1064 nm nanosecond laser induced plasmas in ambient air. For UV pulses at the conditions studied, energy absorption by the plasmas increases relatively gradually with laser pulse energy starting at delivered energy of E~8 mJ. Corresponding measurements of plasma luminosity show that the absorption of UV radiation does not necessarily result in visible plasma emission. For the NIR induced plasmas, the energy absorption profile is far more abrupt and begins at ~55 mJ. In contrast with UV, the absorption of NIR radiation is always accompanied by intense optical emission. The temperatures of both types of plasma have been measured with Rayleigh scattering thermometry (at times after the Thomson signal sufficiently diminishes). The UV plasmas can attain a wider range of temperatures, including lower temperatures, depending on the pulse energy (e.g., T ~400–2000K for E ~7–35 mJ at Δt = 10 μs after the pulse) while the NIR plasmas show only hotter temperatures (e.g., T ~ 12 000 K for E = 75 mJ at Δt = 10 μs after the pulse) as is consistent with the literature. Differences in the fluid dynamics for UV versus NIR pulses are shown with Schlieren imaging. In this work, the contrast in the UV and NIR plasma threshold behavior is attributed to differing roles of avalanche ionization and multiphoton ionization as is also illustrated by a simple numerical model. 3) The present contribution examines the impact of plasma dynamics and plasma-driven fluid dynamics on the flame growth of laser ignited mixtures and shows that a new dual-pulse scheme can be used to control the kernel formation process in ways that extend the lean ignition limit. We perform a comparative study between (conventional) single-pulse laser ignition (λ = 1064 nm) and a novel dual-pulse method based on combining an ultraviolet (UV) pre-ionization pulse (λ = 266 nm) with an overlapped near-infrared (NIR) energy addition pulse (λ = 1064 nm). We employ OH* chemiluminescence to visualize the evolution of the early flame kernel. For single-pulse laser ignition at lean conditions, the flame kernel separates through third lobe detachment, corresponding to high strain rates that extinguish the flame. In this work, we investigate the capabilities of the dual-pulse to control the plasma-driven fluid dynamics by adjusting the axial offset of the two focal points. In particular, we find there exists a beam waist offset whereby the resulting vorticity suppresses formation of the third lobe, consequently reducing flame stretch. With this approach, we demonstrate that the dual-pulse method enables reduced flame speeds (at early times), an extended lean limit, increased combustion efficiency, and decreased laser energy requirements.

Mechanical engineering↗

Mini-cavity-dumped laser

Lasers for use in high precision satellite ranging systems consist typically of an oscillator followed by several amplifier stages. While the shortest optical pulses are achieved by using a mode locked oscillator, such an oscillator is incompatible with the compact design needed in future, highly mobile systems. The laser oscillator achieves pulse lengths approaching those obtainable by mode locking, but in a much more compact and stable design. The oscillator uses two LiNbO3 Pockels cells inside the resonator. One Q-switches the oscillator, and the other is used in a pulse slicing scheme to cavity dump a portion of the circulating optical energy. The length of the optical output pulse measured at 425 + or - 50 picoseconds.

Reed, E.↗

Particle simulation of plasmas on the massively parallel processor

Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.

Gledhill, I. M. A.↗

The effects of spaceflight on open-loop and closed-loop postural control mechanisms: human neurovestibular studies on SLS-2

Stabilogram-diffusion analysis was used to examine how prolonged periods in microgravity affect the open-loop and closed-loop postural control mechanisms. It was hypothesized that following spaceflight: (1) the effective stochastic activity of the open-loop postural control schemes in astronauts is increased; (2) the effective stochastic activity and uncorrelated behavior, respectively, of the closed-loop postural control mechanisms in astronauts are increased; and (3) astronauts utilized open-loop postural controls schemes for shorter time intervals and smaller displacements. Four crew members and two alternates from the 14-day Spacelab Life Sciences 2 Mission were included in the study. Each subject was tested under eyes-open, quiet-standing conditions on multiple preflight and postflight days. The subjects' center-of-pressure trajectories were measured with a force platform and analyzed according to stabilogram-diffusion analysis. It was found that the effective stochastic activity of the open-loop postural control schemes in three of the four crew members was increased following spaceflight. This result is interpreted as an indication that there may be in-flight adaptations to higher-level descending postural control pathways, e.g., a postflight increase in the tonic activation of postural muscles. This change may also be the consequence of a compensatory (e.g., "stiffening") postural control strategy that is adopted by astronauts to account for general feeling of postflight unsteadiness. The crew members, as a group, did not exhibit any consistent preflight/postflight differences in the steady-state behavior of their closed-loop postural control mechanisms or in the functional interaction of their open-loop and closed-loop postural control mechanisms. These results are interpreted as indications that although there may be in-flight adaptations to the vestibular system and/or proprioceptive system, input from the visual system can compensate for such changes during undisturbed stance.

NASA Program Flight↗

Calibration and Data Analysis of the MC-130 Air Balance

Design, calibration, calibration analysis, and intended use of the MC-130 air balance are discussed. The MC-130 balance is an 8.0 inch diameter force balance that has two separate internal air flow systems and one external bellows system. The manual calibration of the balance consisted of a total of 1854 data points with both unpressurized and pressurized air flowing through the balance. A subset of 1160 data points was chosen for the calibration data analysis. The regression analysis of the subset was performed using two fundamentally different analysis approaches. First, the data analysis was performed using a recently developed extension of the Iterative Method. This approach fits gage outputs as a function of both applied balance loads and bellows pressures while still allowing the application of the iteration scheme that is used with the Iterative Method. Then, for comparison, the axial force was also analyzed using the Non-Iterative Method. This alternate approach directly fits loads as a function of measured gage outputs and bellows pressures and does not require a load iteration. The regression models used by both the extended Iterative and Non-Iterative Method were constructed such that they met a set of widely accepted statistical quality requirements. These requirements lead to reliable regression models and prevent overfitting of data because they ensure that no hidden near-linear dependencies between regression model terms exist and that only statistically significant terms are included. Finally, a comparison of the axial force residuals was performed. Overall, axial force estimates obtained from both methods show excellent agreement as the differences of the standard deviation of the axial force residuals are on the order of 0.001 % of the axial force capacity.

Booth, Dennis↗

Hybrid nanocapsules for in situ TEM imaging of gas evolution reactions in confined liquids

Liquid cell transmission electron microscopy (TEM) enables the direct observation of dynamic physical and chemical processes in liquids at the nanoscale. Quantitative investigations into reactions with fast kinetics and/or multiple reagents will benefit from further advances in liquid cell design that facilitate rapid in situ mixing and precise control over reagent volumes and concentrations. This work reports the development of inorganic–organic nanocapsules for high-resolution TEM imaging of nanoscale reactions in liquids with well-defined zeptoliter volumes. These hybrid nanocapsules, with 48 nm average diameter, consist of a thin layer of gold coating a lipid vesicle. As a model reaction, the nucleation, growth, and diffusion of nanobubbles generated by the radiolysis of water is investigated inside the nanocapsules. When the nanobubbles are sufficiently small (10–25 nm diameter), they are mobile in the nanocapsules, but their movement deviates from Brownian motion, which may result from geometric confinement by the nanocapsules. Gases and fluids can be transported between two nanocapsules when they fuse, demonstrating in situ mixing without using complex microfluidic schemes. The ability to synthesize nanocapsules with controlled sizes and to monitor dynamics simultaneously inside multiple nanocapsules provides opportunities to investigate nanoscale processes such as single nanoparticle synthesis in confined volumes and biological processes such as biomineralization and membrane dynamics.

36 MATERIALS SCIENCE↗

Seeded stimulated X-ray emission at 5.9 keV

X-ray free-electron lasers (XFELs) provide intense pulses that can generate stimulated X-ray emission, a phenomenon that has been observed and studied in materials ranging from neon to copper. Two schemes have been employed: amplified spontaneous emission (ASE) and seeded stimulated emission (SSE), where a second color XFEL pulse provides the seed. Both phenomena are currently explored for coherent X-ray laser sources and spectroscopy. Here, we report measurements of ASE and SSE of the 5.9 keV Mn K α 1 fluorescence line from a 3.9 molar NaMnO 4 solution, pumped with 7 femtosecond FWHM XFEL pulses at 6.6 keV. We observed ASE at a pump pulse intensity of 1.7×10 19 W/cm 2 , consistent with earlier findings. We observed SSE at dramatically reduced pump pulse intensities down to 1.1×10 17 W/cm 2 . These intensities are well within the range of many existing XFEL instruments, which supports the experimental feasibility of SSE as a tool to generate coherent X-ray pulses, spectroscopic studies of transition metal complexes, and other applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulation of turbulent mixing layers

A hybrid numerical scheme, combined with a scalar convection scheme, is used to compute turbulent mixing layers. The hybrid scheme is used to solve the vorticity equation by dividing the domain into two regions, applying a Lagrangian method in one region and an Eulerian method in the other. The scalar convection scheme is used to track an interface which separates a dyed fluid from a dye-free fluid, thereby distinguishing the fluid that starts on the high speed side from the fluid that starts on the low speed side. A two-dimensional mixing layer has been simulated using this new method; the roll up of a vortex sheet into structures that interact and pair is observed. The growth of the momentum thickness and the development of the mean turbulence statistics are found to be consistent with experimental observations. Forcing of the layer at moderate and low frequencies shows a dramatic effect on the growth of the momentum thickness and the entrainment ratio. Low frequency forcing enhances the growth of the layer; moderate frequency forcing inhibits the growth.

Mansour, N. N.↗

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗

Revolutionary Vertical Lift Technology Lift + Cruise Concept Vehicle Powertrain Dynamic Model

The Lift + Cruise (L+C) is one of several Revolutionary Vertical Lift Technology (RVLT) concept aircraft identified by NASA to investigate Urban Air Mobility (UAM) requirements. This six-passenger vehicle has a conventional tube-and-wing configuration, but with eight lift rotors (four on each wing) and a pusher rotor on the tail. This arrangement allows it to take off and land vertically as well as fly like an airplane. This report presents a novel dynamic model of the L+C powertrain built using the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS) and the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), both of which are NASA-developed toolboxes, and together are capable of modeling hybrid-electric powertrain architectures. This research effort models the complex behavior of the powertrain at the mechanical timescale (15 to 20 ms). The model consists of a turboshaft engine, generator, direct current bus, energy storage management system, gearboxes, and electric propulsor systems coupled together to create a dynamic model that mimics the operation of the multidomain powertrain. Two control algorithms, designed to equalize the power consumed to the power produced, are compared to assess their practicality: (1) independent proportional-integral (PI) control and (2) linear quadratic integral regulator (LQIR) control. Simulation results highlight the key features present in the model and demonstrate its operation under the two control schemes.

revolutionary vertical lift technology↗

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗

Substrate Effects on Electronic Properties of Atomic Chains

Atomic chains, precise structures of adatoms created on an atomically regulated surface, are candidates for constituent elements in future electronics. It was predicted that Si chains were metallic and Mg chains were semiconducting, and a doping method was also discussed. The substrate was assumed to work as a non-interacting template holding the adatoms. However, this scheme requires a low-temperature environment so that the adatoms will not displace from their ideal positions due to unwanted thermal agitation. For better structural stability, we may seek a scheme to allow the adatoms to form chemical bonding with the substrate atoms and secure their positions. The chemical bonding has two major effects on the chain electronic properties. First, only the remaining s and p orbitals in an adatom not used for the chemical bonding can decide chain band structures, rather than the full set of orbitals previously assumed. Second, because of a possible HOMO energy difference between adatom and substrate atom, semiconducting chains are unintentionally doped. These effects are studied with a self-consistent tight-binding method with universal parameters. With one adatom per unit cell, adatom chains are semiconducting: (1) if adatoms are of group III and form one chemical bond per adatom, or (2) if adatoms are of group IV and form two chemical bonds. 'Me previous result [I I of realizing semiconducting chains by group II adatoms without chemical bonds is consistent with this picture. When the HOMO energy in an adatom is shallower than that in a substrate atom, the entire chain is positively polarized, unintentionally achieving p-type doping. When deeper, the chain is negatively charged, achieving n-type doping.

Yamada, Toshishige↗

The Role and Quality of Software Safety in the NASA Constellation Program

In this study, we examine software safety risk in the early design phase of the NASA Constellation spaceflight program. Obtaining an accurate, program-wide picture of software safety risk is difficult across multiple, independently-developing systems. We leverage one source of safety information, hazard analysis, to provide NASA quality assurance managers with information regarding the ongoing state of software safety across the program. The goal of this research is two-fold: 1) to quantify the relative importance of software with respect to system safety; and 2) to quantify the level of risk presented by software in the hazard analysis. We examined 154 hazard reports created during the preliminary design phase of three major flight hardware systems within the Constellation program. To quantify the importance of software, we collected metrics based on the number of software-related causes and controls of hazardous conditions. To quantify the level of risk presented by software, we created a metric scheme to measure the specificity of these software causes. We found that from 49-70% of hazardous conditions in the three systems could be caused by software or software was involved in the prevention of the hazardous condition. We also found that 12-17% of the 2013 hazard causes involved software, and that 23-29% of all causes had a software control. Furthermore, 10-12% of all controls were software-based. There is potential for inaccuracy in these counts, however, as software causes are not consistently scoped, and the presence of software in a cause or control is not always clear. The application of our software specificity metrics also identified risks in the hazard reporting process. In particular, we found a number of traceability risks in the hazard reports may impede verification of software and system safety.

Layman, Lucas↗

LightForce Photon-pressure Collision Avoidance: Efficiency Analysis in the Current Debris Environment and Long-Term Simulation Perspective

This work provides an efficiency analysis of the LightForce space debris collision avoidance scheme in the current debris environment and describes a simulation approach to assess its impact on the long-term evolution of the space debris environment. LightForce aims to provide just-in-time collision avoidance by utilizing photon pressure from ground-based industrial lasers. These ground stations impart minimal accelerations to increase the miss distance for a predicted conjunction between two objects. In the first part of this paper we will present research that investigates the short-term effect of a few systems consisting of 10kW class lasers directed by 1.5 m diameter telescopes using adaptive optics. The results found such a network of ground stations to mitigate more than 85 percent of conjunctions and could lower the expected number of collisions in Low Earth Orbit (LEO) by an order of magnitude. While these are impressive numbers that indicate LightForce's utility in the short-term, the remaining 15 percent of possible collisions contain (among others) conjunctions between two massive objects that would add large amount of debris if they collide. Still, conjunctions between massive objects and smaller objects can be mitigated. Hence we choose to expand the capabilities of the simulation software to investigate the overall effect of a network of LightForce stations on the long-term debris evolution. In the second part of this paper, we will present the planed simulation approach for that effort.

conjunction analysis↗

Lightforce Photon-Pressure Collision Avoidance: Efficiency Analysis in the Current Debris Environment and Long-Term Simulation Perspective

This work provides an efficiency analysis of the LightForce space debris collision avoidance scheme in the current debris environment and describes a simulation approach to assess its impact on the long-term evolution of the space debris environment. LightForce aims to provide just-in-time collision avoidance by utilizing photon pressure from ground-based industrial lasers. These ground stations impart minimal accelerations to increase the miss distance for a predicted conjunction between two objects. In the first part of this paper we will present research that investigates the short-term effect of a few systems consisting of 20-kilowatt-class lasers directed by 1.5-meter-diameter telescopes using adaptive optics. The results found such a network of ground stations to mitigate more than 85 percent of conjunctions and could lower the expected number of collisions in Low Earth Orbit (LEO) by an order of magnitude. While these are impressive numbers that indicate LightForce's utility in the short-term, the remaining 15 percent of possible collisions contain (among others) conjunctions between two massive objects that would add large amount of debris if they collide. Still, conjunctions between massive objects and smaller objects can be mitigated. Hence, we choose to expand the capabilities of the simulation software to investigate the overall effect of a network of LightForce stations on the long-term debris evolution. In the second part of this paper, we will present the planned simulation approach for that effort. For the efficiency analysis of collision avoidance in the current debris environment, we utilize a simulation approach that uses the entire Two Line Element (TLE) catalog in LEO for a given day as initial input. These objects are propagated for one year and an all-on-all conjunction analysis is performed. For conjunctions that fall below a range threshold, we calculate the probability of collision and record those values. To assess efficiency, we compare a baseline (without collision avoidance) conjunction analysis with an analysis where LightForce is active. Using that approach, we take into account that collision avoidance maneuvers could have effects on third objects. Performing all-on-all conjunction analyses for extended period of time requires significant computer resources; hence we implemented this simulation utilizing a highly parallel approach on the NASA Pleiades supercomputer.

space debris mitigation↗