Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “small system collisions”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Non-conformal attractor in boost-invariant plasmas

We study the dissipative evolution of (0+1)-dimensionally expanding media with Bjorken symmetry using the Boltzmann equation for massive particles in relaxation-time approximation. Breaking conformal symmetry by a mass induces a non-zero bulk viscous pressure in the medium. It is shown that even a small mass (in units of the local temperature) drastically modifies the well-known attractor for the shear Reynolds number previously observed in massless systems. For generic nonzero particle mass, neither the shear nor the bulk viscous pressure relax quickly to a non-equilibrium attractor; they approach the hydrodynamic limit only late, at small values of the inverse Reynolds numbers. Only the longitudinal pressure, which is a combination of thermal, shear and bulk viscous pressures, continues to show early approach to a far-off-equilibrium attractor, driven by the rapid longitudinal expansion at early times. Second-order dissipative hydrodynamics based on a gradient expansion around locally isotropic thermal equilibrium fails to reproduce this attractor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

The fall, recovery, classification, and initial characterization of the Hamburg, Michigan H4 chondrite

Abstract The Hamburg meteorite fell on January 16, 2018, near Hamburg, Michigan, after a fireball event widely observed in the U.S. Midwest and in Ontario, Canada. Several fragments fell onto frozen surfaces of lakes and, thanks to weather radar data, were recovered days after the fall. The studied rock fragments show no or little signs of terrestrial weathering. Here, we present the initial results from an international consortium study to describe the fall, characterize the meteorite, and probe the collision history of Hamburg. About 1 kg of recovered meteorites was initially reported. Petrology, mineral chemistry, trace element and organic chemistry, and O and Cr isotopic compositions are characteristic of H4 chondrites. Cosmic ray exposure ages based on cosmogenic 3 He, 21 Ne, and 38 Ar are ~12 Ma, and roughly agree with each other. Noble gas data as well as the cosmogenic 10 Be concentration point to a small 40–60 cm diameter meteoroid. An 40 Ar‐ 39 Ar age of 4532 ± 24 Ma indicates no major impact event occurring later in its evolutionary history, consistent with data of other H4 chondrites. Microanalyses of phosphates with LA‐ICPMS give an average Pb‐Pb age of 4549 ± 36 Ma. This is in good agreement with the average SIMS Pb‐Pb phosphate age of 4535.3 ± 9.5 Ma and U‐Pb Concordia age of 4535 ± 10 Ma. The weighted average age of 4541.6 ± 9.5 Ma reflects the metamorphic phosphate crystallization age after parent body formation in the early solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Blue Rock CropSpanPV - Low Cost Racking for Agricultural Solar Photovoltaics

In this work, a total of four solar panel layouts, applied to Agrivoltaics, are investigated. This study presents two hypotheses: 1. Prefabrication - Lower total cost can be achieved by constructing a pre-assembled, pre-wired solar array in a factory, with automation, then rapidly deploying it in the field. 2. Tension Structures - Using tension structures in a way to suspend an overhead solar array also reduces costs when compared to conventional support methods. Designs are presented for frames to hold the panels. Configurations are advanced to connect the frames into arrays. Supports and foundations to hold the arrays in the desired layout are calculated, selected and presented. Cost models are prepared for the four layouts, which include definition of the business enterprise, the manufacturing operation, and the required facility for production at scale. Cost models include, detail material and task/labor take-offs for both manufacturing and installation. Comparisons with conventional system adapted to Agrivoltaics reveal that the two hypotheses do not hold true. A prefabricated array requires additional material that is not needed in the base system and the cost of this additional material is not outweighed by the cost savings of reduced field installation. The cost of wire rope for a tension structure does result in an economy of material, however, the cost of end connections and tensioning devices adds significantly to the cost of the overall tension structure. Foundation loads with tension structures, especially the large wind uplift seen with a fixed tilt solar array, impose significant constraints and high costs. The following conclusions are recommended for further consideration: 1. Best Case Agrivoltaic Configuration - The best value/lowest cost configuration for agrivoltaics is the installation of solar panels mounted on a single axis tracker in either a 2-in-portrait or 2-in landscape arrangement. This is useful for grazing lands and for staple crops which require high light levels. These applications represent the vast majority of the potential agrivoltaic market. Important elements in this configuration are the specific geometric layout coordinated with the various field operations and a control system linked between the farm equipment and single axis tracker. The control system is to monitor location and orientation, then actively tilt the solar panels or brake the field equipment to provide clearance and avoid collisions. 2. Long Span Applications - The use of tension structures for an overhead fixed tilt array is beneficial where long spans are a necessity, especially where the base support in the project location has consolidated rock near the surface or some other solid structure. One application is spanning wide irrigation canals. Another may be as an installation at the top deck of parking structures. 3. Direct Ground Mount Configuration – Having skids at the base of a pre-fabricated array that extend when the array is deployed, could be a successful design. This could be a useful strategy for dense, flat, ground-mount configurations like that employed by a couple of successful commercial operations (5G Maverick and Erthos). Arrays would be orientated north and south, deployed very low on prepared grade. Arrays would be aligned closely next to each other and anchored with small ground screws or simply with ballast. This same product arrangement could be very useful for rapid deployment and set-up of solar arrays for temporary deployments, disaster response and the military operations

