Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetically constrained models”

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 37 records · Page 2

Investigating Global Ion and Neutral Atom Populations with IBEX and Voyager

The main objective of this project was to investigate pickup ion (PUI) production in the solar wind and heliosheath (the region between the termination shock and the heliopause) and compute the distributed energetic neutral atom fluxes throughout the helioshpere. The simulations were constrained by comparing the model output against observations from Ulysses, New Horizons, Voyager 1 and 2, and IBEX space probes. As evidenced by the number of peer reviewed journal publications resulting from the project (13 plus three submitted) and their citation rate (156 citations over three years), the project has made a lasting contribution to the field. The outcome is a significant improvement of our understanding of the pickup ion production and distribution in the distant heliosphere. The team has accomplished the entire set of tasks A-H set forth in the proposal. Namely, the transport modeling framework has been augmented with two populations of pickup ions (PUIs), the boundary conditions for the plasma and interstellar neutral hydrogen were verified against Ulysses and New Horizons PUI and an optimal set of velocity diffusion parameters established. The multi-component fluxes of PUIs were computed and isotropic velocity distributions generated for each cell in the computer simulation that covered the heliosphere from 1.5 AU to the heliopause. The distributions were carefully compared with in situ measurements at 3 AU (Ulysses), 12 AU (New Horizons), and 80-90 AU (Voyager 1 and 2) as well as those inferred from ENA fluxes measured by Cassini and IBEX (Wu et al., 2016). Some examples of modeldata comparison are shown in Figure 1. We have used coupled MHD-plasma and kinetic-neutral code to investigate the likely range of plasma and magnetic field parameters in the local interstellar medium (LISM), based on the assumption that the shape of the IBEX ribbon could be used to determine the orientation of the interstellar magnetic field. While the magnetic field is believed to be oriented toward the center of the ribbon, constraining its strength requires comparing the model-predicted angular diameter and circularity of the ribbon with the observations. The study, published in Heerikhuisen et al. (2014), found that the most likely range for the LISM magnetic field strength is between 0.2 and 0.3 nT, which is less than previously thought. Figure 2 shows the IBEX data (left) and compares it to the simulation with a 0.2 nT interstellar magnetic field (center) and a 0.4 nT (right).

Florinski, Vladimir↗

An axisymmetric magnetohydrodynamic model for the Crab pulsar wind bubble

We extend Kennel and Coroniti's (1984) spherical magnetohydrodynamic models for the Crab Nebula to include the pinching effect of the toroidal magnetic field. Since the bulk nebular flow is likely to be very submagnetosonic, a quasi-static treatment is possible. We show that the pinching effect can be responsible for the observed elongation of the pulsar wind bubble, as indicated by the surface brightness contours of optical synchrotron radiation. From the observed elongation we estimate a value for sigma, the ratio of Poynting flux to plasma kinetic energy flux in the free pulsar wind, which is consistent with previous results from spherical models. Using the inferred magnetic field configuration inside the pulsar wind bubble, combined with the observed dimensions of the X-ray nebula, we are able to constrain the particle distribution function. We conclude that, for a power-law injection function, the maximum energy has to be much larger in the pulsar equatorial region than in the polar region.

Begelman, Mitchell C.↗

Transition from Gaseous Compounds to Aerosols in Titan's Atmosphere

We investigate the chemical transition of simple molecules like C2H2 and HCN into aerosol particles in the context of Titan's atmosphere. Experiments that synthesize analogs (tholins) for these aerosols can help understand and constrain these polymerization mechanisms. Using information available from these experiments, we suggest chemical pathways that can link simple molecules to macromolecules, that will be the precursors to aerosol particles: polymers of acetylene and cyanoacetylene, polycyclic aromatics (PAHs), polymers of HCN and other nitriles, and polynes. Although our goal here is not to build a detailed kinetic model for this transition, we propose parameterizations to estimate the production rates of these macromolecules, their C/N and C/H ratios, and the loss of parent molecules (C2H2, HCN, HC3N and other nitriles, C6H6) from the gas phase to the haze. We use a 1-dimensional photochemical model of Titan's atmosphere to estimate the formation rate of precursors macromolecules. We find a production zone slightly lower than 200 km altitude with a total production rate of 4 x 10(exp -14) g/ sq cm s and a C/N approx. = 4. These results are compared with experimental data, and to microphysical models requirements. The Cassini/Huygens mission will bring a detailed picture of the haze distribution and properties, that will be a great challenge for our understanding of those chemical processes.

