Search NASA⌕ Search

SEARCH · Search NASA

Results for “Collision rate coefficient”

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.

124 records · Page 7

Evolution of inclined galactic gas disks. I - A cloud-fluid approach

In this paper, numerical models are developed to follow the time-dependent behavior of inclined gas disks in nonspherical galaxies. Differential equations are developed which describe the time rate of change of the orbits of mass elements composing a thin, smoothly warping and twisting gas disk. These equations include multiple moments due to the nonsphericity of the galactic mass distribution, Coriolis effects arising from determination of orbits in a reference frame which tumbles with the surface figure of the galaxy, and viscous forces arising from cloud-cloud collisions. A cloud-fluid approach is used in determining the viscous forces. Physical arguments are used to determine the maximum permissible local value of the coefficient of kinematic viscosity. For small initial inclinations, analytic expressions are obtained for the inclination as a function of time and for the settling time, and the features of the solution are discussed.

Steiman-Cameron, Thomas Y.↗

Analysis of the upper atmosphere CO2(nu-2) vibrational temperatures retrieved from ATMOS/Spacelab 3 observations

An analysis of the upper atmospheric (80-116 km) CO2(nu-2) vibrational temperatures retrieved from atmospheric trace molecule spectroscopy (ATMOS) experiment Spacelab 3 spectra by using a nonlocal thermodynamic equilibrium (non-LTE) radiative transfer model is presented. Thermal collisions with atmospheric atomic oxygen keep this vibrational state very close to LTE up to around 100 km. Above this height, the different deviations from LTE retrieved from ATMOS/Spacelab 3 spectra for the Northern and the Southern Hemispheres are explained in terms of this collisional process and in terms of the different kinetic temperature profiles measured at those locations. From these simultaneous observations of the kinetic and CO2(nu-2) vibrational temperatures, a deactivation rate of CO2(nu-2) by O(3P) has been derived which leads to a rate coefficient value between 3 and 6 x 10 exp -12 cu cm/s and favors an independent or negative temperature dependence rate constant for the atmospheric temperature range. Cooling rates induced by the CO2 15-micron fundamental band in the upper mesosphere and lower thermosphere were derived from the simultaneous kinetic temperature, CO2(nu-2) vibrational temperature, and CO2 concentration, as measured by ATMOS/Spacelab 3, and found to be a factor of between 5 and 10 times larger than those generally accepted until very recently.

Lopez-Puertas, M.↗

Vibrational nonequilibrium effects on diatomic dissociation rates

The collision-induced dissociation rate of diatomic molecules from a ladder of rotational and anharmonic vibrational states is developed, and the correction for vibrational nonequilibrium is considered. The result is similar to an analytic correction derived by Hammerling et al. (1959) for harmonic oscillators. An empirical correction algorithm suggested by Park (1987, 1990) gives similar results when vibrational temperature is comparable to kinetic temperature but underestimates the dissociation rate when vibrational temperature is small compared with the kinetic temperature. This algorithm uses an effective temperature in the experimentally determined Arrhenius expression for the rate coefficient, which is a weighted average of the vibrational and kinetic temperature, whereas theory indicates that kinetic temperature should appear only in the exponential term of the Arrhenius expression. Nevertheless, an effective temperature can always be found that will numerically duplicate the proper rate coefficient at any given condition, but a constant weighting factor cannot be expected to provide this. However, the algorithm can he adjusted to give reasonable results over a range of conditions if the geometric weighting factor is taken to be a simple linear function of the ratio of vibrational to kinetic temperature in the gas.

Hansen, C. F.↗

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗

Estimated Probability of a Cervical Spine Injury During an ISS Mission