14 SOLAR ENERGY↗
NASA NTRS1990

Mathematical modeling of the flow field and particle motion in a rotating bioreactor at unit gravity and microgravity

The biotechnology group at NASA Johnson Space Center is developing systems for culturing mammalian cells that stimulate some aspect of microgravity and provide a low shear environment for microgravity-based studies on suspension and anchorage dependent cells. The design of these vessels for culturing cells is based on the need to suspend cells and aggregates of cells and microcarrier beads continually in the culturing medium. The design must also provide sufficient circulation for adequate mass transfer of nutrients to the cells and minimize the total force on the cells. Forces, resulting from sources such as hydrodynamic fluid shear and collisions of cells and walls of the vessels, may damage delicate cells and degrade the formation of three dimensional structures. This study examines one particular design in both unit gravity and microgravity based on two concentric cylinders rotating in the same direction at different speeds to create a Couette flow between them. A numerical simulation for the flow field and the trajectories of particles in the vessel. The flow field for the circulation of the culturing medium is modeled by the Navier-Stokes equations. The forces on a particle are assumed to be drag from the fluid's circulation, buoyancy from the gravitational force and centrifugal force from the rotation of the vessel. The problem requires first solving the system of partial differential equations for the fluid flow by a finite difference method and then solving the system of ordinary differential equations for the trajectories by Gear's stiff method. Results of the study indicate that the trajectories in unit gravity and microgravity are very similar except for small spatial deviations on the fast time scale in unit gravity. The total force per unit cross sectional area on a particle in microgravity, however, is significantly smaller than the corresponding value in unit gravity, which is also smaller than anticipated. Hence, this study indicates that this design for a bioreactor with optimal rates of rotation can provide a good environment for culturing cells in microgravity with adequate circulation and minimal force on the cells.

Boyd, Ernest J.↗
DOE OSTI2020

Laser-induced fluorescence measurements of ion fluctuations in electron and ion presheaths

Electrostatic flow-driven instabilities have been predicted in the presheath region of low-temperature plasmas for both ion and electron rich sheaths. The ion-acoustic fluctuations that arise from these instabilities can influence transport through an enhanced collision rate. In this work, ion fluctuation spectra were measured using laser-induced fluorescence (LIF) in the presheath region for both sheath polarities. Furthermore, the non-invasive nature of LIF compared to that of probes is particularly valuable for these measurements because the sheath and presheath around a probe may generate their own flow-driven instabilities, which obscure the measurement. Measurements were made using a recently implemented field programmable gate array-based system able to measure ion fluctuation spectra up to 1 MHz using a two-point correlation function technique. Narrow bandwidth ion fluctuations were observed near 0:46 f pi (460 kHz) for both sheath polarities, where f pi is the ion plasma frequency. The observed fluctuations were significantly stronger in the electron sheath case and were observed over a larger volume of plasma. Fluctuations were measured in the ion sheath case at locations far enough from the sheath that linear theory predicts stability, suggesting that ion-acoustic instabilities generated in the relatively small unstable region near the boundary reflect from the sheath and perturb a much larger volume of the plasma. This is expected to affect probes with both ion and electron rich sheaths, meaning any dc biased probe may effectively act as an ion-acoustic wave antenna. These measurements are consistent with the recent theory and particle-in-cell results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2014