Lebonnois, Sebastien↗

Cluster/Peace Electrons Velocity Distribution Function: Modeling the Strahl in the Solar Wind

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the electron heat-flux and temperature anisotropy. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by fitting the VDF's to a spectral spherical harmonic model from which the moments are derived directly from the model coefficients.

Figueroa-Vinas, Adolfo↗

Modeling the Strahl in the Solar Wind Using CLUSTER/PEACE Electron Distribution Functions

We present a study of kinetic properties of the strahl electron velocity distribution functions (VDF's) in the solar wind. These are used to investigate the pitch-angle scattering and stability of the population to interactions with electromagnetic (whistler) fluctuations. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer. Our study focuses on the mechanisms that control and regulate the pitch-angle and stability of strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. The goal is to check whether the strahl electrons are constrained by some instability (e.g., the whistler instability), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling the 3D-VDFs to a spectral spherical harmonic model from which the moments are obtained directly from the spectral coefficients.

Vinas, Adolfo F.↗

Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars

On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosityin Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, τ, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past Martian climate and environments.

Gale crater↗

Reactive Transport and Mass Balance Modeling of the Stimson Sedimentary Formation and Altered Fracture Zones Constrain Diagenetic Conditions at Gale Crater, Mars

On a planet as cold and dry as present-day Mars, evidence of multiple aqueous episodes offers an intriguing view into very different past environments. Fluvial, lacustrine, and eolian depositional environments are being investigated by the Mars Science Laboratory Curiosity in Gale crater, Mars. Geochemical and mineralogical observations of these sedimentary rocks suggest diagenetic processes affected the sediments. Here, we analyze diagenesis of the Stimson formation eolian parent material, which caused loss of olivine and formation of magnetite. Additional, later alteration in fracture zones resulted in preferential dissolution of pyroxene and precipitation of secondary amorphous silica and Ca sulfate. The ability to compare the unaltered parent material with the reacted material allows constraints to be placed on the characteristics of the altering solutions. In this work we use a combination of a mass balance approach calculating the fraction of a mobile element lost or gained, tau, with fundamental geochemical kinetics and thermodynamics in the reactive transport code CrunchFlow to examine the characteristics of multiple stages of aqueous alteration at Gale crater, Mars. Our model results indicate that early diagenesis of the Stimson sedimentary formation is consistent with leaching of an eolian deposit by a near-neutral solution, and that formation of the altered fracture zones is consistent with a very acidic, high sulfate solution containing Ca, P and Si. These results indicate a range of past aqueous conditions occurring at Gale crater, Mars, with important implications for past martian climate and environments.

Gale↗

Asteroidal Space Weathering: The Major Role of FeS

Space weathering (SW) effects on the lunar surface are reasonably well-understood from sample analyses, remote-sensing data, and experiments, yet our knowledge of asteroidal SW effects are far less constrained. While the same SW processes are operating on asteroids and the Moon, namely solar wind irradiation, impact vaporization and condensation, and impact melting, their relative rates and efficiencies are poorly known, as are their effects on such vastly different parent materials. Asteroidal SW models based on remote-sensing data and experiments are in wide disagreement over the dominant mechanisms involved and their kinetics. Lunar space weathering effects observed in UVVIS-NIR spectra result from surface- and volume-correlated nanophase Fe metal (npFe(sup 0)) particles. In the lunar case, it is the tiny vapor-deposited npFe(sup 0) that provides much of the spectral reddening, while the coarser (largely melt-derived) npFe(sup 0) produce lowered albedos. Nanophase FeS (npFeS) particles are expected to modify reflectance spectra in much the same way as npFe(sup 0) particles. Here we report the results of experiments designed to explore the efficiency of npFeS production via the main space weathering processes operating in the asteroid belt.

Keller, L. P.↗

Triple Oxygen Isotope Kinetics of Sulfate From Pyrite Oxidation

Sulfates from low temperature terrestrial and marine environments are out of equilibrium with their oxygen sources. Instead, oxygen isotope compositions of sulfate imply they are dominated by kinetic signatures of oxidative and reductive processes (Killingsworth et al., 2022) that are under-constrained compared to equilibrium effects (Hayles and Killingsworth, 2022). Consequently, the specific role of pyrite oxidation in setting the oxygen isotope composition of sulfate is also under-constrained. This scenario inhibits a more nuanced understanding of small variations in the oxygen isotopes of sulfate and their implications for tracing oxygen sources and sulfur cycle transformations, particularly in low-temperature relatively sulfate-limited freshwater systems such as rivers, lakes, aquifers, and mining-affected environments. To alleviate these knowledge gaps, we conduct pyrite oxidation experiments at time scales spanning minutes to years, analyze the triple oxygen isotope compositions (∆ 17 O and δ 18 O) of sulfate products, and interpret results using molecular modeling. Our preliminary results show isotopic relationships of 1) pyrite-derived sulfate ∆ 17 O reaching above the meteoric water line as its δ 18 O becomes lighter than water, and 2) sulfate ∆ 17 O and δ 18 O extending between water and O 2 . From these results we interpret a kinetic isotope effect between sulfate and water and mixing between sulfate end members containing most of their oxygen atoms from water and from O 2 .