Introduction: The Integrated Medical Model (IMM) utilizes historical data, cohort data, and external simulations as input factors to provide estimates of crew health, resource utilization and mission outcomes. The Cervical Spine Injury Module (CSIM) is an external simulation designed to provide the IMM with parameter estimates for 1) a probability distribution function (PDF) of the incidence rate, 2) the mean incidence rate, and 3) the standard deviation associated with the mean resulting from injury/trauma of the neck. Methods: An injury mechanism based on an idealized low-velocity blunt impact to the superior posterior thorax of an ISS crewmember was used as the simulated mission environment. As a result of this impact, the cervical spine is inertially loaded from the mass of the head producing an extension-flexion motion deforming the soft tissues of the neck. A multibody biomechanical model was developed to estimate the kinematic and dynamic response of the head-neck system from a prescribed acceleration profile. Logistic regression was performed on a dataset containing AIS1 soft tissue neck injuries from rear-end automobile collisions with published Neck Injury Criterion values producing an injury transfer function (ITF). An injury event scenario (IES) was constructed such that crew 1 is moving through a primary or standard translation path transferring large volume equipment impacting stationary crew 2. The incidence rate for this IES was estimated from in-flight data and used to calculate the probability of occurrence. The uncertainty in the model input factors were estimated from representative datasets and expressed in terms of probability distributions. A Monte Carlo Method utilizing simple random sampling was employed to propagate both aleatory and epistemic uncertain factors. Scatterplots and partial correlation coefficients (PCC) were generated to determine input factor sensitivity. CSIM was developed in the SimMechanics/Simulink environment with a Monte Carlo wrapper (MATLAB) used to integrate the components of the module. Results: The probability of generating an AIS1 soft tissue neck injury from the extension/flexion motion induced by a low-velocity blunt impact to the superior posterior thorax was fitted with a lognormal PDF with mean 0.26409, standard deviation 0.11353, standard error of mean 0.00114, and 95% confidence interval [0.26186, 0.26631]. Combining the probability of an AIS1 injury with the probability of IES occurrence was fitted with a Johnson SI PDF with mean 0.02772, standard deviation 0.02012, standard error of mean 0.00020, and 95% confidence interval [0.02733, 0.02812]. The input factor sensitivity analysis in descending order was IES incidence rate, ITF regression coefficient 1, impactor initial velocity, ITF regression coefficient 2, and all others (equipment mass, crew 1 body mass, crew 2 body mass) insignificant. Verification and Validation (V&V): The IMM V&V, based upon NASA STD 7009, was implemented which included an assessment of the data sets used to build CSIM. The documentation maintained includes source code comments and a technical report. The software code and documentation is under Subversion configuration management. Kinematic validation was performed by comparing the biomechanical model output to established corridors.

Brooker, John E.↗

Charge Transfer Calculations and Database for Astrophysics

A variety of theoretical approaches, having different but overlapping energy ranges of applicability, are applied to investigate charge transfer processes for collisions of atomic ions with atoms and molecules. The methods include quantal molecular-orbital close-coupling, classical trajectory Monte Carlo, and continuum distorted wave methods. Recent collision systems studied include S(+4) + H, S(+4) + He, N(+7) + He, H2O, CO, and CO2, O(+q)(q = 1 - 8) + H, H2, and S(+q)(q = 1 - 16) + H2. The database effort is concentrating on astrophysically important reactions of atomic ions X(+q)(X=H-Zn, q=1-4, and selected higher charges) with H, He, various metal atoms, H2, and other selected molecular targets. Existing data, much of it produced by us, has been compiled and critically evaluated. Data for many reactions missing in the literature are estimated using the multichannel Landau-Zener approximation. Fits to cross sections and rate coefficients using standard functions are provided as well as tabulations of the raw data. The database is available on the World Wide Web at cfadc.phy.ornl.gov/astro/ps/data.

Wang, J. G.↗

Solar observations and atomic data for the 3s2 1s0-3s3p 3P1 transition in S v

High resolution solar observations of the S v intersystem line at 1199.18 A are available from Skylab. This line is potentially useful as a density diagnostic for high density plasmas expected in solar flares. S v lines are also prominent in solar spectra at wavelengths below 1000 A. Collision strengths and radiative decay rates are calculated for levels of the configurations 3 s2, 3s3p, 3p2, and 3s3d. Level populations for the five lowest energy levels have been calculated as a function of electron density. These calculations are carried out assuming the temperature at which S v is most abundant in solar plasmas. The calculated population of the 3s3p 3P1 level reaches a pseudo-Boltzmann equilibrium at a density which is about five times higher than is derived from solar spectra and previous density determinations. It is believed that the discrepancy is due to a significant contribution from resonances to the 3s2-3s3p 3P excitation rate coefficient. It is suggested that the contribution of the resonances to the S v excitation rate coefficient can be estimated from a comparison of the experimental results and the theoretical calculations.