Path Toward a Unifid Geometry for Radiation Transport

The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex computer-aided design (CAD) models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN [high charge and energy transport code developed by NASA Langley Research Center (LaRC)], are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their widespread use for analysis of complex CAD models, leading to the creation and maintenance of toolkit-specific simplistic geometry models. The work presented here builds on the Direct Accelerated Geometry Monte Carlo (DAGMC) toolkit developed for use with the Monte Carlo N-Particle (MCNP) transport code. The workflow for achieving radiation transport on CAD models using MCNP and FLUKA has been demonstrated and the results of analyses on realistic spacecraft/habitats will be presented. Future work is planned that will further automate this process and enable the use of multiple radiation transport codes on identical geometry models imported from CAD. This effort will enhance the modeling tools used by NASA to accurately evaluate the astronaut space radiation risk and accurately determine the protection provided by as-designed exploration mission vehicles and habitats

Lee, Kerry↗
NASA NTRS2014

Evaluation of Cloud-Resolving and Limited Area Model Intercomparison Simulations Using TWP-ICE Observations: Precipitation Microphysics - Part 2

Ten 3-D cloud-resolving model (CRM) simulations and four 3-D limited area model (LAM) simulations of an intense mesoscale convective system observed on 23-24 January 2006 during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) are compared with each other and with observations and retrievals from a scanning polarimetric radar, colocated UHF and VHF vertical profilers, and a Joss-Waldvogel disdrometer in an attempt to explain a low bias in simulated stratiform rainfall. Despite different forcing methodologies, similar precipitation microphysics errors appear in CRMs and LAMs with differences that depend on the details of the bulk microphysics scheme used. One-moment schemes produce too many small raindrops, which biases Doppler velocities low, but produces rainwater contents (RWCs) that are similar to observed. Two-moment rain schemes with a gamma shape parameter (mu) of 0 produce excessive size sorting, which leads to larger Doppler velocities than those produced in one-moment schemes but lower RWCs. Two-moment schemes also produce a convective median volume diameter distribution that is too broad relative to observations and, thus, may have issues balancing raindrop formation, collision-coalescence, and raindrop breakup. Assuming a mu of 2.5 rather than 0 for the raindrop size distribution improves one-moment scheme biases, and allowing mu to have values greater than 0 may improve excessive size sorting in two-moment schemes. Underpredicted stratiform rain rates are associated with underpredicted ice water contents at the melting level rather than excessive rain evaporation, in turn likely associated with convective detrainment that is too high in the troposphere and mesoscale circulations that are too weak. A limited domain size also prevents a large, well-developed stratiform region like the one observed from developing in CRMs, although LAMs also fail to produce such a region.

three dimensional models↗
DOE OSTI2022

Light-front wavefunctions of mesons by design

Abstract We develop a mechanism to build the light-front wavefunctions (LFWFs) of meson bound states on a small-sized basis function representation. Unlike in a standard Hamiltonian formalism, the Hamiltonian in this method is implicit, and the information of the system is carried directly by the functional form and adjustable parameters of the LFWFs. In this work, we model the LFWFs for four charmonium states, $$\eta _c$$ η c , $$J/\psi $$ J / ψ , $$\psi '$$ ψ ′ , and $$\psi (3770)$$ ψ ( 3770 ) as superpositions of orthonormal basis functions. We choose the basis functions as eigenfunctions of an effective Hamiltonian, which has a longitudinal confining potential in addition to the transverse confining potential from light-front holographic QCD. We determine the basis function parameters and superposition coefficients by employing both guidance from the nonrelativistic description of the meson states and the experimental measurements of the meson decay widths. With the obtained wavefunctions, we study the features of those meson states, including charge radii and parton distribution functions. We use the $$J/\psi $$ J / ψ LFWF to calculate the meson production in diffractive deep inelastic scattering and ultra-peripheral heavy-ion collisions, and the $$\eta _c$$ η c LFWF to calculate its diphoton transition form factor. Both results show good agreement with experiments. The obtained LFWFs have simple-functional forms and can be readily used to predict additional experimental observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2024

Survey of Methods for Assessing Safety Compliance of sUAS BVLOS Operations

