SEARCH · Engineering Papers

Results for “parameter space exploration”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Pilot Sensorimotor and Cardiovascular Results from the Joint Russian/U.S. Field Test

The primary goal of this research is to determine functional abilities associated with long-duration space flight crews beginning as soon after landing as possible (< 2 hours) with an additional two follow-up measurements sessions on the day of landing. This goal has both sensorimotor and cardiovascular elements, including evaluations of NASA's new anti-orthostatic compression garment and the Russian Kentavr garment. Functional sensorimotor measurements will include, but are not limited to, assessing hand/eye coordination, standing from a seated position (sit-to-stand), walking normally without falling, measurement of dynamic visual acuity, discriminating different forces generated with both the hands and legs, recovering from a fall (standing from a prone position), coordinated walking involving tandem heel-to-toe placement, and determining postural ataxia while standing. The cardiovascular portion of the investigation includes measuring blood pressure and heart rate during a timed stand test in conjunction with postural ataxia testing (quiet stance sway) as well as cardiovascular responses during the other functional tasks. In addition to the immediate post-landing collection of data for the full FT, postflight data is being acquired twice more within the 24 hours after landing and will continue over the subsequent weeks until functional sensorimotor and cardiovascular responses have returned to preflight normative values. The PFT represents a initial evaluation of the feasibility of testing in the field, and is comprised of a jointly agreed upon subset of tests from the full FT and relies heavily on Russia's Institute of Biomedical Problems Sensory-Motor and Countermeasures Laboratories for content and implementation. The PFT has been collected on several ISS missions. Testing on the U.S. side has included: (1) a sit-to-stand test, (2) recovery from a fall where the crewmember began in the prone position on the ground and then stood for 3 minutes while cardiovascular stability was determined and postural ataxia data were acquired, and (3) a tandem heel-to-toe walk test to determine changes in the central locomotor program. Video, cardiovascular parameters (heart rate and blood pressure), data from bodyworn inertial sensors, and severity of postflight motion sickness were collected during each test session. Our Russian investigators have added measurements associated with: (a) obstacle avoidance, (b) muscle compliance and (c) postural adjustments to perturbations (push) applied to the subject's chest area. The level of functional deficit observed in the crew tested to date is typically beyond what was expected and is clearly triggered by the acquisition of gravity loads immediately after landing when the demands for crew intervention in response to emergency operations will be greatest. Clearly measureable performance parameters such as ability to perform a seat egress, recover from a fall or the ability to see clearly when walking, and related physiologic data (orthostatic responses) are required to provide an evidence base for characterizing programmatic risks and the degree of variability among crewmembers for exploration missions where the crew will be unassisted after landing. Overall, these early functional and related physiologic measurements will allow the estimation of nonlinear sensorimotor and cardiovascular recovery trends that have not been previously captured

Reschke, M. F.↗

NASA NTRS1999

Microgravity Transport Phenomena Experiment (MTPE) Overview

The Microgravity Transport Phenomena Experiment (MTPE) is a fluids experiment supported by the Fundamentals in Biotechnology program in association with the Human Exploration and Development of Space (BEDS) initiative. The MTP Experiment will investigate fluid transport phenomena both in ground based experiments and in the microgravity environment. Many fluid transport processes are affected by gravity. Osmotic flux kinetics in planar membrane systems have been shown to be influenced by gravimetric orientation, either through convective mixing caused by unstably stratified fluid layers, or through a stable fluid boundary layer structure that forms in association with the membrane. Coupled transport phenomena also show gravity related effects. Coefficients associated with coupled transport processes are defined in terms of a steady state condition. Buoyancy (gravity) driven convection interferes with the attainment of steady state, and the measurement of coupled processes. The MTP Experiment measures the kinetics of molecular migration that occurs in fluids, in response to the application of various driving potentials. Three separate driving potentials may be applied to the MTP Experiment fluids, either singly or in combination. The driving potentials include chemical potential, thermal potential, and electrical potential. Two separate fluid arrangements are used to study membrane mediated and bulk fluid transport phenomena. Transport processes of interest in membrane mediated systems include diffusion, osmosis, and streaming potential. Bulk fluid processes of interest include coupled phenomena such as the Soret Effect, Dufour Effect, Donnan Effect, and thermal diffusion potential. MTP Experiments are performed in the Microgravity Transport Apparatus (MTA), an instrument that has been developed specifically for precision measurement of transport processes. Experiment fluids are contained within the MTA fluid cells, designed to create a one dimensional flow geometry of constant cross sectional area, and to facilitate fluid filling and draining operations in microgravity. The fluid cells may be used singly for bulk solutions, or in a Stokes diaphragm configuration to investigate membrane mediated phenomena. Thermal and electrical driving potentials are applied to the experiment fluids through boundary plates located at the ends of the fluid cells. In the ground based instrument, two constant temperature baths circulate through reservoirs adjacent to the boundary plates, and establish the thermal environment within the fluid cells. The boundary plates also serve as electrodes for measurement and application of electrical potentials. The Fluid Manipulation System associated with the MTA is a computer controlled system that enables storage and transfer of experiment fluids during on orbit operations. The system is used to automatically initiate experiments and manipulate fluids by orchestrating pump and valve operations through scripted sequences. Unique technologies are incorporated in the MTA for measurement of fluid properties. Volumetric Flow Sensors have been developed for precision measurement of total fluid volume contained within the fluid cells over time. This data is most useful for measuring the kinetics of osmosis, where fluid is transported from one fluid cell to another through a semipermeable membrane. The MicroSensor Array has been designed to perform in situ measurement of several important fluid parameters, providing simultaneous measurement of solution composition at multiple locations within the experiment fluids. Micromachined sensors and interface electronics have been developed to measure temperature, electrical conductivity, pH, cation activity, and anion activity. The Profile Refractometer uses a laser optical system to directly image the fluid Index of Refraction profile that exists along the MTA fluid cell axis. A video system acquires images of the RI profile over time, and records the transport kinetics that occur upon application of chemical, thermal, or electrical driving potentials. Image processing algorithms have been developed to analyze the refractometer images on a pixel by pixel basis, calibrating and scaling the measured Index of Refraction profile to correlated solution properties of interest such as density, concentration, and temperature. Additional software has been developed to compile the processed images into a three dimensional matrix that contains fluid composition data as a function of experiment time and position in the fluid cell. These data are combined with data from the other sensor systems, and analyzed in the context of transport coefficients associated with the various transport phenomena. Analysis protocols have been developed to measure the transient kinetics, and steady state distribution of fluid components that occur in response to the applied driving potentials. The results are expressed in terms of effective transport coefficients. Experiments have been performed using a variety of solutes, and results generated are that are in agreement with published transport coefficient values.