Feldman, U.↗

Lunar and Planetary Science XXXV: Terrestrial Planets: Building Blocks and Differentiation

The session "Terrestrial Planets: Building Blocks and Differentiation: included the following topics:Magnesium Isotopes in the Earth, Moon, Mars, and Pallasite Parent Body: High-Precision Analysis of Olivine by Laser-Ablation Multi-Collector ICPMS; Meteoritic Constraints on Collision Rates in the Primordial Asteroid Belt and Its Origin; New Constraints on the Origin of the Highly Siderophile Elements in the Earth's Upper Mantle; Further Lu-Hf and Sm-Nd Isotopic Data on Planetary Materials and Consequences for Planetary Differentiation; A Deep Lunar Magma Ocean Based on Neodymium, Strontium and Hafnium Isotope Mass Balance Partial Resetting on Hf-W System by Giant Impacts; On the Problem of Metal-Silicate Equilibration During Planet Formation: Significance for Hf-W Chronometry ; Solid Metal-Liquid Metal Partitioning of Pt, Re, and Os: The Effect of Carbon; Siderophile Element Abundances in Fe-S-Ni-O Melts Segregated from Partially Molten Ordinary Chondrite Under Dynamic Conditions; Activity Coefficients of Silicon in Iron-Nickel Alloys: Experimental Determination and Relevance for Planetary Differentiation; Reinvestigation of the Ni and Co Metal-Silicate Partitioning; Metal/Silicate Paritioning of P, Ga, and W at High Pressures and Temperatures: Dependence on Silicate Melt Composition; and Closure of the Fe-S-Si Liquid Miscibility Gap at High Pressure and Its Implications for Planetary Core Formation.

Source record↗

Dissociation and Internal Excitation of Molecular Nitrogen Due to N + N2 Collisions Using Direct Molecular Simulation

In this work we present a molecular level study of N2+N collisions, focusing on excitation of internal energy modes and non-equilibrium dissociation. The computation technique used here is the direct molecular simulation (DMS) method and the molecular interactions have been modeled using an ab−initio potential energy surface (PES) developed at NASA's Ames Research Center. We carried out vibrational excitation calculations between 5000K and 30000K and found that the characteristic vibrational excitation time for the N + N2 process was an order of magnitude lower than that predicted by the Millikan and White correlation. It is observed that during vibrational excitation the high energy tail of the vibrational energy distribution gets over populated first and the lower energy levels get populated as the system evolves. It is found that the non-equilibrium dissociation rate coefficients for the N + N2 process are larger than those for the N2 + N2 process. This is attributed to the non-equilibrium vibrational energy distributions for the N + N2 process being less depleted than that for the N2 +N2 process. For an isothermal simulation we find that the probability of dissociation goes as 1/T(sub tr) for molecules with internal energy (epsilon(sub int)) less than approximately 9.9eV, while for molecules with epsilon (sub int) greater than 9.9eV the dissociation probability was weakly dependent on translational temperature of the system. We compared non-equilibrium dissociation rate coefficients and characteristic vibrational excitation times obtained by using the ab-initio PES developed at NASA's Ames Research Center to those obtained by using an ab-initio PES developed at the University of Minnesota. Good agreement was found between the macroscopic properties and molecular level description of the system obtained by using the two PESs.

Grover, Maninder S.↗

Development of a Detailed Surface Chemistry Framework in DSMC