Oxygen isotopes↗

Opaque Minerals in the Matrix of the Bishunpur (LL3.1) Chondrite: Constraints on the Chondrule Formation Environment

The chemistry and mineralogy of a group of opaque mineral assemblages in the matrix of the Bishunpur LL3.1 ordinary chondrite provide insight into the nebular environment in which they formed. The assemblages consist of a kamacite (Fe,Ni) core that is rimmed by troilite (FeS) and fayalite (Fe2Si04). Accessory phases in the rims include silica (Si02), chromite (FeCr204), whitlockite (Ca3(P04)2), maricite (FeNaP04), magnetite (Fe304), and tetrataenite (FeNi). We suggest that the metal melted in and equilibrated with an igneous chondrule under high- temperature, reducing conditions. In this environment the molten alloys incorporated varied amounts of Si, Ni, P, Cr, and Co, depending on the oxygen fugacity and temperature of the melt. Some of the metal was subsequently expelled from the chondrule interiors into the surrounding nebular gas. As the temperature dropped, the alloy solidified and volatile elements corroded the metal. The main reaction products were troilite and fayalite. Thermodynamic equilibrium calculations are used to constrain the conditions under which these two phases can form simultaneously in the solar nebula. Kinetic factors are used to place a lower limit on the formation temperature. We determine that the metal corroded between 1173 and 1261 K at a total pressure in the range of 1025.0 to 1024.1 bars and a dust/gas ratio of 302 to 355 x relative to solar composition. These conditions are consistent with our model that the metal corroded in a dust- rich region of the solar nebula that was cooling after a chondrule formation event.

Lauretta, D. S.↗

Optical measurements of winds and kinetic temperatures in the upper atmosphere

A development history is presented for the optical probing of the structure of the upper atmosphere, noting the contributions made by both ground and spaceborne optical sensors, and the role of advancements in thermospheric theory. The significance of optical techniques to current understanding of the upper atmosphere is discussed in the context of the new COSPAR International Reference Atmosphere-1986. While much of the work performed to date has involved the analysis of data obtained by individual observatories, the ability to effectively constrain theoretical models with such single-station data has become more limited with increasing model sophistication.