Mason, Larry W.↗

DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2011

Z-Pinch Magneto-Inertial Fusion Propulsion Engine Design Concept

Fusion-based nuclear propulsion has the potential to enable fast interplanetary transportation. Due to the great distances between the planets of our solar system and the harmful radiation environment of interplanetary space, high specific impulse (Isp) propulsion in vehicles with high payload mass fractions must be developed to provide practical and safe vehicles for human spaceflight missions. Magneto-Inertial Fusion (MIF) is an approach which has been shown to potentially lead to a low cost, small fusion reactor/engine assembly (1). The Z-Pinch dense plasma focus method is an MIF concept in which a column of gas is compressed to thermonuclear conditions by an estimated axial current of approximately 100 MA. Recent advancements in experiments and the theoretical understanding of this concept suggest favorable scaling of fusion power output yield as I(sup 4) (2). The magnetic field resulting from the large current compresses the plasma to fusion conditions, and this is repeated over short timescales (10(exp -6) sec). This plasma formation is widely used in the field of Nuclear Weapons Effects (NWE) testing in the defense industry, as well as in fusion energy research. There is a wealth of literature characterizing Z-Pinch physics and existing models (3-5). In order to be useful in engineering analysis, a simplified Z-Pinch fusion thermodynamic model was developed to determine the quantity of plasma, plasma temperature, rate of expansion, energy production, etc. to calculate the parameters that characterize a propulsion system. The amount of nuclear fuel per pulse, mixture ratio of the D-T and nozzle liner propellant, and assumptions about the efficiency of the engine, enabled the sizing of the propulsion system and resulted in an estimate of the thrust and Isp of a Z-Pinch fusion propulsion system for the concept vehicle. MIF requires a magnetic nozzle to contain and direct the nuclear pulses, as well as a robust structure and radiation shielding. The structure, configuration, and materials of the nozzle must meet many severe requirements. The configuration would focus, in a conical manner, the Deuterium-Tritium (D-T) fuel and Lithium-6/7 liner fluid to meet at a specific point that acts as a cathode so the Li-6 can serve as a current return path to complete the circuit. In addition to serving as a current return path, the Li liner also serves as a radiation shield. The advantage to this configuration is the reaction between neutrons and Li-6 results in the production of additional Tritium, thus adding further fuel to the fusion reaction and boosting the energy output. To understand the applicability of Z-Pinch propulsion to interplanetary travel, it is necessary to design a concept vehicle that uses it. The propulsion system significantly impacts the design of the electrical, thermal control, avionics, radiation shielding, and structural subsystems of a vehicle. The design reference mission is the transport of crew and cargo to Mars and back, with the intention that the vehicle be reused for other missions. Several aspects of this vehicle are based on a previous crewed fusion vehicle study called Human Outer Planet Exploration (HOPE), which employed a Magnetized Target Fusion (MTF) propulsion concept. Analysis of this propulsion system concludes that a 40-fold increase of Isp over chemical propulsion is predicted. This along with a greater than 30% predicted payload mass fraction certainly warrants further development of enabling technologies. The vehicle is designed for multiple interplanetary missions and conceivably may be suited for an automated one-way interstellar voyage.

Miernik, Janie H.↗

NASA NTRS2002

A Chandra X-ray Study of Cygnus A: The Nucleus - 2