Many of the current direct simulation Monte Carlo (DSMC) codes still employ only simple surface catalysis models. These include only basic mechanisms such as dissociation, recombination, and exchange reactions, without any provision for adsorption and finite rate kinetics. Incorporating finite rate chemistry at the surface is increasingly becoming a necessity for various applications such as high speed re-entry flows over thermal protection systems (TPS), micro-electro-mechanical systems (MEMS), surface catalysis, etc. In the recent years, relatively few works have examined finite-rate surface reaction modeling using the DSMC method.In this work, a generalized finite-rate surface chemistry framework incorporating a comprehensive list of reaction mechanisms is developed and implemented into the DSMC solver SPARTA. The various mechanisms include adsorption, desorption, Langmuir-Hinshelwood (LH), Eley-Rideal (ER), Collision Induced (CI), condensation, sublimation, etc. The approach is to stochastically model the various competing reactions occurring on a set of active sites. Both gas-surface (e.g., ER, CI) and pure-surface (e.g., LH, desorption) reaction mechanisms are incorporated. The reaction mechanisms could also be catalytic or surface altering based on the participation of the bulk-phase species (e.g., bulk carbon atoms). Marschall and MacLean developed a general formulation in which multiple phases and surface sites are used and we adopt a similar convention in the current work. Microscopic parameters of reaction probabilities (for gas-surface reactions) and frequencies (for pure-surface reactions) that are require for DSMC are computed from the surface properties and macroscopic parameters such as rate constants, sticking coefficients, etc. The energy and angular distributions of the products are decided based on the reaction type and input parameters. Thus, the user has the capability to model various surface reactions via user-specified reaction rate constants, surface properties and parameters.

Gas↗

High-Sensitivity, Broad-Range Vacuum Gauge Using Nanotubes for Micromachined Cavities

A broad-range vacuum gauge has been created by suspending a single-walled carbon nanotube (SWNT) (metallic or semiconducting) in a Schottky diode format or in a bridge conductor format, between two electrically charged mesas. SWNTs are highly sensitive to molecular collisions because of their extremely small diameters in the range of 1 to 3 nanometers. The measurement parameter will be the change in conductivity of SWNT due to decreasing rate of molecular collisions as the pressure inside a chamber decreases. The rate of heat removal approaches a saturation limit as the mean free path (m.f.p.) lengths of molecules increase due to decreasing pressure. Only those sensing elements that have a long relaxation time can produce a measureable response when m.f.p. of molecules increases (or time between two consecutive collisions increases). A suspended SWNT offers such a capability because of its one-dimensional nature and ultrasmall diameter. In the initial approach, similar architecture was used as that of a SWNT-Schottky diode that has been developed at JPL, and has its changing conductivity measured as the test chamber is pumped down from atmospheric pressure to high vacuum (10(exp -7) Torr). Continuous response of decreasing conductivity has been measured as a function of decreasing pressure (SWNT is a negative thermal coefficient material) from atmosphere to less than 10(exp -6) Torr. A measureable current change in the hundreds of nA range has been recorded in the 10(exp -6) Torr regime.

Manohara, Harish↗

First measurement of helium on Mars: Implications for the problem of radiogenic gases on the terrestrial planets

108 +/- 11 photons of the martian He 584-A airglow detected by the Extreme Ultraviolet Explorer (EUVE) satellite during a 2-day exposure (January 22-23, 1993) correspond to the effective disk average intensity of 43 +/- 10 Rayleigh (Ra). Radiative transfer calculations, using a model atmosphere appropriate to the conditions of the observation and having an exospheric temperature of 210 +/- 20 K, result in a He mixing ratio of 1.1 +/- 0.4 ppm in the lower atmosphere. Nonthermal escape of helium is due to electron impact ionization and pickup of He(+) by the solar wind, to collisions with hot oxygen atoms, and to charge exchange with molecular species with corresponding column loss rates of 1.4 x 10(exp 5), 3 x 10(exp 4), and 7 x 10(exp 3)/sq cm/s, respectively. The lifetime of helium on Mars is 5 x 10(exp 4) years. the He outgassing rate, coupled with the Ar-40 atmospheric abundance and with the K:U:Th ratio measured in the surface rocks, is used as input to a single two-reservoir degassing model which is applied to Mars and then to Venus. A similar model with known abundances if K, U, and Th is applied to Earth. The models for Earth and Mars presume loss of all argon accumulated in the atmospheres during the first billion years by large-scale meteorite and planetesimal impacts. The models show that the degassing coefficients for all three planets may be approximated by function delta = delta(sub 0) x (t(sub 0)/t)(exp 1/2) with delta(sub 0) = 0.1, 0.04, and 0.0125 Byr for Earth, Venus, and Mars, respectively. After a R(exp 2) correction this means that outgassing processes on Venus and Mars are weaker than on Earth by factors of 3 and 30, respectively. Mass ratios of U and Th are almost the same for all three planets, while potassiumis depleted by a factor of 2 in Venus and Mars. Mass ratio of helium and argon are close to 5 x 10(exp -9) and 2 x 10(exp -8) g/g in the interiors of all three planets. The implications of these results are discussed.