To ensure the safety of small Unmanned Aircraft System (sUAS) Beyond Line of Sight (BVLOS) operations in the National Airspace System (NAS), the FAA offers a guideline on how sUAS operators can demonstrate compliance with FAA rules, including regulations, advisory circulars, policy statements, and acceptable means of compliance (AMOC), using a formal and top-down approach. According to FAA Advisory Circular 23.2010-1, the AMOC refers to “one method, but not the only method, to show compliance with a regulatory requirement.” It is common for regulations to include a performance and safety standard rather than a detailed design or operation requirement, which allows for flexibility in fulfilling the regulatory requirements while still achieving a predetermined level of safety. This technical report explores existing frameworks for sUAS BVLOS safety analysis. It begins by discussing how sUAS BVLOS operations can be strategically deconflicted and then discusses their hazards. The next step is to establish safety assessment methods commonly used by the aviation industry and the FAA, illustrate how these methods can be applied to several safety-critical air traffic systems, and list collision models that the FAA and industry use. This investigation documents the processes for assessing safety risks in sUAS BVLOS operations and will allow sUAS BVLOS operations to be assessed for safety compliance more efficiently.

sUAS BVLOS Operations↗
NASA NTRS1991

Cosmological gamma-ray bursts

The distribution in angle and flux of gamma-ray bursts indicates that the majority of gamma-ray bursters are at cosmological distances, i.e., at z of about 1. The rate is then about 10 exp -8/yr in a galaxy like the Milky Way, i.e., orders of magnitude lower than the estimated rate for collisions between neutron stars in close binary systems. The energy per burst is about 10 exp 51 ergs, assuming isotropic emission. The events appear to be less energetic and more frequent if their emission is strongly beamed. Some tests for the distance scale are discussed: a correlation between the burst's strength and its spectrum; the absorption by the Galactic gas below about 2 keV; the X-ray tails caused by forward scattering by the Galactic dust; about 1 month recurrence of some bursts caused by gravitational lensing by foreground galaxies; and a search for gamma-ray bursts in M31. The bursts appear to be a manifestation of something exotic, but conventional compact objects can provide an explanation. The best possibility is offered by a decay of a bindary composed of a spinning-stellar-mass black-hole primary and a neutron or a strange-quark star secondary. In the final phase the secondary is tidally disrupted, forms an accretion disk, and up to 10 exp 54 ergs are released. A very small fraction of this energy powers the gamma-ray burst.

Paczynski, Bohdan↗
DOE OSTI2022

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2011

Small Orbital Debris Mitigation Mission Architecture

Small orbital debris in LEO (1-10 cm in size) presents a clear and present danger to operational LEO spacecraft. This debris field has dramatically increased (nearly doubled) in recent years following the Chinese ASAT Test in 2007 and the Iridium/Cosmos collision in 2009. Estimates of the number of small debris have grown to 500,000 objects after these two events; previously the population was 300,000 objects. These small, untracked debris objects (appproximately 500,000) outnumber the larger and tracked objects (appproximately 20,000) by a factor 25 to 1. Therefore, the risk of the small untracked debris objects to operational spacecraft is much greater than the risk posed by the larger and tracked LEO debris objects. A recent study by The Aerospace Corporation found that the debris environment will increase the costs of maintaining a constellation of government satellites by 5%, a constellation of large commercial satellites by 11%, and a constellation of factory built satellites by 26% from $7.6 billion to $9.57 billion. Based upon these facts, the NASA Marshall Space Flight Center (MSFC) Advanced Concepts Office (ACO) performed an architecture study on Small Orbital Debris Active Removal (SODAR) using a space-based nonweapons- class laser satellite for LEO debris removal. The goal of the SODAR study was to determine the ability of a space-based laser system to remove the most pieces of debris (1 cm to 10 cm, locations unknown), in the shortest amount of time, with the fewest number of spacecraft. The ESA developed MASTER2005 orbital debris model was used to probabilistically classify the future debris environment including impact velocity, magnitude, and directionality. The study ground rules and assumptions placed the spacecraft into a high inclination Low Earth Orbit at 800 km as an initial reference point. The architecture study results found that a spacecraft with an integrated forward-firing laser is capable of reducing the small orbital debris flux within a 60 to 100 km orbital shell by a significant amount within the one spacecraft's operational lifetime. The technology developments required for such an architecture to be successfully employed are: 1) A pico-pulsed, space qualified laser and 2) The ability to detect and track a small LEO small debris object from a space based platform. Therefore, a conceptual design for a demonstration satellite to showcase the ability to detect and track small orbital debris was completed