We report Chandra Advanced CCD Imaging Spectrometer and quasi-simultaneous Rossi X-Ray Timing Explorer (RXTE) observations of the nearby, powerful radio galaxy Cygnus A, with the present paper focusing on the properties of the active nucleus. In the Chandra observation, the hard (less than a few keV) X-ray emission is spatially unresolved with a size is approximately 1" (1.5 kpc, H(sub 0) = 50 km/s/Mpc) and coincides with the radio and near-infrared nuclei. In contrast, the soft (less than 2 keV) emission exhibits a bipolar nebulosity that aligns with the optical bipolar continuum and emission-line structures and approximately with the radio jet. In particular, the soft X-ray emission corresponds very well with the [O III] (lambda)5007 and H(alpha) + [N II] lambda(lambda)6548, 6583 nebulosity imaged with Hubble Space Telescope. At the location of the nucleus, there is only weak soft X-ray emission, an effect that may be intrinsic or result from a dust lane that crosses the nucleus perpendicular to the source axis. The spectra of the various X-ray components have been obtained by simultaneous fits to the six detectors. The compact nucleus is detected to 100 keV and is well described by a heavily absorbed power-law spectrum with Gamma(sub h) = 1.52(sup + 0.12, sub -0.12) (similar to other 0.12 narrow-line radio galaxies) and equivalent hydrogen column N(sub H)(nuc) = 2.0(sup +0.1, sub -0.1) x 10(exp 23)/sq cm. This 0.2 column is compatible with the dust obscuration to the near-infrared source for a normal gas-to-dust ratio. The soft (less than 2 keV) emission from the nucleus may be described by a power-law spectrum with the same index (i.e., Gamma(sub l) = Gamma(sub h), although direct fits suggest a slightly larger value for Gamma(sub l). Narrow emission lines from highly ionized neon and silicon, as well as a "neutral" Fe K(alpha) line, are detected in the nucleus and its vicinity (r approximately less than 2 kpc). The equivalent width (EW) of the Fe K(alpha) line (182(sup +40, sub -54) eV) is in good agreement with theoretical predictions for the EW versus N(sub H)(nuc) relationship in various geometries. An Fe K edge is also seen. The RXTE observations indicate a temperature of kT = 6.9(sup +0., sub -1.0) keV for the cluster gas (discussed in Paper III of this series) and cluster emission lines of Fe K(alpha) and Fe K(beta) and/or Ni K(alpha). We consider the possibility that the extended soft X-ray emission is electron-scattered nuclear radiation. Given that 1% of the unabsorbed 2 - 10 keV nuclear radiation would have to be scattered, the necessary gas column [N(sub H)(Scattering) approx. = 3.5 x 10(exp 22)/sq cm] would absorb the X-rays rather than scatter them if the gas is cold. Thus, the scattering plasma must be highly ionized. If this ionization is achieved through photoionization by the nucleus, the ionization parameter zeta greater than 1 ergs cm/s and the electron density n(sub e) approx. = 6 cc given the observed distance of the soft X-ray emission from the nucleus. The electron column density inferred from the X-ray observations is much too low to account for the extended optical scattered light, strongly suggesting that the polarized optical light is scattered by dust. The presence of highly ionized Ne lines in the soft X-ray spectrum requires 20 ergs cm/s approximately less than zeta approximately less than 300 ergs cm/s these lines may originate closer to the nucleus than the extended soft continuum or in a lower density gas. A collisionally ionized thermal model of the extended soft X-rays cannot be ruled out but is unattractive in view of the low metal abundance required (Z = 0.03 Z(mass)). The hard X-ray to far-infrared ratio for the nucleus of Cygnus A is similar to that seen in Seyfert 1 and unobscured radio galaxies. By means of the correlation between hard X-ray luminosity and nuclear optical absolute magnitude for these classes of object, we estimate M(sub B) = -22.4 for Cygnus A, near the .borderline between Seyfert galaxies and QSOs.

Young, Andrew J.↗

DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗

DOE OSTI2023

Polyhedral Distortions and Unusual Magnetic Order in Spinel FeMn 2 O 4

Spinel compounds AB 2 X 4 consist of both tetrahedral (AX 4 ) and octahedral (BX 6) environments with the former forming a diamond lattice and the latter a geometrically frustrated pyrochlore lattice. Exploring the fascinating physical properties and their correlations with structural features is critical in understanding these materials. FeMn 2 O 4 has been reported to exhibit one structural transition and two successive magnetic transitions. In this work, we report the polyhedral distortions and their correlations to the structural and two magnetic transitions in FeMn 2 O 4 by employing the high-resolution neutron powder diffraction. The cation distribution is found to be ($Mn_{0.9}^{2+}Fe_{0.1}^{3+}$)$_{A}$($Mn^{3+}Fe_{0.9}^{3+}Mn_{0.1}^{2+}$)$_{B}$O 4 . While large trigonal distortion is found even in the high-temperature cubic phase, the first-order cubic-tetragonal structural transition associated with the elongation of both tetrahedra and octahedra with shared oxygen atoms along the c axis occurs at T S ≈ 750 K, driven by the Jahn–Teller effect of the orbital active B-site Mn 3+ cation. Strong magnetoelastic coupling is unveiled at T N1 ≈ 400 K as manifested by the appearance of Néel-type collinear ferrimagnetic order, an anomaly in both tetrahedral and octahedral distortions, as well as an anomalous decrease of the lattice constants c and a weak anomaly of a. Upon cooling to T N2 ≈ 65 K, it evolves to a noncollinear ferrimagnetic order accompanied by the different moments at the split magnetic sites B1 and B2. Only one-half of the B-site Mn 3+ /Fe 3+ spins, i.e., the B2-site spins in the pyrochlore lattice, are canted, which is a unique magnetic order among spinels. The canting angle between A-site and B2-site moments is ~25°, but the B1-site moment stays antiparallel to the A-site moment even at 10 K. This noncollinear order is accompanied by a modification of the O–B–O bond angles in the octahedra without significant change in lattice constants or tetrahedral/octahedral distortion parameters, indicating a distinct magnetoelastic coupling. We demonstrate distinct roles of the A-site and B-site magnetic cations in the structural and magnetic properties of FeMn 2 O 4 . Our study indicates that FeMn 2 O 4 is a wonderful platform to unveil interesting magnetic order and to investigate their correlations with polyhedral distortions and lattice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Backlighting extended gas halos around luminous red galaxies: Kinematic Sunyaev-Zel’dovich effect from DESI Y1 and ACT data

The gas density profile around galaxies, shaped by feedback and affecting the galaxy lensing signal, is imprinted on the cosmic microwave background (CMB) by the kinematic Sunyaev-Zel’dovich effect (kSZ). We precisely measure this effect (𝑆/𝑁 ≈ 10) via velocity stacking with 825,283 spectroscopically confirmed luminous red galaxies (LRG) from the Dark Energy Spectroscopic Instrument Year 1 (DESI Y1) survey, which overlap with the Atacama Cosmology Telescope (ACT) Data Release 6 temperature maps over ≥ 4,000 deg 2 . We explore the kSZ dependence with various galaxy parameters and find no significant trend with redshift but clear trends with stellar mass and absolute magnitude in 𝑔, 𝑟, and 𝑧 bands. Our analysis suggests that the gas extends beyond the dark matter halo (99.5% confidence level, i.e., probability to exceed (PTE)⁢ = 0.005). We find a tentative preference for hydrodynamical simulation models with stronger feedback that drives gas further out (Illustris 𝑧 = 0.5, PTE⁢ = 0.37) over weaker-feedback cases (IllustrisTNG 𝑧 = 0.8, PTE⁢ = 0.045), though with limited statistical significance. In all cases, a free multiplicative amplitude was fit to the simulated profiles, and further modeling work is required to firm up these conclusions. We find consistency between kSZ profiles around spectroscopic and photometric LRG, with comparable statistical power, thus increasing our confidence in the photometric analysis. Additionally, we present the first kSZ measurement around the DESI Y1 bright galaxy sample (BGS) and the emission-line galaxies (ELG) whose features match qualitative expectations. Finally, we forecast 𝑆/𝑁 ∼ 50 for future stacked kSZ measurements using data from the ACT, the DESI Y3, and the Rubin Observatory. Furthermore, these measurements will serve as an input for galaxy formation models and baryonic uncertainties in galaxy lensing.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Dark energy survey year 3 results: cosmology from galaxy clustering and galaxy–galaxy lensing in harmonic space