Krasnopolsky, V. A.↗

Modeling of Sheath Ion-Molecule Reactions in Plasma Enhanced Chemical Vapor Deposition of Carbon Nanotubes

In many plasma simulations, ion-molecule reactions are modeled using ion energy independent reaction rate coefficients that are taken from low temperature selected-ion flow tube experiments. Only exothermic or nearly thermoneutral reactions are considered. This is appropriate for plasma applications such as high-density plasma sources in which sheaths are collisionless and ion temperatures 111 the bulk p!asma do not deviate significantly from the gas temperature. However, for applications at high pressure and large sheath voltages, this assumption does not hold as the sheaths are collisional and ions gain significant energy in the sheaths from Joule heating. Ion temperatures and thus reaction rates vary significantly across the discharge, and endothermic reactions become important in the sheaths. One such application is plasma enhanced chemical vapor deposition of carbon nanotubes in which dc discharges are struck at pressures between 1-20 Torr with applied voltages in the range of 500-700 V. The present work investigates The importance of the inclusion of ion energy dependent ion-molecule reaction rates and the role of collision induced dissociation in generating radicals from the feedstock used in carbon nanotube growth.

Hash, David B.↗

Decision Aid for Conjunction Risk Mitigation by Differential Drag

In the previous five years, the rate of conjunctions that the NASA Conjunction Assessment Risk Analysis (CARA) team processed and analyzed has more than tripled. (NASA CARA, 2024) New missions in the early development phases are now required to plan for dealing with conjunctions under the present space environment, and also projecting forward into a future likely with even further increased utilization of the space environment. Some missions are investigating the possibility of using differential drag to remediate conjunctions without expending limited fuel or for missions without propulsive capabilities. The NASA CARA team studied the historical record of conjunctions to evaluate the circumstances under which differential drag may be successfully applied and have developed a series of tables to use as a decision aid for missions considering differential drag. Currently, if a CARA-protected mission with maneuvering capabilities is predicted to have a conjunction with probability of collision (Pc) greater than 7E-5 (the default value of the ‘yellow threshold’, which may have some other value agreed by CARA and the mission during the Orbital Collison Avoidance Planning (OCAP) process), CARA will use its Maneuver Trade Space (MTS) tool to evaluate and recommend options for the timing and magnitude of a risk mitigation maneuver (RMM), based on the mission’s capabilities. If the conjunction’s Pc is greater than 1E-4, the ‘red threshold’, then an RMM must be executed per NASA Procedural Requirements (NPR) 8079.1 (NASA, 2023), although missions may execute an RMM even if the Pc is lower. The magnitude of the maneuver is typically a few cm/s, and CARA estimates how many will be required for the mission’s nominal lifetime – typically a few per year – during the OCAP process, to inform the mission’s delta-V requirement. For traditional satellites, this is usually smaller than other requirements for orbit insertion, maintenance, and disposal, but for CubeSats or other small satellite missions, a propulsion system may not provide much more than a few cm/s of delta-V or may not fit at all within the available budget of money, time, size, weight, and/or power (SWaP). Conversely, CubeSats often have deployable solar panels, which offer the capacity to have much higher areas facing some directions than others. Such a mission can instead use ‘differential drag’ to remediate a conjunction -- in other words, change its drag area (usually increasing) to deviate from the predicted collision course. This is how Planet’s Dove spacecraft maintain their formations and remediate conjunction risks without having on-board propulsion (Foster, et al., 2017) (Griffith, et al., 2021). CARA has been developing improvements to MTS to support differential-drag for NASA's missions -- where it is effective. CARA records all conjunctions of their protected payloads, with historical records starting in 2005 (with significant conjunction events starting to occur on or after 2013). From this record, approximately 7,300 had a Pc greater than 1E-4 at 3 days prior to the time of closest approach, the time analyzed for differential drag efficacy. Of those, approximately 4,300 had fully-defined covariance matrices stored for both the primary and secondary objects; this set of conjunctions is the basis for the analysis of this work. CARA’s MTS tool was used to propagate the primary satellite forward from that decision point with varying degrees of increase to ballistic coefficient (BC). Because these conjunctions came from multiple missions, the nondimensional ‘delta-BC’ factor was used to quantify and normalize the increase in ballistic coefficient, defined as follows: Delta-BC = BC_new / BC_old - 1 Positive delta-BC factors represent an increase in drag compared to the nominal attitude, while negative delta-BC factors (to a minimum of -1) represent a decrease in drag. At the conclusion of the differential-drag ‘maneuver’, the Pc was recalculated to evaluate whether or not the conjunction was mitigated (Pc < 3E-6). These results were then binned and sorted along several axes, including altitude, amount of delta-BC, and (pre-maneuver) rate of energy dissipation (EDR), to identify underlying patterns. The altitude plot is shown in Figure 1. To validate this analysis, we consulted the record of a NASA mission which uses differential drag to maintain its orbit and remediate conjunction risk. CARA’s empirical record of the mission’s orbit history suggests it achieves a delta-BC of 2.2. Of the twenty-one RMM plans that were submitted by this mission, eighteen were matched with conjunctions in the historical record; of those, twelve were successfully remediated (final measured Pc < 3E-6), and six were not. This is consistent with the expected efficacy for missions orbiting at that altitude. We are presently simulating this mission’s RMMs with MTS; this work is ongoing, but so far, the MTS results are qualitatively in agreement with the empirical results -- correctly predicting that a maneuver would or would not remediate a conjunction, if not exactly matching the final post-remediation Pc value. We found that differential drag was most successful for satellites with perigees below 560 km, and which could adopt an average delta-BC of 2 or greater (that is, increasing their ballistic coefficient by a factor of 3). However, this is a difficult threshold for a mission to clear; very few spacecraft are capable of adopting a high-drag configuration for 72 hours continuously. Planet’s Dove spacecraft use differential drag to remediate conjunctions (Griffith, et al., 2021), and they have a maximum delta-BC factor of 9, but in practice (with mission and charging constraints) they achieve a time-averaged delta-BC that is closer to 2 (Foster, et al., 2017). A mission’s differential drag utility strongly depends on the operational constraints that has the capacity to limit the time-averaged delta-BC. A mission with a high maximum delta-BC of 5 or more can have an effective delta-BC of less than 1 due to operational constraints such as instrument and solar panel pointing, especially if this constraint results in holding an intermediate drag value for most of its orbit. CARA has developed tables that can be used as decision aids to advise missions-in-development about the best way to utilize their differential drag capabilities. For missions below 560 km with the operational flexibility to devote multiple days to holding a high-drag configuration (or a sufficiently high drag ratio to compensate for limitations on that time), they are -- more likely than not -- able to successfully remediate high-risk conjunctions. Conversely, missions that do not meet these exacting criteria -- most missions -- can instead be advised to use on-board propulsion systems to perform RMMs, or to turn their minimum-area face towards the approach vector, thereby reducing Pc at the moment of conjunction due to the decreased Hard-Body Radius (HBR), that is a strongly correlated variable in the Pc calculations. (NASA, 2023)