Wiegmann, Bruce M.↗
DOE OSTI2024

Adiabatic and post-adiabatic hyperspherical treatment of the huge ungerade proton-hydrogen scattering length

While the hydrogen molecular ion is the simplest molecule in nature and very well studied in all of its properties, it remains an interesting system to use for explorations of fundamental questions. One such question treated in this study relates to finding an optimal adiabatic representation of the physics, i.e., the best adiabatic description that minimizes the role of nonadiabatic effects. As a test case explored here in detail, we consider the ungerade symmetry of $H$$^{+}_{2}$, which is known to have a huge scattering length of order 750 Bohr radii, and an incredibly weakly bound excited state. We show that a hyperspherical adiabatic description does an excellent job of capturing the main physics. Furthermore, our calculation yields a competitive scattering length and shows that nonadiabatic corrections are small and can even be adequately captured using the post-adiabatic theory of Klar and Fano.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1972

Model for radon diffusion through the lunar regolith.

Description of a model for radon diffusion through the lunar regolith in which the atom migrates by random walk. The regolith is represented by a system of randomly oriented baffles in which the mean distance which the atom travels between two collisions takes on the role of a mean free path. The effective mean time between two collisions depends on two entities: the actual mean time-of-flight and the mean sticking time on grain surfaces for one collision. The latter depends strongly on the temperature and the heat of adsorption of radon on regolith materials. Both the mean free path as well as the heat of adsorption are either poorly known or unknown for the lunar regolith; hence these quantities are treated as free parameters. Because of the greatly different mean lifetimes against radioactive decay of Rn219, Rn220, and Rn222, the regolith acts as a powerful 'filter' for these species. Rn222 escape is significant (32%) even for a mean free path of 1 micron, a heat of adsorption of 7.0 kcal/mole and a regolith depth of 4 m. Calculations of radon escape from a 4 m thick regolith, using mean free paths of 1, 10, and 80 microns and heats of adsorption of 4.0, 5.2, and 7.0 kcal/mole show that the Rn222/Rn220 escape ratio can be as small as 7.7 and as large as, or larger than 47. The small value of 7.7 is of particular interest, because it is nearly equal to the escape ratio inferred by Turkevich et al. (1970) from their Surveyor 5 results.

Friesen, L. J.↗
NASA NTRS2019

The NASA Orbital Debris Engineering Model 3.1: Development, Verification, and Validation

The NASA Orbital Debris Program Office has developed the Orbital Debris Engineering Model (ORDEM) primarily as a tool for spacecraft designers and other users to understand the long-term risk of collisions with orbital debris. The newest version, ORDEM 3.1, incorporates the latest and highest fidelity datasets available to build and validate representative orbital debris populations encompassing low Earth orbit (LEO) to geosynchronous orbit (GEO) altitudes for the years 2016-2050. ORDEM 3.1 models fluxes for object sizes > 10 μm within or transiting LEO and > 10 cm in GEO. The deterministic portion of the populations in ORDEM 3.1 is based on the U.S. Space Surveillance Network (SSN) catalog, which provides coverage down to approximately 10 cm in LEO and 1 m in GEO. Observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller sizes and orbit regions that are not well-covered by the SSN catalog, yet may pose the greatest threat to operational spacecraft. Objects in LEO ranging from approximately 5 mm to 10 cm are modeled using observational data from ground-based radar, namely the Haystack Ultrawideband Satellite Imaging Radar (HUSIR – formerly known as Haystack). The LEO population smaller than approximately 3 mm in size is characterized based on a reanalysis of in situ data from impacts to the windows and radiators of the U.S. Space Transportation System orbiter vehicle, i.e., the Space Shuttle. Data from impacts on the Hubble Space Telescope are also used to validate the sub-millimeter model populations in LEO. Debris in GEO with sizes ranging from 10 cm to 1 m is modeled using optical measurement data from the Michigan Orbital DEbris Survey Telescope (MODEST). Specific, major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos 2251 debris clouds, and unique populations, such as sodium-potassium droplets, have been re-examined and are modeled and added to the ORDEM environment separately. The debris environment greater than 1 mm is forecast using NASA’s LEO-to- GEO ENvironment Debris model (LEGEND). Future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically, and the NASA Standard Satellite Breakup Model is used to generate fragments from these events. Fragments smaller than 10 cm are further differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to spacecraft. The future projection of the sub-millimeter environment is computed using a special small-particle degradation model where small particles are created from intact spacecraft and rocket bodies. This work discusses the development, features, and capabilities of the ORDEM 3.1 model; the ne new data analyses used to build the model populations; and sample verification and validation results.