We present the joint tomographic analysis of galaxy-galaxy lensing and galaxy clustering in harmonic space (HS), using galaxy catalogues from the first three years of observations by the Dark Energy Survey (DES Y3). We utilize the redMaGiC and MagLim catalogues as lens galaxies and the metacalibration catalogue as source galaxies. The measurements of angular power spectra are performed using the pseudo-$C_\ell$ method, and our theoretical modelling follows the fiducial analyses performed by DES Y3 in configuration space, accounting for galaxy bias, intrinsic alignments, magnification bias, shear magnification bias and photometric redshift uncertainties. We explore different approaches for scale cuts based on non-linear galaxy bias and baryonic effects contamination. Our fiducial covariance matrix is computed analytically, accounting for mask geometry in the Gaussian term, and including non-Gaussian contributions and super-sample covariance terms. To validate our HS pipelines and covariance matrix, we used a suite of 1800 log-normal simulations. We also perform a series of stress tests to gauge the robustness of our HS analysis. In the $\Lambda$CDM model, the clustering amplitude $S_8 =\sigma _8(\Omega _m/0.3)^{0.5}$ is constrained to $S_8 = 0.704\pm 0.029$ and $S_8 = 0.753\pm 0.024$ (68 per cent C.L.) for the redMaGiC and MagLim catalogues, respectively. For the wCDM, the dark energy equation of state is constrained to $w = -1.28 \pm 0.29$ and $w = -1.26^{+0.34}_{-0.27}$, for redMaGiC and MagLim catalogues, respectively. These results are compatible with the corresponding DES Y3 results in configuration space and pave the way for HS analyses using the DES Y6 data.

(cosmology:) cosmological parameters↗

DOE OSTI2025

Exploring HOD-dependent systematics for the DESI 2024 Full-Shape galaxy clustering analysis

We analyze the robustness of the DESI 2024 cosmological inference from the full shape of the galaxy power spectrum to uncertainties in the Halo Occupation Distribution (HOD) model of the galaxy-halo connection and the choice of priors on nuisance parameters. We assess variations in the recovered cosmological parameters across a range of mocks populated with different HOD models and find that shifts are often greater than 20% of the expected statistical uncertainties from the DESI data. We encapsulate the effect of such shifts in terms of a systematic covariance term, C HOD , and an additional diagonal contribution quantifying the impact of our choice of nuisance parameter priors on the ability of the effective field theory (EFT) model to correctly recover the cosmological parameters of the simulations. These two covariance contributions are designed to be added to the usual covariance term, C stat , describing the statistical uncertainty in the power spectrum measurement, in order to fairly represent these sources of systematic uncertainty. This novel approach should be more general and robust to the choice of model or additional external datasets used in cosmological fits than the alternative approach of adding systematic uncertainties to the recovered marginalised parameter posteriors. We compare the approaches within the context of a fixed ΛCDM model and demonstrate that our method gives conservative estimates of the systematic uncertainty that nevertheless have little impact on the final posteriors obtained from DESI data.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

NASA LSII Lunar Simulant Project

Introduction: America has entered a new era of exploration. NASA’s Artemis program will lead humanity forward to the Moon and prepare us for the next giant leap, the exploration of Mars [1]. To champion technologies needed to live on and explore the Moon, NASA’s Space Technology Mission Directorate (STMD) established the Lunar Surface Innovation Initiative (LSII)[2].LSII's technology development portfolio includes: Utilizing the Moon’s resources; Establishing sustainable surface power; Building machinery and electronics that work in extreme environments, like super-chilly permanently shadowed craters; Mitigating lunar dust; Carrying out surface excavation, manufacturing and construction duties; and Extreme access which includes navigating and exploring the surface/subsurface. To support the development and testing of these technologies, LSII created the lunar simulant project, to create and/or acquire low-, medium-and high-fidelity lunar simulants to match the needs of STMD projects at all levels of technology readiness levels (TRL), as well as other NASA programs. There is not one bulk lunar simulant that will satisfy the needs of all projects. NASA's Approach to Simulants: Just as LSII's activities are being implemented through a combination of unique NASA work and public-private partnerships, NASA will work with commercial simulant providers to acquire simulants that meet NASA's needs. If warranted, NASA will develop simulants using government agencies, as was done with the NASA/USGS Lunar Highlands Type (NU-LHT)series of lunar simulants [3].NASA is also collaborating with the Johns Hopkins University Applied Physics Laboratory(JHUAPL) Lunar Surface Innovation Consortium (LSIC) in the development and characterization of lunar simulants [4]. Within NASA, a small team (< 10 people) is coordinating simulant activities across the agency, with team members located at several NASA centers. The overall objective of the project is to procure lunar simulants in sufficient amounts for earth-based testing of subsystems and systems in a variety of environments (i.e., laboratory, high-bay, thermal-vacuum chambers), required for Artemis missions to the Moon, as well as other missions carrying NASA lunar payloads, such as the Commercial Lunar Payload Services (CLPS)program [5]. Lunar Highlands Simulant: NASA's Artemis Program is targeting the lunar south pole region for initial human missions and the Artemis Base Camp. Hence, the LSII lunar simulant project is currently focusing on the mineralogy and properties of lunar highlands regolith [6 and 7]. Plagioclase-rich rocks (e.g., anorthosite, norite) are the dominant constituent in highlands simulants, with Shawmere, Stillwater, and White Mountain anorthosites being used in commercially available simulants. However, because these feedstocks are terrestrial in nature, they include hydrated minerals, carbon-bearing minerals, and other chemical signatures that are not present on the Moon, and these minor mineral assemblages need to be taken into account when trying to understand test procedures and results. Glass Component: While much attention has been placed on the rock/mineral component of lunar simulants, glass is just as important when creating simulants. The glass component in lunar regoliths is often greater than 50% by volume [8]. This component includes impact melt glass, dark matrix breccias, and agglutinates. However, this component is difficult, time-consuming and expensive to make. Most lunar simulants, past and current, have relied on basaltic cinder as a feedstock for glass. Getting better glass components at a lower cost, particularly agglutinates and glass with an anorthositic composition is a near-term objective that the NASA simulant project is trying to address . Characterization: It is extremely important for lunar simulants to be characterized by several analytical methods. Gruener et al. [9] and JHUAPL LSIC [10] conducted initial assessments of some of the commercially available simulants in 2019 and early 2020, before the global pandemic. Further analyses are needed to better quantify important parameters such as, modal mineralogy and glass content, particle shape, and particle size distribution. These quantified results can then be used in determining figures of merit (FOM) that show how well simulants compare to lunar regolith [11 and 12].