conjunction assessment↗

Neotectonics of Asia: Thin-shell finite-element models with faults

As India pushed into and beneath the south margin of Asia in Cenozoic time, it added a great volume of crust, which may have been (1) emplaced locally beneath Tibet, (2) distributed as regional crustal thickening of Asia, (3) converted to mantle eclogite by high-pressure metamorphism, or (4) extruded eastward to increase the area of Asia. The amount of eastward extrusion is especially controversial: plane-stress computer models of finite strain in a continuum lithosphere show minimal escape, while laboratory and theoretical plane-strain models of finite strain in a faulted lithosphere show escape as the dominant mode. We suggest computing the present (or neo)tectonics by use of the known fault network and available data on fault activity, geodesy, and stress to select the best model. We apply a new thin-shell method which can represent a faulted lithosphere of realistic rheology on a sphere, and provided predictions of present velocities, fault slip rates, and stresses for various trial rheologies and boundary conditions. To minimize artificial boundaries, the models include all of Asia east of 40 deg E and span 100 deg on the globe. The primary unknowns are the friction coefficient of faults within Asia and the amounts of shear traction applied to Asia in the Himalayan and oceanic subduction zones at its margins. Data on Quaternary fault activity prove to be most useful in rating the models. Best results are obtained with a very low fault friction of 0.085. This major heterogeneity shows that unfaulted continum models cannot be expected to give accurate simulations of the orogeny. But, even with such weak faults, only a fraction of the internal deformation is expressed as fault slip; this means that rigid microplate models cannot represent the kinematics either. A universal feature of the better models is that eastern China and southeast Asia flow rapidly eastward with respect to Siberia. The rate of escape is very sensitive to the level of shear traction in the Pacific subduction zones, which is below 6 MPa. Because this flow occurs across a wide range of latitudes, the net eastward escape is greater than the rate of crustal addition in the Himalaya. The crustal budget is balanced by extension and thinning, primarily within the Tibetan plateau and the Baikal rift. The low level of deviation stresses in the best models suggests that topographic stress plays a major role in the orogeny; thus, we have to expect that different topography in the past may have been linked with fundamentally different modes of continental collision.