Hernandez, G.↗

Kinetic Simulations of Current-Sheet Formation and Reconnection at a Magnetic X Line

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the plasma physics community, particularly in the context of magnetic reconnection. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high resolution particle-in-cell (PIC) simulations. The initial potential magnetic field is perturbed by thermal pressure introduced into the particle distribution far from the X line. The relaxation of this added stress leads to the development of a current sheet, which reconnects for imposed stress of sufficient strength. We compare the evolution and final state of our PIC simulations with magnetohydrodynamic simulations assuming both uniform and localized resistivities, and with force-free magnetic-field equilibria in which the amount of reconnect ion across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the terrestrial magnetotail and solar corona.

Black, C.↗

Characterization of the High-Albedo NEA 3691 Bede

Characterization of NEAs provides important inputs to models for atmospheric entry, risk assessment and mitigation. Diameter is a key parameter because diameter translates to kinetic energy in atmospheric entry. Diameters can be derived from the absolute magnitude, H(PA=0deg), and from thermal modeling of observed IR fluxes. For both methods, the albedo (pv) is important - high pv surfaces have cooler temperatures, larger diameters for a given Hmag, and shallower phase curves (larger slope parameter G). Thermal model parameters are coupled, however, so that a higher thermal inertia also results in a cooler surface temperature. Multiple parameters contribute to constraining the diameter. Observations made at multiple observing geometries can contribute to understanding the relationships between and potentially breaking some of the degeneracies between parameters. We present data and analyses on NEA 3691 Bede with the aim of best constraining the diameter and pv from a combination of thermal modeling and light curve analyses. We employ our UKIRT+Michelle mid-IR photometric observations of 3691 Bede's thermal emission at 2 phase angles (27&43 deg 2015-03-19 & 04-13), in addition to WISE data (33deg 2010-05-27, Mainzer+2011). Observing geometries differ by solar phase angles and by moderate changes in heliocentric distance (e.g., further distances produce somewhat cooler surface temperatures). With the NEATM model and for a constant IR beaming parameter (eta=constant), there is a family of solutions for (diameter, pv, G, eta) where G is the slope parameter from the H-G Relation. NEATM models employing Pravec+2012's choice of G=0.43, produce D=1.8 km and pv≈0.4, given that G=0.43 is assumed from studies of main belt asteroids (Warner+2009). We present an analysis of the light curve of 3691 Bede to constrain G from observations. We also investigate fitting thermophysical models (TPM, Rozitis+11) to constrain the coupled parameters of thermal inertia (Gamma) and surface roughness, which in turn affect diameter and pv. Surface composition can be related to pv. This study focuses on understanding and characterizing the dependency of parameters with the aim of constraining diameter, pv and thermal inertia for 3691 Bede.

Wooden, Diane H.↗

Characterization of the High-Albedo NEA 3691 Bede