J E Gruener↗

NASA NTRS2009

Developing Lightning Prediction Tools for the CCAFS Dual-Polarimetric Radar

NASA Marshall Space Flight Center and the University of Alabama Huntsville are collaborating with the 45th Weather Squadron (45WS) to develop improved lightning prediction capabilities for the new C-band dual-polarimetric weather radar being acquired for use by 45WS and launch weather forecasters at Cape Canaveral Air Force Station (CCAFS). In particular, these algorithms will focus on lightning onset, cessation and combined lightning-radar applications for convective winds assessment. Research using radar reflectivity (Z) data for prediction of lightning onset has been extensively discussed in the literature and subsequently applied by launch weather forecasters as it pertains to lightning nowcasting. Currently the forecasters apply a relatively straight forward but effective temperature-Z threshold algorithm for assessing the likelihood of lightning onset in a given storm. In addition, a layered VIL above the freezing level product is used as automated guidance for the onset of lightning. Only limited research and field work has been conducted on lightning cessation using Z and vertically-integrated Z for determining cessation. Though not used operationally vertically-integrated Z (basis for VIL) has recently shown promise as a tool for use in nowcasting lightning cessation. The work discussed herein leverages and expands upon these and similar reflectivity-threshold approaches via the application/addition of over two decades of polarimetric radar research focused on distinct multi-parameter radar signatures of ice/mixed-phase initiation and ice-crystal orientation in highly electrified convective clouds. Specifically, our approach is based on numerous previous studies that have observed repeatable patterns in the behavior of the vertical hydrometeor column as it relates to the temporal evolution of differential reflectivity and depolarization (manifested in either LDR or p(sub hv)), development of in-situ mixed and ice phase microphysics, electric fields, and ensuing lightning in the sub-tropical/tropical convection typical of the southeastern U.S., Maritime Continent, and southwestern Amazon. The polarimetric signatures detected in this setting provide a basis for automated 3-D detection of hydrometeor types in fuzzy logic hydrometeor identification algorithms (HID). Our working hypothesis is that improvement in lightning onset warning lead time and specificity for a given storm, relative to application of a Z-threshold algorithm, should arise as a consequence of the ability of dual-polarimetric radar to unambiguously detect and identify (through HID algorithms) the updraft elevation of rain-water cores above the freezing level and subsequent onset of drop freezing, riming, and robust mixed phase processes leading to significant charge separation and lightning. This type of algorithm, though dependent on the quality of the polarimetric data should be less susceptible to variable Z-calibration that can impact a given Z-threshold approach. To facilitate development of the algorithm while the 45WS dual-pol radar is in its current test stages and to evaluate the impact of polarimetric data quality (e.g., modified scan parameters and sampling) on the ensuing algorithms, we are using the ARMOR C-band dual-pol radar in Huntsville combined with N. Alabama LMA data and ARMOR HID algorithms [NCAR algorithm modified for application at C-band] in a testbed fashion. For lightning cessation we are revisiting the application of differential propagation phase variables for the monitoring of ice crystal alignment driven by in-cloud electric fields combined with metrics of ice water path (i.e., vertically integrated reflectivity). Importantly it should be noted that this approach is still very much a research topic and as such, we will explore operational applications that involve radar frequencies other than C-Band by using the UAH MAX X-band dual-pol radar in slow staring modes.

Petersen, W. A.↗

NASA NTRS1996

Grant Proposal for the Continuation of the Voyager Interstellar Mission: LECP Investigation