Matney, M.↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
NASA NTRS1990

A multi-frequency study of the peculiar interacting system Arp 206

Arp 206 is a nearby, relatively large, and bright interacting system comprising unequal members: NGC 3432 and UGC 5983. A third anonymous galaxy, Arp 206c, is visible in the field. The CCD images show a well-developed bridge between NGC 3432 and UGC 5983. On the other hand, the complex H I tails are not visible in the optical. In the total H I map, the bridge is lost in a general envelope encompassing both galaxies. The bridge also appears to have some radio emission. On the Total H I map the system is rather edge-on, far more than it would appear in optical wavelengths. UGC 5983 falls exactly in line with NGC 3432. The velocity of the centers of mass of NGC 3432 and UGC 5983 are 530 km s(exp -1) and 630 km s(exp -1), respectively. In view of the considerable damage sustained by NGC 3432 and the apparent low mass of UGC 5983, it appears that the passage must have been at near parabolic speed, with a small pericentric distance and a very low inclination with rspect to the disk of NGC 3432. The apparent distribution of H I along the z axis of the galaxy could be accounted for by projection effects. The tidal appendage found at higher velocities, which rises at a P.A. approx. equal to 25 degrees west of the main body of the galaxy is probably the tail, the part of the tidal damage away from the perturbing companion. The bridge may be rising north-east from the galaxy and then continue under to the south of the galaxy. The relative sizes of the appendages would indicate that the pericenter was crossed recently. Any further inferences about the collision parameters will need to await the results of detailed computational modelling of the interaction. The authors also summarize the observational characteristics of NGC 3432, UGC 5983, and Arp 206c.

Noreau, Louis↗
NASA NTRS2001

Formation of Planetary Systems

Stars are observed to be forming within cold regions of the galaxy called molecular clouds. A clump of gas and dust within a molecular cloud can collapse into a rotationally- supported disk orbiting the pre s sure -supported star. Such a disk has the same initial elemental composition as the growing star. At sufficient distances from the central star, it is cool enough for approx. 1 - 2% of this material to be in solid form, either remnant interstellar grains or condensates formed within the disk. During the infall stage, the disk is very active and probably highly turbulent. When the infall slows substantially or stops, the disk becomes more quiescent. The dust grows from micron-sized dust to kilometer-sized planetesimals by physical collisions, possibly aided by collective gravity. The dynamics of larger solid bodies within protoplanetary disks are better characterized. The primary perturbations on the orbits of kilometer-sized and larger planetesimals in protoplanetary disks are mutual gravitational interactions and physical collisions. These interactions lead to accretion (and in some cases erosion and fragmentation) of planetesimals. The most massive planets have the largest gravitationally-enhanced collision cross-sections, and accrete almost everything with which they collide. The size distribution of solid bodies becomes quite skewed, with a few large bodies growing much faster than the rest of the swarm, until they have accumulated most of the small bodies within their gravitational reach. Slower growth continues (at least for solid, earth-like planets) as the eccentricities of planetary embryos are pumped up by long-range mutual gravitational perturbations. As planetary masses increase, they become more efficient at stirring random velocities of neighboring bodies. If sufficiently massive and dense planets exist far enough from the star, they can eject material into interstellar space. In most models, giant planets begin their growth like terrestrial planets, but if they become massive enough before the protoplanetary disk dissipates, then they are able to accumulate substantial amounts of gas. Giant planet growth times predicted by current models are similar to estimates of the lifetime of the gaseous protoplanetary disk, so we don't know whether or not giant planets form in most protoplanetary disks. However, planets that become massive while a substantial amount of gas remains in the disk may migrate into the star as a consequence of their gravitational interactions with the disk.

Lissauer, Jack J.↗