Characterization of NEAs provides important inputs to models for atmospheric entry, risk assessment and mitigation. Diameter is a key parameter because diameter translates to kinetic energy in atmospheric entry. Diameters can be derived from the absolute magnitude, H(PA=0deg), and from thermal modeling of observed IR fluxes. For both methods, the albedo (pv) is important - high pv surfaces have cooler temperatures, larger diameters for a given Hmag, and shallower phase curves (larger slope parameter G). Thermal model parameters are coupled, however, so that a higher thermal inertia also results in a cooler surface temperature. Multiple parameters contribute to constraining the diameter. Observations made at multiple observing geometries can contribute to understanding the relationships between and potentially breaking some of the degeneracies between parameters. We present data and analyses on NEA 3691 Bede with the aim of best constraining the diameter and pv from a combination of thermal modeling and light curve analyses. We employ our UKIRT+Michelle mid-IR photometric observations of 3691 Bede's thermal emission at 2 phase angles (27&43 deg 2015-03-19 & 04-13), in addition to WISE data (33deg 2010-05-27, Mainzer+2011). Observing geometries differ by solar phase angles and by moderate changes in heliocentric distance (e.g., further distances produce somewhat cooler surface temperatures). With the NEATM model and for a constant IR beaming parameter (eta=constant), there is a family of solutions for (diameter, pv, G, eta) where G is the slope parameter from the H-G Relation. NEATM models employing Pravec+2012's choice of G=0.43, produce D=1.8 km and pv≈0.4, given that G=0.43 is assumed from studies of main belt asteroids (Warner+2009). We present an analysis of the light curve of 3691 Bede to constrain G from observations. We also investigate fitting thermophysical models (TPM, Rozitis+11) to constrain the coupled parameters of thermal inertia (Gamma) and surface roughness, which in turn affect diameter and pv. Surface composition can be related to pv. This study focuses on understanding and characterizing the dependency of parameters with the aim of constraining diameter, pv and thermal inertia for 3691 Bede.

magnitude↗

Kinetics of the BrO + NO2 Association Reaction. Temperature and Pressure Dependence in the Falloff Regime

A laser flash photolysis-long path absorption technique has been employed to study the kinetics of the reaction BrO + NO2 + M yields (k1) products as a function of temperature (248-346 K), pressure (16-800 torr), and buffer gas identity (N2,CF4) The reaction is found to be in the falloff regime between third and second-order over the entire range of conditions investigated This is the first study where temperature-dependent measurements of k1(P,T) have been reported at pressures greater than 12 torr; hence, our results help constrain choices of k1(P,T) for use in models of lower stratospheric BrO(x) chemistry. Approximate falloff parameters in a convenient form for atmospheric modeling are derived.

Thron, R. P.↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗

Current-Sheet Formation and Reconnection at a Magnetic X Line in Particle-in-Cell Simulations

The integration of kinetic effects into macroscopic numerical models is currently of great interest to the heliophysics community, particularly in the context of magnetic reconnection. Reconnection governs the large-scale energy release and topological rearrangement of magnetic fields in a wide variety of laboratory, heliophysical, and astrophysical systems. We are examining the formation and reconnection of current sheets in a simple, two-dimensional X-line configuration using high-resolution particle-in-cell (PIC) simulations. The initial minimum-energy, potential magnetic field is perturbed by excess thermal pressure introduced into the particle distribution function far from the X line. Subsequently, the relaxation of this added stress leads self-consistently to the development of a current sheet that reconnects for imposed stress of sufficient strength. We compare the time-dependent evolution and final state of our PIC simulations with macroscopic magnetohydrodynamic simulations assuming both uniform and localized electrical resistivities (C. R. DeVore et al., this meeting), as well as with force-free magnetic-field equilibria in which the amount of reconnection across the X line can be constrained to be zero (ideal evolution) or optimal (minimum final magnetic energy). We will discuss implications of our results for understanding magnetic-reconnection onset and cessation at kinetic scales in dynamically formed current sheets, such as those occurring in the solar corona and terrestrial magnetotail.

Black, C.↗

The physics and chemistry of the L134N molecular core

The dark cloud L134N is studied in detail via millimeter- and centimeter-wavelength emission-line spectra. A high-density core of molecular gas exists in L134N which has a kinetic temperature of about 12 K, a peak molecular hydrogen density of about 10 exp 4.5/cu cm, and a mass of about 23 solar. The core may be the site of future star formation. Maps of emission from (C-18)O, CS, H(C-13)O(+), SO, NH3, and C3H2 reveal morphologically different distributions resulting in part from both varying physical conditions within the cloud and optical depth effects. Significant differences also exist which are probably due to chemical abundance variations. A consistent set of LTE chemical abundances has been estimated at as many as seven positions, which can be used to constrain chemical models of dark clouds.

Swade, Daryl A.↗