This proposal documents the plans of the Low Energy Charged Particle (LECP) investigation team for participation in NASA's Voyager Interstellar Mission (VIM) as the Voyager 1 and 2 spacecraft explore the outer reaches of the heliosphere and search for the termination shock and the heliopause. The proposal covers the four year period from 1 January 1997 to 31 December 2000. The LECP instruments on Voyager 1 and 2 measure in situ intensities of charged particles with energies from about 30 keV to 100 MeV for ions, and about 20 keV to greater than 10 MeV for electrons. The instruments provide detailed spectral, angular, and compositional information about the particles. Composition is available for greater than 200 keV/nuc using multi-parameter measurements. Angular information is obtained by a mechanically scanned platform that rotates at various commanded rates. Measurements of low energy ion and electron intensities versus time and spatial location within the heliosphere contain an abundance of information regarding various transport and acceleration processes on both local (approx. 1 hr, approx. 0.01 AU) and global (approx. 11 yrs, approx. 100 AU) scales. The LECP instruments provide unique observations of such dynamical processes, and we anticipate that it will return critical information regarding the boundaries of the heliosphere. Several recent and exciting discoveries based on LECP measurements emphasize the important role that low energy charged particle distributions play in physical processes in the interplanetary medium. Yet, at the same time, these discoveries also underscore the fact that our understanding of processes in the outer heliosphere is, in most cases, incomplete, and in others, only rudimentary at best. Among the discoveries referred to above are the following: (1) Shocks: Examination of greater than 30 keV ion intensities have revealed: (a) a total absence of acceleration beyond only -100-200 keV at a strong transient shock in May 1991 at 35 AU, despite an enhanced level of seed particles; (b) a large transient shock in September 1991 of global scale, with intensities of shock-accelerated ions greater than or equal to 30 keV to approx. 30 MeV showing complex, highly energy-dependent spatial evolution, and small-scale (approx. few gyroradii), often anisotropic, micro-structures; (c) recurrent intensity increases in greater than or equal to 30 keV to -few MeV ions, with structures that, in some cases, show no correlation with the associated corotating shock. (2) Superthermal ion pressure: A global merged interaction region with a leading shock, downstream of which the superthermal ion (greater than or equal to 30 keV to approx. 4 MeV) pressure is comparable to that of the thermal plasma, and the total particle pressure yields a plasma beta of order unity. (3) Pickup ions: Measurements of the C/O ratio within transient structures at 35-45 AU showing the first clear evidence that transient shocks can pre-accelerate interstellar pickup ions from approx. 1 keV/nuc to at least 1 MeV/nuc. (4) Seed particles: Injection of ions for acceleration to high energies at the termination shock is unlikely to be a problem, since interplanetary transient and recurrent shocks are continually accelerating ions, of solar wind or interstellar origin, to highly superthermal energies. (5) Precursor electrons: Ambient solar electrons (greater than or equal to few tens of keV) that exist in the outer heliosphere ca form a broad precursor, several days wide, that is upstream of the termination shock and potentially observable a few months prior to the shock crossing. (6) Solar wind velocity at Voyager 1: We can use LECP ion data to obtain the solar wind velocity at Voyager 1, enabling us to provide critical measurement of the plasma flow as we approach and encounter the termination shock and other regions (necessary due to the partial failure of the Voyager 1 PLS experiment). The work of the LECP investigator team during the VIM will include: (1) Continuing operations with regard to the receipt, processing, verification, cataloging, display, and distribution of the data from the LECP instruments on Voyager 1 and 2, (2) Monitoring the health and performance of the LECP instruments, and evaluating and characterizing the response of the LECP instruments to various energetic particle and plasma environments, (3) Participating in, and supporting Voyager Project planning exercises and other coordinated activities relevant to exploration of the outer heliosphere, (4) Developing analysis techniques and operational procedures suitable for searching for and characterizing the boundaries and unique regions of the outher heliosphere, (5) Continuing the preparation of data sets appropriate for submission to the National Space Sciences Data Center (NSSDC) and, where appropriate, the Planetary Data System (PDS), (6) Maintaining direct Web access to online LECP data through the JHU/APL Voyager LECP home page, (7) Performing scientific evaluations of the Voyager 1 and 2 LECP data sets in conjunction with other data sets and other investigators, with particular focus on the outer regions of the heliosphere, and (8) Publishing the results of these evaluations in the scientific literature and presenting the results in scientific conferences.

Krimigis, Stamatios M.↗

DOE OSTI2026

Cosmological constraints from the DESI DR1 joint power spectrum and bispectrum analysis

We derive cosmological parameter constraints from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) galaxy clustering data, based on a joint full-shape analysis of the power spectrum multipoles and the bispectrum monopole using the ShapeFit framework. This is the follow-up of our previous work, in which we obtained for the first time constraints on the ShapeFit parameters using the bispectrum of DESI DR1. Here we present the first ShapeFit cosmological inference results using the bispectrum of DESI DR1. We recover values for the matter density parameter and Hubble constant of respectively $Ω_m=0.310\pm0.012$ and $H_0=[68.92\pm0.97]\,\mathrm{km\, s^{-1} Mpc^{-1}}$, consistent with previous results from the full DESI DR1 dataset that did not use the bispectrum signal. The inclusion of the bispectrum significantly tightens the constraints on the amplitude of fluctuations, reducing the error-bars in $\ln(A_s\times10^{10})$ by approximately 20%, compared to using the power spectrum alone. We also explore extended cosmological models by performing fits for the evolving dark energy equation of state $w_0w_a$, and the sum of neutrino masses $\sum m_ν$. In these cases, we obtain constraints slightly larger than the ones from previous works from the DESI collaboration, due to not combining the full-shape results with other probes in all tracers. We find no strong evidence of deviations from standard $Λ$CDM, with the dark energy equation-of-state remaining within 2$σ$ from a cosmological constant $Λ$, and the neutrino mass being consistent with the normal hierarchy, $\sum m_ν<0.1\,[eV]$ at 95% confidence limit. These constraints are broadly consistent with other DESI DR1 analyses, thus validating the robustness of the ShapeFit compression approach and the inclusion of the bispectrum for cosmological inference.