Kong, Xianghong↗

Sheet Flows, Avalanches, and Dune Evolution on Earth and Mars

This investigation is a collaboration between researchers at Cornell University, the University of Florida, and the University of Rennes 1, France. Flow modeling at Cornell University focused on mechanisms for the suspension and transport of wind-blown sand that are important in both terrestrial and Martian environments. These mechanisms include the saltation (or jumping) of grains, collisions between grains, and the interaction of grains with the velocity fluctuations of the turbulent wind. Of particular interest are sheet flows; these are relatively thin, highly concentrated regions of grains flowing near the ground under the influence of a strong turbulent wind. In them, the grains are suspended by interparticle collisions. Sheet flows may be relatively rare events, but they have the capacity to move great amounts of sand. In order to describe sheet flows, a turbulent mixture theory was formulated for particles in a fluid in which fluctuations in the volume fiaction of the particles take place on the scale of the turbulent eddies. Ensemble averaged equations for particle and fluid mass, momentum, and energy and fluid rate of dissipation were expressed in terms of Farve (concentration) averaged velocities and the associated velocity fluctuations. Correlations that describe the turbulent suspension of particles and dissipation of turbulent energy of both phases due to fluid particle interactions were modeled and boundary conditions at the bed and at the upper surface of the collisional flow were formulated. The boundary conditions at the upper surface were tested in a numerical simulation developed at the University of Florida. Steady and unsteady solutions for steady and unsteady fully-developed flows were obtained over a range of wind speeds fiom the lowest for which collisional between particles occurred to at which turbulent suspension is found to dominate collisional suspension. Below the value of the wind speed at which collisions between particles were unimportant, numerical solutions were obtained for the velocity distribution function and the resulting fields of concentration, particle and gas mean velocity, and particle shear stress for the steady two-dimensional saltation of spherical sand particles driven by a turbulent wind over a bed characterized by a simple relationship (the splash function) between the properties of incoming particles and those of the rebounding particles and other particles ejected fiom the bed. At the University of Rennes 1, experiments devoted to the characterization of the splash function for beds consisting of either random or ordered arrays of spheres in two- dimensions were completed. These indicated the role played by the packing geometry in the rebound and ejection of grains. Preliminary experiments on response of a three- dimensional collision bed to a collision with a single particle were performed. Data was taken with a single camera focused on the plane of collision. Here, for example, the decrease of the effective coefficient of restitution of the bed with an increase of the angle of incidence of the incoming particle has been measured. Other experiments on avalanches at Rennes studied the properties of the flows of particles that are responsible for the motion of the leeward side of a dune. In these, the dependence of the initiation of avalanches on the packing and depth of the particles was measured. Particle migration was studied in inclined flows of a binary mixture of disks and the mechanisms of diffision and segregation were isolated and characterized. The influence of side wall on dense, rapid inclined flows was measured and shown to be the reason why the angle of the free surface in such flows can exceed the static angle of repose. Future research will be devoted to a better understanding the transition between saltating (collisionless) and collisional flows as the wind speed the increases. This will involve the understanding of the evolution of the splash function as clisions with the bed become more numerous, more frequent, and more violent.

Source record↗