Novell-Masot, S. [ICC, Barcelona U.; Geneva U., De↗

DOE OSTI2026

Cosmological constraints from the DESI DR1 joint power spectrum and bispectrum analysis

Novell-Masot, S. [ICC, Barcelona U.; Geneva U., De↗

DOE OSTI2024

GeoThermalCloud: Cloud Fusion of Big Data and Multi-Physics Models using Machine Learning for Discovery, Exploration, and Development of Hidden Geothermal Resources

The primary goals of this project are exploring hidden geothermal resources in the U.S.A. and designing profitable enhanced geothermal systems (EGS). Many processes and parameters control geothermal exploration and energy production from geothermal fields. Diverse datasets (e.g., geology, geochemistry, geophysics, satellite, airborne geophysics) are available to help characterize subsurface geothermal conditions. Sparse and multi-scale characteristics of these datasets prohibit properly leveraging these datasets for geothermal exploration and profitable EGS design. Recent advancements in machine learning (ML) promise to resolve these issues. The tremendous challenges and risks of geothermal exploration and production bring the demand for novel ML methods and tools that can (1) analyze large field datasets, (2) assimilate model simulations (large inputs and outputs), (3) process sparse datasets, (4) perform transfer learning (between sites with different exploratory levels), (5) extract hidden geothermal signatures in the field and simulation data, (6) label geothermal resources and processes, (7) identify high-value data acquisition targets, and (8) guide geothermal exploration and production by selecting optimal exploration, production, and drilling strategies. To address these necessities, ML-based geothermal resources exploration and enhanced geothermal systems (EGS) design tools have been developed. The exploration tool is called GeoThermalCloud and EGS design tool is called GeoDT-ML. GeoThermalCloud (https://github.com/SmartTensors/GeoThermalCloud.jl) utilizes a LANL unsupervised ML platform called SmartTensors (https://tensors.lanl.gov/) to automate data analyses and interpretations by extracting hidden signatures to identify geothermal prospects. Also, it enables the identification of critical measurements needed to identify geothermal resource signatures. Alternatively, GeoDT-ML (https://github.com/SmartTensors/GeoThermalCloud.jl/tree/master/EGS) is an ML-based alternative to GeoDT (https://github.com/GeoDesignTool/GeoDT.git), a fast, simplified multi-physics solver to evaluate EGS project designs in uncertain geologic systems. GeoDT-ML leverages recent advances in deep learning and high-performance computing. It is a faster and simpler version of GeoDT. To make this project a success, we used capabilities of LANL, PNNL, Google, Stanford, and Julia Computing. We analyzed eight datasets of the U.S.A. using GeothermalCloud and demonstrated potential highly prospective geothermal resources and identified key factors defining highly prospective sites. The first data set includes 44 locations in southwest New Mexico and 18 geological, hydrogeological, geophysical, geothermal, geochemical attributes. We defined low- and medium-temperature hydrothermal systems and discovered a new highly prospective site. The second data set analyzed 18 shallow water chemistry attributes at 14,342 locations in the Great Basin. It demarcated modestly, moderately, and highly prospective sites including key attributes for each type of prospectivity. The third data set analyzed Utah FORGE data including satellite (InSAR), geophysical (gravity, seismic), geochemical, and geothermal attributes. Here, we performed prospectivity analysis to identify future drilling locations using geological, geochemical, and geophysical attributes. Maps of temperature at depth and heat flow are constructed based on the available data. Prospectivity maps were generated, and drilling locations were proposed for future geothermal field exploration. The fourth data set analyzed 21 attributes at 120 locations in Tularosa Basin, New Mexico; data comes from past play fairway analyses in this region. ML analyses identified geothermal signatures associated with modestly, moderately, and highly hydrothermal systems. We also defined dominant attributes and spatial distribution of the geothermal signatures. The fifth, sixth, seventh, and eighth datasets include Tohatchi Springs, New Mexico, Hawaii, Brady site, Nevada, and EGS Collab, respectively. Moreover, we coupled GeothermalCloud and magnetotellurics data to pinpoint drilling locations for developing geothermal projects in the Tularosa Basin, New Mexico. GeothermalCloud found potential prospective locations for geothermal resources near White Sands Missile Range and McGregor Range at Fort Bliss. Magnetotellurics data determined the potential depth (~1800m) of geothermal prospects at McGregor Range based on apparent resistivity structures/layers in the subsurface. The McGregor Range consists of three resistivity layers and two resistivity structures. Magnetotellurics data also helps identify that the western portion of the McGregor Range has thick and low-resistivity earth materials. The low resistivity to the west is most likely for a fault system. Assuming temperature is consistent with a geothermal reservoir, the west-central part of the McGregor Range has the highest geothermal potential because of the increase in porosity and associated permeability attributed to the interpreted fault system. Also, we devised a coupling strategy between a process model and GeothermalCloud to characterize hydrogeological conditions and geothermal conditions, respectively. The process model characterizes hydrogeological and geothermal conditions on highly prospective geothermal sites provided by GeothermalCloud. We developed a physics-informed neural network (PINN) version of the Burns equation that can be easily coupled with GeothermalCloud. Furthermore, we performed an optimal design decision maximizing the economic value of an EGS power plant. This study optimized the range of well spacing between injection and production wells maximizing net present value in dollars (NPV). For this task, we used the GeoDT to simulate the Utah FORGE EGS development cycle from the initial well design to the end of production. Next, we accomplished another crucial task, which is predicting permeability of geothermal reservoirs. Predicting permeability of geothermal reservoirs is a non-trivial task because of huge computational runtime of simulation and lack of measurements. To avoid these limitations, we used easy-to-measure chemical concentrations in the subsurface as measurement data and convolutional neural network based ML model of a high-fidelity model. Next, we predicted permeability using Markov chain Monte Carlo simulation. We found that Markov chain Monte Carlo simulation predicts permeability with a high certainty if the prediction zone in the simulation area has chemical concentration data. Finally, we analyzed the DOE funded INGENIOUS and GeoDAWN projects data. For discovering hidden geothermal systems in the Great Basin, the INGENIOUS project accumulated old data, collected new data, and released them in 2022. The dataset includes a total of 24 geological, geophysical, and geochemical attributes. Data resolution and scale significantly vary prohibiting an appropriate usage. To avoid such limitations, we brought all data in the same resolution and scale by applying the inverse distance weighting interpolation technique for predicting data in unsampled locations. Subsequently, we analyzed LiDAR data of the GeoDAWN project. We received data in tiles format. The DOE’s overarching goal is to use ML on LiDAR data for finding favorable geological structures (e.g., step up faults in Brady, Nevada). To serve the purpose, we need to label favorable geologic structures that correspond to LiDAR data. We wrote an algorithm to label the LiDAR data with the favorable geologic structures.

15 GEOTHERMAL ENERGY↗

NASA NTRSDate not supplied

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

NASA NTRS2013

Dynamics and Control of a Disordered System in Space

In this paper, we present some ideas regarding the modeling, dynamics and control aspects of granular spacecraft. Granular spacecraft are complex multibody systems composed of a spatially disordered distribution of a large number of elements, for instance a cloud of N grains in orbit, with N greater than 10(exp 3). These grains can be large (Cubesat-size) or small (mm-size), and can be active, i.e., a fully equipped vehicle capable sensing their own position and attitude, and enabled with propulsion means, or entirely passive. The ultimate objective would be to study the behavior of the single grains and of large ensembles of grains in orbit and to identify ways to guide and control the shape of a cloud composed of these grains so that it can perform a useful function in space, for instance, as an element of an optical imaging system for astrophysical applications. This concept, in which the aperture does not need to be continuous and monolithic, would increase the aperture size several times compared to large NASA observatories such as ATLAST, allowing for a true Terrestrial Planet Imager that would be able to resolve exo-planet details and do meaningful spectroscopy on distant world. In the paper, we address the modeling and autonomous operation of a distributed assembly (the cloud) of large numbers of highly miniaturized space-borne elements (the grains). A multi-scale, multi-physics model is proposed of the dynamics of the cloud in orbit, as well as a control law for cloud shape maintenance, and preliminary simulation studies yield an estimate of the computational effort, indicating a scale factor of approximately N(exp 1.4) as a function of the number of grains. A granular spacecraft can be defined as a collection of a large number of space-borne elements (in the 1000s) designed and controlled such that a desirable collective behavior emerges, either from the interactions among neighboring grains, and/or between the grains and the environment. In this paper, each grain is considered to be a highly miniaturized spacecraft which has limited size and mass, hence it has limited actuation, limited propulsive capability, limited power, limited sensing, limited communication, limited computational resources, limited range of motion, limited lifetime, and may be expendable. The modeling and dynamics of clouds of vehicles is more challenging than with conventional vehicles because we are faced with a probabilistic vehicle composed of a large number of physically disconnected vehicles. First, different scales of motion occur simultaneously in a cloud: translations and rotations of the cloud as a whole (macro-dynamics), relative rotation and translation of one cloud member with respect to another (meso-dynamics), and individual cloud member dynamics (micro-dynamics). Second, the control design needs to be tolerant of the system complexity, of the system architecture (centralized vs. decentralized large scale system control) as well as robust to un-modeled dynamics and noise sources. Figure 1, top left, shows the kinematic parameters of a 1000 element cloud in orbit. The motion of the system is described with respect to a local vertical-local horizontal (LV-LH) orbiting reference frame (x,y,z)=F(sub ORF) of origin O(sub ORF) which rotates with mean motion omega and orbital semi-major axis R(sub 0). The orbital geometry at the initial time is defined in terms of its six orbital elements, and the orbital dynamics equation for point O(sub ORF) is propagated forward in time under the influence of the gravitational field of the primary and other external perturbations, described below. The origin of this frame coincides with the initial position of the center of mass of the system, and the coordinate axes are z along the local vertical, x toward the flight direction, and y in the orbit normal direction. The assumptions we used to model the dynamics are as follows: 1) The inertial frame is fixed at Earth's center. 2) The orbiting Frame ORF follows Keplerian orbit. 3) the cloud system dynamics is referred to ORF. 4) the attitude of each grain uses the principal body frame as body fixed frame. 5) the atmosphere is assumed to be rigidly rotating with the Earth. Regarding the grains forming the cloud: 1) each grain is modeled as a rigid body; 2) a simple attitude estimator provides attitude estimates, 3) a simple guidance logic commands the position and attitude of each grain, 4) a simple local feedback controller based on PD control of local states is used to stabilize the attitude of the vehicle. Regarding the cloud: 1) the cloud as a whole is modeled as an equivalent rigid body in orbit, and 2) an associated graph establishes agent connectivity and enables coupling between modes of motion at the micro and macro scales; 3) a simple guidance and estimation logic is modeled to estimate and command the attitude of this equivalent rigid body; 4) a cloud shape maintenance controller is based on the dynamics of a stable virtual truss in the orbiting frame. Regarding the environmental perturbations acting on the cloud: 1) a non-spherical gravity field including JO (Earth's spherical field) zonal component, J2 (Earth's oblateness) and J3 zonal components is implemented; 2) atmospheric drag is modeled with an exponential model; 3) solar pressure is modeled assuming the Sun is inertially fixed; and 4) the Earth's magnetic field is model using an equivalent dipole model. The equations of motion are written in a referential system with respect to the origin of the orbiting frame and the state is propagated forward in time using an incremental predictor-corrector scheme. A representative cloud with varying number of grains is simulated to identify the limitations in computation time as the number of grains grows. We derive a control law to track a desired surface in the ORF (equivalently to maintain a reference cloud shape) by defining an error from a desired surface shape, and designing a control law that is exponentially stable and reduces the tracking error to zero. Figure 1 (top right) shows a comparison of various requirements for simulation of single spacecraft vs. granular spacecraft, indicating the high degree of complexity that needs to be taken into consideration. The ORF components of control force required by one of the grains is, for this particular case, in the micro-Newton range. However, no attempt has been made yet to reconfigure (or re-orient) the cloud configuration internally, for which forces in the milli-Newton level are expected, depending on the time required to do the reconfiguration. Figure 1, bottom, shows the computation time as a function of the number of grains, indicating an order N(exp 1.43) scaling on a 8 Gb, 1067 MHz RAM MacOSX computer with a 3.06 GHz Intel Core 2 Duo processor. With this metric, the same simulation for a system of N=1000 grains would take 5.4 hours, and 146 hours (i.e., 6 days) for a system with N=10,000 grains. Therefore, efficient ways to simulate this complex system, where not only the time scales of natural system dynamics, but also the sampling times of the Guidance, Navigation, and Control are included, remain to be explored. Additional details on the cloud modeling, dynamics, and control will be described in the paper.

simulation↗